[cig-commits] r14876 - in short/3D/PyLith/trunk: . libsrc/bc

[brad at geodynamics.org]

Author: brad
Date: 2009-05-05 13:39:13 -0700 (Tue, 05 May 2009)
New Revision: 14876

Modified:
   short/3D/PyLith/trunk/TODO
   short/3D/PyLith/trunk/libsrc/bc/AbsorbingDampers.cc
   short/3D/PyLith/trunk/libsrc/bc/Neumann.cc
Log:
Cleaned up use of heightStratum() and depthStratum().

Modified: short/3D/PyLith/trunk/TODO
===================================================================
--- short/3D/PyLith/trunk/TODO	2009-05-05 20:32:30 UTC (rev 14875)
+++ short/3D/PyLith/trunk/TODO	2009-05-05 20:39:13 UTC (rev 14876)
@@ -46,8 +46,6 @@
     Check use of label_sequence.  label_sequence - iterators are
     cached, so use sequence or cache begin/end to maintain access
 
-    Cleanup use of heightStratum() and depthStratum().
-
     Cleanup logging. Constraints and Integrators should log at the C++
     level using the C++ EventLogger. Add finer grain logging at C++
     level as in ElasticityImplicit.

Modified: short/3D/PyLith/trunk/libsrc/bc/AbsorbingDampers.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/bc/AbsorbingDampers.cc	2009-05-05 20:32:30 UTC (rev 14875)
+++ short/3D/PyLith/trunk/libsrc/bc/AbsorbingDampers.cc	2009-05-05 20:39:13 UTC (rev 14876)
@@ -84,13 +84,13 @@
   // Get 'surface' cells (1 dimension lower than top-level cells)
   const ALE::Obj<SieveSubMesh>& submesh = _boundaryMesh->sieveMesh();
   assert(!submesh.isNull());
-  const int boundaryDepth = submesh->depth()-1; // depth of bndry cells
   const ALE::Obj<SieveSubMesh::label_sequence>& cells = 
-    submesh->depthStratum(boundaryDepth);
+    submesh->heightStratum(1);
   assert(!cells.isNull());
   const SieveSubMesh::label_sequence::iterator cellsEnd = cells->end();
 
   // Make sure surface cells are compatible with quadrature.
+  const int boundaryDepth = submesh->depth()-1; // depth of bndry cells
   for (SieveSubMesh::label_sequence::iterator c_iter=cells->begin();
        c_iter != cellsEnd;
        ++c_iter) {
@@ -267,9 +267,8 @@
   // Get cell information
   const ALE::Obj<SieveSubMesh>& submesh = _boundaryMesh->sieveMesh();
   assert(!submesh.isNull());
-  const int boundaryDepth = submesh->depth()-1; // depth of bndry cells
   const ALE::Obj<SieveSubMesh::label_sequence>& cells = 
-    submesh->depthStratum(boundaryDepth);
+    submesh->heightStratum(1);
   assert(!cells.isNull());
   const SieveSubMesh::label_sequence::iterator cellsEnd = cells->end();
 
@@ -374,9 +373,8 @@
   // Get cell information
   const ALE::Obj<SieveSubMesh>& submesh = _boundaryMesh->sieveMesh();
   assert(!submesh.isNull());
-  const int boundaryDepth = submesh->depth()-1; // depth of bndry cells
   const ALE::Obj<SieveSubMesh::label_sequence>& cells = 
-    submesh->depthStratum(boundaryDepth);
+    submesh->heightStratum(1);
   assert(!cells.isNull());
   const SieveSubMesh::label_sequence::iterator cellsEnd = cells->end();
 

Modified: short/3D/PyLith/trunk/libsrc/bc/Neumann.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/bc/Neumann.cc	2009-05-05 20:32:30 UTC (rev 14875)
+++ short/3D/PyLith/trunk/libsrc/bc/Neumann.cc	2009-05-05 20:39:13 UTC (rev 14876)
@@ -81,13 +81,13 @@
   // Get 'surface' cells (1 dimension lower than top-level cells)
   const ALE::Obj<SieveSubMesh>& submesh = _boundaryMesh->sieveMesh();
   assert(!submesh.isNull());
-  const int boundaryDepth = submesh->depth()-1; // depth of bndry cells
   const ALE::Obj<SieveSubMesh::label_sequence>& cells = 
-    submesh->depthStratum(boundaryDepth);
+    submesh->heightStratum(1);
   assert(!cells.isNull());
   const SieveSubMesh::label_sequence::iterator cellsEnd = cells->end();
 
   // Make sure surface cells are compatible with quadrature.
+  const int boundaryDepth = submesh->depth()-1; // depth of bndry cells
   for (SieveSubMesh::label_sequence::iterator c_iter=cells->begin();
        c_iter != cellsEnd;
        ++c_iter) {
@@ -266,9 +266,8 @@
   // Get cell information
   const ALE::Obj<SieveSubMesh>& submesh = _boundaryMesh->sieveMesh();
   assert(!submesh.isNull());
-  const int boundaryDepth = submesh->depth()-1; // depth of bndry cells
   const ALE::Obj<SieveSubMesh::label_sequence>& cells = 
-    submesh->depthStratum(boundaryDepth);
+    submesh->heightStratum(1);
   assert(!cells.isNull());
   const SieveSubMesh::label_sequence::iterator cellsEnd = cells->end();