[cig-commits] r7502 - short/3D/PyLith/trunk/libsrc/faults

[knepley at geodynamics.org]

Author: knepley
Date: 2007-06-26 09:47:23 -0700 (Tue, 26 Jun 2007)
New Revision: 7502

Modified:
   short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc
Log:
Fixed faults in parallel


Modified: short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc	2007-06-26 16:47:09 UTC (rev 7501)
+++ short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc	2007-06-26 16:47:23 UTC (rev 7502)
@@ -135,14 +135,15 @@
   const ALE::Obj<Mesh::label_sequence>& fFaces = (*fault)->heightStratum(1);
   int faultDepth      = (*fault)->depth()-1; // Depth of fault cells
   int numFaultCorners = 0; // The number of vertices in a fault cell
+  int faultFaceSize   = 0; // The number of vertices in a face between fault cells
   PointArray flippedFaces; // Incorrectly oriented fault cells
 
   if (!(*fault)->commRank()) {
     numFaultCorners = faultSieve->nCone(*fFaces->begin(), faultDepth)->size();
     if (debug) std::cout << "  Fault corners " << numFaultCorners << std::endl;
     assert(numFaultCorners == faceSize);
+    faultFaceSize = _numFaceVertices(*fFaces->begin(), (*fault), faultDepth);
   }
-  int faultFaceSize = _numFaceVertices(*fFaces->begin(), (*fault), faultDepth);
   if (debug) std::cout << "  Fault face size " << faultFaceSize << std::endl;
   for(Mesh::label_sequence::iterator e_iter = fFaces->begin(); e_iter != fFaces->end(); ++e_iter) {
     if (debug) std::cout << "  Checking fault face " << *e_iter << std::endl;