Re: Contributing to GSL

[Tuomo Keskitalo <Tuomo.Keskitalo@iki.fi>]

Hello,

the TODO in ode-initval2 directory is slightly outdated. A BDF method 
has already been implemented in msbdf stepper in gsl-1.15. The selection 
of steppers in ode-initval2 is pretty nice, I think, as far as dense 
general-purpose ode-iv methods go. I sometimes thought about trying to 
implement Radau 5 method described in the second book (Hairer, E., 
Wanner, G., Solving Ordinary Differential Equations II: Stiff and 
Differential-Algebraic Problems) but never did. If you want to fight 
with ODEs, you could try that. However, I cannot promise I'd have much 
time to contribute.

DAEs (differential-algebraic equations) are a different story, even 
though they're close to ODEs. There's currently no DAE solver framework 
in GSL. It is something that would be very nice to have, though. If 
anyone is willing to try to design and implement one, feel free to do it!

BR,
Tuomo

On 03/28/2012 07:11 AM, Varad Deshmukh wrote:

> Hi,
> I am Varad Deshmukh, a Master's student at UCSB majoring in Computer
> Science. My research interests are in Scientific Computing and HPC, and
> I am interested in contributing to GSL. I recently took up a course in
> numerical simulation and I figured a good way to start would be looking
> into ODE solvers.
>
> The TODO list mentions work to be done for adding other stepper
> methods - especially the BDF methods for DAE's implemented in DASSL and
> DASPK. Could you please elaborate more on which new methods you want to
> focus on? I have been reading the code for the Runge-Kutta methods - the
> driver, evolve and control formats - and have a rough idea about the
> design.
>
> It would be great if I could get some help regarding GSL development in
> general and ode-initval2 in particular. If you feel that I should work
> with smaller bugs and fix them before going into a new implementation, I
> would like to hear suggestions.
>
> Thank you,
> Varad.
>


-- 
Tuomo.Keskitalo@iki.fi
http://iki.fi/tuomo.keskitalo