DDX_FMM

Use the fast multipole method to accelerate the solver