[cig-commits] r16795 - in short/3D/PyLith/trunk: . libsrc/faults

[brad at geodynamics.org]

Author: brad
Date: 2010-05-26 17:36:29 -0700 (Wed, 26 May 2010)
New Revision: 16795

Modified:
   short/3D/PyLith/trunk/TODO
   short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc
   short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveLagrange.cc
Log:
Added TODO comments.

Modified: short/3D/PyLith/trunk/TODO
===================================================================
--- short/3D/PyLith/trunk/TODO	2010-05-27 00:25:32 UTC (rev 16794)
+++ short/3D/PyLith/trunk/TODO	2010-05-27 00:36:29 UTC (rev 16795)
@@ -109,6 +109,9 @@
 
 CODE
 
+  * Fix parallel bugs
+    + Need to limit MatGetValues to local.
+
   * Better preconditioning [Matt/Brad]
 
     Diagonal preconditioner does not work very well. It looks like the

Modified: short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc	2010-05-27 00:25:32 UTC (rev 16794)
+++ short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveDyn.cc	2010-05-27 00:36:29 UTC (rev 16795)
@@ -1669,6 +1669,8 @@
 
   const int indexOffset = (negativeSide) ? 0 : 3*submatrixSize + subnrows;
 
+  // :TODO: Use globalOrder->isLocal(idx) to make sure indices are local.
+
   for (SieveMesh::label_sequence::iterator c_iter=cellsCohesiveBegin;
        c_iter != cellsCohesiveEnd;
        ++c_iter) {

Modified: short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveLagrange.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveLagrange.cc	2010-05-27 00:25:32 UTC (rev 16794)
+++ short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveLagrange.cc	2010-05-27 00:36:29 UTC (rev 16795)
@@ -690,6 +690,8 @@
     indicesN = indicesRel + globalOrder->getIndex(v_negative);
     indicesP = indicesRel + globalOrder->getIndex(v_positive);
 
+    // Use globalOrder->isLocal(idx) to make sure indices are local.
+
     PetscErrorCode err = 0;
 
     err = MatGetValues(jacobianMatrix,