From: Wolfgang Bangerth <bangerth@colostate.edu>
Date: Thu, 5 Oct 2023 21:01:30 +0000 (-0600)
Subject: Add another assertion.
X-Git-Tag: relicensing~428^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=d22bfd0f9b5ac6ceea1012f67da907d7d70b40a7;p=dealii.git

Add another assertion.
---

diff --git a/include/deal.II/lac/affine_constraints.templates.h b/include/deal.II/lac/affine_constraints.templates.h
index c00d803f46..00a912d874 100644
--- a/include/deal.II/lac/affine_constraints.templates.h
+++ b/include/deal.II/lac/affine_constraints.templates.h
@@ -2671,9 +2671,23 @@ AffineConstraints<number>::distribute(VectorType &vec) const
             // pass as arguments...
 #ifdef DEBUG
           if (needed_elements_for_distribute != IndexSet())
-            for (const auto i : vec_owned_elements)
-              Assert(needed_elements_for_distribute.is_element(i),
-                     ExcInternalError());
+            {
+              Assert(vec_owned_elements.size() ==
+                       needed_elements_for_distribute.size(),
+                     ExcMessage("You have previously initialized this "
+                                "AffineConstraints object with an index set "
+                                "that stated that vectors have size " +
+                                std::to_string(locally_owned_dofs.size()) +
+                                " entries, but you are now calling "
+                                "AffineConstraints::distribute() with a vector "
+                                "of size " +
+                                std::to_string(vec_owned_elements.size()) +
+                                "."));
+
+              for (const auto i : vec_owned_elements)
+                Assert(needed_elements_for_distribute.is_element(i),
+                       ExcInternalError());
+            }
 #endif
 
           VectorType ghosted_vector;