“…WG-FEM has the flexibility of dealing with discontinuous elements while sharing the same weak formulation with the classical conforming finite element methods. Following its first development in [18,19] for second order elliptic equations, the method has gained a lot attention and popularity and has been applied to various PDEs including the Stokes equation, biharmonic equation, elasticity equation, and Maxwell's equations, see [16,21,11] and the references therein for more details.…”