load "Element_P3" load "Element_P4" int k = 1*2*2; mesh Th = square(10*k, 10*k); //plot(Th); func Pk = P4; fespace Vhc ( Th, Pk ); Vhc uc,ucX,ucY,fC,fCh,ucEX,err; //--exact solution, partial derivatives and rhs---- ucEX = exp(x^2/4+y^2/4); ucX = exp(x^2/4+y^2/4)*x/2; ucY = exp(x^2/4+y^2/4)*y/2; …

Incorrect convergence rate for "Element_P4" in 2d

prj December 21, 2021, 8:53pm 2

You need to use an appropriate qforder, see Quadrature formulae or this similar topic Default quadrature rule.

Topic		Replies	Views
Doubt in the variational formulation of LaplaceRT.edp General Discussion	5	248	May 24, 2023
Sorry in optimization element_OP plein 3d General Discussion	0	269	October 27, 2021
Let the exact solution is unknown for the Poisson equation. How to compute the L2 error estimate u_h - u_h/2 General Discussion	4	308	August 14, 2023
Non-Dimensional Equation error General Discussion	0	431	July 4, 2020
Parabolic equation General Discussion	2	108	June 12, 2024