maxwelldg.hh uses incorrect discretization
The parameter class for maxwelldg.hh provides position-dependent eps and mu, but the discretization in maxwelldg.hh cannot deal with that. The problem is that we transform the unknowns into characteristic variables, which are the eigenvectors of some matrix A on each intersection. However, this A depends on eps and mu, and thus if eps and mu differ on the inside and the outside of the intersection the characteristic bases on the inside and the outside differ. The code will run and will seemingly produce results, but there will be spurious reflections.
See my Diss for a better discretization.