From: David Wells Date: Mon, 10 Jan 2022 14:58:59 +0000 (-0500) Subject: Fix a check for zero entries in the lumped mass matrix. X-Git-Tag: v9.4.0-rc1~632^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F13211%2Fhead;p=dealii.git Fix a check for zero entries in the lumped mass matrix. --- diff --git a/include/deal.II/matrix_free/operators.h b/include/deal.II/matrix_free/operators.h index 64f18fae96..9a60978f5b 100644 --- a/include/deal.II/matrix_free/operators.h +++ b/include/deal.II/matrix_free/operators.h @@ -1960,7 +1960,7 @@ namespace MatrixFreeOperators // entries here. for (size_type i = 0; i < locally_owned_size; ++i) { - if (inverse_lumped_diagonal_vector.local_element(i) == Number(0.)) + if (lumped_diagonal_vector.local_element(i) == Number(0.)) inverse_lumped_diagonal_vector.local_element(i) = Number(1.); else inverse_lumped_diagonal_vector.local_element(i) =