A novel numerical method is presented for solving the anisotropic diffusion equation in magnetic fields confined to a periodic box. The method is accurate and provably stable.
