From: Timo Heister Date: Mon, 31 Aug 2015 17:55:47 +0000 (-0400) Subject: remove Exception from Hdiv_seminorm X-Git-Tag: v8.4.0-rc2~504^2~4 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=12e2bcb02a0d2c91c5c509e5843f31c7a562b539;p=dealii.git remove Exception from Hdiv_seminorm also add some comments --- diff --git a/include/deal.II/numerics/vector_tools.templates.h b/include/deal.II/numerics/vector_tools.templates.h index e5821a98d0..8b6d4c272b 100644 --- a/include/deal.II/numerics/vector_tools.templates.h +++ b/include/deal.II/numerics/vector_tools.templates.h @@ -6218,6 +6218,8 @@ namespace VectorTools } double diff = 0; + + // First work on function values: switch (norm) { case mean: @@ -6245,7 +6247,7 @@ namespace VectorTools diff += sum * fe_values.JxW(q); } - // Compute the root only, if no derivative values are added later + // Compute the root only if no derivative values are added later if (!(update_flags & update_gradients)) diff = std::pow(diff, 1./exponent); break; @@ -6275,8 +6277,10 @@ namespace VectorTools break; case H1_seminorm: + case Hdiv_seminorm: case W1p_seminorm: case W1infty_seminorm: + // function values are not used for these norms break; default: @@ -6284,6 +6288,7 @@ namespace VectorTools break; } + // Now compute terms depending on derivatives: switch (norm) { case W1p_seminorm: