From 6652d6743afedbf62b3bdec4321441b958bd1db2 Mon Sep 17 00:00:00 2001 From: David Wells Date: Thu, 19 May 2022 09:53:29 -0400 Subject: [PATCH] Don't set ghost values in diagonal matrices. We never read them. --- doc/news/changes/incompatibilities/20210817DavidWells | 4 ++++ include/deal.II/matrix_free/operators.h | 6 ++---- 2 files changed, 6 insertions(+), 4 deletions(-) create mode 100644 doc/news/changes/incompatibilities/20210817DavidWells diff --git a/doc/news/changes/incompatibilities/20210817DavidWells b/doc/news/changes/incompatibilities/20210817DavidWells new file mode 100644 index 0000000000..aae77d62b1 --- /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. +
+(David Wells, 2021/08/17) diff --git a/include/deal.II/matrix_free/operators.h b/include/deal.II/matrix_free/operators.h index f0843f42a6..44f8bfaede 100644 --- 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 } -- 2.39.5