 Discrete elliptic barrier

For a monotone discrete elliptic operator, a nonnegative mesh function $\varphi$ with $L_h\varphi \ge d>0$ can bound errors through the discrete maximum principle. If $L_{h}e=\tau$ and the boundary error is zero, the subsolutions $\pm e+(\Vert \tau {\Vert }_{\infty }/d)\varphi$ yield $\Vert e{\Vert }_{\infty }\le (\Vert \tau {\Vert }_{\infty }/d){\max }_{\partial \Omega }\varphi$. On different subregions one may use different lower bounds for $L_h\varphi$ and different truncation-error bounds, taking the maximum of their ratios. Nonnegativity, or an explicit subtraction of the minimum when $L_{h}1=0$, is essential for the unshifted estimate.