Thrust and stiff ODE solvers in odeint

I am studying a parameterized system of stiff ODEs and need to solve it for many values of its parameters. I am working in C/C++/CUDA and I want to solve these systems in parallel.

I was wondering if I could use odeint for this. I found this example, but according to this page of the documentation, stiff solvers require ublas types.

Does it mean that their stiff solvers are incompatible with Thrust?