This is to ensure that the tests give the same results on GPUs and on
CPUs. If one changes the Poisson solver to fft on GPUs, but uses isf on
CPUs, the differences are usually larger than the tolerances in the
tests.