Softened 2D gravity solved each frame with an adaptive-quadtree Fast Multipole Method using Cartesian Taylor expansions.