[cig-commits] r20715 - short/3D/PyLith/trunk/libsrc/pylith/problems

[knepley at geodynamics.org]

Author: knepley
Date: 2012-09-15 22:00:17 -0700 (Sat, 15 Sep 2012)
New Revision: 20715

Modified:
   short/3D/PyLith/trunk/libsrc/pylith/problems/SolverLinear.cc
Log:
Converted SolverLinear

Modified: short/3D/PyLith/trunk/libsrc/pylith/problems/SolverLinear.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/pylith/problems/SolverLinear.cc	2012-09-16 04:40:13 UTC (rev 20714)
+++ short/3D/PyLith/trunk/libsrc/pylith/problems/SolverLinear.cc	2012-09-16 05:00:17 UTC (rev 20715)
@@ -119,13 +119,13 @@
   } // else
   jacobian->resetValuesChanged();
 
+#if 0 // OBSOLETE? Replaced by stuff in Solver?
   // Update KSP operators with custom preconditioner if necessary.
   const ALE::Obj<RealSection>& solutionSection = solution->section();
   assert(!solutionSection.isNull());
   const ALE::Obj<SieveMesh>& sieveMesh = solution->mesh().sieveMesh();
   assert(!sieveMesh.isNull());
 
-#if 0 // OBSOLETE? Replaced by stuff in Solver?
   if (solutionSection->getNumSpaces() > sieveMesh->getDimension() &&
       _jacobianPCFault) {
     
@@ -150,8 +150,8 @@
   } // if
 #endif
 
-  const PetscVec residualVec = residual.vector();
-  const PetscVec solutionVec = solution->vector();
+  const PetscVec residualVec = residual.globalVector();
+  const PetscVec solutionVec = solution->globalVector();
 
   _logger->eventEnd(setupEvent);
   _logger->eventBegin(solveEvent);