[Scipy-tickets] [SciPy] #1434: VODE and ZVODE are not reentrant (and the wrappers don't know this)

SciPy Trac scipy-tickets@scipy....
Mon May 2 15:24:21 CDT 2011


#1434: VODE and ZVODE are not reentrant (and the wrappers don't know this)
-----------------------------+----------------------------------------------
 Reporter:  pv               |       Owner:  somebody
     Type:  defect           |      Status:  new     
 Priority:  highest          |   Milestone:  0.10.0  
Component:  scipy.integrate  |     Version:  0.9.0   
 Keywords:                   |  
-----------------------------+----------------------------------------------
 The VODE and ZVODE integrators in `scipy.integrate` use SAVE statements,
 and are therefore not reentrant. (I didn't check how much of there is
 saved state, and whether there are ways to swap it.)

 Here's an example of how things go wrong because of this:
 http://permalink.gmane.org/gmane.comp.python.scientific.user/28822

 The best possible fix with the current object-oriented interface seems to
 me to be that a new integrator would "invalidate" old ones. This then
 makes the object-oriented interface just a slightly more verbose way to
 spell a procedural call.

 (The "best" fix would be that someone crawls through the FORTRAN code and
 makes it re-entrant. However, redesigning the ODE integrator interfaces in
 `scipy.integrate` from the ground up would be more productive, IMHO.)

-- 
Ticket URL: <http://projects.scipy.org/scipy/ticket/1434>
SciPy <http://www.scipy.org>
SciPy is open-source software for mathematics, science, and engineering.


More information about the Scipy-tickets mailing list