From: Timo Heister <timo.heister@gmail.com>
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<dim; ++k)
-                sum += (data.psi_grads[q][k][k] * data.psi_grads[q][k][k]) *
-                       data.weight_vectors[q](k);
-              diff += sum * fe_values.JxW(q);
+                sum += data.psi_grads[q][k][k] * data.weight_vectors[q](k);
+              diff += sum * sum * fe_values.JxW(q);
             }
           diff = std::sqrt(diff);
           break;