Finally working version of laplace_driver with new fmm.
the "p_j" in the potential returned by the new fmm subroutine is just the
density. That was the mistake.