High-performance n-body solver achieving 128× speedup via checkerboard domain decomposition, toroidal neighbor communication, and sweep-and-prune for scalable uniform/gaussian workloads. - View it on GitHub