[cig-commits] r19728 - short/3D/PyLith/trunk/unittests/libtests/materials

[brad at geodynamics.org]

Author: brad
Date: 2012-03-06 10:48:30 -0800 (Tue, 06 Mar 2012)
New Revision: 19728

Modified:
   short/3D/PyLith/trunk/unittests/libtests/materials/TestElasticMaterial.cc
Log:
Improved test condition for elastic constants by accounting for scale.

Modified: short/3D/PyLith/trunk/unittests/libtests/materials/TestElasticMaterial.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/materials/TestElasticMaterial.cc	2012-03-05 22:39:06 UTC (rev 19727)
+++ short/3D/PyLith/trunk/unittests/libtests/materials/TestElasticMaterial.cc	2012-03-06 18:48:30 UTC (rev 19728)
@@ -564,12 +564,14 @@
     
     const PylithScalar tolerance = (sizeof(double) == sizeof(PylithScalar)) ? 1.0e-06 : 1.0e-05;
     for (int i=0; i < numConsts; ++i)
-      if (fabs(elasticConstsE[i]) > tolerance)
+      if (fabs(elasticConstsE[i]) > tolerance) {
 	CPPUNIT_ASSERT_DOUBLES_EQUAL(1.0, elasticConsts[i]/elasticConstsE[i], 
 				     tolerance);
-      else
+      } else {
+	const double stressScale = 1.0e+9;
 	CPPUNIT_ASSERT_DOUBLES_EQUAL(elasticConstsE[i], elasticConsts[i],
-				     tolerance);
+				     tolerance*stressScale);
+      } // if/else
   } // for
 } // _testCalcElasticConsts