Hello, that's an interesting use case. May I ask why you needed to integrate backwards? Did you use driver_apply for integration?
I haven't really considered anyone "hotswapping" the direction of integration using driver functions.. I think that the stepper at least must be reset if integration direction is changed, to be safe with multistep methods. This needs some thinking and testing. I can try to see to this in a week or few. Tuomo On Fri, Mar 2, 2012 at 6:13 PM, Patrick Alken <[email protected]> wrote: > I had a problem recently where I needed to integrate forward along a curve > (positive step size h) and then integrate backward again (negative h). But > there is currently no way to change the sign of the step size h using the > driver routines, so I had to modify driver->h directly. > > I propose adding a routine gsl_odeiv2_driver_set_h() to allow the user to > reset the step size to what they want. Tuomo: do you see any negative issues > with this? I can make the addition if you want. > > Thanks, > Patrick -- [email protected] http://iki.fi/tuomo.keskitalo
