pymaxwellfdfd
Solving Maxwell’s equations via A python implementation of the 3D curlcurl Efield equations. This code contains additional work to engineer the eignspectrum for better convergence with iterative solvers (using the BeltramiLaplace operator). You can control this in the main function through the input parameter $s = {0,1,1}$
There is also a preconditioners to render the system matrix symmetric.
important notes about implementation
 Note that arrays are ordered using columnmajor (or Fortan) ordering whereas numpy is natively rowmajor or C ordering. You will see this in operations like reshape where I specify ordering.
Examples

Plane Wave

Dipole in Vacuum
Future
Expect integration of this with ceviche for autograd