From: Maximilian Bergbauer 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::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) ||