Total mass is sometimes zero. As a temporary workaround, i set it to f64::MIN if it is zero, but really it should never be zero.
I need to refactor apply_velocity((f64,f64,f64), bool)) to a simulation_results struct and allow the implementor to do with the simulation what they want.
Docs need to reflect all of this.
Replace my own write data file functions with serde stuff.
More particles = way more gravity...
average of all accel forces instead of adding?
use nalgebra for points and point math, perhaps. or just abstract into my own matrix math/simd thing
apply_velocity((f64,f64,f64), bool))
to asimulation_results
struct and allow the implementor to do with the simulation what they want.