Dear all experts, I’m a beginner in Freefem. I tried to write the Freefem code for the eq.(28) in the paper " Mixed Discontinuous Galerkin methods for Darcy flow. F. Brezzi, T.J.R. Hughes, L.D. Marini and A. Masud" (2005):
which is,
I did the code but the result is not good.In the paper, the jump for vector valued function is defined as follows:
So I am not sure about the macro jump definition for the vector valued function I wrote in my freefem code is correct based on the jump definition.
macro jumpV(vec1, vec2) (jump(vec1)*N.x + jump(vec2)*N.y) //
I tried to fix but it does’t seem to work. Looking forwards to any ideas.
Here is my code for the eq.(28)
DarcyFlow.edp (5.7 KB)
Thanks in advance.