Home > p2plib > errcheck.m

errcheck

PURPOSE

compute a posteriori error estimate

SYNOPSIS

function err=errcheck(p)

DESCRIPTION

 compute a posteriori error estimate

CROSS-REFERENCE INFORMATION

This function calls:

bgradu2f include b grad u into f for error-est and mesh-adapt.

This function is called by:

bifdetec Check for change of sign(det(A)). If yes, locate bif-point and save.
cont main cont. routine for -div(c grad u)+au-f-b grad u=0,
inistep initial step in cont; yields first 2 points on branch and first tangent
meshcheck check if FEM-solution is converged to PDE-solution by refining mesh
pmnewtonloop parallel newtonloops called from pmcont

SOURCE CODE

0001 function err=errcheck(p)
0002 % compute a posteriori error estimate
0003 alfa=0.15;beta=0.15;mexp=1;
0004 [c,a,f,b]=p.f(p,p.u,p.lam); 
0005 if any(b) f=bgradu2f(p,f,b,p.u); end
0006 errv=pdejmps(p.points,p.tria,c,a,f,p.u,alfa,beta,mexp);
0007 err=max(max(errv));