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

[brad at geodynamics.org]

Author: brad
Date: 2008-02-05 16:45:07 -0800 (Tue, 05 Feb 2008)
New Revision: 9249

Modified:
   short/3D/PyLith/trunk/pylith/problems/Implicit.py
Log:
Fixed bug in forming residual (arg is t+dt, not t).

Modified: short/3D/PyLith/trunk/pylith/problems/Implicit.py
===================================================================
--- short/3D/PyLith/trunk/pylith/problems/Implicit.py	2008-02-05 23:39:17 UTC (rev 9248)
+++ short/3D/PyLith/trunk/pylith/problems/Implicit.py	2008-02-06 00:45:07 UTC (rev 9249)
@@ -162,12 +162,18 @@
     import pylith.topology.topology as bindings
     bindings.zeroRealSection(dispIncr)
 
-    self._reformResidual(t, dt)
+    self._reformResidual(t+dt, dt)
 
     self._info.log("Solving equations.")
     residual = self.fields.getReal("residual")
     self.solver.solve(dispIncr, self.jacobian, residual)
 
+    # BEGIN TEMPORARY
+    import pylith.topology.topology as bindings
+    bindings.sectionView(self.fields.getReal("dispIncr"), "SOLUTION");
+    bindings.sectionView(self.fields.getReal("residual"), "RESIDUAL");
+    # END TEMPORARY
+
     self._logger.eventEnd(logEvent)
     return