[cig-commits] r16037 - short/3D/PyLith/branches/pylith-friction/libsrc/problems

[brad at geodynamics.org]

Author: brad
Date: 2009-11-25 08:50:18 -0800 (Wed, 25 Nov 2009)
New Revision: 16037

Modified:
   short/3D/PyLith/branches/pylith-friction/libsrc/problems/SolverNonlinear.cc
Log:
Added use of customized line search for nonlinear solver.

Modified: short/3D/PyLith/branches/pylith-friction/libsrc/problems/SolverNonlinear.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/libsrc/problems/SolverNonlinear.cc	2009-11-25 16:46:26 UTC (rev 16036)
+++ short/3D/PyLith/branches/pylith-friction/libsrc/problems/SolverNonlinear.cc	2009-11-25 16:50:18 UTC (rev 16037)
@@ -180,7 +180,7 @@
   *flag   = PETSC_TRUE;
 
   // ======================================================================
-  // Code to constraint solution space.
+  // Code to constrain solution space.
   //
   // Matt- It seems like I should adjust both x (current iterate) and
   // the search direction, with the Lagrange multipliers in x modified
@@ -190,7 +190,7 @@
   assert(0 != lsctx);
   Formulation* formulation = (Formulation*) lsctx;
   assert(0 != formulation);
-  formulation->constrainSolnSpace(&y);
+  formulation->constrainSolnSpace(&x); // Adjust x? y? both?
   // ======================================================================
 
   ierr = VecNorm(y,NORM_2,ynorm);CHKERRQ(ierr);