This directory contains code for operating the Maxwell operator on a vector field:
MaxwellOp(V) = curl(1/eps * curl(V))
We store V in a Fourier basis, so that the curl operation is diagonal. To multiply by epsilon-inverse, we use a pair of FFTs to transform to and from position space.
Also included are initialization routines for the auxiliary data and storage required by MaxwellOp, as well as subsidiary routines to compute things like the curl and the energy density.