“…Numerical implementation of WGFEM was discussed in [20]. This method successfully applied to some PDEs such as parabolic equation [21,22], Helmholtz equation [23,24], biharmonic equation [25][26][27][28], Brinkman equation [29,30], elliptic interface problem [31], linear elasticity problems [32,33], Oseen equation [34], Darcy equation [35,36], Darcy-stokes equation [37][38][39], Maxwell equation [40,41], and so on. An important feature of WG methods is applying totally discontinuous piecewise polynomials on the finite element partition.…”