[cig-commits] r22615 - short/3D/PyLith/trunk/unittests/libtests/bc

[knepley at geodynamics.org]

Author: knepley
Date: 2013-07-15 06:20:22 -0700 (Mon, 15 Jul 2013)
New Revision: 22615

Modified:
   short/3D/PyLith/trunk/unittests/libtests/bc/TestAbsorbingDampers.cc
   short/3D/PyLith/trunk/unittests/libtests/bc/TestNeumann.cc
Log:
Prevent cascading test failures

Modified: short/3D/PyLith/trunk/unittests/libtests/bc/TestAbsorbingDampers.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/bc/TestAbsorbingDampers.cc	2013-07-15 13:19:57 UTC (rev 22614)
+++ short/3D/PyLith/trunk/unittests/libtests/bc/TestAbsorbingDampers.cc	2013-07-15 13:20:22 UTC (rev 22615)
@@ -139,6 +139,7 @@
   PetscErrorCode err = 0;
   PetscInt dp = 0;
   for(PetscInt c = cStart; c < cEnd; ++c) {
+    PetscInt  vertices[32];
     PetscInt *closure = NULL;
     PetscInt closureSize, numCorners = 0;
 
@@ -146,14 +147,14 @@
     for(PetscInt p = 0; p < closureSize*2; p += 2) {
       const PetscInt point = closure[p];
       if ((point >= vStart) && (point < vEnd)) {
-        closure[numCorners++] = point;
+        vertices[numCorners++] = point;
       }
     }
+    err = DMPlexRestoreTransitiveClosure(subMesh, c, PETSC_TRUE, &closureSize, &closure);PYLITH_CHECK_ERROR(err);
     CPPUNIT_ASSERT_EQUAL(_data->numCorners, numCorners);
     for(PetscInt p = 0; p < numCorners; ++p, ++dp) {
-      CPPUNIT_ASSERT_EQUAL(_data->cells[dp], closure[p]);
+      CPPUNIT_ASSERT_EQUAL(_data->cells[dp], vertices[p]);
     } // for
-    err = DMPlexRestoreTransitiveClosure(subMesh, c, PETSC_TRUE, &closureSize, &closure);PYLITH_CHECK_ERROR(err);
   } // for
 
   // Check damping constants

Modified: short/3D/PyLith/trunk/unittests/libtests/bc/TestNeumann.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/bc/TestNeumann.cc	2013-07-15 13:19:57 UTC (rev 22614)
+++ short/3D/PyLith/trunk/unittests/libtests/bc/TestNeumann.cc	2013-07-15 13:20:22 UTC (rev 22615)
@@ -191,6 +191,7 @@
   PetscInt dp = 0;
   PetscErrorCode err = 0;
   for (PetscInt c = cStart; c < cEnd; ++c) {
+    PetscInt  vertices[32];
     PetscInt *closure = PETSC_NULL;
     PetscInt  closureSize, numCorners = 0;
 
@@ -198,14 +199,14 @@
     for(PetscInt p = 0; p < closureSize*2; p += 2) {
       const PetscInt point = closure[p];
       if ((point >= vStart) && (point < vEnd)) {
-        closure[numCorners++] = point;
+        vertices[numCorners++] = point;
       } // if
     } // for
+    err = DMPlexRestoreTransitiveClosure(subMesh, c, PETSC_TRUE, &closureSize, &closure);PYLITH_CHECK_ERROR(err);
     CPPUNIT_ASSERT_EQUAL(_data->numCorners, numCorners);
     for(PetscInt p = 0; p < numCorners; ++p, ++dp) {
-      CPPUNIT_ASSERT_EQUAL(_data->cells[dp], closure[p]);
+      CPPUNIT_ASSERT_EQUAL(_data->cells[dp], vertices[p]);
     } // for
-    err = DMPlexRestoreTransitiveClosure(subMesh, c, PETSC_TRUE, &closureSize, &closure);PYLITH_CHECK_ERROR(err);
   } // for
 
   // Check traction values