py-maxwell-fdfd
Solving Maxwell’s equations via A python implementation of the 3D curl-curl E-field equations. This code contains additional work to engineer the eignspectrum for better convergence with iterative solvers (using the Beltrami-Laplace 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 column-major (or Fortan) ordering whereas numpy is natively row-major 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
