Add new function compute_potential in PoissonFFT.