This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Error scaling in ODE-IV
- From: Dave Auslander <dma at me dot berkeley dot edu>
- To: gsl-discuss at sources dot redhat dot com
- Date: Tue, 23 Jul 2002 13:08:31 -0700
- Subject: Error scaling in ODE-IV
I am using the ODE-initial value routines from GSL. In looking at section
24.3 or the GSL manual, I note that there does not appear to be any scaling
array for the error calculation in the adaptive algorithms. eps_rel and
eps_abs are scalars. Most ODE systems, however, are not naturally scaled so
all components have the same scale so it is very difficult to optimize
performance.
Am I missing something or is this the intended usage?
Thanks
Dave Auslander
PS - many thanks to the GSL developers. This is a huge, very useful effort.
+++++++++++++++++++++++++++++++++++++++++++++++
David M. Auslander, Prof. of ME and Assoc. Dean
Mechanical Engineering, Univ of Calif, Berkeley, CA 94720-1740
Dean's office: 510-642-7594 ME: 510-642-4930 Fax: 510-643-5599
dma@me.berkeley.edu www.me.berkeley.edu/~dma
+++++++++++++++++++++++++++++++++++++++++++++++