From 0ee2147070cedc457a32e4779b48256f08b881fb Mon Sep 17 00:00:00 2001 From: David Wells Date: Mon, 10 Jan 2022 09:58:59 -0500 Subject: [PATCH] Fix a check for zero entries in the lumped mass matrix. --- include/deal.II/matrix_free/operators.h | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) 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) = -- 2.39.5