[cig-commits] r7532 - short/3D/PyLith/trunk/libsrc/feassemble

[knepley at geodynamics.org]

Author: knepley
Date: 2007-06-27 09:50:29 -0700 (Wed, 27 Jun 2007)
New Revision: 7532

Modified:
   short/3D/PyLith/trunk/libsrc/feassemble/ElasticityImplicit.cc
Log:
Complete residual at the end of integrateResidual()


Modified: short/3D/PyLith/trunk/libsrc/feassemble/ElasticityImplicit.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/feassemble/ElasticityImplicit.cc	2007-06-27 14:55:41 UTC (rev 7531)
+++ short/3D/PyLith/trunk/libsrc/feassemble/ElasticityImplicit.cc	2007-06-27 16:50:29 UTC (rev 7532)
@@ -223,9 +223,16 @@
       _material->calcStress(totalStrain);
     CALL_MEMBER_FN(*this, elasticityResidualFn)(stress);
 
+#if 0
+    std::cout << "Updating residual for cell " << *c_iter << std::endl;
+    for(int i = 0; i < _quadrature->spaceDim() * _quadrature->numBasis(); ++i) {
+      std::cout << "  v["<<i<<"]: " << _cellVector[i] << std::endl;
+    }
+#endif
     // Assemble cell contribution into field
     mesh->updateAdd(residual, *c_iter, _cellVector);
   } // for
+  ALE::Distribution<Mesh>::completeSection(mesh, residual);
 } // integrateResidual