Non-contiguous indexing

[juliagpu]

Non-contiguous indexing

[MikeInnes]

This actually works if the index is a CuArray:

julia> xs[cu[1,2], :]
2×3 CuArray{Float64,2}:
 0.213428  0.977906  0.12309 
 0.416557  0.371804  0.730723

Should be pretty easy to do that conversion automatically.

[MikeInnes]

Actually, I remember why I didn't add this earlier: Base uses scalar indexing to check bounds, so converting to CuArray will be bad for performance.

Easy to fix by doing the conversion after bounds checking; for now, ironically, cpu-array indexing will be faster, but we can fix that if cuarray indexing is needed.