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

[willic3 at geodynamics.org]

Author: willic3
Date: 2007-05-29 09:25:26 -0700 (Tue, 29 May 2007)
New Revision: 6991

Modified:
   short/3D/PyLith/trunk/pylith/problems/Implicit.py
Log:
Added some comments to describe functionality.



Modified: short/3D/PyLith/trunk/pylith/problems/Implicit.py
===================================================================
--- short/3D/PyLith/trunk/pylith/problems/Implicit.py	2007-05-29 16:24:09 UTC (rev 6990)
+++ short/3D/PyLith/trunk/pylith/problems/Implicit.py	2007-05-29 16:25:26 UTC (rev 6991)
@@ -141,6 +141,8 @@
         integrator.timeStep(dt)
         integrator.integrateJacobian(self.jacobian, self.fields.cppHandle)
       petsc.mat_assemble(self.jacobian)
+    # Put in loop over integrators to see if stiffness needs
+    # reforming.  If so, then reform it.
     return
 
 
@@ -181,6 +183,9 @@
     self._info.log("Updating integrators states.")
     for integrator in self.integrators:
       integrator.updateState(self.fields.getReal("dispT"))
+    # call update state here.
+    # This will recompute stresses and state variables.
+    # Need to loop over integrators, and put this into integrator.
     return