GravityLab: Astrophysical N-body simulator.

A fast astrophysical N-body simulator.

GravityLab is written in Python and Cython, with performance-critical parts implemented in Cython and the main simulation module written in Python. This makes GravityLab flexible and very fast. The core physics solvers are written in Cython and are optimized to make GravityLab fast. The ODEs governing motions of the bodies are solved using the Leapfrog integration scheme. The Python module is flexible and powerful, while keeping it simple for users to use it.

Installation

Dependencies

The following are required to run GravityLab:

• Cython
• NumPy
• Matplotlib

Install

To install this code, run the following commands in your terminal:

> git clone https://github.com/Kushaalkumar-pothula/gravitylab.git
> cd gravitylab
> python setup.py build_ext --inplace
> python setup.py install

Usage

You can launch a Python prompt or IPython shell from the gravitylab/ directory, which allows you to import the modules directly. You should also all the gravitylab to PYTHONPATH to be able to access it from anywhere throughout the system. To run gravitylab in a Python prompt or IPython shell, you can do:

> cd gravitylab
> python

or if you use IPython (recommended):

> cd gravitylab
> ipython

You can run the simulation using the gravitylab.simulate module. This module provides a powerful (yet simple to use) class Simulation. An instance of the Simulation class can be instantiated with the following command (parameters not shown here):

from gravitylab.simulate import Simulation
my_simulation = Simulation()

Or if you are inside the gravitylab directory:

from simulate import Simulation
my_simulation = Simulation()

Set up simulations

Add initial conditions to a simulation instance:

my_simulation.initial_conditions()

And then run a simulation using the .run() method:

my_simulation.run()

You can also get the calculated positions returned by the above command:

positions = my_simulation.run()

Next, you can create a plot of the evolution positions throughout the simulation:

my_simulation.plot()

The .plot() method takes a many parameters, including 3D plotting.

An preset example (two equally massive bodies with random positions ans velocities) for running GravityLab simulations can be found in gravitylab/example.py:

> python gravitylab/example.py

Viewing help

In [2]: Simulation?
Init signature: Simulation(N, dt, t, tEnd)
Docstring:
Main simulation class for GravityLab

Attributes
----------
N    : Number of bodies
dt   : Timestep
t    : Start time
tEnd : End time

Methods
-------
initial_conditions : Initial conditions for simulation
run                : Run simulation
plot               : Plot positions
Init docstring:
Initialize a simulation

Parameters
----------
N    : Number of bodies
dt   : Timestep
t    : Start time
tEnd : End time

>>> help(Simulation)
Help on class Simulation in module simulate:

class Simulation(builtins.object)
|  Simulation(N, dt, t, tEnd)
|
|  Main simulation class for GravityLab
|
|  Attributes
|  ----------
|  N    : Number of bodies
|  dt   : Timestep
|  t    : Start time
|  tEnd : End time
|
|  Methods
|  -------
|  initial_conditions : Initial conditions for simulation
|  run                : Run simulation
|  plot               : Plot positions
|
|  Methods defined here:
|
|  __init__(self, N, dt, t, tEnd)
|      Initialize a simulation
|
|      Parameters
|      ----------
|      N    : Number of bodies
|      dt   : Timestep
|      t    : Start time
|      tEnd : End time
|
|  animate(self)
|      Animate
|
|  animate_func(self, i)
|      Animate positions
|
|  init_animation(self)
|      Initialize animation
|
|  initial_conditions(self, pos, vel, mass)
|      Initial conditions for a simulation
|
|      Parameters
|      ----------
|      pos : N x 3 array of positions
|      vel : N x 3 array of velocities
|      mass : N x 1 array of masses
|
|  plot(self, three_dimensional=False, start_pos=True, color='blue', alpha=0.6)
|      Plot positions of bodies
|
|      Parameters
|      ----------
|      three_dimensional (optional) : 3D plotting
|      color (optional)             : Desired colour of markers
|      alpha (optional)             : Opacity of markers
|
|  run(self, verbose=False)
|      Run the simulation
|
|      Parameters
|      ----------
|      verbose (optional) : Give more verbose output
|
|      Returns
|      -------
|      pos_arr : N x 3 array of positions
|
|  ----------------------------------------------------------------------

The same step is applicable for methods of Simulation().

Author

This code was written by Kushaal Kumar Pothula. I wrote this code during my 10th grade while learning about numerical simulations during the summer as a learning project.