[Top][All Lists]

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

Re: [Help-gsl] [ODE] GSL and Octave comparison

From: Brian Gough
Subject: Re: [Help-gsl] [ODE] GSL and Octave comparison
Date: Sat, 26 Aug 2006 10:56:11 +0100
User-agent: Wanderlust/2.14.0 (Africa) Emacs/21.3 Mule/5.0 (SAKAKI)

At Thu, 24 Aug 2006 21:24:58 +0200,
Marco Maggi wrote:
> and then recomputing with GSL step function rk2
> control function 'y' eps_abs = eps_rel = 0.001,

Those are large tolerances -- they are applied to each step, so the
global error is O(Nsteps * tolerance).  Also rk2 is the least accurate
method, try one of the others, e.g. rk4 at least.

Brian Gough

Network Theory Ltd,
Publishing the GSL Manual -

reply via email to

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