]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fix a check for zero entries in the lumped mass matrix. 13211/head
authorDavid Wells <drwells@email.unc.edu>
Mon, 10 Jan 2022 14:58:59 +0000 (09:58 -0500)
committerDavid Wells <drwells@email.unc.edu>
Mon, 10 Jan 2022 14:58:59 +0000 (09:58 -0500)
include/deal.II/matrix_free/operators.h

index 64f18fae96a94da7a64a4ec2bdf1ac541949c154..9a60978f5bdbd1577b53e7705d73b36de229510e 100644 (file)
@@ -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) =

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.