[cig-commits] r12817 - short/3D/PyLith/trunk/libsrc/meshio

[knepley at geodynamics.org]

Author: knepley
Date: 2008-09-04 20:54:06 -0700 (Thu, 04 Sep 2008)
New Revision: 12817

Modified:
   short/3D/PyLith/trunk/libsrc/meshio/MeshIOLagrit.cc
Log:
Kludgey fix for multiple fault problem


Modified: short/3D/PyLith/trunk/libsrc/meshio/MeshIOLagrit.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/meshio/MeshIOLagrit.cc	2008-09-05 03:16:39 UTC (rev 12816)
+++ short/3D/PyLith/trunk/libsrc/meshio/MeshIOLagrit.cc	2008-09-05 03:54:06 UTC (rev 12817)
@@ -157,10 +157,15 @@
       fin >> faceCells[c*2+1];
     }
 
+    // Determine the number of cells
+    //   Only do this once since we add cohesive cells after that
+    static int numCells = -1;
+
+    if (numCells == -1) {numCells = mesh->heightStratum(0)->size();}
+
     // Renumber vertices and use zero based indices
     // UCD file has one-based indices for both vertexIDs and fCells
     //   Also, vertex numbers are offset by the number of cells
-    const int numCells = mesh->heightStratum(0)->size();
     for(int c = 0; c < numFCells; ++c)
       for(int corner = 0; corner < numFCorners; ++corner)
         fCells[c*numFCorners+corner] =