[cig-commits] r19174 - short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/faults

[brad at geodynamics.org]

Author: brad
Date: 2011-11-10 09:53:28 -0800 (Thu, 10 Nov 2011)
New Revision: 19174

Modified:
   short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/faults/CohesiveTopology.cc
Log:
Added some asserts.

Modified: short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/faults/CohesiveTopology.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/faults/CohesiveTopology.cc	2011-11-10 17:53:11 UTC (rev 19173)
+++ short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/faults/CohesiveTopology.cc	2011-11-10 17:53:28 UTC (rev 19174)
@@ -661,11 +661,16 @@
   for (SieveSubMesh::label_sequence::iterator v_iter = fVertices2Begin;
        v_iter != fVertices2EndNew;
        ++v_iter) {
+    assert(coordinates->getFiberDimension(*v_iter) ==
+	   coordinates->getFiberDimension(vertexRenumber[*v_iter]));
     coordinates->updatePoint(vertexRenumber[*v_iter], 
 			     coordinates->restrictPoint(*v_iter));
-    if (constraintCell)
+    if (constraintCell) {
+      assert(coordinates->getFiberDimension(*v_iter) ==
+	     coordinates->getFiberDimension(vertexLagrangeRenumber[*v_iter]));
       coordinates->updatePoint(vertexLagrangeRenumber[*v_iter],
 			       coordinates->restrictPoint(*v_iter));
+    } // if
   } // for
   if (debug)
     coordinates->view("Coordinates with shadow vertices");