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

[knepley at geodynamics.org]

Author: knepley
Date: 2007-06-14 10:05:06 -0700 (Thu, 14 Jun 2007)
New Revision: 7232

Modified:
   short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc
Log:
Fixed small bug with groups in CohesiveTopology


Modified: short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc	2007-06-14 16:26:48 UTC (rev 7231)
+++ short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc	2007-06-14 17:05:06 UTC (rev 7232)
@@ -150,14 +150,18 @@
       const ALE::Obj<int_section_type>& group = mesh->getIntSection(*name);
 
       if (group->hasPoint(*v_iter)) {
-        group->setFiberDimension(newPoint, 1);
+        group->addPoint(newPoint, 1);
         if (constraintCell) {
-          group->setFiberDimension(newPoint+1, 1);
+          group->addPoint(newPoint+1, 1);
         }
       }
     } // for
     if (constraintCell) newPoint++;
   } // for
+  for(std::set<std::string>::const_iterator name = groupNames->begin();
+      name != groupNames->end(); ++name) {
+    mesh->reallocate(mesh->getIntSection(*name));
+  } // for
 
   // Split the mesh along the fault sieve and create cohesive elements
   const ALE::Obj<Mesh::label_sequence>& faces = (*fault)->depthStratum(1);
@@ -261,7 +265,6 @@
   const ALE::Obj<Mesh::label_type>&           label          = mesh->createLabel(std::string("censored depth"));
   const ALE::Obj<std::set<Mesh::point_type> > modifiedPoints = new std::set<Mesh::point_type>();
   _computeCensoredDepth(mesh, label, mesh->getSieve(), mesh->getSieve()->roots(), firstCohesiveCell, modifiedPoints);
-  label->view("");
   if (debug)
     mesh->view("Mesh with Cohesive Elements");