Re: [Bug-gsl] Bug in odeiv integrator

[Brian Gough]

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