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

brad at geodynamics.org brad at geodynamics.org
Tue Oct 14 17:04:20 PDT 2008


Author: brad
Date: 2008-10-14 17:04:20 -0700 (Tue, 14 Oct 2008)
New Revision: 13031

Modified:
   short/3D/PyLith/trunk/pylith/problems/Formulation.py
Log:
Fixed bug in where completeResidual is called in integrateResidual().

Modified: short/3D/PyLith/trunk/pylith/problems/Formulation.py
===================================================================
--- short/3D/PyLith/trunk/pylith/problems/Formulation.py	2008-10-14 22:59:43 UTC (rev 13030)
+++ short/3D/PyLith/trunk/pylith/problems/Formulation.py	2008-10-15 00:04:20 UTC (rev 13031)
@@ -470,6 +470,8 @@
     for integrator in self.integrators:
       integrator.timeStep(dt)
       integrator.integrateResidual(residual, t, self.fields)
+    self._info.log("Completing residual.")
+    bindings.completeSection(self.mesh.cppHandle, residual)
 
     # Add in contributions that do not require assembly
     self._info.log("Integrating assembled residual term in operator.")
@@ -477,9 +479,6 @@
       integrator.timeStep(dt)
       integrator.integrateResidualAssembled(residual, t, self.fields)
 
-    self._info.log("Completing residual.")
-    bindings.completeSection(self.mesh.cppHandle, residual)
-
     return
 
 



More information about the cig-commits mailing list