Generic nbody ODE solver

So, I wrote a generic nbody ODE solver, in pure python. It’s likely quite slow, but it works – seemed like too much work at the time to interface to scipy’s ODE solvers, which seemed designed for smaller, less coupled systems.

I may have misjudged, though which is why I’m noting this link http://www.scipy.org/Developer_Zone/Recipes/ODE_Integration_on_Banded_Systems ,which has an example of what seems to be a large dimensional ODE system. Thing is, we still need to do things like update adjacency lists and distances, and recompute summation densities between steps. Worth a look, though.

Advertisements

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s

%d bloggers like this: