Partial Differential Equation Toolbox Previous page   Next Page
pdejmps

Error estimates for adaption

Syntax

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. For details, see initmesh.

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

u is the solution vector. For details, see assempde.

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

where ntau is the unit normal of edge tau and the braced term is the jump in flux across the element edge, where alpha and beta are weight indices and m is an order parameter. The norm is an L2 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.

See Also

adaptmesh, pdeadgsc, pdeadworst


Previous page  pdeintrp pdemdlcv Next page