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

[knepley at geodynamics.org]

Author: knepley
Date: 2007-03-31 16:54:27 -0700 (Sat, 31 Mar 2007)
New Revision: 6496

Modified:
   short/3D/PyLith/trunk/libsrc/feassemble/ExplicitElasticity.cc
Log:
Fixed mass matrix assembly


Modified: short/3D/PyLith/trunk/libsrc/feassemble/ExplicitElasticity.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/feassemble/ExplicitElasticity.cc	2007-03-31 23:48:14 UTC (rev 6495)
+++ short/3D/PyLith/trunk/libsrc/feassemble/ExplicitElasticity.cc	2007-03-31 23:54:27 UTC (rev 6496)
@@ -338,13 +338,14 @@
       const double wt = 
 	quadWts[iQuad] * jacobianDet[iQuad] * density[iQuad] / dt2;
       for (int iBasis=0, iQ=iQuad*numBasis; iBasis < numBasis; ++iBasis) {
-        const int iBlock = (iBasis * spaceDim) * (spaceDim * numBasis);
         const double valI = wt*basis[iQ+iBasis];
         for (int jBasis=0; jBasis < numBasis; ++jBasis) {
-          const int jBlock = (jBasis * spaceDim);
           const double valIJ = valI * basis[iQ+jBasis];
-          for (int iDim=0; iDim < spaceDim; ++iDim)
-            _cellMatrix[iBlock+jBlock+iDim*spaceDim+iDim] += valIJ;
+          for (int iDim=0; iDim < spaceDim; ++iDim) {
+            const int iBlock = (iBasis * spaceDim + iDim) * (spaceDim * numBasis);
+            const int jBlock = (jBasis * spaceDim + iDim);
+            _cellMatrix[iBlock+jBlock] += valIJ;
+          } // for
         } // for
       } // for
     } // for