[SciPy-User] Field Liner Tracing
Fri Oct 29 07:51:37 CDT 2010
I have a vector field defined on a 3D grid and would like to trace field
What's the best method to do this?
I suppose I could use scipy.integrate.odeint, but then I would have to
define my own dy/dt method that interpolates the values on the grid. It
seems to me that this is a task common enough for having a standard
solution somewhere in scipy already.
Also, I am bothered by the cost of having to define dy/dt as a Python
function. It seems to me that this function is probably called quite a
lot by the integrator and that Python overhead may easily become a
It would be wonderful to have a method like trace_field_lines((Vx, Vy, Vz), X,
Y, Z) where (Vx, Vy, Vz) are 3d numpy arrays containing the components of
the vector field and X, Y, Z 1d grids containing the grid coordinates.
Is there anything like that around already?
»Time flies like an arrow, fruit flies like a Banana.«
PGP fingerprint: 5B93 61F8 4EA2 E279 ABF6 02CF A9AD B7F8 AE4E 425C
More information about the SciPy-User