[cig-commits] r15293 - in short/3D/PyLith/trunk: libsrc/problems pylith/apps

[knepley at geodynamics.org]

Author: knepley
Date: 2009-06-15 19:30:06 -0700 (Mon, 15 Jun 2009)
New Revision: 15293

Modified:
   short/3D/PyLith/trunk/libsrc/problems/SolverNonlinear.cc
   short/3D/PyLith/trunk/pylith/apps/PyLithApp.py
Log:
Fixed SNES init


Modified: short/3D/PyLith/trunk/libsrc/problems/SolverNonlinear.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/problems/SolverNonlinear.cc	2009-06-16 01:22:10 UTC (rev 15292)
+++ short/3D/PyLith/trunk/libsrc/problems/SolverNonlinear.cc	2009-06-16 02:30:06 UTC (rev 15293)
@@ -67,7 +67,6 @@
     CHECK_PETSC_ERROR(err);
   } // if    
   err = SNESCreate(fields.mesh().comm(), &_snes); CHECK_PETSC_ERROR(err);
-  err = SNESSetFromOptions(_snes); CHECK_PETSC_ERROR(err);
 
   const topology::Field<topology::Mesh>& residual = fields.get("residual");
   const PetscVec residualVec = residual.vector();
@@ -79,6 +78,8 @@
   err = SNESSetJacobian(_snes, jacobianMat, jacobianMat, reformJacobian,
 			(void*) formulation);
   CHECK_PETSC_ERROR(err);
+
+  err = SNESSetFromOptions(_snes); CHECK_PETSC_ERROR(err);
 } // initialize
 
 // ----------------------------------------------------------------------

Modified: short/3D/PyLith/trunk/pylith/apps/PyLithApp.py
===================================================================
--- short/3D/PyLith/trunk/pylith/apps/PyLithApp.py	2009-06-16 01:22:10 UTC (rev 15292)
+++ short/3D/PyLith/trunk/pylith/apps/PyLithApp.py	2009-06-16 02:30:06 UTC (rev 15293)
@@ -112,11 +112,11 @@
     self.problem.finalize()
     self._eventLogger.stagePop()
 
+    self.compilePerformanceLog()
+    if self.perfLogger.verbose: self.perfLogger.show()
+
     del mesh
     del self.problem
-
-    self.compilePerformanceLog()
-    if self.perfLogger.verbose: self.perfLogger.show()
     return