From: Wolfgang Bangerth Date: Thu, 29 Mar 2018 22:25:20 +0000 (-0600) Subject: Add a changelog entry. X-Git-Tag: v9.0.0-rc1~200^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=b6d3b26c8d7599689b430950749699d2e3ab2d0a;p=dealii.git Add a changelog entry. --- diff --git a/doc/news/changes/minor/20180329Bangerth b/doc/news/changes/minor/20180329Bangerth new file mode 100644 index 0000000000..776fd97740 --- /dev/null +++ b/doc/news/changes/minor/20180329Bangerth @@ -0,0 +1,25 @@ +Fixed: The VectorTools::integrate_difference() function allows users +to provide a weight function that can also serve as a component mask +to select individual components of the solution vector for error +computation. For components not selected, such a mask would then +simply be zero. +
+In some cases, the solution vector contains NaN numbers, for example +when one uses the FE_FaceQ element for certain components of the +solution vector and uses a quadrature formula for error evaluation +that has quadrature points in the interior of the cell. For any +"regular" solution component for which the component mask has a zero +weight, the value of that component will be multiplied by zero and +consequently does not add anything to the error computation. However, +if the NaNs of a FE_FaceQ are multiplied with zero weights, the result +is still a NaN, and adding it to the values times weights of the other +components results in NaNs -- in effect rendering it impossible to get +any information out of the VectorTools::integrate_difference() +function if one of the finite elements involved is FE_FaceQ. +
+This is now fixed by simply skipping vector components for which the +weight vector is zero. This has the same result as before for all +"normal" situations, but also properly skips the NaN case outlined +above. +
+(Wolfgang Bangerth, 2018/03/29)