Open ElliotB256 opened 3 years ago
Turns out the tree method has a name, it's called Barnes-Hut https://en.wikipedia.org/wiki/Barnes%E2%80%93Hut_simulation
There is a recent repo for implementing this algorithmn, wouldn't hurt to try it and see how it performs. I expect it may need a bit of profiling https://docs.rs/nbody_barnes_hut/0.1.0/nbody_barnes_hut/
Good overview in Leyuan Wu, Geophys. J. Int.(2018)215,1865–1886
Long-ranged forces are important for 3D MOTs with lots of atoms, as the repulsive forces reduce the phase space density.
Some thoughts:
MacroAtom
approach, where relevant quantities like repulsive forces are scaled so that a small number of atoms can simulate the PSD of an (intractable) larger number of atoms.Possible implementation:
leaf.force
+ force calculated pair-wise between other atoms in the leaf.Building tree will be O(L) O(N).