[SciPy-User] weired results with ode solver

Rob Clewley rob.clewley@gmail....
Sat Feb 6 13:30:49 CST 2010


Hi,

> The behavior of array int devision is a bit funny, but after you get
> used to it, it's fine.

You can always use

from __future__ import division

at the top of all your files and you'll be safe to do this division
(at an almost negligible performance hit). The mathematically
intuitive version of division is standard in Python 3000.

>>>> timerange=arange(1,12)
>>>> timerange/2
> array([0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5]) #this is why i got multiple values
>>>> timerange/2.0
> array([ 0.5,  1. ,  1.5,  2. ,  2.5,  3. ,  3.5,  4. ,  4.5,  5. ,  5.5])

You may wish to check out numpy's linspace command for generating
float arrays with arbitrary step sizes in a numerically safe manner
(and with endpoints included by default).

-Rob


More information about the SciPy-User mailing list