From: Timo Heister Date: Mon, 31 Aug 2015 20:17:42 +0000 (-0400) Subject: do not use square of weight for Hdiv_seminorm X-Git-Tag: v8.4.0-rc2~504^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F1491%2Fhead;p=dealii.git do not use square of weight for Hdiv_seminorm The interpretation of weight in integrate_difference is that it scales the norm not the function being integrated. We therefore need to scale the components of div u with the square root of the weight. --- diff --git a/include/deal.II/numerics/vector_tools.templates.h b/include/deal.II/numerics/vector_tools.templates.h index 454c44881c..7737b179ba 100644 --- a/include/deal.II/numerics/vector_tools.templates.h +++ b/include/deal.II/numerics/vector_tools.templates.h @@ -6326,10 +6326,9 @@ namespace VectorTools "with at least 'dim' components. In that case, this function " "will take the divergence of the first 'dim' components.")); double sum = 0; - // take the trace of the derivatives, square it, multiply it - // with the weight function + // take the trace of the derivatives scaled by the weight and square it for (unsigned int k=0; k