[Bug-gsl] ode-initval2 redundant user function calls

[Tuomo Keskitalo]

Hi,

On 12/22/2011 09:55 PM, Farkas, Illes wrote:

> One more question. According to test.out.txt the r.h.s. function is
> accessed 3 or 2 times at each time point. Can this be reduced to accessing
> only once (in order to speed up rk2imp) ?

I finally managed to look into this. I think the culprit is evolve_apply 
routine, which does not reuse derivative data from previous accepted 
step. The attached patch should remedy the situation, and yield some 
gain in efficiency. Thanks for the bug report!

Brian, could you please apply the patch, thanks!

--
address@hidden
http://iki.fi/tuomo.keskitalo

odeiv2-evolve.bzr.patch
Description: Text Data