[cig-commits] r18347 - short/3D/PyLith/branches/v1.5-stable/libsrc/problems

[brad at geodynamics.org]

Author: brad
Date: 2011-05-11 16:48:12 -0700 (Wed, 11 May 2011)
New Revision: 18347

Modified:
   short/3D/PyLith/branches/v1.5-stable/libsrc/problems/SolverNonlinear.cc
Log:
Added call to constrainSoln() before integrating residual in nonlinear solve so residual reflects constrained solution. Otherwise, residual for unconstrained solution may be zero and SNES thinks solution has converged.

Modified: short/3D/PyLith/branches/v1.5-stable/libsrc/problems/SolverNonlinear.cc
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/libsrc/problems/SolverNonlinear.cc	2011-05-11 23:45:49 UTC (rev 18346)
+++ short/3D/PyLith/branches/v1.5-stable/libsrc/problems/SolverNonlinear.cc	2011-05-11 23:48:12 UTC (rev 18347)
@@ -164,6 +164,10 @@
   Formulation* formulation = (Formulation*) context;
   assert(0 != formulation);
 
+
+  // TEMPORARY - ASK MATT
+  formulation->constrainSolnSpace(&tmpSolutionVec);
+
   // Reform residual
   formulation->reformResidual(&tmpResidualVec, &tmpSolutionVec);