From 92426414d23f1f5a1619d904884523a6a5148e8d Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Tue, 28 Apr 1998 13:23:13 +0000 Subject: [PATCH] Doc update. git-svn-id: https://svn.dealii.org/trunk@219 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/include/numerics/error_estimator.h | 8 ++++++-- 1 file changed, 6 insertions(+), 2 deletions(-) diff --git a/deal.II/deal.II/include/numerics/error_estimator.h b/deal.II/deal.II/include/numerics/error_estimator.h index 8fa092328f..42dbd01404 100644 --- a/deal.II/deal.II/include/numerics/error_estimator.h +++ b/deal.II/deal.II/include/numerics/error_estimator.h @@ -37,10 +37,14 @@ class dVector; In principle, the implementation of the error estimation is simple: let $$ \eta_K^2 = - \frac{h}{24} \int_{\partial K} \left[\frac{\partial u_h}{\partial n}\right]^2 do + h \int_{\partial K} \left[\frac{\partial u_h}{\partial n}\right]^2 do $$ be the error estimator for cell $K$. $[\cdot]$ denotes the jump of the - argument at the face. + argument at the face. In the paper of Ainsworth, $h$ is divided by $24$, + but this factor is a bit esoteric, stemming from interpolation estimates + and stability constants which may hold for the Poisson problem, but may + not hold for more general situations. In the implementation, this factor + is dropped for these reasons. To perform the integration, use is made of the #FEFaceValues# class and the integration is performed for each cell, i.e. no use is made of the fact, that -- 2.39.5