Partial Differential Equation Toolbox |

pdejmps

Error estimates for adaption

Syntax

-     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. 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 n is the unit normal of edge and the braced term is the jump in flux across the element edge, where and 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`

| pdeintrp |  | pdemdlcv |