[cig-commits] r7032 - in short/3D/PyLith/trunk: . pylith/problems

[brad at geodynamics.org]

Author: brad
Date: 2007-06-01 10:10:08 -0700 (Fri, 01 Jun 2007)
New Revision: 7032

Modified:
   short/3D/PyLith/trunk/TODO
   short/3D/PyLith/trunk/pylith/problems/Implicit.py
Log:
Removed redundant solver.initialize() from Implicit time integration formulation.

Modified: short/3D/PyLith/trunk/TODO
===================================================================
--- short/3D/PyLith/trunk/TODO	2007-06-01 16:40:51 UTC (rev 7031)
+++ short/3D/PyLith/trunk/TODO	2007-06-01 17:10:08 UTC (rev 7032)
@@ -2,11 +2,25 @@
 MAIN PRIORITIES (Brad)
 ======================================================================
 
-1. Unit tests for Integrator stuff.
+0. Unit tests for Integrator stuff.
 
-  a. Add unit test for IntegratorElasticity::calcTotalStrain
+   a. Add unit test for IntegratorElasticity::calcTotalStrain
+   b. Incorporate Quad/Hex basis stuff
 
-2. Add dualBasis to Quadrature.
+1. Simple tests with analytical solutions
+
+   a. 1-D: axial compression
+   b. 2-D: axial compression of a "box"
+   c. 2-DD: shearing of a "box"
+   d. 3-D: axial compression of a "box"
+   3. 3-D: shearing of a "box"
+
+2. Allow use of all elasticity constants (9 for 2-D, 36 for 3-D).
+   a. Materials C++ code
+   b. Integrator C++ code
+   b. Material C++ unit tests
+
+3. Add dualBasis to Quadrature.
    a. Python
      ReferenceCell
      FIATSimplex
@@ -16,7 +30,7 @@
    c. C++ unit tests
    d. Python unit tests
 
-3. Implement faults for kinematic source
+4. Implement faults for kinematic source
    a. Creation of cohesive cells
      i. Add tests for interpolated meshes.
         Double check consistency in ordering of vertices (positive/negative).
@@ -63,8 +77,6 @@
 
 4. Implement absorbing boundary conditions
 
-5. Create suite of simple full test cases
-
 ======================================================================
 SECONDARY PRIORITIES
 ======================================================================

Modified: short/3D/PyLith/trunk/pylith/problems/Implicit.py
===================================================================
--- short/3D/PyLith/trunk/pylith/problems/Implicit.py	2007-06-01 16:40:51 UTC (rev 7031)
+++ short/3D/PyLith/trunk/pylith/problems/Implicit.py	2007-06-01 17:10:08 UTC (rev 7032)
@@ -100,10 +100,10 @@
     self.fields.copyLayout("dispTBctpdt")
     self.jacobian = mesh.createMatrix(self.fields.getReal("dispTBctpdt"))
 
+    self.solver.initialize(mesh, self.fields.getReal("dispIncr"))
+
     from pyre.units.time import s
     self._solveElastic(mesh, materials, t=0.0*s, dt=dt)
-
-    self.solver.initialize(mesh, self.fields.getReal("dispIncr"))
     return
 
 
@@ -219,8 +219,6 @@
     import pylith.utils.petsc as petsc
     petsc.mat_assemble(self.jacobian)
 
-    self.solver.initialize(mesh, dispIncr)
-
     self._info.log("Solving equations.")
     self.solver.solve(dispIncr, self.jacobian, residual)