[cig-commits] [commit] knepley/upgrade-petsc-interface: Added topology consistency checks to adjustTopology() test. (7bbdc39)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/pylith

On branch  : knepley/upgrade-petsc-interface
Link       : https://github.com/geodynamics/pylith/compare/fa51e996b6e488bee3553868d0bcf40e778caa3b...7bbdc391b50f931fe7d30225d5e3d6a3ac8b4e6f

>---------------------------------------------------------------

commit 7bbdc391b50f931fe7d30225d5e3d6a3ac8b4e6f
Author: Brad Aagaard <baagaard at usgs.gov>
Date:   Wed Oct 30 12:11:48 2013 -0700

    Added topology consistency checks to adjustTopology() test.


>---------------------------------------------------------------

7bbdc391b50f931fe7d30225d5e3d6a3ac8b4e6f
 unittests/libtests/faults/TestFaultCohesive.cc | 18 ++++++++++++++++++
 1 file changed, 18 insertions(+)

diff --git a/unittests/libtests/faults/TestFaultCohesive.cc b/unittests/libtests/faults/TestFaultCohesive.cc
index 3466971..b3d3a94 100644
--- a/unittests/libtests/faults/TestFaultCohesive.cc
+++ b/unittests/libtests/faults/TestFaultCohesive.cc
@@ -683,6 +683,15 @@ pylith::faults::TestFaultCohesive::_testAdjustTopology(Fault* fault,
   CPPUNIT_ASSERT_EQUAL(data.cellDim, mesh.dimension());
   dmMesh = mesh.dmMesh();CPPUNIT_ASSERT(dmMesh);
 
+  // Check consistency
+  PetscBool isSimplexMesh = PETSC_TRUE;
+  if ((data.cellDim == 2 && data.numCorners[0] == 4) ||
+      (data.cellDim == 3 && data.numCorners[0] == 8)) {
+    isSimplexMesh = PETSC_FALSE;
+  } // if
+  err = DMPlexCheckSymmetry(dmMesh);CPPUNIT_ASSERT(!err);
+  err = DMPlexCheckSkeleton(dmMesh, isSimplexMesh);CPPUNIT_ASSERT(!err);
+
   // Check vertices
   topology::Stratum verticesStratum(dmMesh, topology::Stratum::DEPTH, 0);
   const PetscInt vStart = verticesStratum.begin();
@@ -827,6 +836,15 @@ pylith::faults::TestFaultCohesive::_testAdjustTopology(Fault* faultA,
   CPPUNIT_ASSERT_EQUAL(data.cellDim, mesh.dimension());
   dmMesh = mesh.dmMesh();CPPUNIT_ASSERT(dmMesh);
 
+  // Check consistency
+  PetscBool isSimplexMesh = PETSC_TRUE;
+  if ((data.cellDim == 2 && data.numCorners[0] == 4) ||
+      (data.cellDim == 3 && data.numCorners[0] == 8)) {
+    isSimplexMesh = PETSC_FALSE;
+  } // if
+  err = DMPlexCheckSymmetry(dmMesh);CPPUNIT_ASSERT(!err);
+  err = DMPlexCheckSkeleton(dmMesh, isSimplexMesh);CPPUNIT_ASSERT(!err);
+
   // Check vertices
   topology::Stratum verticesStratum(dmMesh, topology::Stratum::DEPTH, 0);
   const PetscInt vStart = verticesStratum.begin();