New Version of hmatrix (more ODE solvers)

Previous Topic Next Topic
classic Classic list List threaded Threaded
1 message Options
Reply | Threaded
Open this post in threaded view

New Version of hmatrix (more ODE solvers)

Dominic Steinitz-2

I am pleased to announce a new version of hmatrix:

  • This is not intended to be a breaking change but a lot of modules have been modified to ensure that continuous integration (which has now been set up) is green.
  • Support for SUNDIALS has been added. It should be possible to replace Numeric.GSL.ODE with Numeric.Sundials.ARKode.ODE and have your program work as before bearing in mind that the methods and error control might differ (even for those with the same names!).

The packages that comprise hmatrix are:

  • hmatrix
  • hmatrix-gsl
  • hmatrix-special
  • hmatrix-sparse
  • hmatrix-glpk
  • hmatrix-tests
  • hmatrix-sundials

I am planning to split these up into separate projects to make it easier for contributors. To this end, I have created an hmatrix organisation in github. Please contact me if you wish to become a member.

I strongly recommend anyone solving ODEs to look as hmatrix-sundials as the underlying package provides many more methods, more documentation and more diagnostic information (great when your 200 variable model fails). It is a truism that any sufficiently sophisticated ODE solver package is a domain specific language (I include the ODE solving part of GSL here). The current interface both to GSL and to SUNDIALS is somewhat rudimentary but will hopefully evovle to a nicer DSL (something at which Haskell excels).

Dominic Steinitz
Twitter: @idontgetoutmuch

Haskell mailing list
[hidden email]