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

[willic3 at geodynamics.org]

Author: willic3
Date: 2009-06-12 19:09:26 -0700 (Fri, 12 Jun 2009)
New Revision: 15225

Modified:
   short/3D/PyLith/trunk/unittests/libtests/materials/data/PowerLaw3DTimeDep.py
   short/3D/PyLith/trunk/unittests/libtests/materials/data/PowerLaw3DTimeDepData.cc
Log:
Fixed minor probem to make sure correct strains are being used when
computing elastic constants.



Modified: short/3D/PyLith/trunk/unittests/libtests/materials/data/PowerLaw3DTimeDep.py
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/materials/data/PowerLaw3DTimeDep.py	2009-06-13 01:52:36 UTC (rev 15224)
+++ short/3D/PyLith/trunk/unittests/libtests/materials/data/PowerLaw3DTimeDep.py	2009-06-13 02:09:26 UTC (rev 15225)
@@ -294,8 +294,9 @@
     Function to compute a particular stress component as a function of a
     strain component.
     """
-    strainTpdt[strainComp] = strainVal
-    stressTpdt, visStrainTpdt = self._computeStress(strainTpdt, muV, lambdaV,
+    strainTest = numpy.array(strainTpdt, dtype=numpy.float64)
+    strainTest[strainComp] = strainVal
+    stressTpdt, visStrainTpdt = self._computeStress(strainTest, muV, lambdaV,
                                                     viscosityCoeffV,
                                                     powerLawExpV, visStrainT,
                                                     stressT,

Modified: short/3D/PyLith/trunk/unittests/libtests/materials/data/PowerLaw3DTimeDepData.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/materials/data/PowerLaw3DTimeDepData.cc	2009-06-13 01:52:36 UTC (rev 15224)
+++ short/3D/PyLith/trunk/unittests/libtests/materials/data/PowerLaw3DTimeDepData.cc	2009-06-13 02:09:26 UTC (rev 15225)
@@ -223,28 +223,28 @@
 const double pylith::materials::PowerLaw3DTimeDepData::_elasticConsts[] = {
   6.74326513e+10,
   2.25336738e+10,
-  2.25336747e+10,
+  2.25336738e+10,
   0.00000000e+00,
   0.00000000e+00,
   0.00000000e+00,
   6.74326513e+10,
-  2.25336747e+10,
+  2.25336738e+10,
   0.00000000e+00,
   0.00000000e+00,
   0.00000000e+00,
-  6.74326513e+10,
+  6.74326522e+10,
   0.00000000e+00,
   0.00000000e+00,
   0.00000000e+00,
   4.48989770e+10,
   0.00000000e+00,
   0.00000000e+00,
-  4.48989775e+10,
+  4.48989770e+10,
   0.00000000e+00,
   4.48989770e+10,
   8.64000013e+09,
   2.88000004e+09,
-  2.87999958e+09,
+  2.88000004e+09,
   0.00000000e+00,
   0.00000000e+00,
   0.00000000e+00,
@@ -257,12 +257,12 @@
   0.00000000e+00,
   0.00000000e+00,
   0.00000000e+00,
-  5.75999985e+09,
+  5.76000009e+09,
   0.00000000e+00,
   0.00000000e+00,
   5.75999985e+09,
   0.00000000e+00,
-  5.76000009e+09,
+  5.75999985e+09,
 };
 
 const double pylith::materials::PowerLaw3DTimeDepData::_initialStress[] = {