|
From: | Tuomo Keskitalo |
Subject: | Re: [Bug-gsl] Re: [bug #30947] Please, include fixed step size control object for ode suite |
Date: | Fri, 03 Sep 2010 12:33:20 +0300 |
User-agent: | Mozilla/5.0 (X11; U; Linux x86_64; en-US; rv:1.9.1.11) Gecko/20100713 Thunderbird/3.0.6 |
Hi, On 09/03/2010 11:36 AM, Sergey B Kirpichev wrote:
Also, I'm afraid that using your cfxd may not guarantee the use of a fixed step size; also evolve_apply can change the step length, so you might have to have a function something like "evolve_apply_fixed_step_size" to do the job.Why? fc_control_hadjust() always returns GSL_ODEIV_HADJ_NIL - there is no reason to change step size.
Please see function gsl_odeiv2_evolve_apply in evolve.c. If step_apply returns with a failure, then a smaller step size is attempted. evolve_apply tries hard to make a step forward, even if it has to use a smaller step size than the user or previous step suggested.
-- address@hidden http://iki.fi/tuomo.keskitalo
[Prev in Thread] | Current Thread | [Next in Thread] |