[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] Bug in odeiv integrator
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] Bug in odeiv integrator |
Date: |
Thu, 16 Apr 2009 20:31:53 +0100 |
User-agent: |
Wanderlust/2.14.0 (Africa) Emacs/22.2 Mule/5.0 (SAKAKI) |
At Wed, 15 Apr 2009 22:47:08 +0300,
Taneli Kalvas wrote:
> The behaviour of gsl_odeiv_evolve_apply() differs from documentation.
> The documentation says:
>
> "If the user-supplied functions defined in the system dydt return a
> status other than GSL_SUCCESS the step will be aborted. In this case, t
> and y will be restored to their pre-step values and the error code from
> the user-supplied function will be returned."
>
> The y values are restored, but the t value is not. This happens in case,
> where the control function has first rejected a step because of too
> large error and when taking a new step, the user function returns an
> error. See the following print for a simple example:
Hello,
Thanks for the bug report and example program.
I will take a look at that.
--
Brian Gough
(GSL Maintainer)
Support freedom by joining the FSF!
http://www.fsf.org/associate/support_freedom/join_fsf?referrer=37