Test prepare_for_use() with constraint matrices

tests/systems/constraint_operator_test.C

@@ -8,6 +8,7 @@
 #include <libmesh/mesh.h>
 #include <libmesh/mesh_generation.h>
 #include <libmesh/mesh_function.h>
+#include <libmesh/mesh_tools.h>
 #include <libmesh/numeric_vector.h>
 #include <libmesh/parallel.h>
 #include <libmesh/quadrature.h>
@@ -178,6 +179,13 @@ public:
 
     mesh.copy_constraint_rows(*matrix);
 
+    // We should be prepared at this point
+    CPPUNIT_ASSERT(mesh.is_prepared());
+    CPPUNIT_ASSERT(MeshTools::valid_is_prepared(mesh));
+
+    // Do a redundant prepare-for-use to catch any issues there
+    mesh.prepare_for_use();
+
     es.init ();
     sys.solve();
 
@@ -244,6 +252,8 @@ public:
 
     mesh.read(meshfile);
 
+    CPPUNIT_ASSERT(mesh.is_prepared());
+
     // For these matrices Coreform has been experimenting with PETSc
     // solvers which take the transpose of what we expect, so we'll
     // un-transpose here.
@@ -253,6 +263,13 @@ public:
 
     mesh.copy_constraint_rows(*matrix);
 
+    // We should be prepared at this point
+    CPPUNIT_ASSERT(mesh.is_prepared());
+    CPPUNIT_ASSERT(MeshTools::valid_is_prepared(mesh));
+
+    // Do a redundant prepare-for-use to catch any issues there
+    mesh.prepare_for_use();
+
     es.init ();
     sys.solve();
 
@@ -316,6 +333,13 @@ public:
 
     mesh.copy_constraint_rows(*matrix);
 
+    // We should be prepared at this point
+    CPPUNIT_ASSERT(mesh.is_prepared());
+    CPPUNIT_ASSERT(MeshTools::valid_is_prepared(mesh));
+
+    // Do a redundant prepare-for-use to catch any issues there
+    mesh.prepare_for_use();
+
     es.init ();
     sys.solve();