Mauro, I would be curious to learn details of what you did.
How did you define the reference solution? I presume it is the value of the solution at some point in time? Did you use extrapolation to the limit? How did you measure the number of significant digits? In particular, for the fixed--step integrators one would not expect the increase in accuracy measured in significant digits to be a nice and smooth function of the number of steps (as a proxy for the wall-clock time). I know I could probably find out all of this by rooting around in your code, but I think it's going to be more efficient, if you have the time, to simply ask for an explanation. Thanks a lot, Petr On Friday, January 30, 2015 at 6:03:24 AM UTC-8, Mauro wrote: > > Hi all, > > those of you who solve initial value problems (IVP) of ordinary and > algebraic differential equations (ODE/DAE) might be interested in: > https://github.com/mauro3/IVPTestSuite.jl It provides several test cases > for ODE and DAE solvers based on previous, well-known test sets. It can > easily be adapted to run with your solvers. > > I ran the tests for the three ODE-solver packages ODE.jl, DASSL.jl and > Sundials.jl. Results are here: > https://github.com/mauro3/IVPTestSuite.jl/blob/master/results/results.md > > I found: > - DASSL.jl seems to be as capable as the Sundials solvers if sometimes > slower but also sometimes faster. > - For large(-ish) systems, e.g. stemming from PDEs, DASSL.jl seems to > be the best and fastest choice at the moment because of the patchy > support of > sparse Jacobians in Sundials. (please correct me if wrong). > - ODE.ode23s does ok too but is generally a lot slower than DASSL, > presumably because of its lower order. > > However, take those results with a grain of salt as I haven't spent too > much time optimising yet. > > -- Mauro >
