Provide own (non)-linear solvers

[hrittich]

Hello,

is it possible to use user defined linear and non-linear solvers for stiff systems, especially sparse solvers? I would like to use odeint for solving partial differential equations. Hence, the Jacobians that I consider are large and sparse matrices.

Regards

[bjodah]

I think it's on their roadmap (cf. #169).

Until odeint supports user provided linear solvers, I can personally recommend SUNDIALS
(versions sundials >= 3 even support user supplied direct solvers, and sundials<=2.7.0 accept user defined preconditioners).

[hrittich]

Thank you for the hint.