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

[brad at geodynamics.org]

Author: brad
Date: 2009-05-06 08:51:00 -0700 (Wed, 06 May 2009)
New Revision: 14884

Modified:
   short/3D/PyLith/trunk/libsrc/feassemble/QuadratureRefCell.cc
Log:
Fixed bug for 0-D quadrature. Can't use dimensions of arrays to set cellDim because cellDim is 0; assume if number of basis functions is 1, then cellDim is 0.

Modified: short/3D/PyLith/trunk/libsrc/feassemble/QuadratureRefCell.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/feassemble/QuadratureRefCell.cc	2009-05-06 06:10:56 UTC (rev 14883)
+++ short/3D/PyLith/trunk/libsrc/feassemble/QuadratureRefCell.cc	2009-05-06 15:51:00 UTC (rev 14884)
@@ -78,13 +78,13 @@
 { // initialize
   const int numQuadPts = numQuadPts1;
   const int numBasis = numBasis1;
-  const int cellDim = cellDim2;
+  const int cellDim = (numBasis != 1) ? cellDim2 : 0;
 
   assert(numQuadPts == numQuadPts2);
   assert(numQuadPts == numQuadPts3);
   assert(numQuadPts == numQuadPts4);
   assert(numBasis == numBasis2);
-  assert(cellDim == cellDim3);
+  assert(cellDim2 == cellDim3);
 
   if (0 == basis ||
       0 == basisDerivRef ||