[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Bug-gsl] ode-initval/bsimp is always high order

From: Tom Quinn
Subject: [Bug-gsl] ode-initval/bsimp is always high order
Date: Mon, 20 Sep 2010 13:18:28 -0700 (PDT)
User-agent: Alpine 1.10 (LRH 962 2008-03-14)

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:       address@hidden
Phone:          206-685-9009

reply via email to

[Prev in Thread] Current Thread [Next in Thread]