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.
If you would like to refer to this comment somewhere else in this project, copy and paste the following link:
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?
If you would like to refer to this comment somewhere else in this project, copy and paste the following link:
Hi,
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?