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

[brad at geodynamics.org]

Author: brad
Date: 2009-08-11 10:59:40 -0700 (Tue, 11 Aug 2009)
New Revision: 15531

Modified:
   short/3D/PyLith/trunk/libsrc/meshio/UCDFaultFile.cc
Log:
Fixed bug in creating fault using faces. Still not working.

Modified: short/3D/PyLith/trunk/libsrc/meshio/UCDFaultFile.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/meshio/UCDFaultFile.cc	2009-08-11 17:48:33 UTC (rev 15530)
+++ short/3D/PyLith/trunk/libsrc/meshio/UCDFaultFile.cc	2009-08-11 17:59:40 UTC (rev 15531)
@@ -180,11 +180,18 @@
       for (int i=0; i < 2; ++i)
 	faceCells[c*2+i] -= 1;
   } // if
+
+  // Create Sieve mesh for fault
+  ALE::Obj<pylith::topology::Mesh::SieveSubMesh>& faultSieveMesh =
+    faultMesh->sieveMesh();
+  faultSieveMesh =
+    new pylith::topology::Mesh::SieveSubMesh(mesh.comm(), mesh.dimension()-1,
+					     mesh.debug());
   
   assert(!sieveMesh->getSieve().isNull());
   const int firstFaultCell = 
     sieveMesh->getSieve()->getBaseSize() + sieveMesh->getSieve()->getCapSize();
-  MeshBuilder::buildFaultMesh(faultMesh->sieveMesh(), faultBoundary, 
+  MeshBuilder::buildFaultMesh(faultSieveMesh, faultBoundary, 
 			      fCoordinates, numFVertices, fSpaceDim, fCells, 
 			      numFCells, numFCorners, firstFaultCell, 
 			      faceCells, faultDim);