Comment by wcrossbow

They might have been using the fast multipole expansion method