[cig-commits] r19038 - in short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults: . data
brad at geodynamics.org
brad at geodynamics.org
Thu Oct 6 17:32:19 PDT 2011
Author: brad
Date: 2011-10-06 17:32:19 -0700 (Thu, 06 Oct 2011)
New Revision: 19038
Modified:
short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/TestFaultCohesiveKin.cc
short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataQuad4.cc
short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataQuad4e.cc
short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTet4.cc
short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTet4e.cc
short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTri3.cc
short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTri3d.cc
Log:
Updated test data to use collocated quadrature/vertices.
Modified: short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/TestFaultCohesiveKin.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/TestFaultCohesiveKin.cc 2011-10-07 00:28:44 UTC (rev 19037)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/TestFaultCohesiveKin.cc 2011-10-07 00:32:19 UTC (rev 19038)
@@ -226,6 +226,23 @@
} // for
} // for
+ // Check area
+ const ALE::Obj<RealSection>& areaSection =
+ fault._fields->get("area").section();
+ CPPUNIT_ASSERT(!areaSection.isNull());
+ iVertex = 0;
+ for (SieveSubMesh::label_sequence::iterator v_iter=verticesBegin;
+ v_iter != verticesEnd;
+ ++v_iter, ++iVertex) {
+ const int fiberDim = areaSection->getFiberDimension(*v_iter);
+ CPPUNIT_ASSERT_EQUAL(1, fiberDim);
+ const double* areaVertex = areaSection->restrictPoint(*v_iter);
+ CPPUNIT_ASSERT(0 != areaVertex);
+
+ const double tolerance = 1.0e-06;
+ CPPUNIT_ASSERT_DOUBLES_EQUAL(_data->area[iVertex], areaVertex[0],
+ tolerance);
+ } // for
} // testInitialize
// ----------------------------------------------------------------------
@@ -271,7 +288,7 @@
fault.useSolnIncr(false);
fault.integrateResidual(residual, t, &fields);
- //residual.view("RESIDUAL"); // DEBUGGING
+ residual.view("RESIDUAL"); // DEBUGGING
// Check values
const double* valsE = _data->residual;
@@ -399,7 +416,7 @@
for (int iCol=0; iCol < ncols; ++iCol) {
const int index = ncols*iRow+iCol;
const double valE = valsE[index];
-#if 0 // DEBUGGING
+#if 1 // DEBUGGING
if (fabs(valE-vals[index]) > tolerance)
std::cout << "ERROR: iRow: " << iRow << ", iCol: " << iCol
<< "valE: " << valE
@@ -883,6 +900,8 @@
residual.newSection(topology::FieldBase::VERTICES_FIELD, spaceDim);
residual.allocate();
fields->copyLayout("residual");
+
+ fault->verifyConfiguration(*mesh);
} // _initialize
// ----------------------------------------------------------------------
Modified: short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataQuad4.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataQuad4.cc 2011-10-07 00:28:44 UTC (rev 19037)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataQuad4.cc 2011-10-07 00:32:19 UTC (rev 19038)
@@ -58,24 +58,24 @@
const int pylith::faults::CohesiveKinDataQuad4::_numBasis = 2;
-const int pylith::faults::CohesiveKinDataQuad4::_numQuadPts = 1;
+const int pylith::faults::CohesiveKinDataQuad4::_numQuadPts = 2;
const double pylith::faults::CohesiveKinDataQuad4::_quadPts[] = {
- 0.0,
+ -1.0, 1.0,
};
const double pylith::faults::CohesiveKinDataQuad4::_quadWts[] = {
- 2.0,
+ 1.0, 1.0
};
const double pylith::faults::CohesiveKinDataQuad4::_basis[] = {
- 0.5,
- 0.5
+ 1.0, 0.0,
+ 0.0, 1.0,
};
const double pylith::faults::CohesiveKinDataQuad4::_basisDeriv[] = {
- -0.5,
- 0.5
+ -0.5, 0.5,
+ -0.5, 0.5,
};
const double pylith::faults::CohesiveKinDataQuad4::_verticesRef[] = {
@@ -171,31 +171,31 @@
const double pylith::faults::CohesiveKinDataQuad4::_residual[] = {
0.0, 0.0,
0.0, 0.0,
- +8.4, +9.4, // 4
- +8.4, +9.4, // 5
+ +8.8, +9.8, // 4
+ +8.0, +9.0, // 5
0.0, 0.0,
0.0, 0.0,
- -8.4, -9.4, // 8
- -8.4, -9.4, // 9
- -0.5*(8.7-8.3 + 8.9-8.4) + 0.5*(0.14794836271+0.08241148423),
- -0.5*(9.7-9.3 + 9.9-9.4) + 0.5*(1.77538035254+1.89546413727), // 10
- -0.5*(8.7-8.3 + 8.9-8.4) + 0.5*(0.14794836271+0.08241148423),
- -0.5*(9.7-9.3 + 9.9-9.4) + 0.5*(1.77538035254+1.89546413727), // 11
+ -8.8, -9.8, // 8
+ -8.0, -9.0, // 9
+ -(8.7-8.3) + 0.14794836271,
+ -(9.7-9.3) + 1.77538035254, // 10
+ -(8.9-8.4) + 0.08241148423,
+ -(9.9-9.4) + 1.89546413727, // 11
};
const double pylith::faults::CohesiveKinDataQuad4::_residualIncr[] = {
0.0, 0.0,
0.0, 0.0,
- +8.4, +9.4, // 4
- +8.4, +9.4, // 5
+ +8.8, +9.8, // 4
+ +8.0, +9.0, // 5
0.0, 0.0,
0.0, 0.0,
- -8.4, -9.4, // 8
- -8.4, -9.4, // 9
- -0.5*(8.7-8.3 + 8.9-8.4) + 0.5*(0.14794836271+0.08241148423),
- -0.5*(9.7-9.3 + 9.9-9.4) + 0.5*(1.77538035254+1.89546413727), // 10
- -0.5*(8.7-8.3 + 8.9-8.4) + 0.5*(0.14794836271+0.08241148423),
- -0.5*(9.7-9.3 + 9.9-9.4) + 0.5*(1.77538035254+1.89546413727), // 11
+ -8.8, -9.8, // 8
+ -8.0, -9.0, // 9
+ -(8.7-8.3) + 0.14794836271,
+ -(9.7-9.3) + 1.77538035254, // 10
+ -(8.9-8.4) + 0.08241148423,
+ -(9.9-9.4) + 1.89546413727, // 11
};
const double pylith::faults::CohesiveKinDataQuad4::_jacobian[] = {
@@ -247,8 +247,8 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- -0.5, 0.0, // 10
- -0.5, 0.0, // 11
+ -1.0, 0.0, // 10
+ 0.0, 0.0,
0.0, 0.0, // 4y
0.0, 0.0,
0.0, 0.0,
@@ -257,8 +257,8 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,-0.5, // 10
- 0.0,-0.5, // 11
+ 0.0,-1.0, // 10
+ 0.0, 0.0,
0.0, 0.0, // 5x
0.0, 0.0,
0.0, 0.0,
@@ -267,8 +267,8 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- -0.5, 0.0, // 10
- -0.5, 0.0, // 11
+ 0.0, 0.0,
+ -1.0, 0.0, // 11
0.0, 0.0, // 5y
0.0, 0.0,
0.0, 0.0,
@@ -277,8 +277,8 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,-0.5, // 10
- 0.0,-0.5, // 11
+ 0.0, 0.0,
+ 0.0,-1.0, // 11
0.0, 0.0, // 6x
0.0, 0.0,
0.0, 0.0,
@@ -327,8 +327,8 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- +0.5, 0.0, // 10
- +0.5, 0.0, // 11
+ +1.0, 0.0, // 10
+ 0.0, 0.0,
0.0, 0.0, // 8y
0.0, 0.0,
0.0, 0.0,
@@ -337,8 +337,8 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,+0.5, // 10
- 0.0,+0.5, // 11
+ 0.0,+1.0, // 10
+ 0.0, 0.0,
0.0, 0.0, // 9x
0.0, 0.0,
0.0, 0.0,
@@ -347,8 +347,8 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- +0.5, 0.0, // 10
- +0.5, 0.0, // 11
+ 0.0, 0.0,
+ +1.0, 0.0, // 11
0.0, 0.0, // 9y
0.0, 0.0,
0.0, 0.0,
@@ -357,48 +357,48 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,+0.5, // 10
- 0.0,+0.5, // 11
+ 0.0, 0.0,
+ 0.0,+1.0, // 11
0.0, 0.0, // 10x
0.0, 0.0,
- -0.5, 0.0, // 4
- -0.5, 0.0, // 5
+ -1.0, 0.0, // 4
0.0, 0.0,
0.0, 0.0,
- +0.5, 0.0, // 8
- +0.5, 0.0, // 9
0.0, 0.0,
+ +1.0, 0.0, // 8
0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 10y
0.0, 0.0,
- 0.0,-0.5, // 4
- 0.0,-0.5, // 5
+ 0.0,-1.0, // 4
0.0, 0.0,
0.0, 0.0,
- 0.0,+0.5, // 8
- 0.0,+0.5, // 9
0.0, 0.0,
+ 0.0,+1.0, // 8
0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 11x
0.0, 0.0,
- -0.5, 0.0, // 4
- -0.5, 0.0, // 5
0.0, 0.0,
+ -1.0, 0.0, // 5
0.0, 0.0,
- +0.5, 0.0, // 8
- +0.5, 0.0, // 9
0.0, 0.0,
0.0, 0.0,
+ +1.0, 0.0, // 9
+ 0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 11y
0.0, 0.0,
- 0.0,-0.5, // 4
- 0.0,-0.5, // 5
0.0, 0.0,
+ 0.0,-1.0, // 5
0.0, 0.0,
- 0.0,+0.5, // 8
- 0.0,+0.5, // 9
0.0, 0.0,
0.0, 0.0,
+ 0.0,+1.0, // 9
+ 0.0, 0.0,
+ 0.0, 0.0,
};
const double pylith::faults::CohesiveKinDataQuad4::_fieldIncrAdjusted[] = {
Modified: short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataQuad4e.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataQuad4e.cc 2011-10-07 00:28:44 UTC (rev 19037)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataQuad4e.cc 2011-10-07 00:32:19 UTC (rev 19038)
@@ -77,24 +77,24 @@
const int pylith::faults::CohesiveKinDataQuad4e::_numBasis = 2;
-const int pylith::faults::CohesiveKinDataQuad4e::_numQuadPts = 1;
+const int pylith::faults::CohesiveKinDataQuad4e::_numQuadPts = 2;
const double pylith::faults::CohesiveKinDataQuad4e::_quadPts[] = {
- 0.0,
+ -1.0, 1.0,
};
const double pylith::faults::CohesiveKinDataQuad4e::_quadWts[] = {
- 2.0,
+ 1.0, 1.0,
};
const double pylith::faults::CohesiveKinDataQuad4e::_basis[] = {
- 0.5,
- 0.5
+ 1.0, 0.0,
+ 0.0, 1.0,
};
const double pylith::faults::CohesiveKinDataQuad4e::_basisDeriv[] = {
- -0.5,
- 0.5
+ -0.5, 0.5,
+ -0.5, 0.5,
};
const double pylith::faults::CohesiveKinDataQuad4e::_verticesRef[] = {
@@ -167,7 +167,7 @@
2.2, 3.2, // 14
2.4, 3.4, // 15
1.0, 1.0, // 16
- 2.0, 2.0, // 17
+ 1.0, 1.0, // 17
1.0, 1.0, // 18
};
@@ -211,47 +211,43 @@
const double pylith::faults::CohesiveKinDataQuad4e::_residual[] = {
0.0, 0.0,
0.0, 0.0,
- +4.2, +6.2, // 6
- +8.6,+12.6, // 7
+ +1.0*4.1, +1.0*6.1, // 6
+ +2.0*4.3, +2.0*6.3, // 7
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- +4.4, +6.4, // 11
+ +1.0*4.5, +1.0*6.5, // 11
0.0, 0.0,
- -4.2, -6.2, // 13
- -8.6,-12.6, // 15
- -4.4, -6.4, // 17
- -0.5*(3.0-3.3 + 4.2-3.4) - 0.5*(0.14794836271 + 0.08241148423),
- -0.5*(5.0-5.3 + 6.2-5.4) - 0.5*(1.77538035254 + 1.89546413727), // 16
- -0.5*(3.0-3.3 + 4.2-3.4) - 0.5*(0.14794836271 + 0.08241148423) +
- -0.5*(4.2-3.4 + 4.4-3.8) - 0.5*(0.08241148423 + 0.19186497837),
- -0.5*(5.0-5.3 + 6.2-5.4) - 0.5*(1.77538035254 + 1.89546413727) +
- -0.5*(6.2-5.4 + 6.4-5.8) - 0.5*(1.89546413727 + 1.59887481971), // 17
- -0.5*(4.2-3.4 + 4.4-3.8) - 0.5*(0.08241148423 + 0.19186497837),
- -0.5*(6.2-5.4 + 6.4-5.8) - 0.5*(1.89546413727 + 1.59887481971), // 18
+ -1.0*4.1, -1.0*6.1, // 13
+ -2.0*4.3, -2.0*6.3, // 14
+ -1.0*4.5, -1.0*6.5, // 15
+ -1.0*(3.0-3.3 + 0.14794836271),
+ -1.0*(5.0-5.3 + 1.77538035254), // 16
+ -2.0*(4.2-3.4 + 0.08241148423),
+ -2.0*(6.2-5.4 + 1.89546413727), // 17
+ -1.0*(4.4-3.8 + 0.19186497837),
+ -1.0*(6.4-5.8 + 1.59887481971), // 18
};
const double pylith::faults::CohesiveKinDataQuad4e::_residualIncr[] = {
0.0, 0.0,
0.0, 0.0,
- +4.2, +6.2, // 6
- +8.6,+12.6, // 7
+ +1.0*4.1, +1.0*6.1, // 6
+ +2.0*4.3, +2.0*6.3, // 7
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- +4.4, +6.4, // 11
+ +1.0*4.5, +1.0*6.5, // 11
0.0, 0.0,
- -4.2, -6.2, // 13
- -8.6,-12.6, // 15
- -4.4, -6.4, // 17
- -0.5*(3.0-3.3 + 4.2-3.4) - 0.5*(0.14794836271 + 0.08241148423),
- -0.5*(5.0-5.3 + 6.2-5.4) - 0.5*(1.77538035254 + 1.89546413727), // 16
- -0.5*(3.0-3.3 + 4.2-3.4) - 0.5*(0.14794836271 + 0.08241148423) +
- -0.5*(4.2-3.4 + 4.4-3.8) - 0.5*(0.08241148423 + 0.19186497837),
- -0.5*(5.0-5.3 + 6.2-5.4) - 0.5*(1.77538035254 + 1.89546413727) +
- -0.5*(6.2-5.4 + 6.4-5.8) - 0.5*(1.89546413727 + 1.59887481971), // 17
- -0.5*(4.2-3.4 + 4.4-3.8) - 0.5*(0.08241148423 + 0.19186497837),
- -0.5*(6.2-5.4 + 6.4-5.8) - 0.5*(1.89546413727 + 1.59887481971), // 18
+ -1.0*4.1, -1.0*6.1, // 13
+ -2.0*4.3, -2.0*6.3, // 14
+ -1.0*4.5, -1.0*6.5, // 15
+ -1.0*(3.0-3.3 + 0.14794836271),
+ -1.0*(5.0-5.3 + 1.77538035254), // 16
+ -2.0*(4.2-3.4 + 0.08241148423),
+ -2.0*(6.2-5.4 + 1.89546413727), // 17
+ -1.0*(4.4-3.8 + 0.19186497837),
+ -1.0*(6.4-5.8 + 1.59887481971), // 18
};
const double pylith::faults::CohesiveKinDataQuad4e::_jacobian[] = {
@@ -327,9 +323,9 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- -0.5, 0.0, // 16
- -0.5, 0.0, // 17
+ -1.0, 0.0, // 16
0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 6y
0.0, 0.0,
0.0, 0.0,
@@ -342,9 +338,9 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,-0.5, // 16
- 0.0,-0.5, // 17
+ 0.0,-1.0, // 16
0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 7x
0.0, 0.0,
0.0, 0.0,
@@ -357,9 +353,9 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- -0.5, 0.0, // 16
- -1.0, 0.0, // 17
- -0.5, 0.0, // 18
+ 0.0, 0.0,
+ -2.0, 0.0, // 17
+ 0.0, 0.0,
0.0, 0.0, // 7y
0.0, 0.0,
0.0, 0.0,
@@ -372,9 +368,9 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,-0.5, // 16
- 0.0,-1.0, // 17
- 0.0,-0.5, // 18
+ 0.0, 0.0,
+ 0.0,-2.0, // 17
+ 0.0, 0.0,
0.0, 0.0, // 8x
0.0, 0.0,
0.0, 0.0,
@@ -478,8 +474,8 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- -0.5, 0.0, // 17
- -0.5, 0.0, // 18
+ 0.0, 0.0,
+ -1.0, 0.0, // 18
0.0, 0.0, // 11y
0.0, 0.0,
0.0, 0.0,
@@ -493,8 +489,8 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,-0.5, // 17
- 0.0,-0.5, // 18
+ 0.0, 0.0,
+ 0.0,-1.0, // 18
0.0, 0.0, // 12x
0.0, 0.0,
0.0, 0.0,
@@ -537,9 +533,9 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- +0.5, 0.0, // 16
- +0.5, 0.0, // 17
+ +1.0, 0.0, // 16
0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 13y
0.0, 0.0,
0.0, 0.0,
@@ -552,9 +548,9 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,+0.5, // 16
- 0.0,+0.5, // 17
+ 0.0,+1.0, // 16
0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 14x
0.0, 0.0,
0.0, 0.0,
@@ -567,9 +563,9 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- +0.5, 0.0, // 16
- +1.0, 0.0, // 17
- +0.5, 0.0, // 18
+ 0.0, 0.0,
+ +2.0, 0.0, // 17
+ 0.0, 0.0,
0.0, 0.0, // 14y
0.0, 0.0,
0.0, 0.0,
@@ -582,9 +578,9 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,+0.5, // 16
- 0.0,+1.0, // 17
- 0.0,+0.5, // 18
+ 0.0, 0.0,
+ 0.0,+2.0, // 17
+ 0.0, 0.0,
0.0, 0.0, // 15x
0.0, 0.0,
0.0, 0.0,
@@ -598,8 +594,8 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- +0.5, 0.0, // 17
- +0.5, 0.0, // 18
+ 0.0, 0.0,
+ +1.0, 0.0, // 18
0.0, 0.0, // 15y
0.0, 0.0,
0.0, 0.0,
@@ -613,98 +609,98 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,+0.5, // 17
- 0.0,+0.5, // 18
+ 0.0, 0.0,
+ 0.0,+1.0, // 18
0.0, 0.0, // 16x
0.0, 0.0,
- -0.5, 0.0, // 6
- -0.5, 0.0, // 7
+ -1.0, 0.0, // 6
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- +0.5, 0.0, // 13
- +0.5, 0.0, // 14
0.0, 0.0,
+ +1.0, 0.0, // 13
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 16y
0.0, 0.0,
- 0.0,-0.5, // 6
- 0.0,-0.5, // 7
+ 0.0,-1.0, // 6
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,+0.5, // 13
- 0.0,+0.5, // 14
0.0, 0.0,
+ 0.0,+1.0, // 13
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 17x
0.0, 0.0,
- -0.5, 0.0, // 6
- -1.0, 0.0, // 7
0.0, 0.0,
+ -2.0, 0.0, // 7
0.0, 0.0,
0.0, 0.0,
- -0.5, 0.0, // 11
0.0, 0.0,
- +0.5, 0.0, // 13
- +1.0, 0.0, // 14
- +0.5, 0.0, // 15
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
+ +2.0, 0.0, // 14
+ 0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 17y
0.0, 0.0,
- 0.0,-0.5, // 6
- 0.0,-1.0, // 7
0.0, 0.0,
+ 0.0,-2.0, // 7
0.0, 0.0,
0.0, 0.0,
- 0.0,-0.5, // 11
0.0, 0.0,
- 0.0,+0.5, // 13
- 0.0,+1.0, // 14
- 0.0,+0.5, // 15
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
+ 0.0,+2.0, // 14
+ 0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 18x
0.0, 0.0,
0.0, 0.0,
- -0.5, 0.0, // 7
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- -0.5, 0.0, // 11
0.0, 0.0,
+ -1.0, 0.0, // 11
0.0, 0.0,
- +0.5, 0.0, // 14
- +0.5, 0.0, // 15
0.0, 0.0,
0.0, 0.0,
+ +1.0, 0.0, // 15
0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 18y
0.0, 0.0,
0.0, 0.0,
- 0.0,-0.5, // 7
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,-0.5, // 11
0.0, 0.0,
+ 0.0,-1.0, // 11
0.0, 0.0,
- 0.0,+0.5, // 14
- 0.0,+0.5, // 15
0.0, 0.0,
0.0, 0.0,
+ 0.0,+1.0, // 15
0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
};
pylith::faults::CohesiveKinDataQuad4e::CohesiveKinDataQuad4e(void)
Modified: short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTet4.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTet4.cc 2011-10-07 00:28:44 UTC (rev 19037)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTet4.cc 2011-10-07 00:32:19 UTC (rev 19038)
@@ -44,24 +44,34 @@
const int pylith::faults::CohesiveKinDataTet4::_numBasis = 3;
-const int pylith::faults::CohesiveKinDataTet4::_numQuadPts = 1;
+const int pylith::faults::CohesiveKinDataTet4::_numQuadPts = 3;
const double pylith::faults::CohesiveKinDataTet4::_quadPts[] = {
- -3.33333333e-01, -3.33333333e-01,
+ -1.00000000e+00, -1.00000000e+00,
+ 1.00000000e+00, -1.00000000e+00,
+ -1.00000000e+00, 1.00000000e+00,
};
const double pylith::faults::CohesiveKinDataTet4::_quadWts[] = {
- 2.0,
+ 2.0/3.0, 2.0/3.0, 2.0/3.0,
};
const double pylith::faults::CohesiveKinDataTet4::_basis[] = {
- 3.33333333e-01, 3.33333333e-01,
- 3.33333333e-01,};
+ 1.0, 0.0, 0.0,
+ 0.0, 1.0, 0.0,
+ 0.0, 0.0, 1.0,
+};
const double pylith::faults::CohesiveKinDataTet4::_basisDeriv[] = {
-0.50000000e+00, -0.50000000e+00,
0.50000000e+00, 0.00000000e+00,
0.00000000e+00, 0.50000000e+00,
+ -0.50000000e+00, -0.50000000e+00,
+ 0.50000000e+00, 0.00000000e+00,
+ 0.00000000e+00, 0.50000000e+00,
+ -0.50000000e+00, -0.50000000e+00,
+ 0.50000000e+00, 0.00000000e+00,
+ 0.00000000e+00, 0.50000000e+00,
};
const double pylith::faults::CohesiveKinDataTet4::_verticesRef[] = {
@@ -120,9 +130,9 @@
1.6, 1.6, 1.6, // 7
1.8, 1.8, 1.8, // 8
1.0, 1.0, 1.0, // 9
- 1.0/3.0, 1.0/3.0, 1.0/3.0, // 10
- 1.0/3.0, 1.0/3.0, 1.0/3.0, // 11
- 1.0/3.0, 1.0/3.0, 1.0/3.0, // 12
+ 1.0, 1.0, 1.0, // 10
+ 1.0, 1.0, 1.0, // 11
+ 1.0, 1.0, 1.0, // 12
};
const int pylith::faults::CohesiveKinDataTet4::_numFaultVertices = 3;
@@ -163,60 +173,42 @@
const double pylith::faults::CohesiveKinDataTet4::_residual[] = {
0.0, 0.0, 0.0,
- +(7.7+7.9+7.1)/9.0, +(8.7+8.9+8.1)/9.0, +(9.7+9.9+9.1)/9.0, // 3
- +(7.7+7.9+7.1)/9.0, +(8.7+8.9+8.1)/9.0, +(9.7+9.9+9.1)/9.0, // 4
- +(7.7+7.9+7.1)/9.0, +(8.7+8.9+8.1)/9.0, +(9.7+9.9+9.1)/9.0, // 5
+ +7.7/3.0, +8.7/3.0, +9.7/3.0, // 3
+ +7.9/3.0, +8.9/3.0, +9.9/3.0, // 4
+ +7.1/3.0, +8.1/3.0, +9.1/3.0, // 5
0.0, 0.0, 0.0,
- -(7.7+7.9+7.1)/9.0, -(8.7+8.9+8.1)/9.0, -(9.7+9.9+9.1)/9.0, // 7
- -(7.7+7.9+7.1)/9.0, -(8.7+8.9+8.1)/9.0, -(9.7+9.9+9.1)/9.0, // 8
- -(7.7+7.9+7.1)/9.0, -(8.7+8.9+8.1)/9.0, -(9.7+9.9+9.1)/9.0, // 9
- -(7.6-7.2 + 7.8-7.3 + 7.0-7.4)/9.0 +
- -(-0.07938069066 + -0.14140241667 + -0.18205179147)/9.0,
- -(8.6-8.2 + 8.8-8.3 + 9.0-9.4)/9.0 +
- -(-1.69682900001 + -1.82575588523 + -1.51709826228)/9.0,
- -(9.6-9.2 + 9.8-9.3 + 9.0-9.4)/9.0 +
- -(0.55566483464 + 0.56560966667 + 0.54615537442)/9.0, // 10
- -(7.6-7.2 + 7.8-7.3 + 7.0-7.4)/9.0 +
- -(-0.07938069066 + -0.14140241667 + -0.18205179147)/9.0,
- -(8.6-8.2 + 8.8-8.3 + 9.0-9.4)/9.0 +
- -(-1.69682900001 + -1.82575588523 + -1.51709826228)/9.0,
- -(9.6-9.2 + 9.8-9.3 + 9.0-9.4)/9.0 +
- -(0.55566483464 + 0.56560966667 + 0.54615537442)/9.0, // 11
- -(7.6-7.2 + 7.8-7.3 + 7.0-7.4)/9.0 +
- -(-0.07938069066 + -0.14140241667 + -0.18205179147)/9.0,
- -(8.6-8.2 + 8.8-8.3 + 9.0-9.4)/9.0 +
- -(-1.69682900001 + -1.82575588523 + -1.51709826228)/9.0,
- -(9.6-9.2 + 9.8-9.3 + 9.0-9.4)/9.0 +
- -(0.55566483464 + 0.56560966667 + 0.54615537442)/9.0, // 12
+ -7.7/3.0, -8.7/3.0, -9.7/3.0, // 7
+ -7.9/3.0, -8.9/3.0, -9.9/3.0, // 8
+ -7.1/3.0, -8.1/3.0, -9.1/3.0, // 9
+ -1.0/3.0*(7.6-7.2 + -0.07938069066),
+ -1.0/3.0*(8.6-8.2 + -1.82575588523),
+ -1.0/3.0*(9.6-9.2 + 0.55566483464), // 10
+ -1.0/3.0*(7.8-7.3 + -0.14140241667),
+ -1.0/3.0*(8.8-8.3 + -1.69682900001),
+ -1.0/3.0*(9.8-9.3 + 0.56560966667), // 11
+ -1.0/3.0*(7.0-7.4 + -0.18205179147),
+ -1.0/3.0*(8.0-8.4 + -1.51709826228),
+ -1.0/3.0*(9.0-9.4 + 0.54615537442), // 12
};
const double pylith::faults::CohesiveKinDataTet4::_residualIncr[] = {
0.0, 0.0, 0.0,
- +(7.7+7.9+7.1)/9.0, +(8.7+8.9+8.1)/9.0, +(9.7+9.9+9.1)/9.0, // 3
- +(7.7+7.9+7.1)/9.0, +(8.7+8.9+8.1)/9.0, +(9.7+9.9+9.1)/9.0, // 4
- +(7.7+7.9+7.1)/9.0, +(8.7+8.9+8.1)/9.0, +(9.7+9.9+9.1)/9.0, // 5
+ +7.7/3.0, +8.7/3.0, +9.7/3.0, // 3
+ +7.9/3.0, +8.9/3.0, +9.9/3.0, // 4
+ +7.1/3.0, +8.1/3.0, +9.1/3.0, // 5
0.0, 0.0, 0.0,
- -(7.7+7.9+7.1)/9.0, -(8.7+8.9+8.1)/9.0, -(9.7+9.9+9.1)/9.0, // 7
- -(7.7+7.9+7.1)/9.0, -(8.7+8.9+8.1)/9.0, -(9.7+9.9+9.1)/9.0, // 8
- -(7.7+7.9+7.1)/9.0, -(8.7+8.9+8.1)/9.0, -(9.7+9.9+9.1)/9.0, // 9
- -(7.6-7.2 + 7.8-7.3 + 7.0-7.4)/9.0 +
- -(-0.07938069066 + -0.14140241667 + -0.18205179147)/9.0,
- -(8.6-8.2 + 8.8-8.3 + 9.0-9.4)/9.0 +
- -(-1.69682900001 + -1.82575588523 + -1.51709826228)/9.0,
- -(9.6-9.2 + 9.8-9.3 + 9.0-9.4)/9.0 +
- -(0.55566483464 + 0.56560966667 + 0.54615537442)/9.0, // 10
- -(7.6-7.2 + 7.8-7.3 + 7.0-7.4)/9.0 +
- -(-0.07938069066 + -0.14140241667 + -0.18205179147)/9.0,
- -(8.6-8.2 + 8.8-8.3 + 9.0-9.4)/9.0 +
- -(-1.69682900001 + -1.82575588523 + -1.51709826228)/9.0,
- -(9.6-9.2 + 9.8-9.3 + 9.0-9.4)/9.0 +
- -(0.55566483464 + 0.56560966667 + 0.54615537442)/9.0, // 11
- -(7.6-7.2 + 7.8-7.3 + 7.0-7.4)/9.0 +
- -(-0.07938069066 + -0.14140241667 + -0.18205179147)/9.0,
- -(8.6-8.2 + 8.8-8.3 + 9.0-9.4)/9.0 +
- -(-1.69682900001 + -1.82575588523 + -1.51709826228)/9.0,
- -(9.6-9.2 + 9.8-9.3 + 9.0-9.4)/9.0 +
- -(0.55566483464 + 0.56560966667 + 0.54615537442)/9.0, // 12
+ -7.7/3.0, -8.7/3.0, -9.7/3.0, // 7
+ -7.9/3.0, -8.9/3.0, -9.9/3.0, // 8
+ -7.1/3.0, -8.1/3.0, -9.1/3.0, // 9
+ -1.0/3.0*(7.6-7.2 + -0.07938069066),
+ -1.0/3.0*(8.6-8.2 + -1.82575588523),
+ -1.0/3.0*(9.6-9.2 + 0.55566483464), // 10
+ -1.0/3.0*(7.8-7.3 + -0.14140241667),
+ -1.0/3.0*(8.8-8.3 + -1.69682900001),
+ -1.0/3.0*(9.8-9.3 + 0.56560966667), // 11
+ -1.0/3.0*(7.0-7.4 + -0.18205179147),
+ -1.0/3.0*(8.0-8.4 + -1.51709826228),
+ -1.0/3.0*(9.0-9.4 + 0.54615537442), // 12
};
const double pylith::faults::CohesiveKinDataTet4::_jacobian[] = {
@@ -261,9 +253,9 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0/9.0, 0.0, 0.0, // 10
- -1.0/9.0, 0.0, 0.0, // 11
- -1.0/9.0, 0.0, 0.0, // 12
+ -1.0/3.0, 0.0, 0.0, // 10
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 3y
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -272,9 +264,9 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0/9.0, 0.0, // 10
- 0.0,-1.0/9.0, 0.0, // 11
- 0.0,-1.0/9.0, 0.0, // 12
+ 0.0,-1.0/3.0, 0.0, // 10
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 3z
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -283,9 +275,9 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0/9.0, // 10
- 0.0, 0.0,-1.0/9.0, // 11
- 0.0, 0.0,-1.0/9.0, // 12
+ 0.0, 0.0,-1.0/3.0, // 10
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 4x
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -294,9 +286,9 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0/9.0, 0.0, 0.0, // 10
- -1.0/9.0, 0.0, 0.0, // 11
- -1.0/9.0, 0.0, 0.0, // 12
+ 0.0, 0.0, 0.0,
+ -1.0/3.0, 0.0, 0.0, // 11
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 4y
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -305,9 +297,9 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0/9.0, 0.0, // 10
- 0.0,-1.0/9.0, 0.0, // 11
- 0.0,-1.0/9.0, 0.0, // 12
+ 0.0, 0.0, 0.0,
+ 0.0,-1.0/3.0, 0.0, // 11
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 4z
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -316,9 +308,9 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0/9.0, // 10
- 0.0, 0.0,-1.0/9.0, // 11
- 0.0, 0.0,-1.0/9.0, // 12
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0,-1.0/3.0, // 11
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 5x
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -327,9 +319,9 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0/9.0, 0.0, 0.0, // 10
- -1.0/9.0, 0.0, 0.0, // 11
- -1.0/9.0, 0.0, 0.0, // 12
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ -1.0/3.0, 0.0, 0.0, // 12
0.0, 0.0, 0.0, // 5y
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -338,9 +330,9 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0/9.0, 0.0, // 10
- 0.0,-1.0/9.0, 0.0, // 11
- 0.0,-1.0/9.0, 0.0, // 12
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0,-1.0/3.0, 0.0, // 12
0.0, 0.0, 0.0, // 5z
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -349,9 +341,9 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0/9.0, // 10
- 0.0, 0.0,-1.0/9.0, // 11
- 0.0, 0.0,-1.0/9.0, // 12
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0,-1.0/3.0, // 12
0.0, 0.0, 0.0, // 6x
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -393,9 +385,9 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +1.0/9.0, 0.0, 0.0, // 10
- +1.0/9.0, 0.0, 0.0, // 11
- +1.0/9.0, 0.0, 0.0, // 12
+ +1.0/3.0, 0.0, 0.0, // 10
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 7y
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -404,9 +396,9 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0/9.0, 0.0, // 10
- 0.0,+1.0/9.0, 0.0, // 11
- 0.0,+1.0/9.0, 0.0, // 12
+ 0.0,+1.0/3.0, 0.0, // 10
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 7z
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -415,9 +407,9 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0/9.0, // 10
- 0.0, 0.0,+1.0/9.0, // 11
- 0.0, 0.0,+1.0/9.0, // 12
+ 0.0, 0.0,+1.0/3.0, // 10
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 8x
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -426,9 +418,9 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +1.0/9.0, 0.0, 0.0, // 10
- +1.0/9.0, 0.0, 0.0, // 11
- +1.0/9.0, 0.0, 0.0, // 12
+ 0.0, 0.0, 0.0,
+ +1.0/3.0, 0.0, 0.0, // 11
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 8y
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -437,9 +429,9 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0/9.0, 0.0, // 10
- 0.0,+1.0/9.0, 0.0, // 11
- 0.0,+1.0/9.0, 0.0, // 12
+ 0.0, 0.0, 0.0,
+ 0.0,+1.0/3.0, 0.0, // 11
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 8z
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -448,9 +440,9 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0/9.0, // 10
- 0.0, 0.0,+1.0/9.0, // 11
- 0.0, 0.0,+1.0/9.0, // 12
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0,+1.0/3.0, // 11
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 9x
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -459,9 +451,9 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +1.0/9.0, 0.0, 0.0, // 10
- +1.0/9.0, 0.0, 0.0, // 11
- +1.0/9.0, 0.0, 0.0, // 12
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ +1.0/3.0, 0.0, 0.0, // 12
0.0, 0.0, 0.0, // 9y
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -470,9 +462,9 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0/9.0, 0.0, // 10
- 0.0,+1.0/9.0, 0.0, // 11
- 0.0,+1.0/9.0, 0.0, // 12
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0,+1.0/3.0, 0.0, // 12
0.0, 0.0, 0.0, // 9z
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -481,108 +473,108 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0/9.0, // 10
- 0.0, 0.0,+1.0/9.0, // 11
- 0.0, 0.0,+1.0/9.0, // 12
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0,+1.0/3.0, // 12
0.0, 0.0, 0.0, // 10x
- -1.0/9.0, 0.0, 0.0, // 3
- -1.0/9.0, 0.0, 0.0, // 4
- -1.0/9.0, 0.0, 0.0, // 5
+ -1.0/3.0, 0.0, 0.0, // 3
0.0, 0.0, 0.0,
- +1.0/9.0, 0.0, 0.0, // 7
- +1.0/9.0, 0.0, 0.0, // 8
- +1.0/9.0, 0.0, 0.0, // 9
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ +1.0/3.0, 0.0, 0.0, // 7
0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 10y
- 0.0,-1.0/9.0, 0.0, // 3
- 0.0,-1.0/9.0, 0.0, // 4
- 0.0,-1.0/9.0, 0.0, // 5
+ 0.0,-1.0/3.0, 0.0, // 3
0.0, 0.0, 0.0,
- 0.0,+1.0/9.0, 0.0, // 7
- 0.0,+1.0/9.0, 0.0, // 8
- 0.0,+1.0/9.0, 0.0, // 9
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ 0.0,+1.0/3.0, 0.0, // 7
0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 10z
- 0.0, 0.0,-1.0/9.0, // 3
- 0.0, 0.0,-1.0/9.0, // 4
- 0.0, 0.0,-1.0/9.0, // 5
+ 0.0, 0.0,-1.0/3.0, // 3
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0/9.0, // 7
- 0.0, 0.0,+1.0/9.0, // 8
- 0.0, 0.0,+1.0/9.0, // 9
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ 0.0, 0.0,+1.0/3.0, // 7
0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 11x
- -1.0/9.0, 0.0, 0.0, // 3
- -1.0/9.0, 0.0, 0.0, // 4
- -1.0/9.0, 0.0, 0.0, // 5
0.0, 0.0, 0.0,
- +1.0/9.0, 0.0, 0.0, // 7
- +1.0/9.0, 0.0, 0.0, // 8
- +1.0/9.0, 0.0, 0.0, // 9
+ -1.0/3.0, 0.0, 0.0, // 4
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ +1.0/3.0, 0.0, 0.0, // 8
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 11y
- 0.0,-1.0/9.0, 0.0, // 3
- 0.0,-1.0/9.0, 0.0, // 4
- 0.0,-1.0/9.0, 0.0, // 5
0.0, 0.0, 0.0,
- 0.0,+1.0/9.0, 0.0, // 7
- 0.0,+1.0/9.0, 0.0, // 8
- 0.0,+1.0/9.0, 0.0, // 9
+ 0.0,-1.0/3.0, 0.0, // 4
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ 0.0,+1.0/3.0, 0.0, // 8
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 11z
- 0.0, 0.0,-1.0/9.0, // 3
- 0.0, 0.0,-1.0/9.0, // 4
- 0.0, 0.0,-1.0/9.0, // 5
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0/9.0, // 7
- 0.0, 0.0,+1.0/9.0, // 8
- 0.0, 0.0,+1.0/9.0, // 9
+ 0.0, 0.0,-1.0/3.0, // 4
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ 0.0, 0.0,+1.0/3.0, // 8
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 12x
- -1.0/9.0, 0.0, 0.0, // 3
- -1.0/9.0, 0.0, 0.0, // 4
- -1.0/9.0, 0.0, 0.0, // 5
0.0, 0.0, 0.0,
- +1.0/9.0, 0.0, 0.0, // 7
- +1.0/9.0, 0.0, 0.0, // 8
- +1.0/9.0, 0.0, 0.0, // 9
0.0, 0.0, 0.0,
+ -1.0/3.0, 0.0, 0.0, // 5
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ +1.0/3.0, 0.0, 0.0, // 9
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 12y
- 0.0,-1.0/9.0, 0.0, // 3
- 0.0,-1.0/9.0, 0.0, // 4
- 0.0,-1.0/9.0, 0.0, // 5
0.0, 0.0, 0.0,
- 0.0,+1.0/9.0, 0.0, // 7
- 0.0,+1.0/9.0, 0.0, // 8
- 0.0,+1.0/9.0, 0.0, // 9
0.0, 0.0, 0.0,
+ 0.0,-1.0/3.0, 0.0, // 5
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0,+1.0/3.0, 0.0, // 9
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 12z
- 0.0, 0.0,-1.0/9.0, // 3
- 0.0, 0.0,-1.0/9.0, // 4
- 0.0, 0.0,-1.0/9.0, // 5
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0/9.0, // 7
- 0.0, 0.0,+1.0/9.0, // 8
- 0.0, 0.0,+1.0/9.0, // 9
0.0, 0.0, 0.0,
+ 0.0, 0.0,-1.0/3.0, // 5
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0,+1.0/3.0, // 9
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
};
const double pylith::faults::CohesiveKinDataTet4::_fieldIncrAdjusted[] = {
Modified: short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTet4e.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTet4e.cc 2011-10-07 00:28:44 UTC (rev 19037)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTet4e.cc 2011-10-07 00:32:19 UTC (rev 19038)
@@ -44,24 +44,34 @@
const int pylith::faults::CohesiveKinDataTet4e::_numBasis = 3;
-const int pylith::faults::CohesiveKinDataTet4e::_numQuadPts = 1;
+const int pylith::faults::CohesiveKinDataTet4e::_numQuadPts = 3;
const double pylith::faults::CohesiveKinDataTet4e::_quadPts[] = {
- -3.33333333e-01, -3.33333333e-01,
+ -1.00000000e+00, -1.00000000e+00,
+ 1.00000000e+00, -1.00000000e+00,
+ -1.00000000e+00, 1.00000000e+00,
};
const double pylith::faults::CohesiveKinDataTet4e::_quadWts[] = {
- 2.0,
+ 2.0/3.0, 2.0/3.0, 2.0/3.0,
};
const double pylith::faults::CohesiveKinDataTet4e::_basis[] = {
- 3.33333333e-01, 3.33333333e-01,
- 3.33333333e-01,};
+ 1.0, 0.0, 0.0,
+ 0.0, 1.0, 0.0,
+ 0.0, 0.0, 1.0,
+};
const double pylith::faults::CohesiveKinDataTet4e::_basisDeriv[] = {
-0.50000000e+00, -0.50000000e+00,
0.50000000e+00, 0.00000000e+00,
0.00000000e+00, 0.50000000e+00,
+ -0.50000000e+00, -0.50000000e+00,
+ 0.50000000e+00, 0.00000000e+00,
+ 0.00000000e+00, 0.50000000e+00,
+ -0.50000000e+00, -0.50000000e+00,
+ 0.50000000e+00, 0.00000000e+00,
+ 0.00000000e+00, 0.50000000e+00,
};
const double pylith::faults::CohesiveKinDataTet4e::_verticesRef[] = {
@@ -128,10 +138,10 @@
4.9, 6.9, 8.9, // 11
3.1, 5.1, 7.1, // 12
3.3, 5.3, 7.3, // 13
- 2.0/3.0, 2.0/3.0, 2.0/3.0, // 14
- 1.0/3.0, 1.0/3.0, 1.0/3.0, // 15
- 2.0/3.0, 2.0/3.0, 2.0/3.0, // 16
- 1.0/3.0, 1.0/3.0, 1.0/3.0, // 17
+ 1.0, 1.0, 1.0, // 14
+ 1.0, 1.0, 1.0, // 15
+ 1.0, 1.0, 1.0, // 16
+ 1.0, 1.0, 1.0, // 17
};
@@ -175,124 +185,60 @@
const double pylith::faults::CohesiveKinDataTet4e::_residual[] = {
0.0, 0.0, 0.0,
- +(3.8+3.0+4.2 + 3.8+4.2+4.4)/9.0,
- +(5.8+5.0+6.2 + 5.8+6.2+6.4)/9.0,
- +(7.8+7.0+8.2 + 7.8+8.2+8.4)/9.0, // 5
- +(3.8+3.0+4.2)/9.0, +(5.8+5.0+6.2)/9.0, +(7.8+7.0+8.2)/9.0, // 6
- +(3.8+3.0+4.2 + 3.8+4.2+4.4)/9.0,
- +(5.8+5.0+6.2 + 5.8+6.2+6.4)/9.0,
- +(7.8+7.0+8.2 + 7.8+8.2+8.4)/9.0, // 5
- +(3.8+4.2+4.4)/9.0, +(5.8+6.2+6.4)/9.0, +(7.8+8.2+8.4)/9.0, // 8
+ +2.0/3.0*3.8, +2.0/3.0*5.8, +2.0/3.0*7.8, // 5
+ +1.0/3.0*3.0, +1.0/3.0*5.0, +1.0/3.0*7.0, // 6
+ +2.0/3.0*4.2, +2.0/3.0*6.2, +2.0/3.0*8.2, // 7
+ +1.0/3.0*4.4, +1.0/3.0*6.4, +1.0/3.0*8.4, // 8
0.0, 0.0, 0.0,
- -(3.8+3.0+4.2 + 3.8+4.2+4.4)/9.0,
- -(5.8+5.0+6.2 + 5.8+6.2+6.4)/9.0,
- -(7.8+7.0+8.2 + 7.8+8.2+8.4)/9.0, // 10
- -(3.8+3.0+4.2)/9.0, -(5.8+5.0+6.2)/9.0, -(7.8+7.0+8.2)/9.0, // 11
- -(3.8+3.0+4.2 + 3.8+4.2+4.4)/9.0,
- -(5.8+5.0+6.2 + 5.8+6.2+6.4)/9.0,
- -(7.8+7.0+8.2 + 7.8+8.2+8.4)/9.0, // 12
- -(3.8+4.2+4.4)/9.0, -(5.8+6.2+6.4)/9.0, -(7.8+8.2+8.4)/9.0, // 13
+ -2.0/3.0*3.8, -2.0/3.0*5.8, -2.0/3.0*7.8, // 10
+ -1.0/3.0*3.0, -1.0/3.0*5.0, -1.0/3.0*7.0, // 11
+ -2.0/3.0*4.2, -2.0/3.0*6.2, -2.0/3.0*8.2, // 12
+ -1.0/3.0*4.4, -1.0/3.0*6.4, -1.0/3.0*8.4, // 13
- -(3.7-3.2 + 3.9-3.3 + 4.1-3.4)/9.0 +
- -(3.7-3.2 + 4.1-3.4 + 4.3-3.5)/9.0 +
- +(0.07938069066 + 0.14140241667 + 0.18205179147)/9.0 +
- +(0.07938069066 + 0.18205179147 + 0.19904410828)/9.0,
- -(5.7-5.2 + 5.9-5.3 + 6.1-5.4)/9.0 +
- -(5.7-5.2 + 6.1-5.4 + 6.3-5.5)/9.0 +
- +(1.82575588523 + 1.69682900001 + 1.51709826228)/9.0 +
- +(1.82575588523 + 1.51709826228 + 1.29378670385)/9.0,
- -(7.7-7.2 + 7.9-7.3 + 8.1-7.4)/9.0 +
- -(7.7-7.2 + 8.1-7.4 + 8.3-7.5)/9.0 +
- +(-0.55566483464 + -0.56560966667 + -0.54615537442)/9.0 +
- +(-0.55566483464 + -0.54615537442 + -0.49761027071)/9.0, // 14
+ -2.0/3.0*(3.7-3.2 + -0.07938069066),
+ -2.0/3.0*(5.7-5.2 + -1.82575588523),
+ -2.0/3.0*(7.7-7.2 + 0.55566483464), // 14
- -(3.7-3.2 + 3.9-3.3 + 4.1-3.4)/9.0 +
- +(0.07938069066 + 0.14140241667 + 0.18205179147)/9.0,
- -(5.7-5.2 + 5.9-5.3 + 6.1-5.4)/9.0 +
- +(1.82575588523 + 1.69682900001 + 1.51709826228)/9.0,
- -(7.7-7.2 + 7.9-7.3 + 8.1-7.4)/9.0 +
- +(-0.55566483464 + -0.56560966667 + -0.54615537442)/9.0, // 15
+ -1.0/3.0*(3.9-3.3 + -0.14140241667),
+ -1.0/3.0*(5.9-5.3 + -1.69682900001),
+ -1.0/3.0*(7.9-7.3 + 0.56560966667), // 15
- -(3.7-3.2 + 3.9-3.3 + 4.1-3.4)/9.0 +
- -(3.7-3.2 + 4.1-3.4 + 4.3-3.5)/9.0 +
- +(0.07938069066 + 0.14140241667 + 0.18205179147)/9.0 +
- +(0.07938069066 + 0.18205179147 + 0.19904410828)/9.0,
- -(5.7-5.2 + 5.9-5.3 + 6.1-5.4)/9.0 +
- -(5.7-5.2 + 6.1-5.4 + 6.3-5.5)/9.0 +
- +(1.82575588523 + 1.69682900001 + 1.51709826228)/9.0 +
- +(1.82575588523 + 1.51709826228 + 1.29378670385)/9.0,
- -(7.7-7.2 + 7.9-7.3 + 8.1-7.4)/9.0 +
- -(7.7-7.2 + 8.1-7.4 + 8.3-7.5)/9.0 +
- +(-0.55566483464 + -0.56560966667 + -0.54615537442)/9.0 +
- +(-0.55566483464 + -0.54615537442 + -0.49761027071)/9.0, // 16
+ -2.0/3.0*(4.1-3.4 + -0.18205179147),
+ -2.0/3.0*(6.1-5.4 + -1.51709826228),
+ -2.0/3.0*(8.1-7.4 + 0.54615537442), // 16
- -(3.7-3.2 + 4.1-3.4 + 4.3-3.5)/9.0 +
- +(0.07938069066 + 0.18205179147 + 0.19904410828)/9.0,
- -(5.7-5.2 + 6.1-5.4 + 6.3-5.5)/9.0 +
- +(1.82575588523 + 1.51709826228 + 1.29378670385)/9.0,
- -(7.7-7.2 + 8.1-7.4 + 8.3-7.5)/9.0 +
- +(-0.55566483464 + -0.54615537442 + -0.49761027071)/9.0, // 17
+ -1.0/3.0*(4.3-3.5 + -0.19904410828),
+ -1.0/3.0*(6.3-5.5 + -1.29378670385),
+ -1.0/3.0*(8.3-7.5 + 0.49761027071), // 17
};
const double pylith::faults::CohesiveKinDataTet4e::_residualIncr[] = {
0.0, 0.0, 0.0,
- +(3.8+3.0+4.2 + 3.8+4.2+4.4)/9.0,
- +(5.8+5.0+6.2 + 5.8+6.2+6.4)/9.0,
- +(7.8+7.0+8.2 + 7.8+8.2+8.4)/9.0, // 5
- +(3.8+3.0+4.2)/9.0, +(5.8+5.0+6.2)/9.0, +(7.8+7.0+8.2)/9.0, // 6
- +(3.8+3.0+4.2 + 3.8+4.2+4.4)/9.0,
- +(5.8+5.0+6.2 + 5.8+6.2+6.4)/9.0,
- +(7.8+7.0+8.2 + 7.8+8.2+8.4)/9.0, // 5
- +(3.8+4.2+4.4)/9.0, +(5.8+6.2+6.4)/9.0, +(7.8+8.2+8.4)/9.0, // 8
+ +2.0/3.0*3.8, +2.0/3.0*5.8, +2.0/3.0*7.8, // 5
+ +1.0/3.0*3.0, +1.0/3.0*5.0, +1.0/3.0*7.0, // 6
+ +2.0/3.0*4.2, +2.0/3.0*6.2, +2.0/3.0*8.2, // 7
+ +1.0/3.0*4.4, +1.0/3.0*6.4, +1.0/3.0*8.4, // 8
0.0, 0.0, 0.0,
- -(3.8+3.0+4.2 + 3.8+4.2+4.4)/9.0,
- -(5.8+5.0+6.2 + 5.8+6.2+6.4)/9.0,
- -(7.8+7.0+8.2 + 7.8+8.2+8.4)/9.0, // 10
- -(3.8+3.0+4.2)/9.0, -(5.8+5.0+6.2)/9.0, -(7.8+7.0+8.2)/9.0, // 11
- -(3.8+3.0+4.2 + 3.8+4.2+4.4)/9.0,
- -(5.8+5.0+6.2 + 5.8+6.2+6.4)/9.0,
- -(7.8+7.0+8.2 + 7.8+8.2+8.4)/9.0, // 12
- -(3.8+4.2+4.4)/9.0, -(5.8+6.2+6.4)/9.0, -(7.8+8.2+8.4)/9.0, // 13
+ -2.0/3.0*3.8, -2.0/3.0*5.8, -2.0/3.0*7.8, // 10
+ -1.0/3.0*3.0, -1.0/3.0*5.0, -1.0/3.0*7.0, // 11
+ -2.0/3.0*4.2, -2.0/3.0*6.2, -2.0/3.0*8.2, // 12
+ -1.0/3.0*4.4, -1.0/3.0*6.4, -1.0/3.0*8.4, // 13
- -(3.7-3.2 + 3.9-3.3 + 4.1-3.4)/9.0 +
- -(3.7-3.2 + 4.1-3.4 + 4.3-3.5)/9.0 +
- +(0.07938069066 + 0.14140241667 + 0.18205179147)/9.0 +
- +(0.07938069066 + 0.18205179147 + 0.19904410828)/9.0,
- -(5.7-5.2 + 5.9-5.3 + 6.1-5.4)/9.0 +
- -(5.7-5.2 + 6.1-5.4 + 6.3-5.5)/9.0 +
- +(1.82575588523 + 1.69682900001 + 1.51709826228)/9.0 +
- +(1.82575588523 + 1.51709826228 + 1.29378670385)/9.0,
- -(7.7-7.2 + 7.9-7.3 + 8.1-7.4)/9.0 +
- -(7.7-7.2 + 8.1-7.4 + 8.3-7.5)/9.0 +
- +(-0.55566483464 + -0.56560966667 + -0.54615537442)/9.0 +
- +(-0.55566483464 + -0.54615537442 + -0.49761027071)/9.0, // 14
+ -2.0/3.0*(3.7-3.2 + -0.07938069066),
+ -2.0/3.0*(5.7-5.2 + -1.82575588523),
+ -2.0/3.0*(7.7-7.2 + 0.55566483464), // 14
- -(3.7-3.2 + 3.9-3.3 + 4.1-3.4)/9.0 +
- +(0.07938069066 + 0.14140241667 + 0.18205179147)/9.0,
- -(5.7-5.2 + 5.9-5.3 + 6.1-5.4)/9.0 +
- +(1.82575588523 + 1.69682900001 + 1.51709826228)/9.0,
- -(7.7-7.2 + 7.9-7.3 + 8.1-7.4)/9.0 +
- +(-0.55566483464 + -0.56560966667 + -0.54615537442)/9.0, // 15
+ -1.0/3.0*(3.9-3.3 + -0.14140241667),
+ -1.0/3.0*(5.9-5.3 + -1.69682900001),
+ -1.0/3.0*(7.9-7.3 + 0.56560966667), // 15
- -(3.7-3.2 + 3.9-3.3 + 4.1-3.4)/9.0 +
- -(3.7-3.2 + 4.1-3.4 + 4.3-3.5)/9.0 +
- +(0.07938069066 + 0.14140241667 + 0.18205179147)/9.0 +
- +(0.07938069066 + 0.18205179147 + 0.19904410828)/9.0,
- -(5.7-5.2 + 5.9-5.3 + 6.1-5.4)/9.0 +
- -(5.7-5.2 + 6.1-5.4 + 6.3-5.5)/9.0 +
- +(1.82575588523 + 1.69682900001 + 1.51709826228)/9.0 +
- +(1.82575588523 + 1.51709826228 + 1.29378670385)/9.0,
- -(7.7-7.2 + 7.9-7.3 + 8.1-7.4)/9.0 +
- -(7.7-7.2 + 8.1-7.4 + 8.3-7.5)/9.0 +
- +(-0.55566483464 + -0.56560966667 + -0.54615537442)/9.0 +
- +(-0.55566483464 + -0.54615537442 + -0.49761027071)/9.0, // 16
+ -2.0/3.0*(4.1-3.4 + -0.18205179147),
+ -2.0/3.0*(6.1-5.4 + -1.51709826228),
+ -2.0/3.0*(8.1-7.4 + 0.54615537442), // 16
- -(3.7-3.2 + 4.1-3.4 + 4.3-3.5)/9.0 +
- +(0.07938069066 + 0.18205179147 + 0.19904410828)/9.0,
- -(5.7-5.2 + 6.1-5.4 + 6.3-5.5)/9.0 +
- +(1.82575588523 + 1.51709826228 + 1.29378670385)/9.0,
- -(7.7-7.2 + 8.1-7.4 + 8.3-7.5)/9.0 +
- +(-0.55566483464 + -0.54615537442 + -0.49761027071)/9.0, // 17
+ -1.0/3.0*(4.3-3.5 + -0.19904410828),
+ -1.0/3.0*(6.3-5.5 + -1.29378670385),
+ -1.0/3.0*(8.3-7.5 + 0.49761027071), // 17
};
const double pylith::faults::CohesiveKinDataTet4e::_jacobian[] = {
@@ -348,10 +294,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -2.0/9.0, 0.0, 0.0, // 14
- -1.0/9.0, 0.0, 0.0, // 15
- -2.0/9.0, 0.0, 0.0, // 16
- -1.0/9.0, 0.0, 0.0, // 17
+ -2.0/3.0, 0.0, 0.0, // 14
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 5y
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -362,10 +308,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-2.0/9.0, 0.0, // 14
- 0.0,-1.0/9.0, 0.0, // 15
- 0.0,-2.0/9.0, 0.0, // 16
- 0.0,-1.0/9.0, 0.0, // 17
+ 0.0,-2.0/3.0, 0.0, // 14
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 5z
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -376,10 +322,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-2.0/9.0, // 14
- 0.0, 0.0,-1.0/9.0, // 15
- 0.0, 0.0,-2.0/9.0, // 16
- 0.0, 0.0,-1.0/9.0, // 17
+ 0.0, 0.0,-2.0/3.0, // 14
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 6x
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -390,10 +336,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0/9.0, 0.0, 0.0, // 14
- -1.0/9.0, 0.0, 0.0, // 15
- -1.0/9.0, 0.0, 0.0, // 16
0.0, 0.0, 0.0,
+ -1.0/3.0, 0.0, 0.0, // 15
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 6y
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -404,10 +350,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0/9.0, 0.0, // 14
- 0.0,-1.0/9.0, 0.0, // 15
- 0.0,-1.0/9.0, 0.0, // 16
0.0, 0.0, 0.0,
+ 0.0,-1.0/3.0, 0.0, // 15
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 6z
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -418,10 +364,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0/9.0, // 14
- 0.0, 0.0,-1.0/9.0, // 15
- 0.0, 0.0,-1.0/9.0, // 16
0.0, 0.0, 0.0,
+ 0.0, 0.0,-1.0/3.0, // 15
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 7x
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -432,10 +378,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -2.0/9.0, 0.0, 0.0, // 14
- -1.0/9.0, 0.0, 0.0, // 15
- -2.0/9.0, 0.0, 0.0, // 16
- -1.0/9.0, 0.0, 0.0, // 17
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ -2.0/3.0, 0.0, 0.0, // 16
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 7y
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -446,10 +392,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-2.0/9.0, 0.0, // 14
- 0.0,-1.0/9.0, 0.0, // 15
- 0.0,-2.0/9.0, 0.0, // 16
- 0.0,-1.0/9.0, 0.0, // 17
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0,-2.0/3.0, 0.0, // 16
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 7z
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -460,10 +406,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-2.0/9.0, // 14
- 0.0, 0.0,-1.0/9.0, // 15
- 0.0, 0.0,-2.0/9.0, // 16
- 0.0, 0.0,-1.0/9.0, // 17
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0,-2.0/3.0, // 16
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 8x
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -474,10 +420,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0/9.0, 0.0, 0.0, // 14
0.0, 0.0, 0.0,
- -1.0/9.0, 0.0, 0.0, // 16
- -1.0/9.0, 0.0, 0.0, // 17
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ -1.0/3.0, 0.0, 0.0, // 17
0.0, 0.0, 0.0, // 8y
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -488,10 +434,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0/9.0, 0.0, // 14
0.0, 0.0, 0.0,
- 0.0,-1.0/9.0, 0.0, // 16
- 0.0,-1.0/9.0, 0.0, // 17
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0,-1.0/3.0, 0.0, // 17
0.0, 0.0, 0.0, // 8z
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -502,10 +448,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0/9.0, // 14
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0/9.0, // 16
- 0.0, 0.0,-1.0/9.0, // 17
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0,-1.0/3.0, // 17
0.0, 0.0, 0.0, // 9x
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -558,10 +504,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +2.0/9.0, 0.0, 0.0, // 14
- +1.0/9.0, 0.0, 0.0, // 15
- +2.0/9.0, 0.0, 0.0, // 16
- +1.0/9.0, 0.0, 0.0, // 17
+ +2.0/3.0, 0.0, 0.0, // 14
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 10y
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -572,10 +518,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+2.0/9.0, 0.0, // 14
- 0.0,+1.0/9.0, 0.0, // 15
- 0.0,+2.0/9.0, 0.0, // 16
- 0.0,+1.0/9.0, 0.0, // 17
+ 0.0,+2.0/3.0, 0.0, // 14
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 10z
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -586,10 +532,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+2.0/9.0, // 14
- 0.0, 0.0,+1.0/9.0, // 15
- 0.0, 0.0,+2.0/9.0, // 16
- 0.0, 0.0,+1.0/9.0, // 17
+ 0.0, 0.0,+2.0/3.0, // 14
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 11x
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -600,10 +546,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +1.0/9.0, 0.0, 0.0, // 14
- +1.0/9.0, 0.0, 0.0, // 15
- +1.0/9.0, 0.0, 0.0, // 16
0.0, 0.0, 0.0,
+ +1.0/3.0, 0.0, 0.0, // 15
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 11y
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -614,10 +560,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0/9.0, 0.0, // 14
- 0.0,+1.0/9.0, 0.0, // 15
- 0.0,+1.0/9.0, 0.0, // 16
0.0, 0.0, 0.0,
+ 0.0,+1.0/3.0, 0.0, // 15
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 11z
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -628,10 +574,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0/9.0, // 14
- 0.0, 0.0,+1.0/9.0, // 15
- 0.0, 0.0,+1.0/9.0, // 16
0.0, 0.0, 0.0,
+ 0.0, 0.0,+1.0/3.0, // 15
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 12x
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -642,10 +588,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +2.0/9.0, 0.0, 0.0, // 14
- +1.0/9.0, 0.0, 0.0, // 15
- +2.0/9.0, 0.0, 0.0, // 16
- +1.0/9.0, 0.0, 0.0, // 17
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ +2.0/3.0, 0.0, 0.0, // 16
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 12y
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -656,10 +602,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+2.0/9.0, 0.0, // 14
- 0.0,+1.0/9.0, 0.0, // 15
- 0.0,+2.0/9.0, 0.0, // 16
- 0.0,+1.0/9.0, 0.0, // 17
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0,+2.0/3.0, 0.0, // 16
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 12z
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -670,10 +616,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+2.0/9.0, // 14
- 0.0, 0.0,+1.0/9.0, // 15
- 0.0, 0.0,+2.0/9.0, // 16
- 0.0, 0.0,+1.0/9.0, // 17
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0,+2.0/3.0, // 16
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 13x
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -684,10 +630,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +1.0/9.0, 0.0, 0.0, // 14
0.0, 0.0, 0.0,
- +1.0/9.0, 0.0, 0.0, // 16
- +1.0/9.0, 0.0, 0.0, // 17
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ +1.0/3.0, 0.0, 0.0, // 17
0.0, 0.0, 0.0, // 13y
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -698,10 +644,10 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0/9.0, 0.0, // 14
0.0, 0.0, 0.0,
- 0.0,+1.0/9.0, 0.0, // 16
- 0.0,+1.0/9.0, 0.0, // 17
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0,+1.0/3.0, 0.0, // 17
0.0, 0.0, 0.0, // 13z
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -712,178 +658,178 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0/9.0, // 14
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0/9.0, // 16
- 0.0, 0.0,+1.0/9.0, // 17
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0,+1.0/3.0, // 17
0.0, 0.0, 0.0, // 14x
- -2.0/9.0, 0.0, 0.0, // 5
- -1.0/9.0, 0.0, 0.0, // 6
- -2.0/9.0, 0.0, 0.0, // 7
- -1.0/9.0, 0.0, 0.0, // 8
+ -2.0/3.0, 0.0, 0.0, // 5
0.0, 0.0, 0.0,
- +2.0/9.0, 0.0, 0.0, // 10
- +1.0/9.0, 0.0, 0.0, // 11
- +2.0/9.0, 0.0, 0.0, // 12
- +1.0/9.0, 0.0, 0.0, // 13
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ +2.0/3.0, 0.0, 0.0, // 10
0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 14y
- 0.0,-2.0/9.0, 0.0, // 5
- 0.0,-1.0/9.0, 0.0, // 6
- 0.0,-2.0/9.0, 0.0, // 7
- 0.0,-1.0/9.0, 0.0, // 8
+ 0.0,-2.0/3.0, 0.0, // 5
0.0, 0.0, 0.0,
- 0.0,+2.0/9.0, 0.0, // 10
- 0.0,+1.0/9.0, 0.0, // 11
- 0.0,+2.0/9.0, 0.0, // 12
- 0.0,+1.0/9.0, 0.0, // 13
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ 0.0,+2.0/3.0, 0.0, // 10
0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 14z
- 0.0, 0.0,-2.0/9.0, // 5
- 0.0, 0.0,-1.0/9.0, // 6
- 0.0, 0.0,-2.0/9.0, // 7
- 0.0, 0.0,-1.0/9.0, // 8
+ 0.0, 0.0,-2.0/3.0, // 5
0.0, 0.0, 0.0,
- 0.0, 0.0,+2.0/9.0, // 10
- 0.0, 0.0,+1.0/9.0, // 11
- 0.0, 0.0,+2.0/9.0, // 12
- 0.0, 0.0,+1.0/9.0, // 13
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ 0.0, 0.0,+2.0/3.0, // 10
0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 15x
- -1.0/9.0, 0.0, 0.0, // 5
- -1.0/9.0, 0.0, 0.0, // 6
- -1.0/9.0, 0.0, 0.0, // 7
0.0, 0.0, 0.0,
+ -1.0/3.0, 0.0, 0.0, // 6
0.0, 0.0, 0.0,
- +1.0/9.0, 0.0, 0.0, // 10
- +1.0/9.0, 0.0, 0.0, // 11
- +1.0/9.0, 0.0, 0.0, // 12
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ +1.0/3.0, 0.0, 0.0, // 11
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 15y
- 0.0,-1.0/9.0, 0.0, // 5
- 0.0,-1.0/9.0, 0.0, // 6
- 0.0,-1.0/9.0, 0.0, // 7
0.0, 0.0, 0.0,
+ 0.0,-1.0/3.0, 0.0, // 6
0.0, 0.0, 0.0,
- 0.0,+1.0/9.0, 0.0, // 10
- 0.0,+1.0/9.0, 0.0, // 11
- 0.0,+1.0/9.0, 0.0, // 12
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ 0.0,+1.0/3.0, 0.0, // 11
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 15z
- 0.0, 0.0,-1.0/9.0, // 5
- 0.0, 0.0,-1.0/9.0, // 6
- 0.0, 0.0,-1.0/9.0, // 7
0.0, 0.0, 0.0,
+ 0.0, 0.0,-1.0/3.0, // 6
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0/9.0, // 10
- 0.0, 0.0,+1.0/9.0, // 11
- 0.0, 0.0,+1.0/9.0, // 12
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ 0.0, 0.0,+1.0/3.0, // 11
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 16x
- -2.0/9.0, 0.0, 0.0, // 5
- -1.0/9.0, 0.0, 0.0, // 6
- -2.0/9.0, 0.0, 0.0, // 7
- -1.0/9.0, 0.0, 0.0, // 8
0.0, 0.0, 0.0,
- +2.0/9.0, 0.0, 0.0, // 10
- +1.0/9.0, 0.0, 0.0, // 11
- +2.0/9.0, 0.0, 0.0, // 12
- +1.0/9.0, 0.0, 0.0, // 13
0.0, 0.0, 0.0,
+ -2.0/3.0, 0.0, 0.0, // 7
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ +2.0/3.0, 0.0, 0.0, // 12
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 16y
- 0.0,-2.0/9.0, 0.0, // 5
- 0.0,-1.0/9.0, 0.0, // 6
- 0.0,-2.0/9.0, 0.0, // 7
- 0.0,-1.0/9.0, 0.0, // 8
0.0, 0.0, 0.0,
- 0.0,+2.0/9.0, 0.0, // 10
- 0.0,+1.0/9.0, 0.0, // 11
- 0.0,+2.0/9.0, 0.0, // 12
- 0.0,+1.0/9.0, 0.0, // 13
0.0, 0.0, 0.0,
+ 0.0,-2.0/3.0, 0.0, // 7
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0,+2.0/3.0, 0.0, // 12
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 16z
- 0.0, 0.0,-2.0/9.0, // 5
- 0.0, 0.0,-1.0/9.0, // 6
- 0.0, 0.0,-2.0/9.0, // 7
- 0.0, 0.0,-1.0/9.0, // 8
0.0, 0.0, 0.0,
- 0.0, 0.0,+2.0/9.0, // 10
- 0.0, 0.0,+1.0/9.0, // 11
- 0.0, 0.0,+2.0/9.0, // 12
- 0.0, 0.0,+1.0/9.0, // 13
0.0, 0.0, 0.0,
+ 0.0, 0.0,-2.0/3.0, // 7
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0,+2.0/3.0, // 12
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0,// 17x
- -1.0/9.0, 0.0, 0.0, // 5
0.0, 0.0, 0.0,
- -1.0/9.0, 0.0, 0.0, // 7
- -1.0/9.0, 0.0, 0.0, // 8
0.0, 0.0, 0.0,
- +1.0/9.0, 0.0, 0.0, // 10
0.0, 0.0, 0.0,
- +1.0/9.0, 0.0, 0.0, // 12
- +1.0/9.0, 0.0, 0.0, // 13
+ -1.0/3.0, 0.0, 0.0, // 8
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ +1.0/3.0, 0.0, 0.0, // 13
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 17y
- 0.0,-1.0/9.0, 0.0, // 5
0.0, 0.0, 0.0,
- 0.0,-1.0/9.0, 0.0, // 7
- 0.0,-1.0/9.0, 0.0, // 8
0.0, 0.0, 0.0,
- 0.0,+1.0/9.0, 0.0, // 10
0.0, 0.0, 0.0,
- 0.0,+1.0/9.0, 0.0, // 12
- 0.0,+1.0/9.0, 0.0, // 13
+ 0.0,-1.0/3.0, 0.0, // 8
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ 0.0,+1.0/3.0, 0.0, // 13
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 17z
- 0.0, 0.0,-1.0/9.0, // 5
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0/9.0, // 7
- 0.0, 0.0,-1.0/9.0, // 8
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0/9.0, // 10
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0/9.0, // 12
- 0.0, 0.0,+1.0/9.0, // 13
+ 0.0, 0.0,-1.0/3.0, // 8
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
+ 0.0, 0.0,+1.0/3.0, // 13
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0,
};
pylith::faults::CohesiveKinDataTet4e::CohesiveKinDataTet4e(void)
Modified: short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTri3.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTri3.cc 2011-10-07 00:28:44 UTC (rev 19037)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTri3.cc 2011-10-07 00:32:19 UTC (rev 19038)
@@ -61,24 +61,24 @@
const int pylith::faults::CohesiveKinDataTri3::_numBasis = 2;
-const int pylith::faults::CohesiveKinDataTri3::_numQuadPts = 1;
+const int pylith::faults::CohesiveKinDataTri3::_numQuadPts = 2;
const double pylith::faults::CohesiveKinDataTri3::_quadPts[] = {
- 0.0,
+ -1.0, 1.0,
};
const double pylith::faults::CohesiveKinDataTri3::_quadWts[] = {
- 2.0,
+ 1.0, 1.0
};
const double pylith::faults::CohesiveKinDataTri3::_basis[] = {
- 0.5,
- 0.5
+ 1.0, 0.0,
+ 0.0, 1.0,
};
const double pylith::faults::CohesiveKinDataTri3::_basisDeriv[] = {
- -0.5,
- 0.5
+ -0.5, 0.5,
+ -0.5, 0.5,
};
const double pylith::faults::CohesiveKinDataTri3::_verticesRef[] = {
@@ -166,28 +166,28 @@
const double pylith::faults::CohesiveKinDataTri3::_residual[] = {
0.0, 0.0,
- +8.7, +9.7, // 3
- +8.7, +9.7, // 4
+ +8.6, +9.6, // 3
+ +8.8, +9.8, // 4
0.0, 0.0,
- -8.7, -9.7, // 6
- -8.7, -9.7, // 7
- -0.5*(8.5-8.2 + 8.7-8.3) - 0.5*(0.08241148423+0.14794836271),
- -0.5*(9.5-9.2 + 9.7-9.3) - 0.5*(1.89546413727+1.77538035254), // 8
- -0.5*(8.5-8.2 + 8.7-8.3) - 0.5*(0.08241148423+0.14794836271),
- -0.5*(9.5-9.2 + 9.7-9.3) - 0.5*(1.89546413727+1.77538035254), // 9
+ -8.6, -9.6, // 6
+ -8.8, -9.8, // 7
+ -(8.5-8.2) - (0.08241148423),
+ -(9.5-9.2) - (1.89546413727), // 8
+ -(8.7-8.3) - (0.14794836271),
+ -(9.7-9.3) - (1.77538035254), // 9
};
const double pylith::faults::CohesiveKinDataTri3::_residualIncr[] = {
0.0, 0.0,
- +8.7, +9.7, // 3
- +8.7, +9.7, // 4
+ +8.6, +9.6, // 3
+ +8.8, +9.8, // 4
0.0, 0.0,
- -8.7, -9.7, // 6
- -8.7, -9.7, // 7
- -0.5*(8.5-8.2 + 8.7-8.3) - 0.5*(0.08241148423+0.14794836271),
- -0.5*(9.5-9.2 + 9.7-9.3) - 0.5*(1.89546413727+1.77538035254), // 8
- -0.5*(8.5-8.2 + 8.7-8.3) - 0.5*(0.08241148423+0.14794836271),
- -0.5*(9.5-9.2 + 9.7-9.3) - 0.5*(1.89546413727+1.77538035254), // 9
+ -8.6, -9.6, // 6
+ -8.8, -9.8, // 7
+ -(8.5-8.2) - (0.08241148423),
+ -(9.5-9.2) - (1.89546413727), // 8
+ -(8.7-8.3) - (0.14794836271),
+ -(9.7-9.3) - (1.77538035254), // 9
};
const double pylith::faults::CohesiveKinDataTri3::_jacobian[] = {
@@ -213,32 +213,32 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- -0.5, 0.0, // 8
- -0.5, 0.0, // 9
+ -1.0, 0.0, // 8
+ 0.0, 0.0,
0.0, 0.0, // 3y
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,-0.5, // 8
- 0.0,-0.5, // 9
+ 0.0,-1.0, // 8
+ 0.0, 0.0,
0.0, 0.0, // 4x
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- -0.5, 0.0, // 8
- -0.5, 0.0, // 9
+ 0.0, 0.0,
+ -1.0, 0.0, // 9
0.0, 0.0, // 4y
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,-0.5, // 8
- 0.0,-0.5, // 9
+ 0.0, 0.0,
+ 0.0,-1.0, // 9
0.0, 0.0, // 5x
0.0, 0.0,
0.0, 0.0,
@@ -261,66 +261,66 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- +0.5, 0.0, // 8
- +0.5, 0.0, // 9
+ +1.0, 0.0, // 8
+ 0.0, 0.0,
0.0, 0.0, // 6y
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,+0.5, // 8
- 0.0,+0.5, // 9
+ 0.0,+1.0, // 8
+ 0.0, 0.0,
0.0, 0.0, // 7x
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- +0.5, 0.0, // 8
- +0.5, 0.0, // 9
+ 0.0, 0.0,
+ +1.0, 0.0, // 9
0.0, 0.0, // 7y
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,+0.5, // 8
- 0.0,+0.5, // 9
+ 0.0, 0.0,
+ 0.0,+1.0, // 9
0.0, 0.0, // 8x
- -0.5, 0.0, // 3
- -0.5, 0.0, // 4
+ -1.0, 0.0, // 3
0.0, 0.0,
- +0.5, 0.0, // 6
- +0.5, 0.0, // 7
0.0, 0.0,
+ +1.0, 0.0, // 6
0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 8y
- 0.0,-0.5, // 3
- 0.0,-0.5, // 4
+ 0.0,-1.0, // 3
0.0, 0.0,
- 0.0,+0.5, // 6
- 0.0,+0.5, // 7
0.0, 0.0,
+ 0.0,+1.0, // 6
0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 9x
- -0.5, 0.0, // 3
- -0.5, 0.0, // 4
0.0, 0.0,
- +0.5, 0.0, // 6
- +0.5, 0.0, // 7
+ -1.0, 0.0, // 4
0.0, 0.0,
0.0, 0.0,
+ +1.0, 0.0, // 7
+ 0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 9y
- 0.0,-0.5, // 3
- 0.0,-0.5, // 4
0.0, 0.0,
- 0.0,+0.5, // 6
- 0.0,+0.5, // 7
+ 0.0,-1.0, // 4
0.0, 0.0,
0.0, 0.0,
+ 0.0,+1.0, // 7
+ 0.0, 0.0,
+ 0.0, 0.0,
};
const double pylith::faults::CohesiveKinDataTri3::_fieldIncrAdjusted[] = {
Modified: short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTri3d.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTri3d.cc 2011-10-07 00:28:44 UTC (rev 19037)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTri3d.cc 2011-10-07 00:32:19 UTC (rev 19038)
@@ -76,24 +76,24 @@
const int pylith::faults::CohesiveKinDataTri3d::_numBasis = 2;
-const int pylith::faults::CohesiveKinDataTri3d::_numQuadPts = 1;
+const int pylith::faults::CohesiveKinDataTri3d::_numQuadPts = 2;
const double pylith::faults::CohesiveKinDataTri3d::_quadPts[] = {
- 0.0,
+ -1.0, 1.0,
};
const double pylith::faults::CohesiveKinDataTri3d::_quadWts[] = {
- 2.0,
+ 1.0, 1.0
};
const double pylith::faults::CohesiveKinDataTri3d::_basis[] = {
- 0.5,
- 0.5
+ 1.0, 0.0,
+ 0.0, 1.0,
};
const double pylith::faults::CohesiveKinDataTri3d::_basisDeriv[] = {
- -0.5,
- 0.5
+ -0.5, 0.5,
+ -0.5, 0.5,
};
const double pylith::faults::CohesiveKinDataTri3d::_verticesRef[] = {
@@ -153,7 +153,7 @@
6.7, 8.7, // 10
6.9, 8.9, // 11
7.1, 9.1, // 12
- 2.0, 2.0, // 13
+ 1.0, 1.0, // 13
1.0, 1.0, // 14
1.0, 1.0, // 15
};
@@ -196,54 +196,38 @@
const double pylith::faults::CohesiveKinDataTri3d::_residual[] = {
0.0, 0.0,
- +0.5*(6.8+6.0 + 6.8+7.2), +0.5*(8.8+8.0 + 8.8+9.2), // 5
- +0.5*(6.8+6.0), +0.5*(8.8+8.0), // 6
+ +2.0*6.8, +2.0*8.8, // 5
+ +1.0*6.0, +1.0*8.0, // 6
0.0, 0.0,
- +0.5*(6.8+7.2), +0.5*(8.8+9.2), // 8
+ +1.0*7.2, +1.0*9.2, // 8
0.0, 0.0,
- -0.5*(6.8+6.0 + 6.8+7.2), -0.5*(8.8+8.0 + 8.8+9.2), // 10
- -0.5*(6.8+6.0), -0.5*(8.8+8.0), // 6
- -0.5*(6.8+7.2), -0.5*(8.8+9.2), // 8
- -0.5*(6.7-6.2+7.1-6.5 + 6.7-6.2+6.9-6.3
- -1.59887481971-0.70710678118654757*(1.89546413727-0.08241148423) +
- +0.14794836271-0.70710678118654757*(1.89546413727-0.08241148423)),
- -0.5*(8.7-8.2+9.1-8.5 + 8.7-8.2+8.9-8.3
- +0.19186497837-0.70710678118654757*(-1.89546413727-0.08241148423) +
- +1.77538035254-0.70710678118654757*(-1.89546413727-0.08241148423)), // 13
- -0.5*(6.7-6.2+6.9-6.3
- +0.14794836271-0.70710678118654757*(1.89546413727-0.08241148423)),
- -0.5*(8.7-8.2+8.9-8.3
- +1.77538035254-0.70710678118654757*(-1.89546413727-0.08241148423)), // 14
- -0.5*(6.7-6.2+7.1-6.5
- -1.59887481971-0.70710678118654757*(1.89546413727-0.08241148423)),
- -0.5*(8.7-8.2+9.1-8.5
- +0.19186497837-0.70710678118654757*(-1.89546413727-0.08241148423)), // 15
+ -2.0*6.8, -2.0*8.8, // 10
+ -1.0*6.0, -1.0*8.0, // 6
+ -1.0*7.2, -1.0*9.2, // 8
+ -2.0*(6.7-6.2 -0.70710678118654757*(1.89546413727-0.08241148423)),
+ -2.0*(8.7-8.2 -0.70710678118654757*(-1.89546413727-0.08241148423)), // 13
+ -1.0*(6.9-6.3 +0.14794836271),
+ -1.0*(8.9-8.3 +1.77538035254), // 14
+ -1.0*(7.1-6.5 -1.59887481971),
+ -1.0*(9.1-8.5 +0.19186497837), // 15
};
const double pylith::faults::CohesiveKinDataTri3d::_residualIncr[] = {
0.0, 0.0,
- +0.5*(6.8+6.0 + 6.8+7.2), +0.5*(8.8+8.0 + 8.8+9.2), // 5
- +0.5*(6.8+6.0), +0.5*(8.8+8.0), // 6
+ +2.0*6.8, +2.0*8.8, // 5
+ +1.0*6.0, +1.0*8.0, // 6
0.0, 0.0,
- +0.5*(6.8+7.2), +0.5*(8.8+9.2), // 8
+ +1.0*7.2, +1.0*9.2, // 8
0.0, 0.0,
- -0.5*(6.8+6.0 + 6.8+7.2), -0.5*(8.8+8.0 + 8.8+9.2), // 10
- -0.5*(6.8+6.0), -0.5*(8.8+8.0), // 6
- -0.5*(6.8+7.2), -0.5*(8.8+9.2), // 8
- -0.5*(6.7-6.2+7.1-6.5 + 6.7-6.2+6.9-6.3
- -1.59887481971-0.70710678118654757*(1.89546413727-0.08241148423) +
- +0.14794836271-0.70710678118654757*(1.89546413727-0.08241148423)),
- -0.5*(8.7-8.2+9.1-8.5 + 8.7-8.2+8.9-8.3
- +0.19186497837-0.70710678118654757*(-1.89546413727-0.08241148423) +
- +1.77538035254-0.70710678118654757*(-1.89546413727-0.08241148423)), // 13
- -0.5*(6.7-6.2+6.9-6.3
- +0.14794836271-0.70710678118654757*(1.89546413727-0.08241148423)),
- -0.5*(8.7-8.2+8.9-8.3
- +1.77538035254-0.70710678118654757*(-1.89546413727-0.08241148423)), // 14
- -0.5*(6.7-6.2+7.1-6.5
- -1.59887481971-0.70710678118654757*(1.89546413727-0.08241148423)),
- -0.5*(8.7-8.2+9.1-8.5
- +0.19186497837-0.70710678118654757*(-1.89546413727-0.08241148423)), // 15
+ -2.0*6.8, -2.0*8.8, // 10
+ -1.0*6.0, -1.0*8.0, // 6
+ -1.0*7.2, -1.0*9.2, // 8
+ -2.0*(6.7-6.2 -0.70710678118654757*(1.89546413727-0.08241148423)),
+ -2.0*(8.7-8.2 -0.70710678118654757*(-1.89546413727-0.08241148423)), // 13
+ -1.0*(6.9-6.3 +0.14794836271),
+ -1.0*(8.9-8.3 +1.77538035254), // 14
+ -1.0*(7.1-6.5 -1.59887481971),
+ -1.0*(9.1-8.5 +0.19186497837), // 15
};
const double pylith::faults::CohesiveKinDataTri3d::_jacobian[] = {
@@ -280,9 +264,9 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- -1.0, 0.0, // 13
- -0.5, 0.0, // 14
- -0.5, 0.0, // 15
+ -2.0, 0.0, // 13
+ 0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 5y
0.0, 0.0,
0.0, 0.0,
@@ -292,9 +276,9 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,-1.0, // 13
- 0.0,-0.5, // 14
- 0.0,-0.5, // 15
+ 0.0,-2.0, // 13
+ 0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 6x
0.0, 0.0,
0.0, 0.0,
@@ -304,9 +288,9 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- -0.5, 0.0, // 13
- -0.5, 0.0, // 14
0.0, 0.0,
+ -1.0, 0.0, // 14
+ 0.0, 0.0,
0.0, 0.0, // 6y
0.0, 0.0,
0.0, 0.0,
@@ -316,9 +300,9 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,-0.5, // 13
- 0.0,-0.5, // 14
0.0, 0.0,
+ 0.0,-1.0, // 14
+ 0.0, 0.0,
0.0, 0.0, // 7x
0.0, 0.0,
0.0, 0.0,
@@ -352,9 +336,9 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- -0.5, 0.0, // 13
0.0, 0.0,
- -0.5, 0.0, // 15
+ 0.0, 0.0,
+ -1.0, 0.0, // 15
0.0, 0.0, // 8y
0.0, 0.0,
0.0, 0.0,
@@ -364,9 +348,9 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,-0.5, // 13
0.0, 0.0,
- 0.0,-0.5, // 15
+ 0.0, 0.0,
+ 0.0,-1.0, // 15
0.0, 0.0, // 9x
0.0, 0.0,
0.0, 0.0,
@@ -400,9 +384,9 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- +1.0, 0.0, // 13
- +0.5, 0.0, // 14
- +0.5, 0.0, // 15
+ +2.0, 0.0, // 13
+ 0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 10y
0.0, 0.0,
0.0, 0.0,
@@ -412,9 +396,9 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,+1.0, // 13
- 0.0,+0.5, // 14
- 0.0,+0.5, // 15
+ 0.0,+2.0, // 13
+ 0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 11x
0.0, 0.0,
0.0, 0.0,
@@ -424,9 +408,9 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- +0.5, 0.0, // 13
- +0.5, 0.0, // 14
0.0, 0.0,
+ +1.0, 0.0, // 14
+ 0.0, 0.0,
0.0, 0.0, // 11y
0.0, 0.0,
0.0, 0.0,
@@ -436,9 +420,9 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,+0.5, // 13
- 0.0,+0.5, // 14
0.0, 0.0,
+ 0.0,+1.0, // 14
+ 0.0, 0.0,
0.0, 0.0, // 12x
0.0, 0.0,
0.0, 0.0,
@@ -448,9 +432,9 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- +0.5, 0.0, // 13
0.0, 0.0,
- +0.5, 0.0, // 15
+ 0.0, 0.0,
+ +1.0, 0.0, // 15
0.0, 0.0, // 12y
0.0, 0.0,
0.0, 0.0,
@@ -460,81 +444,81 @@
0.0, 0.0,
0.0, 0.0,
0.0, 0.0,
- 0.0,+0.5, // 13
0.0, 0.0,
- 0.0,+0.5, // 15
+ 0.0, 0.0,
+ 0.0,+1.0, // 15
0.0, 0.0, // 13x
- -1.0, 0.0, // 5
- -0.5, 0.0, // 6
+ -2.0, 0.0, // 5
0.0, 0.0,
- -0.5, 0.0, // 8
0.0, 0.0,
- +1.0, 0.0, // 10
- +0.5, 0.0, // 11
- +0.5, 0.0, // 12
0.0, 0.0,
0.0, 0.0,
+ +2.0, 0.0, // 10
0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 13y
- 0.0,-1.0, // 5
- 0.0,-0.5, // 6
+ 0.0,-2.0, // 5
0.0, 0.0,
- 0.0,-0.5, // 8
0.0, 0.0,
- 0.0,+1.0, // 10
- 0.0,+0.5, // 11
- 0.0,+0.5, // 12
0.0, 0.0,
0.0, 0.0,
+ 0.0,+2.0, // 10
0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 14x
- -0.5, 0.0, // 5
- -0.5, 0.0, // 6
0.0, 0.0,
+ -1.0, 0.0, // 6
0.0, 0.0,
0.0, 0.0,
- +0.5, 0.0, // 10
- +0.5, 0.0, // 11
0.0, 0.0,
0.0, 0.0,
+ +1.0, 0.0, // 11
0.0, 0.0,
0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 14y
- 0.0,-0.5, // 5
- 0.0,-0.5, // 6
0.0, 0.0,
+ 0.0,-1.0, // 6
0.0, 0.0,
0.0, 0.0,
- 0.0,+0.5, // 10
- 0.0,+0.5, // 11
0.0, 0.0,
0.0, 0.0,
+ 0.0,+1.0, // 11
0.0, 0.0,
0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 15x
- -0.5, 0.0, // 5
0.0, 0.0,
0.0, 0.0,
- -0.5, 0.0, // 8
0.0, 0.0,
- +0.5, 0.0, // 10
+ -1.0, 0.0, // 8
0.0, 0.0,
- +0.5, 0.0, // 12
0.0, 0.0,
0.0, 0.0,
+ +1.0, 0.0, // 12
0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
0.0, 0.0, // 15y
- 0.0,-0.5, // 5
0.0, 0.0,
0.0, 0.0,
- 0.0,-0.5, // 8
0.0, 0.0,
- 0.0,+0.5, // 10
+ 0.0,-1.0, // 8
0.0, 0.0,
- 0.0,+0.5, // 12
0.0, 0.0,
0.0, 0.0,
+ 0.0,+1.0, // 12
0.0, 0.0,
+ 0.0, 0.0,
+ 0.0, 0.0,
};
pylith::faults::CohesiveKinDataTri3d::CohesiveKinDataTri3d(void)
More information about the CIG-COMMITS
mailing list