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

[brad at geodynamics.org]

Author: brad
Date: 2013-04-04 09:40:46 -0700 (Thu, 04 Apr 2013)
New Revision: 21715

Modified:
   short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveDyn.cc
   short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveLagrange.cc
Log:
Fixed getting fault cells in submesh (height should be 1).

Modified: short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveDyn.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveDyn.cc	2013-04-04 16:32:40 UTC (rev 21714)
+++ short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveDyn.cc	2013-04-04 16:40:46 UTC (rev 21715)
@@ -1811,7 +1811,7 @@
 
   // Get fault cell information
   PetscDM faultDMMesh = _faultMesh->dmMesh();assert(faultDMMesh);
-  topology::Stratum cellsStratum(faultDMMesh, topology::Stratum::HEIGHT, 0);
+  topology::Stratum cellsStratum(faultDMMesh, topology::Stratum::HEIGHT, 1);
   const PetscInt cStart = cellsStratum.begin();
   const PetscInt cEnd = cellsStratum.end();
   const PetscInt numCells = cellsStratum.size();

Modified: short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveLagrange.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveLagrange.cc	2013-04-04 16:32:40 UTC (rev 21714)
+++ short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveLagrange.cc	2013-04-04 16:40:46 UTC (rev 21715)
@@ -1029,7 +1029,7 @@
   topology::Stratum faultVerticesStratum(faultDMMesh, topology::Stratum::DEPTH, 0);
   const PetscInt fvStart = faultVerticesStratum.begin();
   const PetscInt fvEnd = faultVerticesStratum.end();
-  topology::Stratum faultCellsStratum(faultDMMesh, topology::Stratum::HEIGHT, 0);
+  topology::Stratum faultCellsStratum(faultDMMesh, topology::Stratum::HEIGHT, 1);
   const PetscInt fcStart = faultCellsStratum.begin();
   const PetscInt fcEnd = faultCellsStratum.end();
   assert(ncells == fcEnd-fcStart);
@@ -1284,7 +1284,7 @@
   const PetscInt vStart = verticesStratum.begin();
   const PetscInt vEnd = verticesStratum.end();
 
-  topology::Stratum cellsStratum(faultDMMesh, topology::Stratum::HEIGHT, 0);
+  topology::Stratum cellsStratum(faultDMMesh, topology::Stratum::HEIGHT, 1);
   const PetscInt cStart = cellsStratum.begin();
   const PetscInt cEnd = cellsStratum.end();
 
@@ -1549,7 +1549,7 @@
   const PetscInt vStart = verticesStratum.begin();
   const PetscInt vEnd = verticesStratum.end();
 
-  topology::Stratum cellsStratum(faultDMMesh, topology::Stratum::HEIGHT, 0);
+  topology::Stratum cellsStratum(faultDMMesh, topology::Stratum::HEIGHT, 1);
   const PetscInt cStart = cellsStratum.begin();
   const PetscInt cEnd = cellsStratum.end();
 
@@ -1837,7 +1837,7 @@
   if (0 == strcasecmp("partition", name)) {
 
     PetscDM faultDMMesh = _faultMesh->dmMesh();assert(faultDMMesh);
-    topology::Stratum cellsStratum(faultDMMesh, topology::Stratum::HEIGHT, 0);
+    topology::Stratum cellsStratum(faultDMMesh, topology::Stratum::HEIGHT, 1);
     const PetscInt cStart = cellsStratum.begin();
     const PetscInt cEnd = cellsStratum.end();