[SciPy-user] odeint with digital data

[Rob Clewley]

Hi Ryan,

To my knowledge you cannot do this with odeint unless you change a
constant value on the RHS (i.e. technically change your system) to
reflect the changing input value after every time-step, but it is very
easy to do in PyDSTool. There you can also force integration to be
only at the discrete mesh points of your input signal (and would
linearly interpolate in-between otherwise). I am certainly happy to
help you set up your script to do this if you wish to try it in
PyDSTool.

HTH!
Rob

On 04/06/07, Ryan Krauss <ryanlists at gmail.com> wrote:
> I have a vector of experimental data that I need to use as part of a
> system of ode's.  I would like to solve this system using
> integrate.odeint.  Can odeint be forced to solve only at the discrete
> points in time where the experimentally measured signal is available?
> Or do I need to set up some interpolation function to find the
> measured signal at any time?  If that signal is the only thing that
> explictly depends on time and I set up a digital look-up table that
> returns the same constant value for the range from t to t+dt, would I
> effectively force odeint to do what I want?  Am I making any sense?
> Is there a better way?
>
> Thanks,
>
> Ryan
> _______________________________________________
> SciPy-user mailing list
> SciPy-user at scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-user
>