pdejmps (Partial Differential Equations Toolbox)

Error estimates for adaption.

errf=pdejmps(p,t,c,a,f,u,alfa,beta,m)

Description

errf=pdejmps(p,t,c,a,f,u,alfa,beta,m)

calculates the error indication function used for adaption. The columns of errf correspond to triangles, and the rows correspond to the different equations in the PDE system.

p and t

are mesh data. See the entry on initmesh for details.

c, a, and f

are PDE coefficients. See the entry on assempde for details. c, a, and f must be expanded, so that columns correspond to triangles.

is the solution vector. See the entry on assempde for details.

The formula for computing the error indicator E(K) for each triangle K is

where n

is the unit normal of edge

and the braced term is the jump in flux across the element edge. The norm is an L₂ norm computed over the element K. The error indicator is stored in errf as column vectors, one for each triangle in t. More information can be found in the section "Adaptive Mesh Refinement" in the manual.