The bsimp_alloc() routine uses a very small EPSILON to set the order of
the integrator, and there seems to be no interface to change this.
The result is that even for small steps and large error tolerances, a
minimum of 140 function evaluations are taken per timestep. For my
application, this is more than an order of magnitude more than are
necessary.
Tom Quinn Astronomy, University of Washington
Internet: [email protected]
Phone: 206-685-9009
_______________________________________________
Bug-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/bug-gsl
