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
to triangles, and the rows correspond to the different equations in
the PDE system.
mesh data. For details, see
f are PDE coefficients.
f must be expanded, so that columns correspond
The formula for computing the error indicator E(K) for each triangle K is
where 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
errf as column vectors, one for each triangle
t. More information can be found in the section Adaptive Mesh Refinement.