They might have been using the fast multipole expansion method