From: Maximilian Bergbauer <bergbauer@lnm.mw.tum.de>
Date: Tue, 18 May 2021 04:31:30 +0000 (+0200)
Subject: FEPointEvaluation::integrate(): zero out solution array also if no points given
X-Git-Tag: v9.3.0-rc1~47^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=fdc084d4396dbfc7002f72091846d5e31316ed92;p=dealii.git

FEPointEvaluation::integrate(): zero out solution array also if no points given
---

diff --git a/include/deal.II/fe/fe_point_evaluation.h b/include/deal.II/fe/fe_point_evaluation.h
index f945b72b10..552e3bf7e1 100644
--- a/include/deal.II/fe/fe_point_evaluation.h
+++ b/include/deal.II/fe/fe_point_evaluation.h
@@ -783,7 +783,10 @@ FEPointEvaluation<n_components, dim, spacedim, Number>::integrate(
   const EvaluationFlags::EvaluationFlags &                    integration_flags)
 {
   if (unit_points.size() == 0) // no evaluation points provided
-    return;
+    {
+      std::fill(solution_values.begin(), solution_values.end(), 0.0);
+      return;
+    }
 
   AssertDimension(solution_values.size(), fe->dofs_per_cell);
   if (false /*TODO*/ && (((integration_flags & EvaluationFlags::values) ||