Error estimates for adaptation

`pdejmps`

is not recommended.

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

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

calculates
the error indication function used for adaptation. 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. `c`

, `a`

,
and `f`

must be expanded, so that columns correspond
to triangles.

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

$$E\left(K\right)=\alpha {\Vert {h}^{m}\left(f-au\right)\Vert}_{K}+\beta {\left(\frac{1}{2}{\displaystyle \sum _{\tau \in \partial K}{h}_{\tau}^{2m}{[{n}_{\tau}\cdot \text{\hspace{0.17em}}(c\nabla {u}_{h})]}^{2}}\right)}^{1/2}$$

where $${n}_{\tau}$$ is the unit normal of edge $$\tau $$ 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 *L*_{2} 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.

Was this topic helpful?