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

[knepley at geodynamics.org]

Author: knepley
Date: 2007-04-19 15:16:15 -0700 (Thu, 19 Apr 2007)
New Revision: 6613

Modified:
   short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc
Log:
Fixed material ids on fault


Modified: short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc	2007-04-19 21:58:41 UTC (rev 6612)
+++ short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc	2007-04-19 22:16:15 UTC (rev 6613)
@@ -144,6 +144,7 @@
 
   // Split the mesh along the fault sieve and create cohesive elements
   const ALE::Obj<Mesh::label_sequence>& faces = fault->depthStratum(1);
+  const ALE::Obj<Mesh::label_type>& material = mesh->getLabel("material-id");
   PointArray newVertices;
   
   for(Mesh::label_sequence::iterator f_iter = faces->begin();
@@ -199,6 +200,7 @@
         std::cout << "    vertex " << vertexRenumber[*v_iter] << std::endl;
       sieve->addArrow(vertexRenumber[*v_iter], newPoint, color++);
     }
+    mesh->setValue(material, newPoint, materialId);
   } // for
   mesh->stratify();
   if (debug)