Sparse mass matrix for Radau?

  • Juan Pablo Carbajal


    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.

  • Rob Clewley

    Rob Clewley - 2013-03-18

    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?


Get latest updates about Open Source Projects, Conferences and News.

Sign up for the SourceForge newsletter:

No, thanks