Make OPSA-CUDA work for non-sorted indices

[frostedoyster]

At the moment, it only works with sorted indices.
This is not high-priority, but the main issue is that it doesn't error out if the indices are not sorted, it just gives wrong results