[cig-commits] r8378 - short/3D/PyLith/trunk/pylith/feassemble

[luis at geodynamics.org]

Author: luis
Date: 2007-12-04 16:09:48 -0800 (Tue, 04 Dec 2007)
New Revision: 8378

Modified:
   short/3D/PyLith/trunk/pylith/feassemble/FIATLagrange.py
Log:
Found missing index increment in FIATLagrange tabulation loops (for case dim == 3)


Modified: short/3D/PyLith/trunk/pylith/feassemble/FIATLagrange.py
===================================================================
--- short/3D/PyLith/trunk/pylith/feassemble/FIATLagrange.py	2007-12-05 00:06:32 UTC (rev 8377)
+++ short/3D/PyLith/trunk/pylith/feassemble/FIATLagrange.py	2007-12-05 00:09:48 UTC (rev 8378)
@@ -626,7 +626,8 @@
                       self.basisDeriv[n][m][2] = basis[r][g]*basis[q][f]*basisDeriv[s][h][0]
                       m += 1
                 if not m == self.numCorners: raise RuntimeError('Invalid 3D function tabulation')
-          if not n == self.numQuadPts: raise RuntimeError('Invalid 2D quadrature')
+                n += 1
+          if not n == self.numQuadPts: raise RuntimeError('Invalid 3D quadrature')
         self.vertices = numpy.reshape(self.vertices, (self.numCorners, dim))
         self.quadPts = numpy.reshape(self.quadPts, (self.numQuadPts, dim))
         self.quadWts = numpy.reshape(self.quadWts, (self.numQuadPts))