  In this paper a new technique aimed to obtain accurate estimates of the error
in energy norm using a moving least squares (MLS) recovery-based procedure is
presented. We explore the capabilities of a recovery technique based on an
enhanced MLS fitting, which directly provides continuous interpolated fields,
to obtain estimates of the error in energy norm as an alternative to the
superconvergent patch recovery (SPR). Boundary equilibrium is enforced using a
nearest point approach that modifies the MLS functional. Lagrange multipliers
are used to impose a nearly exact satisfaction of the internal equilibrium
equation. The numerical results show the high accuracy of the proposed error
estimator.
