From b6dd4f1ebbf2beb5092a4d032041633703f66cf9 Mon Sep 17 00:00:00 2001 From: Timo Heister Date: Mon, 31 Aug 2015 13:59:05 -0400 Subject: [PATCH] fix incorrect Hdiv computation || u ||_Hdiv = || div u ||_L2 = sqrt( \int_\Omega (div u)^2 ) = sqrt( \int_\Omega (sum d/dx_i u_i)^2 ) and not = sqrt( \int_\Omega (sum (d/dx_i u_i)^2) ) which was the old implementation. --- 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 8b6d4c272b..454c44881c 100644 --- a/include/deal.II/numerics/vector_tools.templates.h +++ b/include/deal.II/numerics/vector_tools.templates.h @@ -6329,9 +6329,8 @@ namespace VectorTools // take the trace of the derivatives, square it, multiply it // with the weight function for (unsigned int k=0; k