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

[brad at geodynamics.org]

Author: brad
Date: 2009-06-15 18:22:10 -0700 (Mon, 15 Jun 2009)
New Revision: 15292

Modified:
   short/3D/PyLith/trunk/TODO
   short/3D/PyLith/trunk/libsrc/bc/TimeDependentPoints.cc
   short/3D/PyLith/trunk/pylith/problems/Implicit.py
Log:
Fixed bug that set reset displacement increment to zero after setting constraints (should be before setting constraints).

Modified: short/3D/PyLith/trunk/TODO
===================================================================
--- short/3D/PyLith/trunk/TODO	2009-06-16 01:17:57 UTC (rev 15291)
+++ short/3D/PyLith/trunk/TODO	2009-06-16 01:22:10 UTC (rev 15292)
@@ -13,8 +13,17 @@
 
     Constitutive model
 
+  Time dependent BC
+    fixed_dof -> bc_dof (Python only)
+    bc_dof in TimeDependentPoints, not TimeDependent
+
   Time dependent Neumann BC
 
+
+
+  TimeDependent.py
+    fixed_
+
 BONUS
 
   pylithinfo (Python application, send extra args to PyLith)

Modified: short/3D/PyLith/trunk/libsrc/bc/TimeDependentPoints.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/bc/TimeDependentPoints.cc	2009-06-16 01:17:57 UTC (rev 15291)
+++ short/3D/PyLith/trunk/libsrc/bc/TimeDependentPoints.cc	2009-06-16 01:22:10 UTC (rev 15292)
@@ -326,7 +326,7 @@
       } // if
     } // if
 
-    valueSection->updateAddPoint(p_bc, &valuesVertex[0]);
+    valueSection->updatePoint(p_bc, &valuesVertex[0]);
   } // for
 }  // _calculateValue
 
@@ -442,7 +442,7 @@
       } // if/else
     } // if
 
-    valueSection->updateAddPoint(p_bc, &valuesVertex[0]);
+    valueSection->updatePoint(p_bc, &valuesVertex[0]);
   } // for
 }  // _calculateValueIncr
 

Modified: short/3D/PyLith/trunk/pylith/problems/Implicit.py
===================================================================
--- short/3D/PyLith/trunk/pylith/problems/Implicit.py	2009-06-16 01:17:57 UTC (rev 15291)
+++ short/3D/PyLith/trunk/pylith/problems/Implicit.py	2009-06-16 01:22:10 UTC (rev 15292)
@@ -161,6 +161,7 @@
 
     self._info.log("Setting constraints.")
     dispIncr = self.fields.get("dispIncr(t->t+dt)")
+    dispIncr.zero()
     if 0 == self._stepCount:
       for constraint in self.constraints:
         constraint.setField(t+dt, dispIncr)
@@ -183,7 +184,6 @@
     Advance to next time step.
     """
     dispIncr = self.fields.get("dispIncr(t->t+dt)")
-    dispIncr.zero()
 
     self._reformResidual(t+dt, dt)