[cig-commits] r8023 - in short/3D/PyLith/trunk: . libsrc/feassemble

[knepley at geodynamics.org]

Author: knepley
Date: 2007-09-25 10:19:45 -0700 (Tue, 25 Sep 2007)
New Revision: 8023

Modified:
   short/3D/PyLith/trunk/TODO
   short/3D/PyLith/trunk/libsrc/feassemble/Quadrature.cc
Log:
Share layout for precomputation sections


Modified: short/3D/PyLith/trunk/TODO
===================================================================
--- short/3D/PyLith/trunk/TODO	2007-09-25 15:37:00 UTC (rev 8022)
+++ short/3D/PyLith/trunk/TODO	2007-09-25 17:19:45 UTC (rev 8023)
@@ -205,6 +205,12 @@
 
   d. Same as a) for update
 
+  e. Put in flop counting for IR
+
+8. Memory fixes
+
+  a. Share memory in precomputation sections
+
 ======================================================================
 QUESTIONS FOR LEIF
 ======================================================================

Modified: short/3D/PyLith/trunk/libsrc/feassemble/Quadrature.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/feassemble/Quadrature.cc	2007-09-25 15:37:00 UTC (rev 8022)
+++ short/3D/PyLith/trunk/libsrc/feassemble/Quadrature.cc	2007-09-25 17:19:45 UTC (rev 8023)
@@ -278,12 +278,16 @@
 
   _quadPtsPre->setFiberDimension(cells, _numQuadPts * _spaceDim);
   _quadPtsPre->allocatePoint();
+  _jacobianPre->getAtlas()->setAtlas(_quadPtsPre->getAtlas()->getAtlas());
   _jacobianPre->setFiberDimension(cells, _numQuadPts * _cellDim * _spaceDim);
   _jacobianPre->allocatePoint();
+  _jacobianDetPre->getAtlas()->setAtlas(_quadPtsPre->getAtlas()->getAtlas());
   _jacobianDetPre->setFiberDimension(cells, _numQuadPts);
   _jacobianDetPre->allocatePoint();
+  _jacobianInvPre->setAtlas(_jacobianPre->getAtlas());
   _jacobianInvPre->setFiberDimension(cells, _numQuadPts * _cellDim * _spaceDim);
   _jacobianInvPre->allocatePoint();
+  _basisDerivPre->getAtlas()->setAtlas(_quadPtsPre->getAtlas()->getAtlas());
   _basisDerivPre->setFiberDimension(cells, _numQuadPts * _numBasis * _spaceDim);
   _basisDerivPre->allocatePoint();
   for(Mesh::label_sequence::iterator c_iter = cells->begin(); c_iter != end; ++c_iter) {