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
