A FreeFEM++ code uses the Modified Weak Galerkin Finite Element Method to solve the Poisson Problem

Hello everyone,
I wrote a FreeFEM++ code to solve the Poisson Problem -∆u=f in Ω u=g on ∂Ω using the Modified Weak Galerkin Finite Element Method based on the paper:

[1] “A modified weak Galerkin finite element method”. X. Wang, N.S. Malluwawadu, F. Gao, T.C. McMillan.

https://www.sciencedirect.com/science/article/pii/S0377042714002131 3.

I modified a FreeFEM++ code for the Discontinuous Galerkin FEM to solve this problem to fit the Modified Weak Galerkin FEM, but the results are not the same as in the paper. I am wondering what’s wrong with my code.
MWGFEM.edp (1.4 KB)

1 Like