commons-issues mailing list archives

Site index · List index
Message view « Date » · « Thread »
Top « Date » · « Thread »
From "Luc Maisonobe (JIRA)" <>
Subject [jira] [Commented] (MATH-172) implementation of implicit Runge-Kutta method for stiff differential equations
Date Tue, 30 Aug 2011 18:47:46 GMT


Luc Maisonobe commented on MATH-172:

Yes, there have been some work on it. Nothing has been committed yet.
The work was based on the patch attached to this issue, which is an implementation of BDF
(Backward Differential Formula). There are some stability problems with Adams methods and
the BDF formula shares (or could share) some code with Adams.
I first thought simply computing the new point by applying the formulas from high order to
low order (i.e. from small to large numbers) would be sufficient to solve these numerical
instabilities. I was wrong. Another attempt could be made now that we have MathUtils.linearCombination,
which can compute accurate linear combinations even when cancellations occurs, which by construction
do happen in these algorithms.

I would really much like to have stiff ODE solver in 3.0.

If you want to help solving this issue, do not hesitate to propose something.

> implementation of implicit Runge-Kutta method for stiff differential equations
> ------------------------------------------------------------------------------
>                 Key: MATH-172
>                 URL:
>             Project: Commons Math
>          Issue Type: Improvement
>    Affects Versions: 1.2
>            Reporter: James Housden
>            Assignee: Luc Maisonobe
>             Fix For: 3.0
>         Attachments: bdf-2009-07-07.patch
> In version 1.2 several explicit Runge-Kutta methods have been implemented which work
fine for non-stiff problems. However these methods are not suitable for stiff equations. The
implementation of an implicit Runge-Kutta method such as Radau5 would extend the functionality
of the ode package to cover the solution of stiff ordinary differential equations.

This message is automatically generated by JIRA.
For more information on JIRA, see:


View raw message