Is there a way to have the c-code for the
void massMatrix(unsigned n_, unsigned np_, double t, double *Y_, double *p_, double **f_, unsigned wkn_, double *wk_, unsigned xvn_, double *xv_)
to use a sparse representation?
I can't find this on the manual, I am sorry if the answer is there.
I don't believe there is any support for sparse matrices with Radau. But you can search the full Fortran code and docs on Hairer & Wanner's website. I think it's linked from the PyDSTool documentation online. Sparse solvers are of interest to us in future development, though. Maybe you can tell me more about your use case and alternative sparse solving tools that you would like to use with PyDSTool?
Sign up for the SourceForge newsletter:
You seem to have CSS turned off.
Please don't fill out this field.