A parallelized and high-performance C++ implementation of the Barnes-Hut algorithm for simulating an N-body system and visualizing it.

CUDA GPU acceleration, interactive simulation viewer, lossy simulation data compression, and optimized support for other compilers (non-GCC) coming soon. See Roadmap

The above is a fast-forwarded and reduced resolution render from an older project.

Full-size 3D demo coming soon. Simulation in progress on Google Compute Engine.

1. Checkout the repo from github, and build the project with cmake

$ git clone https://github.com/oliverzh2000/parallel-barnes-hut.git
$ cd parallel-barnes-hut
$ cmake ./ DCMAKE_BUILD_TYPE=RelWithDebinfo
$ make
$ ./nbody --help

See Command Line Options for detailed usage information.

2. Initialize the simulation with between ~10^7 and ~10^8 stars with random positions and velocities and then leave the system running for a few days, weeks, or even months. See Recommended Parameter Values for some general tips.

• Important! Make sure you have enough disk space to store the simulation results. (Coming soon) You can use the --estimate command line option to get an estimate of the disk space needed to store the simulation results at your specified level of detail as well as an estimate of the time needed to run the simulation.

3. (Coming soon) Sign up for an email notification when your simulation completes.
4. (Coming soon) Come back when finished and view your spectacular results on the simulation viewer.

Run the simulation driver (the nbody executable) with --help for more information.

Use --sim-directory to give the simulation driver a path to a directory. This directory must contain an init.txt that is used to specify all the parameters of the simulation.

integratorType: 'IntegratorLeapfrog' | 'IntegratorEuler'
    timestep: positive-decimal

forceCalcType: 'ForceCalcAllPairs' | 'ForceCalcBarnesHut' | 'ForceCalcBarnesHutParallel'
    gravConst: positive-decimal
    softening: positive-decimal
    theta: positive-decimal-value  *1.

starsInitMode: 'readStarsInline' | 'readStarsFromLatestBinaryFrame' | 'readStarsFromLatestShortFrame' |createSpiralGalaxy'
    n: number-of-stars
    pos-x pos-y pos-z vel-x vel-y vel-z mass
    ... 
    ... n rows total, one per star  *2.
    
    n: number-of-stars
    centerPos:
        x: decimal
        y: decimal
        z: decimal
    centerVel:
        x: decimal
        y: decimal
        z: decimal
    radialStdDev: positive-decimal
    avgMass: positive-decimal
    avgMassStdDev: positive-decimal
    seed: integer  *3.

Notes:

1. Use this field only with a Barnes-Hut force calculator type.
2. Use this section only with readStarsInline star initialization mode.
3. Use this section only with createSpiralGalaxy star initialization mode.

All units are specified in SI: meters, seconds, kilograms, radian, etc.

integratorType: IntegratorLeapfrog
    timestep: 0.01

forceCalcType: ForceCalcBarnesHutParallel
    gravConst: 1.0
    softening: 0.01
    theta: 0.5

starsInitMode: readStarsFromLatestBinaryFrame

integratorType: IntegratorLeapfrog
    timestep: 0.01

forceCalcType: ForceCalcBarnesHutParallel
    gravConst: 1.0
    softening: 0.01
    theta: 0.5

starsInitMode: createSpiralGalaxy
    n: 100000
    centerPos:
        x: 0
        y: 0
        z: 0
    centerVel:
        x: 0
        y: 0
        z: 0
    radialStdDev: 100
    avgMass: 1
    avgMassStdDev: 0.1
    seed: 0

integratorType: IntegratorLeapfrog
    timestep: 0.01

forceCalcType: ForceCalcBarnesHutParallel
    gravConst: 1.0
    softening: 0.01
    theta: 0.5

starsInitMode: readStarsInline
    n: 5
    0 8 0 0 0 0 1
    5 7 0 0 0 0 1
    5 5 0 0 0 0 1
    7 5 0 0 0 0 1
    8 0 0 0 0 0 1