Pointer to CuArray as a parameter for DifferentialEquations right-hand-side function

[ChrisRackauckas]

@fedoroff

Now I would like to solve this kind of equation on GPU with f being a CuArray. In this situation func will be converted to a GPU kernel function (though I do not know how it actually done in DifferentialEquations, but I try to write my own implementation where it happens). The first problem is the fact that CuArray f is not isbits and therefore it can not be accessed inside func . I can solve this issue by converting f into StaticArray. However, in this case I will loose the advantage of changing f without remaking the prob .

You might want to look at DiffEqGPU.jl where this is done by KernelAbstractions.jl