From fdc084d4396dbfc7002f72091846d5e31316ed92 Mon Sep 17 00:00:00 2001 From: Maximilian Bergbauer Date: Tue, 18 May 2021 06:31:30 +0200 Subject: [PATCH] FEPointEvaluation::integrate(): zero out solution array also if no points given --- include/deal.II/fe/fe_point_evaluation.h | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) 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::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) || -- 2.39.5