[SciPy-User] Trajectory Integration via scipy.integrate or PyDSTool?
Mon May 14 12:22:37 CDT 2012
I think you should use a tool that's appropriate for the task. With an
hour of work and help on this list or the PyDSTool forum at
should be able to get your problem working quickly with PyDSTool. I
expect your problem is easy to solve in PyDSTool, and difficult to do
You can still use scipy's VODE in PyDSTool and you won't need any
fancy installation. Check out the file vode_event_test1.py in the
tests directory. Writing the vector field is simple, but in case you
have specific problems setting up yours then you should send me your
code on list so that I can help. I'd be interested to see what the
equations look like.
It's very simple to define a terminal event, just copy and adapt the
code in the file. You can write a loop so that whenever integration
has stopped because of an event, you reverse the velocity IC and
On Mon, May 14, 2012 at 12:59 PM, Sara Gallian <firstname.lastname@example.org> wrote:
> Hi all,
> I'm trying to integrate particle trajectories in a finite domain, and I need to make them "reflect" at the boundaries (i.e. reverse the velocity perpendicular to the plane they collided with). I started by using scipy.integrate.odeint , but since the integration steps are variable, simply checking the position and reversing the velocity won't work!
> Can anyone suggest the quickest way to obtain this? would Vode be able to handle this, or should I try to learn to use PyDSTool?
> I'm running pretty late on a deadline, so any suggestion is more than appreciated :)
> SciPy-User mailing list
Robert Clewley, Ph.D.
Neuroscience Institute and
Department of Mathematics and Statistics
Georgia State University
PO Box 5030
Atlanta, GA 30302, USA
tel: 404-413-6420 fax: 404-413-5446
More information about the SciPy-User