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

[brad at geodynamics.org]

Author: brad
Date: 2007-06-21 08:39:19 -0700 (Thu, 21 Jun 2007)
New Revision: 7342

Modified:
   short/3D/PyLith/trunk/pylith/problems/Explicit.py
   short/3D/PyLith/trunk/pylith/problems/Implicit.py
Log:
Fixed call to reformJacobian in Implicit.py. Added use of reformJacobian to Explicit.py

Modified: short/3D/PyLith/trunk/pylith/problems/Explicit.py
===================================================================
--- short/3D/PyLith/trunk/pylith/problems/Explicit.py	2007-06-21 06:54:54 UTC (rev 7341)
+++ short/3D/PyLith/trunk/pylith/problems/Explicit.py	2007-06-21 15:39:19 UTC (rev 7342)
@@ -117,13 +117,7 @@
       if integrator.needNewJacobian():
         needNewJacobian = True
     if needNewJacobian:
-      self._info.log("Reforming Jacobian of operator.")
-      import pylith.utils.petsc as petsc
-      petsc.mat_setzero(self.jacobian)
-      for integrator in self.integrators:
-        integrator.timeStep(dt)
-        integrator.integrateJacobian(self.jacobian, t, self.fields)
-      petsc.mat_assemble(self.jacobian)
+      self._reformJacobian(t, dt)
     return
 
 

Modified: short/3D/PyLith/trunk/pylith/problems/Implicit.py
===================================================================
--- short/3D/PyLith/trunk/pylith/problems/Implicit.py	2007-06-21 06:54:54 UTC (rev 7341)
+++ short/3D/PyLith/trunk/pylith/problems/Implicit.py	2007-06-21 15:39:19 UTC (rev 7342)
@@ -137,7 +137,7 @@
       if integrator.needNewJacobian():
         needNewJacobian = True
     if needNewJacobian:
-      self._reformJacobian()
+      self._reformJacobian(t, dt)
     return