Don't set ghost values in diagonal matrices.

We never read them.

diff --git a/doc/news/changes/incompatibilities/20210817DavidWells b/doc/news/changes/incompatibilities/20210817DavidWells
new file mode 100644 (file)
index 0000000..aae77d6
--- /dev/null
+++ b/doc/news/changes/incompatibilities/20210817DavidWells
@@ -0,0 +1,4 @@
+Changed: The DiagonalMatrix objects set up by MatrixFreeOperators::Base and
+its inheriting classes no set ghost entries.
+<br>
+(David Wells, 2021/08/17)

diff --git a/include/deal.II/matrix_free/operators.h b/include/deal.II/matrix_free/operators.h
index f0843f42a6603f553536787f617c2ef762db83ec..44f8bfaede591febefe5ccfb89126ba6407c9da6 100644 (file)
--- a/include/deal.II/matrix_free/operators.h
+++ b/include/deal.II/matrix_free/operators.h
@@ -1909,8 +1909,7 @@ namespace MatrixFreeOperators
           1. / inverse_diagonal_vector.local_element(i);
       }
 
-    inverse_diagonal_vector.update_ghost_values();
-    diagonal_vector.update_ghost_values();
+    // We never need ghost values so don't update them
   }
 
 
@@ -2206,8 +2205,7 @@ namespace MatrixFreeOperators
       else
         inverse_diagonal_vector.local_element(i) = 1.;
 
-    inverse_diagonal_vector.update_ghost_values();
-    diagonal_vector.update_ghost_values();
+    // We never need ghost values so don't update them
   }