Writing extendable and hardware-agnostic GPU libraries

Very interesting. How does this deal with data regions? Do you have an array type that explicitly recides on GPU, i.e. to pass it to a CPU routine like file i/o you need to wrap CUDAmemcopy? or does it manage the device state?

How does the current Julia GPU support compare to other high level GPU languages like Futhark and Halide?