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

[brad at geodynamics.org]

Author: brad
Date: 2008-10-15 13:01:58 -0700 (Wed, 15 Oct 2008)
New Revision: 13091

Modified:
   short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveKin.cc
Log:
Replaced calls to ALE::Completion::completeSectionAdd() with calls to ALE::Distribution<pylith::Mesh>::completeSection() due to bug in completeSectionAdd().

Modified: short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveKin.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveKin.cc	2008-10-15 19:21:35 UTC (rev 13090)
+++ short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveKin.cc	2008-10-15 20:01:58 UTC (rev 13091)
@@ -627,7 +627,11 @@
   _orientation->setDebug(2);
 #endif
   // Assemble orientation information
+#if 1
+  ALE::Distribution<pylith::Mesh>::completeSection(*_faultMesh, _orientation);
+#else
   ALE::Completion::completeSectionAdd(_faultMesh->getSendOverlap(), _faultMesh->getRecvOverlap(), _orientation, _orientation);
+#endif
 
 #if 0
   _orientation->view("ORIENTATION After complete");
@@ -833,10 +837,22 @@
     PetscLogFlops( numQuadPts*(1+numBasis*2) );
   } // for
 
+#if 0
+  _area->view("AREA");
+#endif
+
   // Assemble area information
+#if 1
+  ALE::Distribution<pylith::Mesh>::completeSection(*_faultMesh, _area);
+#else
   ALE::Completion::completeSectionAdd(_faultMesh->getSendOverlap(), _faultMesh->getRecvOverlap(), _area, _area);
+#endif
 
-  //_area->view("AREA");
+#if 0
+  _area->view("AREA");
+  _faultMesh->getSendOverlap()->view("Send fault overlap");
+  _faultMesh->getRecvOverlap()->view("Receive fault overlap");
+#endif
 } // _calcArea
 
 // ----------------------------------------------------------------------
@@ -854,11 +870,6 @@
   assert(!_pseudoStiffness.isNull());
   assert(!_area.isNull());
 
-#if 0
-  _pseudoStiffness->view("PSEUDOSTIFFNESS");
-  _area->view("AREA");
-#endif
-
   const ALE::Obj<Mesh::label_sequence>& vertices    = mesh->depthStratum(0);
   assert(!vertices.isNull());
   const Mesh::label_sequence::iterator  verticesEnd = vertices->end();