<br><br><div class="gmail_quote">On Wed, Aug 15, 2012 at 10:42 AM, Fabrice Silva <span dir="ltr">&lt;<a href="mailto:silva@lma.cnrs-mrs.fr" target="_blank">silva@lma.cnrs-mrs.fr</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Le mardi 14 août 2012 à 21:21 +0200, Ralf Gommers a écrit :<br>
<div class="im">&gt; On Sun, Aug 12, 2012, Fabrice Silva &lt;<a href="mailto:silva@lma.cnrs-mrs.fr">silva@lma.cnrs-mrs.fr</a>&gt; wrote:<br>
&gt;         I made a pull request [1] to integrate the LSODA solver that<br>
&gt;         is used in odeint into the modular scipy.integrate.ode generic<br>
&gt;         class. In a similar way as for vode, it just wraps the already<br>
&gt;         present lsoda.f file (see .pyf file) and exposes it within an<br>
&gt;         IntegratorBase subclass adjusting the coefficients before<br>
&gt;         calling lsoda.<br>
&gt;<br>
&gt; Does that mean that odeint can be made a wrapper around lsoda and that<br>
&gt; the odepack static extension can be completely removed?<br>
<br>
</div>Hi Ralf,<br>
The pull request allows to run the integration using the object-oriented<br>
interface ode, with the same solver than the odeint interface uses, i.e.<br>
lsoda, extending the integrators available for the object-oriented<br>
interface.<br>
<br>
As I understand the scipy.integrate architecture, we are by now<br>
building:<br>
* the odepack library, which has all the fortran sources required by<br>
lsoda and vode at least.<br>
* the _odepack extension, which defines the _odepack module needed by<br>
odeint.<br>
<br>
This latter would be removable, and odeint a wrapper around the lsoda<br>
pyf&#39;ed function. I suppose you are talking about the _odepack extension,<br>
am I wrong?<br></blockquote><div><br>I was mixing it up a bit, but yes: the _odepack extension and the C source for it. Not necessary to do that at once I guess, but wrapping the same function twice is once too many.<br>
</div></div><br>And forgot in my first email: nice PR, looks good to me.<br><br>Ralf<br><br>