Dynamic Sparsity Pattern must be initialized with locally_relevant_rows. Without...

diff --git a/source/lac/sparsity_tools.cc b/source/lac/sparsity_tools.cc
index 348358f4ff4620ec6067f9d9383c7bc42a566e19..a57bada4a10ba7b1e6d229c7ff12de3dbdde6154 100644 (file)
--- a/source/lac/sparsity_tools.cc
+++ b/source/lac/sparsity_tools.cc
@@ -1022,6 +1022,13 @@ namespace SparsityTools
                               const MPI_Comm          mpi_comm,
                               const IndexSet         &locally_relevant_rows)
   {
+    IndexSet rows_dsp = dsp.row_index_set();
+    rows_dsp.subtract_set(locally_relevant_rows);
+    AssertThrow(
+      rows_dsp.n_elements() == 0,
+      ExcMessage(
+        "Dynamic Sparsity Pattern must be initialized with locally_relevant_rows."));
+
     IndexSet requested_rows(locally_relevant_rows);
     requested_rows.subtract_set(locally_owned_rows);