From: David Wells <drwells@email.unc.edu>
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=0ee2147070cedc457a32e4779b48256f08b881fb;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) =