Hi,
this is a nice project. Thank you for sharing.
One nice feature of CVODE/SUNDIALS is the efficient calculation of the sensitivities (gradients with respect to parameters using a set of extra equations). It would be great if this functionality could also be available via the python interface.
(Or if it is already available the documentation is missing)
best,
lohox