[cig-commits] r16929 - short/3D/PyLith/trunk/libsrc/problems

[brad at geodynamics.org]

Author: brad
Date: 2010-06-08 07:48:18 -0700 (Tue, 08 Jun 2010)
New Revision: 16929

Modified:
   short/3D/PyLith/trunk/libsrc/problems/Explicit.cc
   short/3D/PyLith/trunk/libsrc/problems/Implicit.cc
Log:
Fixed bug in which velocity at constrained DOF were not updated (updatePoint -> updatePointAll).

Modified: short/3D/PyLith/trunk/libsrc/problems/Explicit.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/problems/Explicit.cc	2010-06-08 09:26:09 UTC (rev 16928)
+++ short/3D/PyLith/trunk/libsrc/problems/Explicit.cc	2010-06-08 14:48:18 UTC (rev 16929)
@@ -105,10 +105,10 @@
     } // for
     
     assert(velSection->getFiberDimension(*v_iter) == spaceDim);
-    velSection->updatePoint(*v_iter, &velVertex[0]);
+    velSection->updatePointAll(*v_iter, &velVertex[0]);
 
     assert(accSection->getFiberDimension(*v_iter) == spaceDim);
-    accSection->updatePoint(*v_iter, &accVertex[0]);
+    accSection->updatePointAll(*v_iter, &accVertex[0]);
   } // for
 
   PetscLogFlops(vertices->size() * 6*spaceDim);

Modified: short/3D/PyLith/trunk/libsrc/problems/Implicit.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/problems/Implicit.cc	2010-06-08 09:26:09 UTC (rev 16928)
+++ short/3D/PyLith/trunk/libsrc/problems/Implicit.cc	2010-06-08 14:48:18 UTC (rev 16929)
@@ -76,7 +76,7 @@
     velVertex = dispIncrVertex / dt;
     
     assert(velSection->getFiberDimension(*v_iter) == spaceDim);
-    velSection->updatePoint(*v_iter, &velVertex[0]);
+    velSection->updatePointAll(*v_iter, &velVertex[0]);
   } // for
   PetscLogFlops(vertices->size() * spaceDim);