Solvers for Maxwell's electrodynamic equations with time-space varying coefficients based on PyClaw