Dear all experts, Respected Francois Sir, and Respected Hecht sir, I have tried to implement modified weak Galerkin in the heat equation 2D FF++ code using the poisson equation modified weak Galerkin code . I am wondering why they both not giving good results sir. Can you please check my code sir??.
Heat equation in 2D with u=0 in boundary. Code> Heat_WGM_Sine.edp (2.4 KB)
Heat equation in 2D with neumann derivative of u=0 in boundary. Code> MWG_Heat_2D_Cosine.edp (2.3 KB)
You miss exp(t) in the last term of definition of f, and *(nTonEdge-1) in the two intalledges involving the mean Heat_WGM_Sine.edp (2.4 KB)
For Neumann BC, do not put the factor (nTonEdge-1) in the two intalledges involving the mean, but put it in the intalledges involving the jumps.
Thank you so much sir. If you get time kindly look the post of integro-partial differential equation code in FF++. I have stored the solution of previous steps and computed the approximate for integral but results not coming good sir.