Re: [Help-gsl] Performance issue with ODE solvers

[Brian Gough]

At Wed, 2 Jan 2008 17:51:57 +0100,
Jochen Küpper wrote:
> 
> Sebastien,
> 
> have you received any replies, comments?
> Do you meanwhile have more information on the issue?
> 
> Have people made similar observations for other problems?

The main thing I'd need to know before looking at it more is how the
user-specified error tolerances are used by each package -- otherwise
it's not clear that the runtimes can be compared on a like-for-like
basis. If that is ok, it may be that GSL routines are overestimating
the error, and hence taking smaller steps.

> Sundials is BSD licensed - are there ways to (legally) incorporate it  
> into GSL? Brian, how would we proceed to ask the original authors for  
> a contribution to GSL? What's your opinion on that?

The license is compatible, so it can be used side-by-side with GSL as
an independent packages.  I think that is the best way to use it.

> On the other hand, there are several articles in scientific journals,  
> including a long one in ACM TOMS. Anybody up to implementing it in  
> GSL based on these descriptions? Actually, if no-one is taking up on  
> it, maybe with a found background in ODEs could mentor this as a  
> Google summer of code project?

If somebody wanted to implement some of the algorithms in the GSL
framework that would be cool.

-- 
Brian Gough