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

[brad at geodynamics.org]

Author: brad
Date: 2008-10-03 09:55:11 -0700 (Fri, 03 Oct 2008)
New Revision: 12990

Modified:
   short/3D/PyLith/trunk/pylith/problems/Formulation.py
Log:
Added missing mat_assemble() after integrateResidualAssembled(). Need to get matrix in correct state.

Modified: short/3D/PyLith/trunk/pylith/problems/Formulation.py
===================================================================
--- short/3D/PyLith/trunk/pylith/problems/Formulation.py	2008-10-03 16:53:00 UTC (rev 12989)
+++ short/3D/PyLith/trunk/pylith/problems/Formulation.py	2008-10-03 16:55:11 UTC (rev 12990)
@@ -431,6 +431,9 @@
       integrator.integrateJacobianAssembled(self.jacobian, t+dt, self.fields)
       self._debug.log(resourceUsageString()) # TEMPORARY
 
+    self._info.log("Assembling Jacobian of operator.")
+    petsc.mat_assemble(self.jacobian)
+
     if self.viewJacobian:
       filename = self._createJacobianFilename(t+dt)
       petsc.mat_view_binary(self.jacobian, filename)