[SciPy-User] Field Liner Tracing

[Nikolaus Rath]

Hello,

I have a vector field defined on a 3D grid and would like to trace field
lines.

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
bottleneck...

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?


Thanks,

   -Nikolaus

-- 
 »Time flies like an arrow, fruit flies like a Banana.«

  PGP fingerprint: 5B93 61F8 4EA2 E279 ABF6  02CF A9AD B7F8 AE4E 425C