We use a hyperbolic divergence cleaning approach to simulate the fully time dependent Maxwell equations in a
purely hyperbolic form (PHM):