[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

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

From: Brian Gough
Subject: Re: [Help-gsl] Performance issue with ODE solvers
Date: Mon, 03 Dec 2007 14:35:30 +0000
User-agent: Wanderlust/2.14.0 (Africa) Emacs/22.1 Mule/5.0 (SAKAKI)

At Fri, 30 Nov 2007 10:58:38 -0500,
Sebastien Maret wrote:
> I am rather confused by these differences. The fact that both programs
> give the same results make me think that it's not a bug in my program,
> but a real performance difference between the two solvers. Maybe the
> GSL solvers are not suited to solve such a large systems of stiff
> equations? I have tried to use other implicit solvers from the GSL
> (except for bsimp, which require the Jacobian), and they all have
> similar performance.


A couple of questions: 

Are the overall error goals set in the same way?

If you increase the gsl error tolerance by a factor of 10 do the
solutions deviate significantly?

Brian Gough
(GSL Maintainer)

Network Theory Ltd,
Publishing the GSL Manual -

reply via email to

[Prev in Thread] Current Thread [Next in Thread]