On May 5, 2:56 am, Guilherme <guito...@gmail.com> wrote:
> IDA is shown as the last on the series DASSL -> DASPK -> IDA. > Sundials is C ++(?) coded. > > Could you handle Mass matrix with DASSL? DASSL handles implicit ODEs of the form F(t, x, x') = 0, so I think it would work for non-constant mass matrices too, but an adapted solver might be better. I quite like your suggestion of IDA through the Sundials package. It is recommended by the experts, the interface is clean (DASSL basically consists of one Fortran function which takes 15 arguments, all of which have to be in a very specific form) and is relatively well updated (the last version dates from 2006). It is also licensed under the BSD license. On the downside, I think that writing a Python wrapper around even a small part of its functionality might prove to be quite an undertaking. By contrast, DASSL is much easier to understand. I will try to put my C++ wrapper on github in a week or so, after I've cleaned everything up, so you can have a look. All the best, J. -- To post to this group, send an email to sage-devel@googlegroups.com To unsubscribe from this group, send an email to sage-devel+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-devel URL: http://www.sagemath.org
