bugGNU Scientific Library - Bugs: bug #26255, gsl_odeiv_evolve_apply fails to...

 
 

You are not allowed to post comments on this tracker with your current authentication level.

bug #26255: gsl_odeiv_evolve_apply fails to restore t on step failure

Submitted by:  -Deleted Account- <bjg>
Submitted on:  Mon 20 Apr 2009 07:29:48 PM UTC  
 
Category: Runtime errorSeverity: 3 - Normal
Operating System: Status: Fixed
Assigned to: -Deleted Account- <bjg>Open/Closed: Closed
Release: 1.12

Tue 21 Apr 2009 01:36:06 PM UTC, comment #1:

Fixed by
9a69804a1b183005230eaa8fa46376e00ba2b7ce ad008ec33f934ada1ecc71432d6ec8181ad1066b

-Deleted Account- <bjg>In charge of this item.
Mon 20 Apr 2009 07:29:48 PM UTC, original submission:

From: Taneli Kalvas <taneli.kalvas@jyu.fi>
To: -unavailable-
Subject: [Bug-gsl] Bug in odeiv integrator
Date: Wed, 15 Apr 2009 22:47:08 +0300

Hi!

The behaviour of gsl_odeiv_evolve_apply() differs from documentation.
The documentation says:

"If the user-supplied functions defined in the system dydt return a
status other than GSL_SUCCESS the step will be aborted. In this case, t
and y will be restored to their pre-step values and the error code from
the user-supplied function will be returned."

The y values are restored, but the t value is not. This happens in case,
where the control function has first rejected a step because of too
large error and when taking a new step, the user function returns an
error.

It is easy to get around this problem by resetting the t in case of user
function errors, so this is easy to fix, but I hope the problem would be
fixed internally in gsl_odeiv_evolve_apply() to prevent further headaches.

Thanks and best regards,
Taneli Kalvas

--
Taneli Kalvas
M.Sc., Researcher
Physics Department, room FL114
P.O. Box 35 (YFL)
40014 University of Jyväskylä, Finland
Phone: +358-44-314-1602
Fax: +358-14-260-2351
Email: -unavailable-

-Deleted Account- <bjg>In charge of this item.

 

Attached Files
file #17977:  ode.c added by bjg (1KiB - text/x-csrc)
file #17978:  out.dat added by bjg (1KiB - application/octet-stream)

 

Depends on the following items: None found

Items that depend on this one: None found

 

CC list is empty

 

Please enter the title of George Orwell's famous dystopian book (it's a date):

 

 

Follow 4 latest changes.

Date Changed By Updated Field Previous Value => Replaced By
Tue 21 Apr 2009 01:36:05 PM UTCbjgStatusConfirmed=>Fixed
  Open/ClosedOpen=>Closed
Mon 20 Apr 2009 07:29:49 PM UTCbjgAttached File-=>Added out.dat, #17978
Mon 20 Apr 2009 07:29:48 PM UTCbjgAttached File-=>Added ode.c, #17977

Back to the top


Powered by Savane 3.1-cleanup1