MPI-distributed bodies

[blackwer]

Body calculations are largely relegated to rank 0, and then the solutions are distributed to all the ranks currently. This works really well for a few bodies, but will likely bottleneck computations with large numbers of bodies.