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

[knepley at geodynamics.org]

Author: knepley
Date: 2008-01-18 15:11:58 -0800 (Fri, 18 Jan 2008)
New Revision: 9092

Modified:
   short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc
Log:
Fixed censored depth by adding new shadow vertices


Modified: short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc	2008-01-18 00:41:42 UTC (rev 9091)
+++ short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc	2008-01-18 23:11:58 UTC (rev 9092)
@@ -601,6 +601,13 @@
     const ALE::Obj<PointSet>          modifiedPoints = new PointSet();
 
     _computeCensoredDepth(mesh, label, mesh->getSieve(), mesh->getSieve()->roots(), firstCohesiveCell-(constraintCell?numFaultVertices:0), modifiedPoints);
+  } else {
+    // Insert new shadow vertices into existing label
+    const ALE::Obj<Mesh::label_type>& label          = mesh->getLabel(labelName);
+
+    for(std::map<int,int>::const_iterator v_iter = vertexRenumber.begin(); v_iter != vertexRenumber.end(); ++v_iter) {
+      mesh->setValue(label, v_iter->second, 0);
+    }
   }
   if (debug) mesh->view("Mesh with Cohesive Elements");