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

[Brian Gough]

URL:
  <http://savannah.gnu.org/bugs/?31109>

                 Summary:  ode-initval/bsimp is always high order
                 Project: GNU Scientific Library
            Submitted by: bjg
            Submitted on: Thu 23 Sep 2010 02:23:29 PM BST
                Category: Performance
                Severity: 3 - Normal
        Operating System: 
                  Status: Confirmed
             Assigned to: None
             Open/Closed: Open
                 Release: 1.14
         Discussion Lock: Any

    _______________________________________________________

Details:

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 to this item at:

  <http://savannah.gnu.org/bugs/?31109>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/