[cig-commits] r7266 - in short/3D/PyLith/trunk: . libsrc/feassemble
unittests/libtests/feassemble/data
brad at geodynamics.org
brad at geodynamics.org
Fri Jun 15 13:11:09 PDT 2007
Author: brad
Date: 2007-06-15 13:11:09 -0700 (Fri, 15 Jun 2007)
New Revision: 7266
Modified:
short/3D/PyLith/trunk/TODO
short/3D/PyLith/trunk/libsrc/feassemble/ElasticityExplicit.cc
short/3D/PyLith/trunk/libsrc/feassemble/ElasticityImplicit.cc
short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityExplicitData2DLinear.cc
short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityExplicitData2DQuadratic.cc
short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityExplicitData3DLinear.cc
short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityExplicitData3DQuadratic.cc
short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityImplicitData2DLinear.cc
short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityImplicitData2DQuadratic.cc
short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityImplicitData3DLinear.cc
short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityImplicitData3DQuadratic.cc
short/3D/PyLith/trunk/unittests/libtests/feassemble/data/IntegratorElasticity.py
Log:
Fixed bug in integrating Jacobian and residual for implicit time stepping and integrating residual for explicit time stepping. B^T D B must have 1/2 account for shear in D (not in B).
Modified: short/3D/PyLith/trunk/TODO
===================================================================
--- short/3D/PyLith/trunk/TODO 2007-06-15 20:10:58 UTC (rev 7265)
+++ short/3D/PyLith/trunk/TODO 2007-06-15 20:11:09 UTC (rev 7266)
@@ -4,6 +4,8 @@
creation of cohesive cells for hex8 meshes
tests/1d/line3/dislocation.cfg [CHACO]
tests/3d/tet4/shearelasticstrain.cfg [SIEVE INDEXING AND CONSTRAINTS]
+ Note: ordering the Dirichlet BC so that the constrainted DOF are in
+ order is a workaround.
VTK output
zero out values at Lagrange constraint vertices
@@ -13,6 +15,8 @@
MAIN PRIORITIES (Brad)
======================================================================
+Add catching ALE exceptions in modules.
+
1. Additional unit tests
Unit tests with multiple cohesive cells.
@@ -29,8 +33,6 @@
containers
verify containers work
-Check trapping of errors in reading spatialdata files.
-
======================================================================
SECONDARY PRIORITIES
======================================================================
Modified: short/3D/PyLith/trunk/libsrc/feassemble/ElasticityExplicit.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/feassemble/ElasticityExplicit.cc 2007-06-15 20:10:58 UTC (rev 7265)
+++ short/3D/PyLith/trunk/libsrc/feassemble/ElasticityExplicit.cc 2007-06-15 20:11:09 UTC (rev 7266)
@@ -243,8 +243,8 @@
++iBasis) {
const double N1 = wt*basisDeriv[iQ+iBasis*cellDim ];
const double N2 = wt*basisDeriv[iQ+iBasis*cellDim+1];
- _cellVector[iBasis*spaceDim ] -= N1*s11 + 0.5*N2*s12;
- _cellVector[iBasis*spaceDim+1] -= 0.5*N1*s12 + N2*s22;
+ _cellVector[iBasis*spaceDim ] -= N1*s11 + N2*s12;
+ _cellVector[iBasis*spaceDim+1] -= N1*s12 + N2*s22;
} // for
} // for
err = PetscLogFlops(numQuadPts*(1+numBasis*(8+2+9)));
@@ -275,9 +275,9 @@
const double N1 = wt*basisDeriv[iQ+iBasis*cellDim+0];
const double N2 = wt*basisDeriv[iQ+iBasis*cellDim+1];
const double N3 = wt*basisDeriv[iQ+iBasis*cellDim+2];
- _cellVector[iBlock ] -= N1*s11 + 0.5*N2*s12 + 0.5*N3*s13;
- _cellVector[iBlock+1] -= 0.5*N1*s12 + N2*s22 + 0.5*N3*s23;
- _cellVector[iBlock+2] -= 0.5*N1*s13 + 0.5*N2*s23 + N3*s33;
+ _cellVector[iBlock ] -= N1*s11 + N2*s12 + N3*s13;
+ _cellVector[iBlock+1] -= N1*s12 + N2*s22 + N3*s23;
+ _cellVector[iBlock+2] -= N1*s13 + N2*s23 + N3*s33;
} // for
} // for
err = PetscLogFlops(numQuadPts*(1+numBasis*(3+12)));
Modified: short/3D/PyLith/trunk/libsrc/feassemble/ElasticityImplicit.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/feassemble/ElasticityImplicit.cc 2007-06-15 20:10:58 UTC (rev 7265)
+++ short/3D/PyLith/trunk/libsrc/feassemble/ElasticityImplicit.cc 2007-06-15 20:11:09 UTC (rev 7266)
@@ -227,8 +227,8 @@
++iBasis) {
const double N1 = wt*basisDeriv[iQ+iBasis*cellDim ];
const double N2 = wt*basisDeriv[iQ+iBasis*cellDim+1];
- _cellVector[iBasis*spaceDim ] -= N1*s11 + 0.5*N2*s12;
- _cellVector[iBasis*spaceDim+1] -= 0.5*N1*s12 + N2*s22;
+ _cellVector[iBasis*spaceDim ] -= N1*s11 + N2*s12;
+ _cellVector[iBasis*spaceDim+1] -= N1*s12 + N2*s22;
} // for
} // for
PetscErrorCode err = PetscLogFlops(numQuadPts*(1+numBasis*(8+2+9)));
@@ -258,9 +258,9 @@
const double N1 = wt*basisDeriv[iQ+iBasis*cellDim+0];
const double N2 = wt*basisDeriv[iQ+iBasis*cellDim+1];
const double N3 = wt*basisDeriv[iQ+iBasis*cellDim+2];
- _cellVector[iBlock ] -= N1*s11 + 0.5*N2*s12 + 0.5*N3*s13;
- _cellVector[iBlock+1] -= 0.5*N1*s12 + N2*s22 + 0.5*N3*s23;
- _cellVector[iBlock+2] -= 0.5*N1*s13 + 0.5*N2*s23 + N3*s33;
+ _cellVector[iBlock ] -= N1*s11 + N2*s12 + N3*s13;
+ _cellVector[iBlock+1] -= N1*s12 + N2*s22 + N3*s23;
+ _cellVector[iBlock+2] -= N1*s13 + N2*s23 + N3*s33;
} // for
} // for
PetscErrorCode err = PetscLogFlops(numQuadPts*(1+numBasis*(3+12)));
@@ -409,7 +409,7 @@
const double C1112 = elasticConsts[iQuad][2]/2.0;
const double C2222 = elasticConsts[iQuad][3];
const double C2212 = elasticConsts[iQuad][4]/2.0;
- const double C1212 = elasticConsts[iQuad][5]/4.0;
+ const double C1212 = elasticConsts[iQuad][5]/2.0;
for (int iBasis=0, iQ=iQuad*numBasis*cellDim;
iBasis < numBasis;
++iBasis) {
@@ -471,12 +471,12 @@
const double C3312 = elasticConsts[iQuad][12]/2.0;
const double C3323 = elasticConsts[iQuad][13]/2.0;
const double C3313 = elasticConsts[iQuad][14]/2.0;
- const double C1212 = elasticConsts[iQuad][15]/4.0;
- const double C1223 = elasticConsts[iQuad][16]/4.0;
- const double C1213 = elasticConsts[iQuad][17]/4.0;
- const double C2323 = elasticConsts[iQuad][18]/4.0;
- const double C2313 = elasticConsts[iQuad][19]/4.0;
- const double C1313 = elasticConsts[iQuad][20]/4.0;
+ const double C1212 = elasticConsts[iQuad][15]/2.0;
+ const double C1223 = elasticConsts[iQuad][16]/2.0;
+ const double C1213 = elasticConsts[iQuad][17]/2.0;
+ const double C2323 = elasticConsts[iQuad][18]/2.0;
+ const double C2313 = elasticConsts[iQuad][19]/2.0;
+ const double C1313 = elasticConsts[iQuad][20]/2.0;
for (int iBasis=0, iQ=iQuad*numBasis*cellDim;
iBasis < numBasis;
++iBasis) {
Modified: short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityExplicitData2DLinear.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityExplicitData2DLinear.cc 2007-06-15 20:10:58 UTC (rev 7265)
+++ short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityExplicitData2DLinear.cc 2007-06-15 20:11:09 UTC (rev 7266)
@@ -90,9 +90,9 @@
};
const double pylith::feassemble::ElasticityExplicitData2DLinear::_valsResidual[] = {
- -6.32656035e+10, 1.22513006e+11,
- -1.32542399e+10, -1.10498357e+11,
- 7.65207601e+10, -1.20233573e+10,
+ -6.72542399e+10, 1.24660734e+11,
+ -9.57242172e+09, -1.09884721e+11,
+ 7.68275783e+10, -1.47847210e+10,
};
const double pylith::feassemble::ElasticityExplicitData2DLinear::_valsJacobian[] = {
Modified: short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityExplicitData2DQuadratic.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityExplicitData2DQuadratic.cc 2007-06-15 20:10:58 UTC (rev 7265)
+++ short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityExplicitData2DQuadratic.cc 2007-06-15 20:11:09 UTC (rev 7266)
@@ -130,12 +130,12 @@
};
const double pylith::feassemble::ElasticityExplicitData2DQuadratic::_valsResidual[] = {
- 8.94394715e+08, 8.85261110e+09,
- -2.12714501e+10, 9.02616504e+09,
- 2.84000005e+08, -3.31752565e+10,
- 2.76103121e+09, 6.57610682e+10,
- -5.08876156e+09, -8.47107565e+10,
- 2.24330805e+10, 3.42541857e+10,
+ 2.04455943e+09, 9.48741748e+09,
+ -1.86265045e+10, 3.65998688e+09,
+ -2.49512114e+09, -3.26311648e+10,
+ 1.39608918e+10, 6.18995066e+10,
+ -1.90442781e+10, -8.28253071e+10,
+ 2.41727473e+10, 4.04175780e+10,
};
const double pylith::feassemble::ElasticityExplicitData2DQuadratic::_valsJacobian[] = {
Modified: short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityExplicitData3DLinear.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityExplicitData3DLinear.cc 2007-06-15 20:10:58 UTC (rev 7265)
+++ short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityExplicitData3DLinear.cc 2007-06-15 20:11:09 UTC (rev 7266)
@@ -96,10 +96,10 @@
};
const double pylith::feassemble::ElasticityExplicitData3DLinear::_valsResidual[] = {
- -3.17161475e+09, 3.76142900e+10, 1.81833517e+09,
- -1.69858191e+09, 6.63279672e+10, 1.45939393e+10,
- 3.16565291e+09, -1.08396398e+11, -2.23578730e+10,
- 1.70891499e+09, 4.45615803e+09, 5.94458982e+09,
+ -6.07523928e+09, 3.62541599e+10, 3.19640209e+09,
+ -4.00838674e+09, 6.66858204e+10, 2.19195434e+10,
+ 6.67438897e+09, -1.05592494e+11, -3.14998804e+10,
+ 3.41360830e+09, 2.65453164e+09, 6.38292625e+09,
};
const double pylith::feassemble::ElasticityExplicitData3DLinear::_valsJacobian[] = {
Modified: short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityExplicitData3DQuadratic.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityExplicitData3DQuadratic.cc 2007-06-15 20:10:58 UTC (rev 7265)
+++ short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityExplicitData3DQuadratic.cc 2007-06-15 20:11:09 UTC (rev 7266)
@@ -177,16 +177,16 @@
};
const double pylith::feassemble::ElasticityExplicitData3DQuadratic::_valsResidual[] = {
- 8.00450011e+10, -9.00472985e+09, 1.61541766e+11,
- -3.54088543e+10, -6.47782824e+10, 2.41769691e+10,
- -1.55663344e+11, -4.18485549e+11, 2.08248421e+11,
- -9.98985902e+09, -4.26541080e+10, 3.19516877e+11,
- -5.76612056e+10, -2.26094787e+11, -1.40838144e+11,
- 1.37516193e+11, 4.40486089e+11, -1.93421106e+11,
- -8.06475029e+11, -3.08666160e+11, -1.91097837e+11,
- 8.66282546e+10, 4.68394314e+11, -2.59144518e+11,
- 1.66627489e+11, 6.19155802e+11, -4.59691032e+11,
- 5.94347823e+11, -4.58318391e+11, 5.30727623e+11,
+ 1.21810568e+11, -3.29478688e+10, 2.43738946e+11,
+ -2.90354241e+10, -7.34518237e+10, 2.34529293e+10,
+ -2.49211185e+11, -4.63822584e+11, 2.54633336e+11,
+ 3.54901735e+10, -5.62151785e+10, 3.13375453e+11,
+ -1.24954972e+11, -1.83164844e+11, -1.74376488e+11,
+ 2.33420251e+11, 5.42368499e+11, -3.17019958e+11,
+ -9.29980081e+11, -4.08320625e+11, -2.12610662e+11,
+ 1.17953752e+11, 4.72743819e+11, -2.44313374e+11,
+ 2.27857411e+11, 6.71949765e+11, -5.10996159e+11,
+ 5.96615974e+11, -4.69104960e+11, 6.24134995e+11,
};
const double pylith::feassemble::ElasticityExplicitData3DQuadratic::_valsJacobian[] = {
Modified: short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityImplicitData2DLinear.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityImplicitData2DLinear.cc 2007-06-15 20:10:58 UTC (rev 7265)
+++ short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityImplicitData2DLinear.cc 2007-06-15 20:11:09 UTC (rev 7266)
@@ -90,30 +90,30 @@
};
const double pylith::feassemble::ElasticityImplicitData2DLinear::_valsResidual[] = {
- -4.40727273e+10, 1.63800000e+11,
- -1.81636364e+10, -1.54200000e+11,
- 6.22363636e+10, -9.60000000e+09,
+ -3.94636364e+10, 1.61318182e+11,
+ -2.24181818e+10, -1.54909091e+11,
+ 6.18818182e+10, -6.40909091e+09,
};
const double pylith::feassemble::ElasticityImplicitData2DLinear::_valsJacobian[] = {
- 3.15681818e+10, -1.86136364e+10,
- -4.90909091e+09, 9.68181818e+09,
- -2.66590909e+10, 8.93181818e+09,
- -1.86136364e+10, 7.24772727e+10,
- 2.18181818e+09, -6.28636364e+10,
- 1.64318182e+10, -9.61363636e+09,
- -4.90909091e+09, 2.18181818e+09,
- 1.14545455e+10, 4.90909091e+09,
- -6.54545455e+09, -7.09090909e+09,
- 9.68181818e+09, -6.28636364e+10,
- 4.90909091e+09, 5.91818182e+10,
- -1.45909091e+10, 3.68181818e+09,
- -2.66590909e+10, 1.64318182e+10,
- -6.54545455e+09, -1.45909091e+10,
- 3.32045455e+10, -1.84090909e+09,
- 8.93181818e+09, -9.61363636e+09,
- -7.09090909e+09, 3.68181818e+09,
- -1.84090909e+09, 5.93181818e+09,
+ 4.30909091e+10, -2.48181818e+10,
+ -1.55454545e+10, 7.90909091e+09,
+ -2.75454545e+10, 1.69090909e+10,
+ -2.48181818e+10, 7.58181818e+10,
+ 7.90909091e+09, -6.19090909e+10,
+ 1.69090909e+10, -1.39090909e+10,
+ -1.55454545e+10, 7.90909091e+09,
+ 2.12727273e+10, 6.54545455e+09,
+ -5.72727273e+09, -1.44545455e+10,
+ 7.90909091e+09, -6.19090909e+10,
+ 6.54545455e+09, 5.94545455e+10,
+ -1.44545455e+10, 2.45454545e+09,
+ -2.75454545e+10, 1.69090909e+10,
+ -5.72727273e+09, -1.44545455e+10,
+ 3.32727273e+10, -2.45454545e+09,
+ 1.69090909e+10, -1.39090909e+10,
+ -1.44545455e+10, 2.45454545e+09,
+ -2.45454545e+09, 1.14545455e+10,
};
pylith::feassemble::ElasticityImplicitData2DLinear::ElasticityImplicitData2DLinear(void)
Modified: short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityImplicitData2DQuadratic.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityImplicitData2DQuadratic.cc 2007-06-15 20:10:58 UTC (rev 7265)
+++ short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityImplicitData2DQuadratic.cc 2007-06-15 20:11:09 UTC (rev 7266)
@@ -130,87 +130,87 @@
};
const double pylith::feassemble::ElasticityImplicitData2DQuadratic::_valsResidual[] = {
- -4.04084094e+09, -4.29933897e+10,
- -4.13127930e+10, 2.79255564e+10,
- 8.98193790e+08, -1.29686671e+11,
- 2.07693369e+10, 1.58734365e+11,
- -2.29487522e+10, -2.12565853e+11,
- 4.66348554e+10, 1.98585991e+11,
+ -7.03759529e+08, -4.10320114e+10,
+ -3.70610450e+10, 2.11260736e+10,
+ -1.31513835e+10, -1.25323831e+11,
+ 4.33964587e+10, 1.47829671e+11,
+ -5.10800854e+10, -2.16829139e+11,
+ 5.85998147e+10, 2.14229238e+11,
};
const double pylith::feassemble::ElasticityImplicitData2DQuadratic::_valsJacobian[] = {
- 1.92724753e+10, 1.13862332e+10,
- 4.42819754e+09, 2.19495425e+09,
- 2.80230421e+09, 3.73249461e+08,
- -1.92262756e+10, -8.64560662e+09,
- -4.36542846e+09, -2.08992316e+09,
- -2.91127292e+09, -3.21890717e+09,
- 1.13862332e+10, 7.19847471e+10,
- 4.69495425e+09, -1.54251684e+09,
- -2.12675054e+09, 2.75320212e+10,
- -1.86456066e+10, -8.12569923e+08,
- -2.08992316e+09, -6.25120411e+09,
- 6.78109283e+09, -9.09104775e+10,
- 4.42819754e+09, 4.69495425e+09,
- 5.38342257e+10, -1.02511036e+10,
- 1.09718676e+10, -6.50919571e+09,
- -1.30974352e+09, -2.10712128e+10,
- -5.55692740e+10, 2.85686411e+10,
- -1.23552734e+10, 4.56791677e+09,
- 2.19495425e+09, -1.54251684e+09,
- -1.02511036e+10, 1.72387553e+10,
- -4.00919571e+09, 5.80929253e+09,
- -1.10712128e+10, 2.74018801e+10,
- 1.85686411e+10, -4.36343981e+10,
- 4.56791677e+09, -5.27301294e+09,
- 2.80230421e+09, -2.12675054e+09,
- 1.09718676e+10, -4.00919571e+09,
- 2.66480946e+10, -1.60701978e+10,
- -2.19813953e+09, -5.62260428e+07,
- -2.26517201e+10, 1.04034074e+10,
- -1.55724068e+10, 1.18589627e+10,
- 3.73249461e+08, 2.75320212e+10,
- -6.50919571e+09, 5.80929253e+09,
- -1.60701978e+10, 8.90068119e+10,
- -5.62260428e+07, -1.36007171e+10,
- 2.04034074e+10, 1.08657906e+10,
- 1.85896275e+09, -1.19613199e+11,
- -1.92262756e+10, -1.86456066e+10,
- -1.30974352e+09, -1.10712128e+10,
- -2.19813953e+09, -5.62260428e+07,
- 7.51332944e+10, -1.31964480e+10,
- -9.80844073e+09, 1.73687559e+10,
- -4.25906950e+10, 2.56007375e+10,
- -8.64560662e+09, -8.12569923e+08,
- -2.10712128e+10, 2.74018801e+10,
- -5.62260428e+07, -1.36007171e+10,
- -1.31964480e+10, 1.19528505e+11,
- 1.73687559e+10, -1.34841265e+11,
- 2.56007375e+10, 2.32416756e+09,
- -4.36542846e+09, -2.08992316e+09,
- -5.55692740e+10, 1.85686411e+10,
- -2.26517201e+10, 2.04034074e+10,
- -9.80844073e+09, 1.73687559e+10,
- 1.00479896e+11, -2.19113697e+10,
- -8.08503301e+09, -3.23395114e+10,
- -2.08992316e+09, -6.25120411e+09,
- 2.85686411e+10, -4.36343981e+10,
- 1.04034074e+10, 1.08657906e+10,
- 1.73687559e+10, -1.34841265e+11,
- -2.19113697e+10, 1.82672151e+11,
- -3.23395114e+10, -8.81107383e+09,
- -2.91127292e+09, 6.78109283e+09,
- -1.23552734e+10, 4.56791677e+09,
- -1.55724068e+10, 1.85896275e+09,
- -4.25906950e+10, 2.56007375e+10,
- -8.08503301e+09, -3.23395114e+10,
- 8.15146812e+10, -6.46919843e+09,
- -3.21890717e+09, -9.09104775e+10,
- 4.56791677e+09, -5.27301294e+09,
- 1.18589627e+10, -1.19613199e+11,
- 2.56007375e+10, 2.32416756e+09,
- -3.23395114e+10, -8.81107383e+09,
- -6.46919843e+09, 2.22283596e+11,
+ 3.10620755e+10, 1.51816443e+10,
+ 4.03724615e+09, 4.59327234e+09,
+ 7.44201344e+09, -1.16900072e+09,
+ -1.88162512e+10, -1.81941422e+10,
+ -5.31233692e+09, -2.78656421e+09,
+ -1.84127470e+10, 2.37479043e+09,
+ 1.51816443e+10, 7.32318930e+10,
+ 4.59327234e+09, -7.39325357e+08,
+ -1.16900072e+09, 2.72257870e+10,
+ -1.81941422e+10, -4.08528660e+09,
+ -2.78656421e+09, -6.82095745e+09,
+ 2.37479043e+09, -8.88121106e+10,
+ 4.03724615e+09, 4.59327234e+09,
+ 5.52513202e+10, -1.36681382e+10,
+ 1.16542644e+10, -7.01226095e+09,
+ 3.42514288e+09, -2.14282837e+10,
+ -6.14617630e+10, 3.14248548e+10,
+ -1.29062107e+10, 6.09055570e+09,
+ 4.59327234e+09, -7.39325357e+08,
+ -1.36681382e+10, 2.59749439e+10,
+ -7.01226095e+09, 7.52420433e+09,
+ -2.14282837e+10, 2.63944418e+10,
+ 3.14248548e+10, -5.19138623e+10,
+ 6.09055570e+09, -7.24040229e+09,
+ 7.44201344e+09, -1.16900072e+09,
+ 1.16542644e+10, -7.01226095e+09,
+ 4.11450311e+10, -2.14269305e+10,
+ -4.46688705e+09, -7.49680571e+07,
+ -2.01418211e+10, 2.05378765e+10,
+ -3.56326007e+10, 9.14528366e+09,
+ -1.16900072e+09, 2.72257870e+10,
+ -7.01226095e+09, 7.52420433e+09,
+ -2.14269305e+10, 9.10320049e+10,
+ -7.49680571e+07, -1.35889491e+10,
+ 2.05378765e+10, 6.67218742e+09,
+ 9.14528366e+09, -1.18865235e+11,
+ -1.88162512e+10, -1.81941422e+10,
+ 3.42514288e+09, -2.14282837e+10,
+ -4.46688705e+09, -7.49680571e+07,
+ 9.34772297e+10, -1.75952640e+10,
+ -3.26438450e+10, 2.31583412e+10,
+ -4.09753893e+10, 3.41343167e+10,
+ -1.81941422e+10, -4.08528660e+09,
+ -2.14282837e+10, 2.63944418e+10,
+ -7.49680571e+07, -1.35889491e+10,
+ -1.75952640e+10, 1.28993398e+11,
+ 2.31583412e+10, -1.32670104e+11,
+ 3.41343167e+10, -5.04349923e+09,
+ -5.31233692e+09, -2.78656421e+09,
+ -6.14617630e+10, 3.14248548e+10,
+ -2.01418211e+10, 2.05378765e+10,
+ -3.26438450e+10, 2.31583412e+10,
+ 1.28924268e+11, -2.92151596e+10,
+ -9.36450187e+09, -4.31193486e+10,
+ -2.78656421e+09, -6.82095745e+09,
+ 3.14248548e+10, -5.19138623e+10,
+ 2.05378765e+10, 6.67218742e+09,
+ 2.31583412e+10, -1.32670104e+11,
+ -2.92151596e+10, 1.94678071e+11,
+ -4.31193486e+10, -9.94533453e+09,
+ -1.84127470e+10, 2.37479043e+09,
+ -1.29062107e+10, 6.09055570e+09,
+ -3.56326007e+10, 9.14528366e+09,
+ -4.09753893e+10, 3.41343167e+10,
+ -9.36450187e+09, -4.31193486e+10,
+ 1.17291450e+11, -8.62559790e+09,
+ 2.37479043e+09, -8.88121106e+10,
+ 6.09055570e+09, -7.24040229e+09,
+ 9.14528366e+09, -1.18865235e+11,
+ 3.41343167e+10, -5.04349923e+09,
+ -4.31193486e+10, -9.94533453e+09,
+ -8.62559790e+09, 2.29906581e+11,
};
pylith::feassemble::ElasticityImplicitData2DQuadratic::ElasticityImplicitData2DQuadratic(void)
Modified: short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityImplicitData3DLinear.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityImplicitData3DLinear.cc 2007-06-15 20:10:58 UTC (rev 7265)
+++ short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityImplicitData3DLinear.cc 2007-06-15 20:11:09 UTC (rev 7266)
@@ -96,61 +96,61 @@
};
const double pylith::feassemble::ElasticityImplicitData3DLinear::_valsResidual[] = {
- -7.70693928e+08, 3.56969021e+10, 2.39826518e+09,
- 9.38355019e+09, 5.66238848e+10, 2.52838755e+10,
- -1.24875155e+10, -9.54987299e+10, -3.89812887e+10,
- 3.87465923e+09, 3.17794300e+09, 1.12991481e+10,
+ 2.78190830e+07, 3.56885378e+10, 3.94225527e+09,
+ 1.52090149e+10, 5.28108736e+10, 3.64722119e+10,
+ -2.29679058e+10, -8.91228005e+10, -5.23343247e+10,
+ 7.73107187e+09, 6.23389095e+08, 1.19198575e+10,
};
const double pylith::feassemble::ElasticityImplicitData3DLinear::_valsJacobian[] = {
- 7.12856258e+09, 4.35594796e+09, -2.39776952e+08,
- -1.37871747e+09, 1.62314126e+09, -2.45631970e+09,
- -6.44021066e+09, -6.42843866e+09, 4.69330855e+09,
- 6.90365551e+08, 4.49349442e+08, -1.99721190e+09,
- 4.35594796e+09, 2.26675960e+10, -6.07806691e+08,
- -4.12685874e+09, 3.62653346e+10, -7.02881041e+09,
- -4.28438662e+08, -6.26112144e+10, 1.27342007e+10,
- 1.99349442e+08, 3.67828377e+09, -5.09758364e+09,
- -2.39776952e+08, -6.07806691e+08, 4.32032218e+09,
- -9.56319703e+08, -4.02881041e+09, 6.03531599e+09,
- 2.19330855e+09, 7.23420074e+09, -1.17763321e+10,
- -9.97211896e+08, -2.59758364e+09, 1.42069393e+09,
- -1.37871747e+09, -4.12685874e+09, -9.56319703e+08,
- 3.12511617e+10, -1.67181691e+10, 8.97072491e+09,
- -3.32634758e+10, 2.19307621e+10, -1.25604089e+10,
- 3.39103160e+09, -1.08573420e+09, 4.54600372e+09,
- 1.62314126e+09, 3.62653346e+10, -4.02881041e+09,
- -1.67181691e+10, 6.92530204e+10, -1.69753717e+10,
- 1.56807621e+10, -1.13835037e+11, 3.01143123e+10,
- -5.85734201e+08, 8.31668216e+09, -9.11013011e+09,
- -2.45631970e+09, -7.02881041e+09, 6.03531599e+09,
- 8.97072491e+09, -1.69753717e+10, 3.17077138e+10,
- -8.81040892e+09, 2.93643123e+10, -5.25594796e+10,
- 2.29600372e+09, -5.36013011e+09, 1.48164498e+10,
- -6.44021066e+09, -4.28438662e+08, 2.19330855e+09,
- -3.32634758e+10, 1.56807621e+10, -8.81040892e+09,
- 4.56229864e+10, -1.58968401e+10, 9.20074349e+09,
- -5.91929988e+09, 6.44516729e+08, -2.58364312e+09,
- -6.42843866e+09, -6.26112144e+10, 7.23420074e+09,
- 2.19307621e+10, -1.13835037e+11, 2.93643123e+10,
- -1.58968401e+10, 1.90753098e+11, -5.20260223e+10,
- 3.94516729e+08, -1.43068463e+10, 1.54275093e+10,
- 4.69330855e+09, 1.27342007e+10, -1.17763321e+10,
- -1.25604089e+10, 3.01143123e+10, -5.25594796e+10,
- 9.20074349e+09, -5.20260223e+10, 9.11232962e+10,
- -1.33364312e+09, 9.17750929e+09, -2.67874845e+10,
- 6.90365551e+08, 1.99349442e+08, -9.97211896e+08,
- 3.39103160e+09, -5.85734201e+08, 2.29600372e+09,
- -5.91929988e+09, 3.94516729e+08, -1.33364312e+09,
- 1.83790273e+09, -8.13197026e+06, 3.48513011e+07,
- 4.49349442e+08, 3.67828377e+09, -2.59758364e+09,
- -1.08573420e+09, 8.31668216e+09, -5.36013011e+09,
- 6.44516729e+08, -1.43068463e+10, 9.17750929e+09,
- -8.13197026e+06, 2.31188042e+09, -1.21979554e+09,
- -1.99721190e+09, -5.09758364e+09, 1.42069393e+09,
- 4.54600372e+09, -9.11013011e+09, 1.48164498e+10,
- -2.58364312e+09, 1.54275093e+10, -2.67874845e+10,
- 3.48513011e+07, -1.21979554e+09, 1.05503408e+10,
+ 1.08203222e+10, 5.80793061e+09, -3.19702602e+08,
+ 5.03531599e+09, -1.66914498e+09, -2.27509294e+09,
+ -1.72763321e+10, -4.57125155e+09, 4.59107807e+09,
+ 1.42069393e+09, 4.32465923e+08, -1.99628253e+09,
+ 5.80793061e+09, 2.32515489e+10, -8.10408922e+08,
+ -1.66914498e+09, 3.51505576e+10, -7.37174721e+09,
+ -4.57125155e+09, -6.22131351e+10, 1.33122677e+10,
+ 4.32465923e+08, 3.81102850e+09, -5.13011152e+09,
+ -3.19702602e+08, -8.10408922e+08, 8.57372986e+09,
+ -2.27509294e+09, -7.37174721e+09, 1.09665428e+10,
+ 4.59107807e+09, 1.33122677e+10, -2.15452292e+10,
+ -1.99628253e+09, -5.13011152e+09, 2.00495663e+09,
+ 5.03531599e+09, -1.66914498e+09, -2.27509294e+09,
+ 4.48327138e+10, -2.22908922e+10, 1.19609665e+10,
+ -5.65594796e+10, 2.50743494e+10, -1.42472119e+10,
+ 6.69144981e+09, -1.11431227e+09, 4.56133829e+09,
+ -1.66914498e+09, 3.51505576e+10, -7.37174721e+09,
+ -2.22908922e+10, 7.52342007e+10, -2.26338290e+10,
+ 2.50743494e+10, -1.21016729e+11, 3.96524164e+10,
+ -1.11431227e+09, 1.06319703e+10, -9.64684015e+09,
+ -2.27509294e+09, -7.37174721e+09, 1.09665428e+10,
+ 1.19609665e+10, -2.26338290e+10, 4.51979554e+10,
+ -1.42472119e+10, 3.96524164e+10, -7.19962825e+10,
+ 4.56133829e+09, -9.64684015e+09, 1.58317844e+10,
+ -1.72763321e+10, -4.57125155e+09, 4.59107807e+09,
+ -5.65594796e+10, 2.50743494e+10, -1.42472119e+10,
+ 8.56232962e+10, -2.11957869e+10, 1.22676580e+10,
+ -1.17874845e+10, 6.92688971e+08, -2.61152416e+09,
+ -4.57125155e+09, -6.22131351e+10, 1.33122677e+10,
+ 2.50743494e+10, -1.21016729e+11, 3.96524164e+10,
+ -2.11957869e+10, 2.01727385e+11, -6.93680297e+10,
+ 6.92688971e+08, -1.84975217e+10, 1.64033457e+10,
+ 4.59107807e+09, 1.33122677e+10, -2.15452292e+10,
+ -1.42472119e+10, 3.96524164e+10, -7.19962825e+10,
+ 1.22676580e+10, -6.93680297e+10, 1.22023544e+11,
+ -2.61152416e+09, 1.64033457e+10, -2.84820322e+10,
+ 1.42069393e+09, 4.32465923e+08, -1.99628253e+09,
+ 6.69144981e+09, -1.11431227e+09, 4.56133829e+09,
+ -1.17874845e+10, 6.92688971e+08, -2.61152416e+09,
+ 3.67534077e+09, -1.08426270e+07, 4.64684015e+07,
+ 4.32465923e+08, 3.81102850e+09, -5.13011152e+09,
+ -1.11431227e+09, 1.06319703e+10, -9.64684015e+09,
+ 6.92688971e+08, -1.84975217e+10, 1.64033457e+10,
+ -1.08426270e+07, 4.05452292e+09, -1.62639405e+09,
+ -1.99628253e+09, -5.13011152e+09, 2.00495663e+09,
+ 4.56133829e+09, -9.64684015e+09, 1.58317844e+10,
+ -2.61152416e+09, 1.64033457e+10, -2.84820322e+10,
+ 4.64684015e+07, -1.62639405e+09, 1.06452912e+10,
};
pylith::feassemble::ElasticityImplicitData3DLinear::ElasticityImplicitData3DLinear(void)
Modified: short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityImplicitData3DQuadratic.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityImplicitData3DQuadratic.cc 2007-06-15 20:10:58 UTC (rev 7265)
+++ short/3D/PyLith/trunk/unittests/libtests/feassemble/data/ElasticityImplicitData3DQuadratic.cc 2007-06-15 20:11:09 UTC (rev 7266)
@@ -177,319 +177,319 @@
};
const double pylith::feassemble::ElasticityImplicitData3DQuadratic::_valsResidual[] = {
- -7.86265156e+10, -6.50453174e+09, 9.66313773e+10,
- -1.66117754e+10, -1.00795779e+11, 2.53173672e+10,
- -1.53180823e+11, -4.94876361e+11, 3.47626004e+11,
- 8.32634822e+09, -6.82511215e+10, 3.63126117e+11,
- 8.82738604e+10, -2.77446201e+11, -9.32014365e+10,
- 2.31344042e+11, 4.99241249e+11, -1.16153668e+11,
- -6.35712433e+11, -3.27240946e+11, -8.95734051e+10,
- 6.25497326e+10, 6.16264999e+11, -3.42992649e+11,
- 1.54879392e+11, 6.99521034e+11, -6.94797606e+11,
- 3.38758171e+11, -5.39912342e+11, 5.04017899e+11,
+ -4.23293588e+10, -2.77322955e+10, 1.60678757e+11,
+ -9.57785654e+09, -1.16546360e+11, 2.08212509e+10,
+ -2.32850642e+11, -5.45737494e+11, 3.99403558e+11,
+ 4.66589928e+10, -1.10817849e+11, 3.59789383e+11,
+ 3.81294856e+10, -2.30875676e+11, -1.04828412e+11,
+ 3.46864294e+11, 5.91207973e+11, -2.10576075e+11,
+ -7.59895172e+11, -4.18997477e+11, -6.23908371e+10,
+ 1.13569956e+11, 6.51755933e+11, -3.83093543e+11,
+ 2.17299147e+11, 7.79749940e+11, -7.81469581e+11,
+ 2.82131154e+11, -5.72006695e+11, 6.01665499e+11,
};
const double pylith::feassemble::ElasticityImplicitData3DQuadratic::_valsJacobian[] = {
- 1.64587725e+11, 9.84163932e+09, 4.44360094e+10,
- 4.31136969e+10, -9.67377317e+09, 9.68730846e+09,
- 1.72433393e+10, 1.68395158e+10, -1.37427697e+10,
- -2.76740098e+09, -3.75545131e+09, 2.10529583e+10,
- -1.92074818e+11, 2.83889103e+10, -4.11542992e+10,
- -7.33195554e+10, -6.04669733e+10, 3.13088959e+10,
- -3.79452598e+10, 6.83256899e+09, -6.75012498e+10,
- 4.32409462e+10, 1.05899903e+10, -7.12250802e+09,
- -7.32183538e+09, 1.35625469e+10, 5.29603146e+08,
- 4.52431626e+10, -1.21589738e+10, 2.25060515e+10,
- 9.84163932e+09, 3.74110321e+10, 5.45366987e+09,
- -3.78488428e+09, 5.08119747e+09, -2.53520789e+09,
- 9.28396020e+09, 4.62369779e+09, 2.86086181e+09,
- -2.08878464e+09, 2.90757333e+09, 2.26767442e+09,
- 8.61113251e+09, -3.11009516e+10, 6.61532462e+09,
- -3.41336399e+10, -1.80693031e+10, -1.56487044e+10,
- 2.77013431e+08, -2.29618143e+10, -4.07780027e+09,
- 8.70110143e+09, 1.11175517e+10, -7.34453128e+08,
- 6.22921354e+09, 1.16588860e+10, 3.09128144e+09,
- -2.93675159e+09, -6.67869453e+08, 2.70735351e+09,
- 4.44360094e+10, 5.45366987e+09, 7.74890906e+10,
- 1.13539751e+10, -5.97965233e+09, 1.60535378e+10,
- -3.85388086e+09, 8.86086181e+09, -2.40454954e+10,
- 9.49740277e+09, -2.87881134e+08, 3.13514179e+10,
- -4.77098547e+10, 1.79486579e+10, -6.96316741e+10,
- 4.30889593e+09, -3.36487044e+10, 5.58752110e+10,
- -3.39456943e+10, 2.58886639e+09, -1.20057080e+11,
- 1.98860309e+09, 2.82110243e+09, -7.90011998e+09,
- -1.35928574e+09, 5.53572589e+09, 7.44113765e+09,
- 1.52838293e+10, -3.29264649e+09, 3.34239741e+10,
- 4.31136969e+10, -3.78488428e+09, 1.13539751e+10,
- 1.42767964e+11, -4.34832278e+10, 2.15301897e+10,
- -5.43347808e+09, -1.42227223e+10, 1.55895434e+10,
- 5.22011338e+09, 3.68879538e+09, -2.04595732e+10,
- -1.63110784e+11, 2.07249856e+10, -3.39736363e+10,
- 3.93560332e+10, -1.56465667e+10, 2.41542882e+10,
- 6.68142494e+10, 1.75610693e+09, -7.48896515e+09,
- -5.33337677e+10, 6.40011381e+10, -5.74254265e+10,
- 8.24390755e+09, -1.55083612e+10, -1.20603962e+09,
- -8.36379340e+10, 2.47473635e+09, 4.79256444e+10,
- -9.67377317e+09, 5.08119747e+09, -5.97965233e+09,
- -4.34832278e+10, 8.23753811e+10, -1.44200455e+10,
- -5.00050007e+09, 3.20162597e+10, -1.38752685e+10,
- 3.55462049e+08, -1.00466279e+10, 1.34585702e+10,
- 3.73916523e+10, -3.39906977e+10, 1.76114051e+10,
- -1.18687890e+10, 3.41593741e+10, -1.54798356e+10,
- -8.91055973e+09, -4.51310595e+09, 8.38583869e+09,
- 3.42233603e+10, -1.41382125e+11, 4.95083091e+10,
- -8.61947236e+09, 2.25958013e+10, 2.10512559e+08,
- 1.55858475e+10, 1.37045424e+10, -3.94198336e+10,
- 9.68730846e+09, -2.53520789e+09, 1.60535378e+10,
- 2.15301897e+10, -1.44200455e+10, 4.58794614e+10,
- 6.47843228e+09, -1.15419352e+10, 2.21170686e+10,
- -9.68179544e+09, 7.68079238e+09, -1.48994108e+10,
- -3.05291919e+10, 9.38918285e+09, -5.26894578e+10,
- 1.48209548e+10, -1.09242801e+10, 3.20598211e+10,
- 3.99923744e+08, 5.83028314e+09, -3.81491876e+09,
- -3.06476488e+10, 4.09527536e+10, -8.33096254e+10,
- 2.38404821e+08, -1.78948744e+09, -4.53249037e+09,
- 1.77034222e+10, -2.26420559e+10, 4.31360142e+10,
- 1.72433393e+10, 9.28396020e+09, -3.85388086e+09,
- -5.43347808e+09, -5.00050007e+09, 6.47843228e+09,
- 7.16359466e+10, 2.59889613e+10, -2.72686901e+10,
- 1.65440527e+10, 3.16312861e+09, -8.65825157e+09,
- -6.52797899e+08, 1.85457182e+09, -5.63350106e+08,
- -6.89736250e+10, -3.22042140e+10, 1.74017376e+10,
- 1.40671395e+10, 1.19621355e+10, -1.22855569e+10,
- 1.51826692e+10, 7.07040096e+09, -1.39583273e+10,
- -7.74138196e+10, -1.62253654e+10, 3.77047625e+10,
- 1.78005733e+10, -5.89307895e+09, 5.00312450e+09,
- 1.68395158e+10, 4.62369779e+09, 8.86086181e+09,
- -1.42227223e+10, 3.20162597e+10, -1.15419352e+10,
- 2.59889613e+10, 1.80686439e+11, -7.74605257e+10,
- 4.82979528e+09, 1.79598853e+10, -2.06060780e+10,
- -3.43170670e+07, 4.44565836e+10, -3.46608469e+09,
- -5.23153251e+10, -9.12139712e+10, 4.30892780e+09,
- 2.21843577e+10, 3.26826875e+10, -1.58640573e+10,
- 3.37370676e+10, -1.55953018e+11, 6.25890730e+10,
- -2.27809210e+10, -1.34125295e+11, 9.06491951e+10,
- -1.42264123e+10, 6.88667317e+10, -3.74693768e+10,
- -1.37427697e+10, 2.86086181e+09, -2.40454954e+10,
- 1.55895434e+10, -1.38752685e+10, 2.21170686e+10,
- -2.72686901e+10, -7.74605257e+10, 1.88543480e+11,
- -7.88047379e+09, -1.22727446e+10, 6.36519849e+10,
- -3.41127883e+08, -1.15771958e+10, -2.78601808e+09,
- 4.75128487e+10, 2.23089278e+10, 1.90153362e+10,
- -2.01744457e+10, -1.37529462e+10, 3.23406830e+10,
- -4.07361051e+10, 6.80335174e+10, -1.26385136e+11,
- 3.43714292e+10, 6.72047507e+10, -2.58048902e+11,
- 1.26697912e+10, -3.14693768e+10, 8.55969985e+10,
- -2.76740098e+09, -2.08878464e+09, 9.49740277e+09,
- 5.22011338e+09, 3.55462049e+08, -9.68179544e+09,
- 1.65440527e+10, 4.82979528e+09, -7.88047379e+09,
- 4.33158550e+10, 3.21234686e+09, -1.76782633e+10,
- -1.47958880e+09, -3.89776204e+09, -1.86657290e+09,
- 1.79397167e+09, 1.65305987e+09, -1.16943436e+09,
- -2.56311069e+09, 4.73379414e+09, -2.09181973e+10,
- 3.14502449e+10, 5.14358097e+09, -1.42013065e+10,
- -5.62211955e+10, -1.10328159e+10, 2.35686692e+10,
- -3.52929416e+10, -2.90867657e+09, 4.03299716e+10,
- -3.75545131e+09, 2.90757333e+09, -2.87881134e+08,
- 3.68879538e+09, -1.00466279e+10, 7.68079238e+09,
- 3.16312861e+09, 1.79598853e+10, -1.22727446e+10,
- 3.21234686e+09, 4.39124124e+10, -1.87394346e+10,
- -2.45331760e+09, -1.23711093e+10, 6.63934952e+09,
- -1.23582902e+09, 2.48793486e+10, -1.51348107e+10,
- 8.17823858e+09, -1.88437554e+10, 4.75168196e+09,
- 6.58802541e+09, 2.82864270e+10, -1.30242825e+10,
- -7.58837146e+09, -7.21985379e+10, 4.70625956e+10,
- -9.79756546e+09, -4.48561615e+09, -6.67526587e+09,
- 2.10529583e+10, 2.26767442e+09, 3.13514179e+10,
- -2.04595732e+10, 1.34585702e+10, -1.48994108e+10,
- -8.65825157e+09, -2.06060780e+10, 6.36519849e+10,
- -1.76782633e+10, -1.87394346e+10, 2.17621308e+11,
- -2.53323957e+09, 1.51949051e+10, 4.16713128e+09,
- 8.60834341e+09, -1.96903663e+10, 8.41013823e+10,
- -5.44737529e+10, -5.02609582e+09, -1.44647291e+11,
- -2.33124176e+10, -1.70242825e+10, 5.40965949e+10,
- 2.37908914e+10, 7.36181511e+10, -2.54016802e+11,
- 7.36633049e+10, -2.34530436e+10, -4.14263159e+10,
- -1.92074818e+11, 8.61113251e+09, -4.77098547e+10,
- -1.63110784e+11, 3.73916523e+10, -3.05291919e+10,
- -6.52797899e+08, -3.43170670e+07, -3.41127883e+08,
- -1.47958880e+09, -2.45331760e+09, -2.53323957e+09,
- 3.53483013e+11, -5.12535904e+10, 6.65092647e+10,
- 1.02041521e+10, 6.79670426e+10, -5.54293711e+10,
- -3.05391652e+10, -1.37080077e+09, 7.40312483e+10,
- -1.19371939e+10, -7.41825181e+10, 5.62333076e+10,
- -2.96048615e+09, 1.57919013e+09, 1.16536228e+09,
- 3.90676692e+10, 1.37455264e+10, -6.13963978e+10,
- 2.83889103e+10, -3.11009516e+10, 1.79486579e+10,
- 2.07249856e+10, -3.39906977e+10, 9.38918285e+09,
- 1.85457182e+09, 4.44565836e+10, -1.15771958e+10,
- -3.89776204e+09, -1.23711093e+10, 1.51949051e+10,
- -5.12535904e+10, 1.96588810e+11, -9.80288892e+09,
- 3.40225982e+10, -1.77359245e+11, 7.23366995e+10,
- 2.21847548e+10, 6.06520314e+10, -6.27037313e+10,
- -5.28491848e+10, -5.60183528e+10, 7.24945231e+09,
- 1.57919013e+09, -4.05453065e+10, -4.63771111e+09,
- -7.54473600e+08, 4.96882379e+10, -3.33973705e+10,
- -4.11542992e+10, 6.61532462e+09, -6.96316741e+10,
- -3.39736363e+10, 1.76114051e+10, -5.26894578e+10,
- -5.63350106e+08, -3.46608469e+09, -2.78601808e+09,
- -1.86657290e+09, 6.63934952e+09, 4.16713128e+09,
- 6.65092647e+10, -9.80288892e+09, 1.35994046e+11,
- -2.19849267e+10, 5.08366995e+10, -6.56295612e+10,
- 3.63090261e+10, -3.40370646e+10, 6.88617703e+10,
- 2.23999743e+10, -2.36394366e+10, 9.03395269e+10,
- 1.16536228e+09, -4.63771111e+09, -5.77577532e+09,
- -2.68408422e+10, -6.11959277e+09, -1.02849988e+11,
- -7.33195554e+10, -3.41336399e+10, 4.30889593e+09,
- 3.93560332e+10, -1.18687890e+10, 1.48209548e+10,
- -6.89736250e+10, -5.23153251e+10, 4.75128487e+10,
- 1.79397167e+09, -1.23582902e+09, 8.60834341e+09,
- 1.02041521e+10, 3.40225982e+10, -2.19849267e+10,
- 3.12159136e+11, 2.92248145e+10, 5.60209313e+10,
- -7.12403900e+09, -1.47307540e+10, 3.47829963e+10,
- -2.00593645e+11, 2.76861351e+10, -4.22801960e+10,
- 5.14586591e+10, 3.95580486e+09, -6.91804987e+10,
- -6.49610880e+10, 1.93949844e+10, -3.26093490e+10,
- -6.04669733e+10, -1.80693031e+10, -3.36487044e+10,
- -1.56465667e+10, 3.41593741e+10, -1.09242801e+10,
- -3.22042140e+10, -9.12139712e+10, 2.23089278e+10,
- 1.65305987e+09, 2.48793486e+10, -1.96903663e+10,
- 6.79670426e+10, -1.77359245e+11, 5.08366995e+10,
- 2.92248145e+10, 3.58326610e+11, -1.16726238e+11,
- -3.40085318e+10, -1.12164051e+11, 1.04718662e+11,
- 1.89639128e+10, 2.35692392e+10, -1.69355813e+10,
- 5.12247153e+09, 3.59820192e+10, -2.07313406e+10,
- 1.93949844e+10, -7.81100215e+10, 4.07922215e+10,
- 3.13088959e+10, -1.56487044e+10, 5.58752110e+10,
- 2.41542882e+10, -1.54798356e+10, 3.20598211e+10,
- 1.74017376e+10, 4.30892780e+09, 1.90153362e+10,
- -1.16943436e+09, -1.51348107e+10, 8.41013823e+10,
- -5.54293711e+10, 7.23366995e+10, -6.56295612e+10,
- 5.60209313e+10, -1.16726238e+11, 4.18312117e+11,
- 5.47274407e+10, 7.42186615e+10, -2.67327985e+11,
- -7.06690849e+10, -7.32447015e+09, 4.62951788e+09,
- -2.37360543e+10, -2.13424517e+10, -1.30151249e+11,
- -3.26093490e+10, 4.07922215e+10, -1.50884590e+11,
- -3.79452598e+10, 2.77013431e+08, -3.39456943e+10,
- 6.68142494e+10, -8.91055973e+09, 3.99923744e+08,
- 1.40671395e+10, 2.21843577e+10, -2.01744457e+10,
- -2.56311069e+09, 8.17823858e+09, -5.44737529e+10,
- -3.05391652e+10, 2.21847548e+10, 3.63090261e+10,
- -7.12403900e+09, -3.40085318e+10, 5.47274407e+10,
- 3.86375132e+11, 5.02292341e+10, -3.06673830e+09,
- -9.83165906e+10, -1.75914346e+10, 2.59536009e+10,
- -1.91748363e+10, -9.32197943e+10, 5.77923363e+10,
- -2.71593519e+11, 5.06767217e+10, -6.35216965e+10,
- 6.83256899e+09, -2.29618143e+10, 2.58886639e+09,
- 1.75610693e+09, -4.51310595e+09, 5.83028314e+09,
- 1.19621355e+10, 3.26826875e+10, -1.37529462e+10,
- 4.73379414e+09, -1.88437554e+10, -5.02609582e+09,
- -1.37080077e+09, 6.06520314e+10, -3.40370646e+10,
- -1.47307540e+10, -1.12164051e+11, 7.42186615e+10,
- 5.02292341e+10, 1.48564361e+11, -2.74420686e+10,
- -1.75914346e+10, -3.81913183e+10, 1.19072739e+10,
- -4.78864609e+10, -6.23542125e+10, 3.36268953e+09,
- 6.06561062e+09, 1.71291766e+10, -1.76495993e+10,
- -6.75012498e+10, -4.07780027e+09, -1.20057080e+11,
- -7.48896515e+09, 8.38583869e+09, -3.81491876e+09,
- -1.22855569e+10, -1.58640573e+10, 3.23406830e+10,
- -2.09181973e+10, 4.75168196e+09, -1.44647291e+11,
- 7.40312483e+10, -6.27037313e+10, 6.88617703e+10,
- 3.47829963e+10, 1.04718662e+11, -2.67327985e+11,
- -3.06673830e+09, -2.74420686e+10, 3.89321192e+11,
- 2.59536009e+10, 1.19072739e+10, -4.75512737e+10,
- 2.49590029e+10, -1.37484216e+10, 8.43609178e+10,
- -4.84661410e+10, -5.92737708e+09, 8.51398448e+09,
- 4.32409462e+10, 8.70110143e+09, 1.98860309e+09,
- -5.33337677e+10, 3.42233603e+10, -3.06476488e+10,
- 1.51826692e+10, 3.37370676e+10, -4.07361051e+10,
- 3.14502449e+10, 6.58802541e+09, -2.33124176e+10,
- -1.19371939e+10, -5.28491848e+10, 2.23999743e+10,
- -2.00593645e+11, 1.89639128e+10, -7.06690849e+10,
- -9.83165906e+10, -1.75914346e+10, 2.59536009e+10,
- 3.60589444e+11, -1.00112592e+10, -1.35534859e+10,
- -5.00245072e+10, 9.32958569e+09, 7.64197916e+10,
- -3.62575999e+10, -3.10911747e+10, 5.21567726e+10,
- 1.05899903e+10, 1.11175517e+10, 2.82110243e+09,
- 6.40011381e+10, -1.41382125e+11, 4.09527536e+10,
- 7.07040096e+09, -1.55953018e+11, 6.80335174e+10,
- 5.14358097e+09, 2.82864270e+10, -1.70242825e+10,
- -7.41825181e+10, -5.60183528e+10, -2.36394366e+10,
- 2.76861351e+10, 2.35692392e+10, -7.32447015e+09,
- -1.75914346e+10, -3.81913183e+10, 1.19072739e+10,
- -1.00112592e+10, 4.51201237e+11, -1.63563227e+11,
- 1.24406968e+10, 5.24733623e+10, -4.94543585e+10,
- -2.51467303e+10, -1.75103003e+11, 1.37291127e+11,
- -7.12250802e+09, -7.34453128e+08, -7.90011998e+09,
- -5.74254265e+10, 4.95083091e+10, -8.33096254e+10,
- -1.39583273e+10, 6.25890730e+10, -1.26385136e+11,
- -1.42013065e+10, -1.30242825e+10, 5.40965949e+10,
- 5.62333076e+10, 7.24945231e+09, 9.03395269e+10,
- -4.22801960e+10, -1.69355813e+10, 4.62951788e+09,
- 2.59536009e+10, 1.19072739e+10, -4.75512737e+10,
- -1.35534859e+10, -1.63563227e+11, 4.08232344e+11,
- 3.41420138e+10, -3.81210252e+10, 6.86665515e+10,
- 3.22123281e+10, 1.01124461e+11, -3.60818380e+11,
- -7.32183538e+09, 6.22921354e+09, -1.35928574e+09,
- 8.24390755e+09, -8.61947236e+09, 2.38404821e+08,
- -7.74138196e+10, -2.27809210e+10, 3.43714292e+10,
- -5.62211955e+10, -7.58837146e+09, 2.37908914e+10,
- -2.96048615e+09, 1.57919013e+09, 1.16536228e+09,
- 5.14586591e+10, 5.12247153e+09, -2.37360543e+10,
- -1.91748363e+10, -4.78864609e+10, 2.49590029e+10,
- -5.00245072e+10, 1.24406968e+10, 3.41420138e+10,
- 1.54619326e+11, 2.14981019e+10, -5.77090025e+10,
- -1.20521220e+09, 4.00055518e+10, -3.58627619e+10,
- 1.35625469e+10, 1.16588860e+10, 5.53572589e+09,
- -1.55083612e+10, 2.25958013e+10, -1.78948744e+09,
- -1.62253654e+10, -1.34125295e+11, 6.72047507e+10,
- -1.10328159e+10, -7.21985379e+10, 7.36181511e+10,
- 1.57919013e+09, -4.05453065e+10, -4.63771111e+09,
- 3.95580486e+09, 3.59820192e+10, -2.13424517e+10,
- -9.32197943e+10, -6.23542125e+10, -1.37484216e+10,
- 9.32958569e+09, 5.24733623e+10, -3.81210252e+10,
- 2.14981019e+10, 3.50744424e+11, -1.23781855e+11,
- 8.60611074e+10, -1.64231141e+11, 5.70623242e+10,
- 5.29603146e+08, 3.09128144e+09, 7.44113765e+09,
- -1.20603962e+09, 2.10512559e+08, -4.53249037e+09,
- 3.77047625e+10, 9.06491951e+10, -2.58048902e+11,
- 2.35686692e+10, 4.70625956e+10, -2.54016802e+11,
- 1.16536228e+09, -4.63771111e+09, -5.77577532e+09,
- -6.91804987e+10, -2.07313406e+10, -1.30151249e+11,
- 5.77923363e+10, 3.36268953e+09, 8.43609178e+10,
- 7.64197916e+10, -4.94543585e+10, 6.86665515e+10,
- -5.77090025e+10, -1.23781855e+11, 5.36960708e+11,
- -6.90849841e+10, 5.42289908e+10, -4.49040962e+10,
- 4.52431626e+10, -2.93675159e+09, 1.52838293e+10,
- -8.36379340e+10, 1.55858475e+10, 1.77034222e+10,
- 1.78005733e+10, -1.42264123e+10, 1.26697912e+10,
- -3.52929416e+10, -9.79756546e+09, 7.36633049e+10,
- 3.90676692e+10, -7.54473600e+08, -2.68408422e+10,
- -6.49610880e+10, 1.93949844e+10, -3.26093490e+10,
- -2.71593519e+11, 6.06561062e+09, -4.84661410e+10,
- -3.62575999e+10, -2.51467303e+10, 3.22123281e+10,
- -1.20521220e+09, 8.60611074e+10, -6.90849841e+10,
- 3.90836890e+11, -7.42456166e+10, 2.54686407e+10,
- -1.21589738e+10, -6.67869453e+08, -3.29264649e+09,
- 2.47473635e+09, 1.37045424e+10, -2.26420559e+10,
- -5.89307895e+09, 6.88667317e+10, -3.14693768e+10,
- -2.90867657e+09, -4.48561615e+09, -2.34530436e+10,
- 1.37455264e+10, 4.96882379e+10, -6.11959277e+09,
- 1.93949844e+10, -7.81100215e+10, 4.07922215e+10,
- 5.06767217e+10, 1.71291766e+10, -5.92737708e+09,
- -3.10911747e+10, -1.75103003e+11, 1.01124461e+11,
- 4.00055518e+10, -1.64231141e+11, 5.42289908e+10,
- -7.42456166e+10, 2.73208963e+11, -1.03241580e+11,
- 2.25060515e+10, 2.70735351e+09, 3.34239741e+10,
- 4.79256444e+10, -3.94198336e+10, 4.31360142e+10,
- 5.00312450e+09, -3.74693768e+10, 8.55969985e+10,
- 4.03299716e+10, -6.67526587e+09, -4.14263159e+10,
- -6.13963978e+10, -3.33973705e+10, -1.02849988e+11,
- -3.26093490e+10, 4.07922215e+10, -1.50884590e+11,
- -6.35216965e+10, -1.76495993e+10, 8.51398448e+09,
- 5.21567726e+10, 1.37291127e+11, -3.60818380e+11,
- -3.58627619e+10, 5.70623242e+10, -4.49040962e+10,
- 2.54686407e+10, -1.03241580e+11, 5.30212399e+11,
+ 1.73593357e+11, 1.31221858e+10, 5.92480125e+10,
+ 4.41282224e+10, -8.97243830e+09, 1.40275224e+10,
+ 1.34679027e+10, 1.74156506e+10, -1.17311004e+10,
+ 2.50981776e+09, -3.89615730e+09, 2.03669074e+10,
+ -1.97580971e+11, 2.46666952e+10, -5.92427692e+10,
+ -6.39826914e+10, -6.30670755e+10, 2.37451946e+10,
+ -5.75008309e+10, 4.73972161e+09, -6.76312960e+10,
+ 4.15615137e+10, 1.28607278e+10, -3.42260328e+09,
+ -4.09074006e+09, 1.31945069e+10, -5.53121731e+08,
+ 4.78944202e+10, -1.00638169e+10, 2.51932539e+10,
+ 1.31221858e+10, 7.18520028e+10, 7.27155983e+09,
+ -8.97243830e+09, 1.37022228e+10, -5.67657348e+09,
+ 1.74156506e+10, 3.37218949e+09, 7.81448242e+09,
+ -3.89615730e+09, 7.04979720e+09, 1.31986219e+09,
+ 2.46666952e+10, -6.88018779e+10, 1.63759884e+10,
+ -6.30670755e+10, -1.97824896e+10, -3.28649392e+10,
+ 4.73972161e+09, -4.55140745e+10, -9.92622589e+08,
+ 1.28607278e+10, 1.58627981e+10, 1.39109953e+09,
+ 1.31945069e+10, 1.10938371e+10, 5.75133822e+09,
+ -1.00638169e+10, 1.11655945e+10, -3.90195315e+08,
+ 5.92480125e+10, 7.27155983e+09, 1.03914450e+11,
+ 1.40275224e+10, -5.67657348e+09, 2.24800951e+10,
+ -1.17311004e+10, 7.81448242e+09, -1.95631650e+10,
+ 2.03669074e+10, 1.31986219e+09, 2.98048728e+10,
+ -5.92427692e+10, 1.63759884e+10, -9.96264559e+10,
+ 2.37451946e+10, -3.28649392e+10, 3.93731217e+10,
+ -6.76312960e+10, -9.92622589e+08, -1.23190287e+11,
+ -3.42260328e+09, 1.39109953e+09, 6.48660713e+08,
+ -5.53121731e+08, 5.75133822e+09, 7.71963836e+09,
+ 2.51932539e+10, -3.90195315e+08, 3.84390694e+10,
+ 4.41282224e+10, -8.97243830e+09, 1.40275224e+10,
+ 1.54867792e+11, -5.79776370e+10, 2.87069197e+10,
+ 2.95819506e+09, -1.28154816e+10, 1.47119838e+10,
+ 1.21720190e+09, 2.69617162e+09, -2.00942458e+10,
+ -1.67957268e+11, 3.87444252e+10, -4.30018854e+10,
+ 4.73211108e+10, -1.83435705e+10, 2.59834953e+10,
+ 6.22243332e+10, -4.76963520e+09, -4.72602760e+09,
+ -8.43708418e+10, 6.54829989e+10, -5.87153835e+10,
+ 1.05412088e+10, -1.60852224e+10, -6.45089868e+08,
+ -7.09299539e+10, 1.20403892e+10, 4.37527110e+10,
+ -8.97243830e+09, 1.37022228e+10, -5.67657348e+09,
+ -5.79776370e+10, 1.06553726e+11, -1.92267274e+10,
+ -1.28154816e+10, 3.29179853e+10, -1.69448024e+10,
+ 2.69617162e+09, -1.09961911e+10, 1.40929084e+10,
+ 3.87444252e+10, -6.46611992e+10, 1.80003919e+10,
+ -1.83435705e+10, 4.31637836e+10, -1.76027438e+10,
+ -4.76963520e+09, 5.16244894e+09, 9.47741456e+09,
+ 6.54829989e+10, -1.54809527e+11, 6.03073752e+10,
+ -1.60852224e+10, 2.20227238e+10, -1.05264992e+09,
+ 1.20403892e+10, 6.94402728e+09, -4.13745930e+10,
+ 1.40275224e+10, -5.67657348e+09, 2.24800951e+10,
+ 2.87069197e+10, -1.92267274e+10, 7.73569901e+10,
+ 1.47119838e+10, -1.69448024e+10, 2.49986324e+10,
+ -2.00942458e+10, 1.40929084e+10, -1.48784174e+10,
+ -4.30018854e+10, 1.80003919e+10, -7.96202073e+10,
+ 2.59834953e+10, -1.76027438e+10, 4.14841411e+10,
+ -4.72602760e+09, 9.47741456e+09, 5.72099869e+09,
+ -5.87153835e+10, 6.03073752e+10, -1.08351528e+11,
+ -6.45089868e+08, -1.05264992e+09, 3.20090476e+08,
+ 4.37527110e+10, -4.13745930e+10, 3.04892048e+10,
+ 1.34679027e+10, 1.74156506e+10, -1.17311004e+10,
+ 2.95819506e+09, -1.28154816e+10, 1.47119838e+10,
+ 1.23438637e+11, 3.46519484e+10, -3.63582535e+10,
+ 2.79586306e+10, 5.32861592e+09, -1.10258169e+10,
+ 5.63042682e+09, 1.21350317e+09, -6.02985326e+08,
+ -7.63547390e+10, -5.63463594e+10, 4.32763908e+10,
+ 2.31172881e+10, 2.27643288e+10, -2.16400017e+10,
+ -2.79271875e+10, 2.72049791e+10, -3.64629550e+10,
+ -1.32369258e+11, -2.60041909e+10, 4.80507945e+10,
+ 4.00801042e+10, -1.34129942e+10, 1.17819438e+10,
+ 1.74156506e+10, 3.37218949e+09, 7.81448242e+09,
+ -1.28154816e+10, 3.29179853e+10, -1.69448024e+10,
+ 3.46519484e+10, 2.10679032e+11, -1.03280701e+11,
+ 5.32861592e+09, 2.90912967e+10, -2.19192151e+10,
+ 1.21350317e+09, 4.17179320e+10, -1.00288537e+10,
+ -5.63463594e+10, -9.41470160e+10, 1.77452371e+10,
+ 2.27643288e+10, 3.80097265e+10, -1.97446690e+10,
+ 2.72049791e+10, -1.64835738e+11, 8.70817270e+10,
+ -2.60041909e+10, -1.77738439e+11, 1.05235964e+11,
+ -1.34129942e+10, 8.09330309e+10, -4.59591691e+10,
+ -1.17311004e+10, 7.81448242e+09, -1.95631650e+10,
+ 1.47119838e+10, -1.69448024e+10, 2.49986324e+10,
+ -3.63582535e+10, -1.03280701e+11, 2.16964664e+11,
+ -1.10258169e+10, -2.19192151e+10, 6.56449763e+10,
+ -6.02985326e+08, -1.00288537e+10, 3.92385068e+09,
+ 4.32763908e+10, 1.77452371e+10, -5.96357004e+09,
+ -2.16400017e+10, -1.97446690e+10, 3.77361229e+10,
+ -3.64629550e+10, 8.70817270e+10, -1.41181432e+11,
+ 4.80507945e+10, 1.05235964e+11, -2.76877324e+11,
+ 1.17819438e+10, -4.59591691e+10, 9.43172444e+10,
+ 2.50981776e+09, -3.89615730e+09, 2.03669074e+10,
+ 1.21720190e+09, 2.69617162e+09, -2.00942458e+10,
+ 2.79586306e+10, 5.32861592e+09, -1.10258169e+10,
+ 8.03801203e+10, 4.28312915e+09, -2.35710177e+10,
+ -2.63620606e+09, -4.23405310e+09, -2.93320831e+09,
+ 1.80513827e+10, 2.78153902e+08, 4.95927270e+09,
+ -2.69586121e+10, 8.60802181e+09, -5.02613001e+10,
+ 4.22351859e+10, 7.82107092e+09, -2.50091494e+10,
+ -1.02342437e+11, -1.24141249e+10, 3.15730404e+10,
+ -4.04150843e+10, -8.47082802e+09, 7.59955177e+10,
+ -3.89615730e+09, 7.04979720e+09, 1.31986219e+09,
+ 2.69617162e+09, -1.09961911e+10, 1.40929084e+10,
+ 5.32861592e+09, 2.90912967e+10, -2.19192151e+10,
+ 4.28312915e+09, 8.08573662e+10, -2.49859128e+10,
+ -4.23405310e+09, -1.13494224e+10, 1.45561697e+10,
+ 2.78153902e+08, 3.65196843e+10, -2.32167847e+10,
+ 8.60802181e+09, -3.99831278e+10, -1.82942578e+08,
+ 7.82107092e+09, 3.97041316e+10, -2.00323766e+10,
+ -1.24141249e+10, -1.15124311e+11, 8.04538311e+10,
+ -8.47082802e+09, -1.57692240e+10, -2.00855397e+10,
+ 2.03669074e+10, 1.31986219e+09, 2.98048728e+10,
+ -2.00942458e+10, 1.40929084e+10, -1.48784174e+10,
+ -1.10258169e+10, -2.19192151e+10, 6.56449763e+10,
+ -2.35710177e+10, -2.49859128e+10, 2.19824483e+11,
+ -2.93320831e+09, 1.45561697e+10, 1.88117001e+09,
+ 4.95927270e+09, -2.32167847e+10, 8.38973112e+10,
+ -5.02613001e+10, -1.82942578e+08, -1.40625956e+11,
+ -2.50091494e+10, -2.00323766e+10, 6.03522660e+10,
+ 3.15730404e+10, 8.04538311e+10, -2.60578922e+11,
+ 7.59955177e+10, -2.00855397e+10, -4.53217837e+10,
+ -1.97580971e+11, 2.46666952e+10, -5.92427692e+10,
+ -1.67957268e+11, 3.87444252e+10, -4.30018854e+10,
+ 5.63042682e+09, 1.21350317e+09, -6.02985326e+08,
+ -2.63620606e+09, -4.23405310e+09, -2.93320831e+09,
+ 3.85696291e+11, -6.83381205e+10, 8.86790196e+10,
+ -2.67543764e+10, 6.79930939e+10, -5.16095318e+10,
+ -9.58513672e+09, 1.38759693e+10, 7.35601830e+10,
+ -6.19215806e+09, -8.46878019e+10, 5.24221879e+10,
+ -9.76062412e+09, 2.10558684e+09, 1.55381637e+09,
+ 2.91400233e+10, 8.66070187e+09, -5.88248267e+10,
+ 2.46666952e+10, -6.88018779e+10, 1.63759884e+10,
+ 3.87444252e+10, -6.46611992e+10, 1.80003919e+10,
+ 1.21350317e+09, 4.17179320e+10, -1.00288537e+10,
+ -4.23405310e+09, -1.13494224e+10, 1.45561697e+10,
+ -6.83381205e+10, 2.60180928e+11, -1.30705186e+10,
+ 6.79930939e+10, -1.76805094e+11, 8.21155993e+10,
+ 1.38759693e+10, 6.33678206e+10, -6.44938640e+10,
+ -8.46878019e+10, -4.14570852e+10, -1.09266562e+10,
+ 2.10558684e+09, -3.98284804e+10, -6.18361482e+09,
+ 8.66070187e+09, 3.76364782e+10, -2.63446422e+10,
+ -5.92427692e+10, 1.63759884e+10, -9.96264559e+10,
+ -4.30018854e+10, 1.80003919e+10, -7.96202073e+10,
+ -6.02985326e+08, -1.00288537e+10, 3.92385068e+09,
+ -2.93320831e+09, 1.45561697e+10, 1.88117001e+09,
+ 8.86790196e+10, -1.30705186e+10, 2.11705117e+11,
+ -5.16095318e+10, 8.21155993e+10, -8.74213471e+10,
+ 7.35601830e+10, -6.44938640e+10, 6.99356117e+10,
+ 5.24221879e+10, -1.09266562e+10, 7.56292185e+10,
+ 1.55381637e+09, -6.18361482e+09, -1.20128555e+10,
+ -5.88248267e+10, -2.63446422e+10, -8.43941023e+10,
+ -6.39826914e+10, -6.30670755e+10, 2.37451946e+10,
+ 4.73211108e+10, -1.83435705e+10, 2.59834953e+10,
+ -7.63547390e+10, -5.63463594e+10, 4.32763908e+10,
+ 1.80513827e+10, 2.78153902e+08, 4.95927270e+09,
+ -2.67543764e+10, 6.79930939e+10, -5.16095318e+10,
+ 4.13046988e+11, 3.89664193e+10, 7.46945751e+10,
+ -6.36916424e+10, -3.24928572e+10, 5.96736247e+10,
+ -1.86334149e+11, 3.11000319e+10, -7.52995207e+10,
+ 3.47603416e+10, 6.05218426e+09, -6.19443687e+10,
+ -9.60622254e+10, 2.58599791e+10, -4.34791320e+10,
+ -6.30670755e+10, -1.97824896e+10, -3.28649392e+10,
+ -1.83435705e+10, 4.31637836e+10, -1.76027438e+10,
+ -5.63463594e+10, -9.41470160e+10, 1.77452371e+10,
+ 2.78153902e+08, 3.65196843e+10, -2.32167847e+10,
+ 6.79930939e+10, -1.76805094e+11, 8.21155993e+10,
+ 3.89664193e+10, 4.49980968e+11, -1.55634984e+11,
+ -3.24928572e+10, -1.47723652e+11, 1.19291549e+11,
+ 3.11000319e+10, -7.00384183e+09, -1.61733676e+10,
+ 6.05218426e+09, 2.23790298e+10, -2.80491949e+10,
+ 2.58599791e+10, -1.06581372e+11, 5.43896286e+10,
+ 2.37451946e+10, -3.28649392e+10, 3.93731217e+10,
+ 2.59834953e+10, -1.76027438e+10, 4.14841411e+10,
+ 4.32763908e+10, 1.77452371e+10, -5.96357004e+09,
+ 4.95927270e+09, -2.32167847e+10, 8.38973112e+10,
+ -5.16095318e+10, 8.21155993e+10, -8.74213471e+10,
+ 7.46945751e+10, -1.55634984e+11, 4.97969373e+11,
+ 5.96736247e+10, 1.19291549e+11, -2.71854799e+11,
+ -7.52995207e+10, -1.61733676e+10, -2.21556189e+10,
+ -6.19443687e+10, -2.80491949e+10, -1.10527585e+11,
+ -4.34791320e+10, 5.43896286e+10, -1.64801027e+11,
+ -5.75008309e+10, 4.73972161e+09, -6.76312960e+10,
+ 6.22243332e+10, -4.76963520e+09, -4.72602760e+09,
+ 2.31172881e+10, 2.27643288e+10, -2.16400017e+10,
+ -2.69586121e+10, 8.60802181e+09, -5.02613001e+10,
+ -9.58513672e+09, 1.38759693e+10, 7.35601830e+10,
+ -6.36916424e+10, -3.24928572e+10, 5.96736247e+10,
+ 4.47739208e+11, 6.69723122e+10, -4.08898440e+09,
+ -1.06262150e+11, -2.34552462e+10, 3.46048012e+10,
+ -1.49150887e+10, -9.40708368e+10, 5.51675595e+10,
+ -2.54167369e+11, 3.78282216e+10, -7.46585584e+10,
+ 4.73972161e+09, -4.55140745e+10, -9.92622589e+08,
+ -4.76963520e+09, 5.16244894e+09, 9.47741456e+09,
+ 2.27643288e+10, 3.80097265e+10, -1.97446690e+10,
+ 8.60802181e+09, -3.99831278e+10, -1.82942578e+08,
+ 1.38759693e+10, 6.33678206e+10, -6.44938640e+10,
+ -3.24928572e+10, -1.47723652e+11, 1.19291549e+11,
+ 6.69723122e+10, 2.57490592e+11, -3.65894247e+10,
+ -2.34552462e+10, -5.81619320e+10, 1.58763652e+10,
+ -9.40708368e+10, -4.94585896e+10, -6.92382137e+09,
+ 3.78282216e+10, -2.31892124e+10, -1.57179843e+10,
+ -6.76312960e+10, -9.92622589e+08, -1.23190287e+11,
+ -4.72602760e+09, 9.47741456e+09, 5.72099869e+09,
+ -2.16400017e+10, -1.97446690e+10, 3.77361229e+10,
+ -5.02613001e+10, -1.82942578e+08, -1.40625956e+11,
+ 7.35601830e+10, -6.44938640e+10, 6.99356117e+10,
+ 5.96736247e+10, 1.19291549e+11, -2.71854799e+11,
+ -4.08898440e+09, -3.65894247e+10, 4.50096056e+11,
+ 3.46048012e+10, 1.58763652e+10, -6.56498964e+10,
+ 5.51675595e+10, -6.92382137e+09, 6.79135146e+10,
+ -7.46585584e+10, -1.57179843e+10, -3.00813662e+10,
+ 4.15615137e+10, 1.28607278e+10, -3.42260328e+09,
+ -8.43708418e+10, 6.54829989e+10, -5.87153835e+10,
+ -2.79271875e+10, 2.72049791e+10, -3.64629550e+10,
+ 4.22351859e+10, 7.82107092e+09, -2.50091494e+10,
+ -6.19215806e+09, -8.46878019e+10, 5.24221879e+10,
+ -1.86334149e+11, 3.11000319e+10, -7.52995207e+10,
+ -1.06262150e+11, -2.34552462e+10, 3.46048012e+10,
+ 4.71475009e+11, -1.33483456e+10, -1.80713146e+10,
+ -2.93522947e+10, 1.45135217e+10, 7.37078702e+10,
+ -1.14832927e+11, -3.74919366e+10, 5.62460671e+10,
+ 1.28607278e+10, 1.58627981e+10, 1.39109953e+09,
+ 6.54829989e+10, -1.54809527e+11, 6.03073752e+10,
+ 2.72049791e+10, -1.64835738e+11, 8.70817270e+10,
+ 7.82107092e+09, 3.97041316e+10, -2.00323766e+10,
+ -8.46878019e+10, -4.14570852e+10, -1.09266562e+10,
+ 3.11000319e+10, -7.00384183e+09, -1.61733676e+10,
+ -2.34552462e+10, -5.81619320e+10, 1.58763652e+10,
+ -1.33483456e+10, 5.43964443e+11, -2.18084302e+11,
+ 1.45135217e+10, 5.26460008e+10, -5.83835892e+10,
+ -3.74919366e+10, -2.25909250e+11, 1.58943725e+11,
+ -3.42260328e+09, 1.39109953e+09, 6.48660713e+08,
+ -5.87153835e+10, 6.03073752e+10, -1.08351528e+11,
+ -3.64629550e+10, 8.70817270e+10, -1.41181432e+11,
+ -2.50091494e+10, -2.00323766e+10, 6.03522660e+10,
+ 5.24221879e+10, -1.09266562e+10, 7.56292185e+10,
+ -7.52995207e+10, -1.61733676e+10, -2.21556189e+10,
+ 3.46048012e+10, 1.58763652e+10, -6.56498964e+10,
+ -1.80713146e+10, -2.18084302e+11, 5.09589329e+11,
+ 7.37078702e+10, -5.83835892e+10, 6.56005522e+10,
+ 5.62460671e+10, 1.58943725e+11, -3.74481552e+11,
+ -4.09074006e+09, 1.31945069e+10, -5.53121731e+08,
+ 1.05412088e+10, -1.60852224e+10, -6.45089868e+08,
+ -1.32369258e+11, -2.60041909e+10, 4.80507945e+10,
+ -1.02342437e+11, -1.24141249e+10, 3.15730404e+10,
+ -9.76062412e+09, 2.10558684e+09, 1.55381637e+09,
+ 3.47603416e+10, 6.05218426e+09, -6.19443687e+10,
+ -1.49150887e+10, -9.40708368e+10, 5.51675595e+10,
+ -2.93522947e+10, 1.45135217e+10, 7.37078702e+10,
+ 2.80044129e+11, 2.86641359e+10, -7.69453367e+10,
+ -3.25152372e+10, 8.40444394e+10, -6.99651640e+10,
+ 1.31945069e+10, 1.10938371e+10, 5.75133822e+09,
+ -1.60852224e+10, 2.20227238e+10, -1.05264992e+09,
+ -2.60041909e+10, -1.77738439e+11, 1.05235964e+11,
+ -1.24141249e+10, -1.15124311e+11, 8.04538311e+10,
+ 2.10558684e+09, -3.98284804e+10, -6.18361482e+09,
+ 6.05218426e+09, 2.23790298e+10, -2.80491949e+10,
+ -9.40708368e+10, -4.94585896e+10, -6.92382137e+09,
+ 1.45135217e+10, 5.26460008e+10, -5.83835892e+10,
+ 2.86641359e+10, 4.36944208e+11, -1.65042473e+11,
+ 8.40444394e+10, -1.62935980e+11, 7.41942100e+10,
+ -5.53121731e+08, 5.75133822e+09, 7.71963836e+09,
+ -6.45089868e+08, -1.05264992e+09, 3.20090476e+08,
+ 4.80507945e+10, 1.05235964e+11, -2.76877324e+11,
+ 3.15730404e+10, 8.04538311e+10, -2.60578922e+11,
+ 1.55381637e+09, -6.18361482e+09, -1.20128555e+10,
+ -6.19443687e+10, -2.80491949e+10, -1.10527585e+11,
+ 5.51675595e+10, -6.92382137e+09, 6.79135146e+10,
+ 7.37078702e+10, -5.83835892e+10, 6.56005522e+10,
+ -7.69453367e+10, -1.65042473e+11, 5.85917235e+11,
+ -6.99651640e+10, 7.41942100e+10, -6.74743444e+10,
+ 4.78944202e+10, -1.00638169e+10, 2.51932539e+10,
+ -7.09299539e+10, 1.20403892e+10, 4.37527110e+10,
+ 4.00801042e+10, -1.34129942e+10, 1.17819438e+10,
+ -4.04150843e+10, -8.47082802e+09, 7.59955177e+10,
+ 2.91400233e+10, 8.66070187e+09, -5.88248267e+10,
+ -9.60622254e+10, 2.58599791e+10, -4.34791320e+10,
+ -2.54167369e+11, 3.78282216e+10, -7.46585584e+10,
+ -1.14832927e+11, -3.74919366e+10, 5.62460671e+10,
+ -3.25152372e+10, 8.40444394e+10, -6.99651640e+10,
+ 4.91808250e+11, -9.89941555e+10, 3.39581876e+10,
+ -1.00638169e+10, 1.11655945e+10, -3.90195315e+08,
+ 1.20403892e+10, 6.94402728e+09, -4.13745930e+10,
+ -1.34129942e+10, 8.09330309e+10, -4.59591691e+10,
+ -8.47082802e+09, -1.57692240e+10, -2.00855397e+10,
+ 8.66070187e+09, 3.76364782e+10, -2.63446422e+10,
+ 2.58599791e+10, -1.06581372e+11, 5.43896286e+10,
+ 3.78282216e+10, -2.31892124e+10, -1.57179843e+10,
+ -3.74919366e+10, -2.25909250e+11, 1.58943725e+11,
+ 8.40444394e+10, -1.62935980e+11, 7.41942100e+10,
+ -9.89941555e+10, 3.97705908e+11, -1.37655440e+11,
+ 2.51932539e+10, -3.90195315e+08, 3.84390694e+10,
+ 4.37527110e+10, -4.13745930e+10, 3.04892048e+10,
+ 1.17819438e+10, -4.59591691e+10, 9.43172444e+10,
+ 7.59955177e+10, -2.00855397e+10, -4.53217837e+10,
+ -5.88248267e+10, -2.63446422e+10, -8.43941023e+10,
+ -4.34791320e+10, 5.43896286e+10, -1.64801027e+11,
+ -7.46585584e+10, -1.57179843e+10, -3.00813662e+10,
+ 5.62460671e+10, 1.58943725e+11, -3.74481552e+11,
+ -6.99651640e+10, 7.41942100e+10, -6.74743444e+10,
+ 3.39581876e+10, -1.37655440e+11, 6.03308657e+11,
};
pylith::feassemble::ElasticityImplicitData3DQuadratic::ElasticityImplicitData3DQuadratic(void)
Modified: short/3D/PyLith/trunk/unittests/libtests/feassemble/data/IntegratorElasticity.py
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/feassemble/data/IntegratorElasticity.py 2007-06-15 20:10:58 UTC (rev 7265)
+++ short/3D/PyLith/trunk/unittests/libtests/feassemble/data/IntegratorElasticity.py 2007-06-15 20:11:09 UTC (rev 7266)
@@ -116,9 +116,9 @@
C2222 = lambda2mu
C2212 = 0.0
C1212 = 2.0*self.lameMu
- D = numpy.array([ [C1111, C1122, C1112],
- [C1122, C2222, C2212],
- [C1112, C2212, C1212] ],
+ D = numpy.array([ [C1111, C1122, 0.5*C1112],
+ [C1122, C2222, 0.5*C2212],
+ [0.5*C1112, 0.5*C2212, 0.5*C1212] ],
dtype=numpy.float64)
elif 3 == self.cellDim:
lambda2mu = self.lameLambda + 2.0*self.lameMu
@@ -143,12 +143,12 @@
C2323 = 2.0*self.lameMu
C2313 = 0.0
C1313 = 2.0*self.lameMu
- D = numpy.array([ [C1111, C1122, C1133, C1112, C1123, C1113],
- [C1122, C2222, C2233, C2212, C2223, C2213],
- [C1133, C2233, C3333, C3312, C3323, C3313],
- [C1112, C2212, C3312, C1212, C1223, C1213],
- [C1123, C2223, C3323, C1223, C2323, C2313],
- [C1113, C2213, C3313, C1213, C2313, C1313] ],
+ D = numpy.array([ [C1111, C1122, C1133, 0.5*C1112, 0.5*C1123, 0.5*C1113],
+ [C1122, C2222, C2233, 0.5*C2212, 0.5*C2223, 0.5*C2213],
+ [C1133, C2233, C3333, 0.5*C3312, 0.5*C3323, 0.5*C3313],
+ [0.5*C1112, 0.5*C2212, 0.5*C3312, 0.5*C1212, 0.5*C1223, 0.5*C1213],
+ [0.5*C1123, 0.5*C2223, 0.5*C3323, 0.5*C1223, 0.5*C2323, 0.5*C2313],
+ [0.5*C1113, 0.5*C2213, 0.5*C3313, 0.5*C1213, 0.5*C2313, 0.5*C1313] ],
dtype=numpy.float64)
return D
@@ -176,20 +176,20 @@
B[0, iBasis*self.spaceDim+0] = basisDeriv[iQuad, iBasis, 0]
B[1, iBasis*self.spaceDim+1] = basisDeriv[iQuad, iBasis, 1]
B[2, iBasis*self.spaceDim+2] = basisDeriv[iQuad, iBasis, 2]
- B[3, iBasis*self.spaceDim+0] = 0.5*basisDeriv[iQuad, iBasis, 1]
- B[3, iBasis*self.spaceDim+1] = 0.5*basisDeriv[iQuad, iBasis, 0]
- B[4, iBasis*self.spaceDim+1] = 0.5*basisDeriv[iQuad, iBasis, 2]
- B[4, iBasis*self.spaceDim+2] = 0.5*basisDeriv[iQuad, iBasis, 1]
- B[5, iBasis*self.spaceDim+0] = 0.5*basisDeriv[iQuad, iBasis, 2]
- B[5, iBasis*self.spaceDim+2] = 0.5*basisDeriv[iQuad, iBasis, 0]
+ B[3, iBasis*self.spaceDim+0] = basisDeriv[iQuad, iBasis, 1]
+ B[3, iBasis*self.spaceDim+1] = basisDeriv[iQuad, iBasis, 0]
+ B[4, iBasis*self.spaceDim+1] = basisDeriv[iQuad, iBasis, 2]
+ B[4, iBasis*self.spaceDim+2] = basisDeriv[iQuad, iBasis, 1]
+ B[5, iBasis*self.spaceDim+0] = basisDeriv[iQuad, iBasis, 2]
+ B[5, iBasis*self.spaceDim+2] = basisDeriv[iQuad, iBasis, 0]
elif 2 == self.spaceDim:
B = numpy.zeros( (3, self.spaceDim*self.numBasis),
dtype=numpy.float64)
for iBasis in xrange(self.numBasis):
B[0, iBasis*self.spaceDim+0] = basisDeriv[iQuad, iBasis, 0]
B[1, iBasis*self.spaceDim+1] = basisDeriv[iQuad, iBasis, 1]
- B[2, iBasis*self.spaceDim+0] = 0.5*basisDeriv[iQuad, iBasis, 1]
- B[2, iBasis*self.spaceDim+1] = 0.5*basisDeriv[iQuad, iBasis, 0]
+ B[2, iBasis*self.spaceDim+0] = basisDeriv[iQuad, iBasis, 1]
+ B[2, iBasis*self.spaceDim+1] = basisDeriv[iQuad, iBasis, 0]
elif 1 == self.spaceDim:
B = numpy.zeros( (1, self.spaceDim*self.numBasis),
dtype=numpy.float64)
More information about the cig-commits
mailing list