Is there a Mathematica CUDA ode solver?

I have a system of ODEs which I’d like to solve with many (thousands) different initial conditions. In the net there are articles that CUDA can speed up such computations, but I haven’t found anything for Mathematica. So is there something around?