From: Timo Heister Date: Mon, 31 Aug 2015 17:59:05 +0000 (-0400) Subject: fix incorrect Hdiv computation X-Git-Tag: v8.4.0-rc2~504^2~3 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=b6dd4f1ebbf2beb5092a4d032041633703f66cf9;p=dealii.git 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. --- 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