help-gsl
[Top][All Lists]
Advanced

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

Re: [Help-gsl] solving ODE in a loop


From: Frank Reininghaus
Subject: Re: [Help-gsl] solving ODE in a loop
Date: Tue, 9 Sep 2008 19:40:23 +0200
User-agent: KMail/1.9.10

Hi,

On Tuesday 09 September 2008 19:04:34 Lado Samushia wrote:
> I am trying to solve am ode system in a double loop (for different values
> of parameters).
> ...
> It works fine for the first time, but brakes down on second (and
> consequent) loop executions (gives "nan" for function values and jumps
> trough interval in one step).

I think you should call

gsl_odeiv_step_reset (s);
gsl_odeiv_evolve_reset (e);

before each loop execution, see

http://www.gnu.org/software/gsl/manual/html_node/Stepping-Functions.html
http://www.gnu.org/software/gsl/manual/html_node/Evolution.html

I hope that helps!

Regards
Frank




reply via email to

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