[cig-commits] r7239 - in short/3D/PyLith/trunk: . unittests/libtests/feassemble

[brad at geodynamics.org]

Author: brad
Date: 2007-06-14 11:44:35 -0700 (Thu, 14 Jun 2007)
New Revision: 7239

Modified:
   short/3D/PyLith/trunk/TODO
   short/3D/PyLith/trunk/unittests/libtests/feassemble/TestQuadrature.cc
Log:
Added testing of basisDeriv (derivatives with respect to global coordinates) to unit test for Quadrature::computeGeometry().

Modified: short/3D/PyLith/trunk/TODO
===================================================================
--- short/3D/PyLith/trunk/TODO	2007-06-14 18:22:44 UTC (rev 7238)
+++ short/3D/PyLith/trunk/TODO	2007-06-14 18:44:35 UTC (rev 7239)
@@ -2,8 +2,6 @@
 MAIN PRIORITIES (Brad)
 ======================================================================
 
-Create unit tests for Quadrature::basisDeriv().
-
 FaultCohesiveKin::integrateJacobian()
   Need to prevent multiple contributions of orientation information for
   each vertex, because we want to use ADD_VALUES in updateOperator().
@@ -24,7 +22,7 @@
 
   a. FaultCohesive
     i. Add checking of faultMesh [not currently used]
-    ii. Tests for cases with multiple cohesive cells
+    ii. Tests cases with multiple cohesive cells
 
 examples/templates
   containers

Modified: short/3D/PyLith/trunk/unittests/libtests/feassemble/TestQuadrature.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/feassemble/TestQuadrature.cc	2007-06-14 18:22:44 UTC (rev 7238)
+++ short/3D/PyLith/trunk/unittests/libtests/feassemble/TestQuadrature.cc	2007-06-14 18:44:35 UTC (rev 7239)
@@ -263,6 +263,7 @@
   const double* verticesRef = data.verticesRef;
   const double* basis = data.basis;
   const double* basisDerivRef = data.basisDerivRef;
+  const double* basisDeriv = data.basisDeriv;
   const double* quadPtsRef = data.quadPtsRef;
   const double* quadWts = data.quadWts;
 
@@ -324,6 +325,12 @@
   for (int i=0; i < size; ++i)
     CPPUNIT_ASSERT_DOUBLES_EQUAL(jacobianDet[i], pQuad->_jacobianDet[i], 
 				 tolerance);
+
+  size = numQuadPts * numBasis * spaceDim;
+  for (int i=0; i < size; ++i)
+    CPPUNIT_ASSERT_DOUBLES_EQUAL(basisDeriv[i], pQuad->_basisDeriv[i], 
+				 tolerance);
+
 } // testComputeGeometry