From a0bdb4c785965a0eef79c80f999dce0896d8a4ce Mon Sep 17 00:00:00 2001 From: Timo Heister Date: Mon, 31 Aug 2015 16:17:42 -0400 Subject: [PATCH] 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. --- include/deal.II/numerics/vector_tools.templates.h | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) 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