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 _______________________________________________ Help-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/help-gsl
