[cig-commits] r4226 - short/3D/PyLith/branches/pylith-0.8/pylith3d/module

[knepley at geodynamics.org]

Author: knepley
Date: 2006-08-04 12:18:23 -0700 (Fri, 04 Aug 2006)
New Revision: 4226

Modified:
   short/3D/PyLith/branches/pylith-0.8/pylith3d/module/scanner.cc
Log:
Fixed ALE::Point refrences to Sieve points


Modified: short/3D/PyLith/branches/pylith-0.8/pylith3d/module/scanner.cc
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/module/scanner.cc	2006-08-04 04:56:28 UTC (rev 4225)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/module/scanner.cc	2006-08-04 19:18:23 UTC (rev 4226)
@@ -166,7 +166,7 @@
     }
 
     if (constraints[0] || constraints[1] || constraints[2]) {
-      fprintf(f, "%7d %4d %4d %4d % 16.8E % 16.8E % 16.8E\n", (*v_iter).index+1-numElements,
+      fprintf(f, "%7d %4d %4d %4d % 16.8E % 16.8E % 16.8E\n", *v_iter+1-numElements,
               constraints[0], constraints[1], constraints[2], values[0], values[1], values[2]);
     }
   }
@@ -224,21 +224,21 @@
   boundaries->getAtlas()->setTopology(mesh->getTopologyNew());
   // Reverse order allows newer conditions to override older, as required by PyLith
   for(int v = numBoundaryVertices-1; v >= 0; v--) {
-    ALE::Mesh::point_type vertex(0, boundaryVertices[v*(numBoundaryComponents+1)] + numElements);
+    ALE::Mesh::point_type vertex(boundaryVertices[v*(numBoundaryComponents+1)] + numElements);
     int size = 0;
 
-    if (seen.find(vertex.index) == seen.end()) {
+    if (seen.find(vertex) == seen.end()) {
       for(int c = 0; c < numBoundaryComponents; c++) {
         size += boundaryVertices[v*(numBoundaryComponents+1)+c+1];
       }
       boundaries->getAtlas()->setFiberDimension(patch, vertex, size);
-      seen.insert(vertex.index);
+      seen.insert(vertex);
     }
   }
   boundaries->getAtlas()->orderPatches();
   boundaries->allocate();
   for(int v = 0; v < numBoundaryVertices; v++) {
-    ALE::Mesh::point_type vertex(0, boundaryVertices[v*(numBoundaryComponents+1)] + numElements);
+    ALE::Mesh::point_type vertex(boundaryVertices[v*(numBoundaryComponents+1)] + numElements);
     ALE::Mesh::foliated_section_type::value_type values[3];
 
     for(int c = 0, i = 0; c < numBoundaryComponents; c++) {