[cig-commits] r7015 - in short/3D/PyLith/trunk: examples/twotri3 pylith/solver

[brad at geodynamics.org]

Author: brad
Date: 2007-05-31 10:52:01 -0700 (Thu, 31 May 2007)
New Revision: 7015

Modified:
   short/3D/PyLith/trunk/examples/twotri3/pylithapp.cfg
   short/3D/PyLith/trunk/pylith/solver/SolverLinear.py
Log:
Added a little more error checking to Python SolverLinear. Added some additional settings to twotri3 example configuration file.

Modified: short/3D/PyLith/trunk/examples/twotri3/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/twotri3/pylithapp.cfg	2007-05-31 17:32:08 UTC (rev 7014)
+++ short/3D/PyLith/trunk/examples/twotri3/pylithapp.cfg	2007-05-31 17:52:01 UTC (rev 7015)
@@ -5,10 +5,11 @@
 # journal
 # ----------------------------------------------------------------------
 [pylithapp.journal.info]
-eqdeformation = 1
+timedependent = 1
 explicit = 1
 implicit = 1
 petsc = 1
+solverlinear = 1
 meshioascii = 1
 homogeneous = 1
 explicitelasticity = 1
@@ -45,3 +46,9 @@
 db.iohandler.filename = matprops.spatialdb
 quadrature = pylith.feassemble.quadrature.Quadrature2D
 quadrature.cell.shape = triangle
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+pc_type = jacobi

Modified: short/3D/PyLith/trunk/pylith/solver/SolverLinear.py
===================================================================
--- short/3D/PyLith/trunk/pylith/solver/SolverLinear.py	2007-05-31 17:32:08 UTC (rev 7014)
+++ short/3D/PyLith/trunk/pylith/solver/SolverLinear.py	2007-05-31 17:52:01 UTC (rev 7015)
@@ -67,6 +67,7 @@
     Solve linear system.
     """
     self._info.log("Solving linear equations.")
+    assert(None != self.cppHandle)
     self.cppHandle.solve(fieldOut, jacobian, fieldIn)
     return