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

[brad at geodynamics.org]

Author: brad
Date: 2008-08-12 15:55:23 -0700 (Tue, 12 Aug 2008)
New Revision: 12618

Modified:
   short/3D/PyLith/trunk/pylith/problems/Implicit.py
Log:
TEMPORATRY: Added zeroing of entries immediately after creating Jacobian matrix to get better estimate of memory usage.

Modified: short/3D/PyLith/trunk/pylith/problems/Implicit.py
===================================================================
--- short/3D/PyLith/trunk/pylith/problems/Implicit.py	2008-08-12 21:46:39 UTC (rev 12617)
+++ short/3D/PyLith/trunk/pylith/problems/Implicit.py	2008-08-12 22:55:23 UTC (rev 12618)
@@ -102,6 +102,13 @@
 
     self._info.log("Creating Jacobian matrix.")
     self.jacobian = self.mesh.createMatrix(self.fields.getSolution())
+
+    # BEGIN TEMPORARY
+    # Access entries in matrix here to get correct memory usage
+    import pylith.utils.petsc as petsc
+    petsc.mat_setzero(self.jacobian)
+    # END TEMPORARY
+
     self._debug.log(resourceUsageString())
 
     self._info.log("Initializing solver.")