[cig-commits] r18964 - in short/3D/PyLith/branches/v1.6-revisedfault: libsrc/pylith/faults unittests/libtests/faults/data
brad at geodynamics.org
brad at geodynamics.org
Thu Sep 22 17:11:07 PDT 2011
Author: brad
Date: 2011-09-22 17:11:07 -0700 (Thu, 22 Sep 2011)
New Revision: 18964
Modified:
short/3D/PyLith/branches/v1.6-revisedfault/libsrc/pylith/faults/FaultCohesiveLagrange.cc
short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataHex8.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/CohesiveKinDataTet4f.cc
Log:
More work on C++ test data for revised fault implementation.
Modified: short/3D/PyLith/branches/v1.6-revisedfault/libsrc/pylith/faults/FaultCohesiveLagrange.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/libsrc/pylith/faults/FaultCohesiveLagrange.cc 2011-09-22 22:14:21 UTC (rev 18963)
+++ short/3D/PyLith/branches/v1.6-revisedfault/libsrc/pylith/faults/FaultCohesiveLagrange.cc 2011-09-23 00:11:07 UTC (rev 18964)
@@ -406,6 +406,7 @@
<< ", slip: " << slipGlobalCell[jBasis*spaceDim+iDim]
<< ", dispP: " << dispTpdtCell[jBP + iDim]
<< ", dispN: " << dispTpdtCell[jBN + iDim]
+ << ", dispL: " << dispTpdtCell[jBL + iDim]
<< ", residualN: " << residualCell[iBN + iDim]
<< ", residualP: " << residualCell[iBP + iDim]
<< ", residualL: " << residualCell[iBL + iDim]
Modified: short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataHex8.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataHex8.cc 2011-09-22 22:14:21 UTC (rev 18963)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataHex8.cc 2011-09-23 00:11:07 UTC (rev 18964)
@@ -59,10 +59,10 @@
};
const double pylith::faults::CohesiveKinDataHex8::_basis[] = {
- 0.62200847, 0.16666667, 0.16666667, 0.0446582,
- 0.16666667, 0.62200847, 0.0446582, 0.16666667,
- 0.16666667, 0.0446582, 0.62200847, 0.16666667,
- 0.0446582, 0.16666667, 0.16666667, 0.62200847,
+ 0.62200847, 0.16666667, 0.0446582, 0.16666667,
+ 0.16666667, 0.62200847, 0.16666667, 0.0446582,
+ 0.0446582, 0.16666667, 0.62200847, 0.16666667,
+ 0.16666667, 0.0446582, 0.16666667, 0.62200847,
};
const double pylith::faults::CohesiveKinDataHex8::_basisDeriv[] = {
@@ -215,22 +215,204 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -9.4,-5.4,+7.4, // 6
- -9.6,-5.6,+7.6, // 7
- -9.8,-5.8,+7.8, // 8
- -9.0,-5.0,+7.0, // 9
+
+ // 6
+ -0.62200847*(0.62200847*5.4+0.16666667*5.6+0.04465820*5.0+0.16666667*5.8) +
+ -0.16666667*(0.16666667*5.4+0.62200847*5.6+0.16666667*5.0+0.04465820*5.8) +
+ -0.04465820*(0.04465820*5.4+0.16666667*5.6+0.62200847*5.0+0.16666667*5.8) +
+ -0.16666667*(0.16666667*5.4+0.04465820*5.6+0.16666667*5.0+0.62200847*5.8),
+
+ -0.62200847*(0.62200847*7.4+0.16666667*7.6+0.04465820*7.0+0.16666667*7.8) +
+ -0.16666667*(0.16666667*7.4+0.62200847*7.6+0.16666667*7.0+0.04465820*7.8) +
+ -0.04465820*(0.04465820*7.4+0.16666667*7.6+0.62200847*7.0+0.16666667*7.8) +
+ -0.16666667*(0.16666667*7.4+0.04465820*7.6+0.16666667*7.0+0.62200847*7.8),
+
+ -0.62200847*(0.62200847*9.4+0.16666667*9.6+0.04465820*9.0+0.16666667*9.8) +
+ -0.16666667*(0.16666667*9.4+0.62200847*9.6+0.16666667*9.0+0.04465820*9.8) +
+ -0.04465820*(0.04465820*9.4+0.16666667*9.6+0.62200847*9.0+0.16666667*9.8) +
+ -0.16666667*(0.16666667*9.4+0.04465820*9.6+0.16666667*9.0+0.62200847*9.8),
+
+ // 7
+ -0.62200847*(0.62200847*5.6+0.16666667*5.4+0.04465820*5.8+0.16666667*5.0) +
+ -0.16666667*(0.16666667*5.6+0.62200847*5.4+0.16666667*5.8+0.04465820*5.0) +
+ -0.04465820*(0.04465820*5.6+0.16666667*5.4+0.62200847*5.8+0.16666667*5.0) +
+ -0.16666667*(0.16666667*5.6+0.04465820*5.4+0.16666667*5.8+0.62200847*5.0),
+
+ -0.62200847*(0.62200847*7.6+0.16666667*7.4+0.04465820*7.8+0.16666667*7.0) +
+ -0.16666667*(0.16666667*7.6+0.62200847*7.4+0.16666667*7.8+0.04465820*7.0) +
+ -0.04465820*(0.04465820*7.6+0.16666667*7.4+0.62200847*7.8+0.16666667*7.0) +
+ -0.16666667*(0.16666667*7.6+0.04465820*7.4+0.16666667*7.8+0.62200847*7.0),
+
+ -0.62200847*(0.62200847*9.6+0.16666667*9.4+0.04465820*9.8+0.16666667*9.0) +
+ -0.16666667*(0.16666667*9.6+0.62200847*9.4+0.16666667*9.8+0.04465820*9.0) +
+ -0.04465820*(0.04465820*9.6+0.16666667*9.4+0.62200847*9.8+0.16666667*9.0) +
+ -0.16666667*(0.16666667*9.6+0.04465820*9.4+0.16666667*9.8+0.62200847*9.0),
+
+ // 8
+ -0.62200847*(0.62200847*5.8+0.16666667*5.0+0.04465820*5.6+0.16666667*5.4) +
+ -0.16666667*(0.16666667*5.8+0.62200847*5.0+0.16666667*5.6+0.04465820*5.4) +
+ -0.04465820*(0.04465820*5.8+0.16666667*5.0+0.62200847*5.6+0.16666667*5.4) +
+ -0.16666667*(0.16666667*5.8+0.04465820*5.0+0.16666667*5.6+0.62200847*5.4),
+
+ -0.62200847*(0.62200847*7.8+0.16666667*7.0+0.04465820*7.6+0.16666667*7.4) +
+ -0.16666667*(0.16666667*7.8+0.62200847*7.0+0.16666667*7.6+0.04465820*7.4) +
+ -0.04465820*(0.04465820*7.8+0.16666667*7.0+0.62200847*7.6+0.16666667*7.4) +
+ -0.16666667*(0.16666667*7.8+0.04465820*7.0+0.16666667*7.6+0.62200847*7.4),
+
+ -0.62200847*(0.62200847*9.8+0.16666667*9.0+0.04465820*9.6+0.16666667*9.4) +
+ -0.16666667*(0.16666667*9.8+0.62200847*9.0+0.16666667*9.6+0.04465820*9.4) +
+ -0.04465820*(0.04465820*9.8+0.16666667*9.0+0.62200847*9.6+0.16666667*9.4) +
+ -0.16666667*(0.16666667*9.8+0.04465820*9.0+0.16666667*9.6+0.62200847*9.4),
+
+ // 9
+ -0.62200847*(0.62200847*5.0+0.16666667*5.6+0.04465820*5.4+0.16666667*5.8) +
+ -0.16666667*(0.16666667*5.0+0.62200847*5.6+0.16666667*5.4+0.04465820*5.8) +
+ -0.04465820*(0.04465820*5.0+0.16666667*5.6+0.62200847*5.4+0.16666667*5.8) +
+ -0.16666667*(0.16666667*5.0+0.04465820*5.6+0.16666667*5.4+0.62200847*5.8),
+
+ -0.62200847*(0.62200847*7.0+0.16666667*7.6+0.04465820*7.4+0.16666667*7.8) +
+ -0.16666667*(0.16666667*7.0+0.62200847*7.6+0.16666667*7.4+0.04465820*7.8) +
+ -0.04465820*(0.04465820*7.0+0.16666667*7.6+0.62200847*7.4+0.16666667*7.8) +
+ -0.16666667*(0.16666667*7.0+0.04465820*7.6+0.16666667*7.4+0.62200847*7.8),
+
+ -0.62200847*(0.62200847*9.0+0.16666667*9.6+0.04465820*9.4+0.16666667*9.8) +
+ -0.16666667*(0.16666667*9.0+0.62200847*9.6+0.16666667*9.4+0.04465820*9.8) +
+ -0.04465820*(0.04465820*9.0+0.16666667*9.6+0.62200847*9.4+0.16666667*9.8) +
+ -0.16666667*(0.16666667*9.0+0.04465820*9.6+0.16666667*9.4+0.62200847*9.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,
- +9.4,+5.4,-7.4, // 14
- +9.6,+5.6,-7.6, // 15
- +9.8,+5.8,-7.8, // 16
- +9.0,+5.0,-7.0, // 17
- 0.8+1.82575588523, -0.8+-0.55566483464, 0.8+0.07938069066, // 18 (constraint)
- 0.9+1.69682900001, -0.9+-0.56560966667, 0.9+0.14140241667, // 19 (constraint)
- 1.0+1.51709826228, -1.0+-0.54615537442, 1.0+0.18205179147, // 20 (constraint)
- 1.1+1.29378670385, -1.1+-0.49761027071, 1.1+0.19904410828, // 21 (constraint)
+
+ // 14
+ +0.62200847*(0.62200847*5.4+0.16666667*5.6+0.04465820*5.0+0.16666667*5.8) +
+ +0.16666667*(0.16666667*5.4+0.62200847*5.6+0.16666667*5.0+0.04465820*5.8) +
+ +0.04465820*(0.04465820*5.4+0.16666667*5.6+0.62200847*5.0+0.16666667*5.8) +
+ +0.16666667*(0.16666667*5.4+0.04465820*5.6+0.16666667*5.0+0.62200847*5.8),
+
+ +0.62200847*(0.62200847*7.4+0.16666667*7.6+0.04465820*7.0+0.16666667*7.8) +
+ +0.16666667*(0.16666667*7.4+0.62200847*7.6+0.16666667*7.0+0.04465820*7.8) +
+ +0.04465820*(0.04465820*7.4+0.16666667*7.6+0.62200847*7.0+0.16666667*7.8) +
+ +0.16666667*(0.16666667*7.4+0.04465820*7.6+0.16666667*7.0+0.62200847*7.8),
+
+ +0.62200847*(0.62200847*9.4+0.16666667*9.6+0.04465820*9.0+0.16666667*9.8) +
+ +0.16666667*(0.16666667*9.4+0.62200847*9.6+0.16666667*9.0+0.04465820*9.8) +
+ +0.04465820*(0.04465820*9.4+0.16666667*9.6+0.62200847*9.0+0.16666667*9.8) +
+ +0.16666667*(0.16666667*9.4+0.04465820*9.6+0.16666667*9.0+0.62200847*9.8),
+
+ // 15
+ +0.62200847*(0.62200847*5.6+0.16666667*5.4+0.04465820*5.8+0.16666667*5.0) +
+ +0.16666667*(0.16666667*5.6+0.62200847*5.4+0.16666667*5.8+0.04465820*5.0) +
+ +0.04465820*(0.04465820*5.6+0.16666667*5.4+0.62200847*5.8+0.16666667*5.0) +
+ +0.16666667*(0.16666667*5.6+0.04465820*5.4+0.16666667*5.8+0.62200847*5.0),
+
+ +0.62200847*(0.62200847*7.6+0.16666667*7.4+0.04465820*7.8+0.16666667*7.0) +
+ +0.16666667*(0.16666667*7.6+0.62200847*7.4+0.16666667*7.8+0.04465820*7.0) +
+ +0.04465820*(0.04465820*7.6+0.16666667*7.4+0.62200847*7.8+0.16666667*7.0) +
+ +0.16666667*(0.16666667*7.6+0.04465820*7.4+0.16666667*7.8+0.62200847*7.0),
+
+ +0.62200847*(0.62200847*9.6+0.16666667*9.4+0.04465820*9.8+0.16666667*9.0) +
+ +0.16666667*(0.16666667*9.6+0.62200847*9.4+0.16666667*9.8+0.04465820*9.0) +
+ +0.04465820*(0.04465820*9.6+0.16666667*9.4+0.62200847*9.8+0.16666667*9.0) +
+ +0.16666667*(0.16666667*9.6+0.04465820*9.4+0.16666667*9.8+0.62200847*9.0),
+
+ // 16
+ +0.62200847*(0.62200847*5.8+0.16666667*5.0+0.04465820*5.6+0.16666667*5.4) +
+ +0.16666667*(0.16666667*5.8+0.62200847*5.0+0.16666667*5.6+0.04465820*5.4) +
+ +0.04465820*(0.04465820*5.8+0.16666667*5.0+0.62200847*5.6+0.16666667*5.4) +
+ +0.16666667*(0.16666667*5.8+0.04465820*5.0+0.16666667*5.6+0.62200847*5.4),
+
+ +0.62200847*(0.62200847*7.8+0.16666667*7.0+0.04465820*7.6+0.16666667*7.4) +
+ +0.16666667*(0.16666667*7.8+0.62200847*7.0+0.16666667*7.6+0.04465820*7.4) +
+ +0.04465820*(0.04465820*7.8+0.16666667*7.0+0.62200847*7.6+0.16666667*7.4) +
+ +0.16666667*(0.16666667*7.8+0.04465820*7.0+0.16666667*7.6+0.62200847*7.4),
+
+ +0.62200847*(0.62200847*9.8+0.16666667*9.0+0.04465820*9.6+0.16666667*9.4) +
+ +0.16666667*(0.16666667*9.8+0.62200847*9.0+0.16666667*9.6+0.04465820*9.4) +
+ +0.04465820*(0.04465820*9.8+0.16666667*9.0+0.62200847*9.6+0.16666667*9.4) +
+ +0.16666667*(0.16666667*9.8+0.04465820*9.0+0.16666667*9.6+0.62200847*9.4),
+
+ // 17
+ +0.62200847*(0.62200847*5.0+0.16666667*5.6+0.04465820*5.4+0.16666667*5.8) +
+ +0.16666667*(0.16666667*5.0+0.62200847*5.6+0.16666667*5.4+0.04465820*5.8) +
+ +0.04465820*(0.04465820*5.0+0.16666667*5.6+0.62200847*5.4+0.16666667*5.8) +
+ +0.16666667*(0.16666667*5.0+0.04465820*5.6+0.16666667*5.4+0.62200847*5.8),
+
+ +0.62200847*(0.62200847*7.0+0.16666667*7.6+0.04465820*7.4+0.16666667*7.8) +
+ +0.16666667*(0.16666667*7.0+0.62200847*7.6+0.16666667*7.4+0.04465820*7.8) +
+ +0.04465820*(0.04465820*7.0+0.16666667*7.6+0.62200847*7.4+0.16666667*7.8) +
+ +0.16666667*(0.16666667*7.0+0.04465820*7.6+0.16666667*7.4+0.62200847*7.8),
+
+ +0.62200847*(0.62200847*9.0+0.16666667*9.6+0.04465820*9.4+0.16666667*9.8) +
+ +0.16666667*(0.16666667*9.0+0.62200847*9.6+0.16666667*9.4+0.04465820*9.8) +
+ +0.04465820*(0.04465820*9.0+0.16666667*9.6+0.62200847*9.4+0.16666667*9.8) +
+ +0.16666667*(0.16666667*9.0+0.04465820*9.6+0.16666667*9.4+0.62200847*9.8),
+
+ // 18 (constraint)
+ 0.62200847*(0.62200847*(5.3-4.5+0.07938069066)+0.16666667*(5.7-4.7+0.18205179147)+0.04465820*(5.9-4.8+0.19904410828)+0.16666667*(5.5-4.6+0.14140241667)) +
+ 0.16666667*(0.16666667*(5.3-4.5+0.07938069066)+0.62200847*(5.7-4.7+0.18205179147)+0.16666667*(5.9-4.8+0.19904410828)+0.04465820*(5.5-4.6+0.14140241667)) +
+ 0.04465820*(0.04465820*(5.3-4.5+0.07938069066)+0.16666667*(5.7-4.7+0.18205179147)+0.62200847*(5.9-4.8+0.19904410828)+0.16666667*(5.5-4.6+0.14140241667)) +
+ 0.16666667*(0.16666667*(5.3-4.5+0.07938069066)+0.04465820*(5.7-4.7+0.18205179147)+0.16666667*(5.9-4.8+0.19904410828)+0.62200847*(5.5-4.6+0.14140241667)),
+
+ 0.62200847*(0.62200847*(7.3-6.5+1.82575588523)+0.16666667*(7.7-6.7+1.51709826228)+0.04465820*(7.9-6.8+1.29378670385)+0.16666667*(7.5-6.6+1.69682900001)) +
+ 0.16666667*(0.16666667*(7.3-6.5+1.82575588523)+0.62200847*(7.7-6.7+1.51709826228)+0.16666667*(7.9-6.8+1.29378670385)+0.04465820*(7.5-6.6+1.69682900001)) +
+ 0.04465820*(0.04465820*(7.3-6.5+1.82575588523)+0.16666667*(7.7-6.7+1.51709826228)+0.62200847*(7.9-6.8+1.29378670385)+0.16666667*(7.5-6.6+1.69682900001)) +
+ 0.16666667*(0.16666667*(7.3-6.5+1.82575588523)+0.04465820*(7.7-6.7+1.51709826228)+0.16666667*(7.9-6.8+1.29378670385)+0.62200847*(7.5-6.6+1.69682900001)),
+
+ 0.62200847*(0.62200847*(9.3-8.5+0.55566483464)+0.16666667*(9.7-8.7+0.54615537442)+0.04465820*(9.9-8.8+0.49761027071)+0.16666667*(9.5-8.6+0.56560966667)) +
+ 0.16666667*(0.16666667*(9.3-8.5+0.55566483464)+0.62200847*(9.7-8.7+0.54615537442)+0.16666667*(9.9-8.8+0.49761027071)+0.04465820*(9.5-8.6+0.56560966667)) +
+ 0.04465820*(0.04465820*(9.3-8.5+0.55566483464)+0.16666667*(9.7-8.7+0.54615537442)+0.62200847*(9.9-8.8+0.49761027071)+0.16666667*(9.5-8.6+0.56560966667)) +
+ 0.16666667*(0.16666667*(9.3-8.5+0.55566483464)+0.04465820*(9.7-8.7+0.54615537442)+0.16666667*(9.9-8.8+0.49761027071)+0.62200847*(9.5-8.6+0.56560966667)),
+
+ // 19 (constraint)
+ 0.62200847*(0.62200847*(5.5-4.6+0.14140241667)+0.16666667*(5.3-4.5+0.07938069066)+0.04465820*(5.7-4.7+0.18205179147)+0.16666667*(5.9-4.8+0.19904410828)) +
+ 0.16666667*(0.16666667*(5.5-4.6+0.14140241667)+0.62200847*(5.3-4.5+0.07938069066)+0.16666667*(5.7-4.7+0.18205179147)+0.04465820*(5.9-4.8+0.19904410828)) +
+ 0.04465820*(0.04465820*(5.5-4.6+0.14140241667)+0.16666667*(5.3-4.5+0.07938069066)+0.62200847*(5.7-4.7+0.18205179147)+0.16666667*(5.9-4.8+0.19904410828)) +
+ 0.16666667*(0.16666667*(5.5-4.6+0.14140241667)+0.04465820*(5.3-4.5+0.07938069066)+0.16666667*(5.7-4.7+0.18205179147)+0.62200847*(5.9-4.8+0.19904410828)),
+
+ 0.62200847*(0.62200847*(7.5-6.6+1.69682900001)+0.16666667*(7.3-6.5+1.82575588523)+0.04465820*(7.7-6.7+1.51709826228)+0.16666667*(7.9-6.8+1.29378670385)) +
+ 0.16666667*(0.16666667*(7.5-6.6+1.69682900001)+0.62200847*(7.3-6.5+1.82575588523)+0.16666667*(7.7-6.7+1.51709826228)+0.04465820*(7.9-6.8+1.29378670385)) +
+ 0.04465820*(0.04465820*(7.5-6.6+1.69682900001)+0.16666667*(7.3-6.5+1.82575588523)+0.62200847*(7.7-6.7+1.51709826228)+0.16666667*(7.9-6.8+1.29378670385)) +
+ 0.16666667*(0.16666667*(7.5-6.6+1.69682900001)+0.04465820*(7.3-6.5+1.82575588523)+0.16666667*(7.7-6.7+1.51709826228)+0.62200847*(7.9-6.8+1.29378670385)),
+
+ 0.62200847*(0.62200847*(9.5-8.6+0.56560966667)+0.16666667*(9.3-8.5+0.55566483464)+0.04465820*(9.7-8.7+0.54615537442)+0.16666667*(9.9-8.8+0.49761027071)) +
+ 0.16666667*(0.16666667*(9.5-8.6+0.56560966667)+0.62200847*(9.3-8.5+0.55566483464)+0.16666667*(9.7-8.7+0.54615537442)+0.04465820*(9.9-8.8+0.49761027071)) +
+ 0.04465820*(0.04465820*(9.5-8.6+0.56560966667)+0.16666667*(9.3-8.5+0.55566483464)+0.62200847*(9.7-8.7+0.54615537442)+0.16666667*(9.9-8.8+0.49761027071)) +
+ 0.16666667*(0.16666667*(9.5-8.6+0.56560966667)+0.04465820*(9.3-8.5+0.55566483464)+0.16666667*(9.7-8.7+0.54615537442)+0.62200847*(9.9-8.8+0.49761027071)),
+
+ // 20 (constraint)
+ 0.62200847*(0.62200847*(5.7-4.7+0.18205179147)+0.16666667*(5.9-4.8+0.19904410828)+0.04465820*(5.5-4.6+0.14140241667)+0.16666667*(5.3-4.5+0.07938069066)) +
+ 0.16666667*(0.16666667*(5.7-4.7+0.18205179147)+0.62200847*(5.9-4.8+0.19904410828)+0.16666667*(5.5-4.6+0.14140241667)+0.04465820*(5.3-4.5+0.07938069066)) +
+ 0.04465820*(0.04465820*(5.7-4.7+0.18205179147)+0.16666667*(5.9-4.8+0.19904410828)+0.62200847*(5.5-4.6+0.14140241667)+0.16666667*(5.3-4.5+0.07938069066)) +
+ 0.16666667*(0.16666667*(5.7-4.7+0.18205179147)+0.04465820*(5.9-4.8+0.19904410828)+0.16666667*(5.5-4.6+0.14140241667)+0.62200847*(5.3-4.5+0.07938069066)),
+
+ 0.62200847*(0.62200847*(7.7-6.7+1.51709826228)+0.16666667*(7.9-6.8+1.29378670385)+0.04465820*(7.5-6.6+1.69682900001)+0.16666667*(7.3-6.5+1.82575588523)) +
+ 0.16666667*(0.16666667*(7.7-6.7+1.51709826228)+0.62200847*(7.9-6.8+1.29378670385)+0.16666667*(7.5-6.6+1.69682900001)+0.04465820*(7.3-6.5+1.82575588523)) +
+ 0.04465820*(0.04465820*(7.7-6.7+1.51709826228)+0.16666667*(7.9-6.8+1.29378670385)+0.62200847*(7.5-6.6+1.69682900001)+0.16666667*(7.3-6.5+1.82575588523)) +
+ 0.16666667*(0.16666667*(7.7-6.7+1.51709826228)+0.04465820*(7.9-6.8+1.29378670385)+0.16666667*(7.5-6.6+1.69682900001)+0.62200847*(7.3-6.5+1.82575588523)),
+
+ 0.62200847*(0.62200847*(9.7-8.7+0.54615537442)+0.16666667*(9.9-8.8+0.49761027071)+0.04465820*(9.5-8.6+0.56560966667)+0.16666667*(9.3-8.5+0.55566483464)) +
+ 0.16666667*(0.16666667*(9.7-8.7+0.54615537442)+0.62200847*(9.9-8.8+0.49761027071)+0.16666667*(9.5-8.6+0.56560966667)+0.04465820*(9.3-8.5+0.55566483464)) +
+ 0.04465820*(0.04465820*(9.7-8.7+0.54615537442)+0.16666667*(9.9-8.8+0.49761027071)+0.62200847*(9.5-8.6+0.56560966667)+0.16666667*(9.3-8.5+0.55566483464)) +
+ 0.16666667*(0.16666667*(9.7-8.7+0.54615537442)+0.04465820*(9.9-8.8+0.49761027071)+0.16666667*(9.5-8.6+0.56560966667)+0.62200847*(9.3-8.5+0.55566483464)),
+
+ // 21 (constraint)
+ 0.62200847*(0.62200847*(5.9-4.8+0.19904410828)+0.16666667*(5.5-4.6+0.14140241667)+0.04465820*(5.3-4.5+0.07938069066)+0.16666667*(5.7-4.7+0.18205179147)) +
+ 0.16666667*(0.16666667*(5.9-4.8+0.19904410828)+0.62200847*(5.5-4.6+0.14140241667)+0.16666667*(5.3-4.5+0.07938069066)+0.04465820*(5.7-4.7+0.18205179147)) +
+ 0.04465820*(0.04465820*(5.9-4.8+0.19904410828)+0.16666667*(5.5-4.6+0.14140241667)+0.62200847*(5.3-4.5+0.07938069066)+0.16666667*(5.7-4.7+0.18205179147)) +
+ 0.16666667*(0.16666667*(5.9-4.8+0.19904410828)+0.04465820*(5.5-4.6+0.14140241667)+0.16666667*(5.3-4.5+0.07938069066)+0.62200847*(5.7-4.7+0.18205179147)),
+
+ 0.62200847*(0.62200847*(7.9-6.8+1.29378670385)+0.16666667*(7.5-6.6+1.69682900001)+0.04465820*(7.3-6.5+1.82575588523)+0.16666667*(7.7-6.7+1.51709826228)) +
+ 0.16666667*(0.16666667*(7.9-6.8+1.29378670385)+0.62200847*(7.5-6.6+1.69682900001)+0.16666667*(7.3-6.5+1.82575588523)+0.04465820*(7.7-6.7+1.51709826228)) +
+ 0.04465820*(0.04465820*(7.9-6.8+1.29378670385)+0.16666667*(7.5-6.6+1.69682900001)+0.62200847*(7.3-6.5+1.82575588523)+0.16666667*(7.7-6.7+1.51709826228)) +
+ 0.16666667*(0.16666667*(7.9-6.8+1.29378670385)+0.04465820*(7.5-6.6+1.69682900001)+0.16666667*(7.3-6.5+1.82575588523)+0.62200847*(7.7-6.7+1.51709826228)),
+
+ 0.62200847*(0.62200847*(9.9-8.8+0.49761027071)+0.16666667*(9.5-8.6+0.56560966667)+0.04465820*(9.3-8.5+0.55566483464)+0.16666667*(9.7-8.7+0.54615537442)) +
+ 0.16666667*(0.16666667*(9.9-8.8+0.49761027071)+0.62200847*(9.5-8.6+0.56560966667)+0.16666667*(9.3-8.5+0.55566483464)+0.04465820*(9.7-8.7+0.54615537442)) +
+ 0.04465820*(0.04465820*(9.9-8.8+0.49761027071)+0.16666667*(9.5-8.6+0.56560966667)+0.62200847*(9.3-8.5+0.55566483464)+0.16666667*(9.7-8.7+0.54615537442)) +
+ 0.16666667*(0.16666667*(9.9-8.8+0.49761027071)+0.04465820*(9.5-8.6+0.56560966667)+0.16666667*(9.3-8.5+0.55566483464)+0.62200847*(9.7-8.7+0.54615537442)),
+
};
const double pylith::faults::CohesiveKinDataHex8::_residualIncr[] = {
@@ -238,22 +420,203 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -9.4,-5.4,+7.4, // 6
- -9.6,-5.6,+7.6, // 7
- -9.8,-5.8,+7.8, // 8
- -9.0,-5.0,+7.0, // 9
+
+ // 6
+ -0.62200847*(0.62200847*5.4+0.16666667*5.6+0.04465820*5.0+0.16666667*5.8) +
+ -0.16666667*(0.16666667*5.4+0.62200847*5.6+0.16666667*5.0+0.04465820*5.8) +
+ -0.04465820*(0.04465820*5.4+0.16666667*5.6+0.62200847*5.0+0.16666667*5.8) +
+ -0.16666667*(0.16666667*5.4+0.04465820*5.6+0.16666667*5.0+0.62200847*5.8),
+
+ -0.62200847*(0.62200847*7.4+0.16666667*7.6+0.04465820*7.0+0.16666667*7.8) +
+ -0.16666667*(0.16666667*7.4+0.62200847*7.6+0.16666667*7.0+0.04465820*7.8) +
+ -0.04465820*(0.04465820*7.4+0.16666667*7.6+0.62200847*7.0+0.16666667*7.8) +
+ -0.16666667*(0.16666667*7.4+0.04465820*7.6+0.16666667*7.0+0.62200847*7.8),
+
+ -0.62200847*(0.62200847*9.4+0.16666667*9.6+0.04465820*9.0+0.16666667*9.8) +
+ -0.16666667*(0.16666667*9.4+0.62200847*9.6+0.16666667*9.0+0.04465820*9.8) +
+ -0.04465820*(0.04465820*9.4+0.16666667*9.6+0.62200847*9.0+0.16666667*9.8) +
+ -0.16666667*(0.16666667*9.4+0.04465820*9.6+0.16666667*9.0+0.62200847*9.8),
+
+ // 7
+ -0.62200847*(0.62200847*5.6+0.16666667*5.4+0.04465820*5.8+0.16666667*5.0) +
+ -0.16666667*(0.16666667*5.6+0.62200847*5.4+0.16666667*5.8+0.04465820*5.0) +
+ -0.04465820*(0.04465820*5.6+0.16666667*5.4+0.62200847*5.8+0.16666667*5.0) +
+ -0.16666667*(0.16666667*5.6+0.04465820*5.4+0.16666667*5.8+0.62200847*5.0),
+
+ -0.62200847*(0.62200847*7.6+0.16666667*7.4+0.04465820*7.8+0.16666667*7.0) +
+ -0.16666667*(0.16666667*7.6+0.62200847*7.4+0.16666667*7.8+0.04465820*7.0) +
+ -0.04465820*(0.04465820*7.6+0.16666667*7.4+0.62200847*7.8+0.16666667*7.0) +
+ -0.16666667*(0.16666667*7.6+0.04465820*7.4+0.16666667*7.8+0.62200847*7.0),
+
+ -0.62200847*(0.62200847*9.6+0.16666667*9.4+0.04465820*9.8+0.16666667*9.0) +
+ -0.16666667*(0.16666667*9.6+0.62200847*9.4+0.16666667*9.8+0.04465820*9.0) +
+ -0.04465820*(0.04465820*9.6+0.16666667*9.4+0.62200847*9.8+0.16666667*9.0) +
+ -0.16666667*(0.16666667*9.6+0.04465820*9.4+0.16666667*9.8+0.62200847*9.0),
+
+ // 8
+ -0.62200847*(0.62200847*5.8+0.16666667*5.0+0.04465820*5.6+0.16666667*5.4) +
+ -0.16666667*(0.16666667*5.8+0.62200847*5.0+0.16666667*5.6+0.04465820*5.4) +
+ -0.04465820*(0.04465820*5.8+0.16666667*5.0+0.62200847*5.6+0.16666667*5.4) +
+ -0.16666667*(0.16666667*5.8+0.04465820*5.0+0.16666667*5.6+0.62200847*5.4),
+
+ -0.62200847*(0.62200847*7.8+0.16666667*7.0+0.04465820*7.6+0.16666667*7.4) +
+ -0.16666667*(0.16666667*7.8+0.62200847*7.0+0.16666667*7.6+0.04465820*7.4) +
+ -0.04465820*(0.04465820*7.8+0.16666667*7.0+0.62200847*7.6+0.16666667*7.4) +
+ -0.16666667*(0.16666667*7.8+0.04465820*7.0+0.16666667*7.6+0.62200847*7.4),
+
+ -0.62200847*(0.62200847*9.8+0.16666667*9.0+0.04465820*9.6+0.16666667*9.4) +
+ -0.16666667*(0.16666667*9.8+0.62200847*9.0+0.16666667*9.6+0.04465820*9.4) +
+ -0.04465820*(0.04465820*9.8+0.16666667*9.0+0.62200847*9.6+0.16666667*9.4) +
+ -0.16666667*(0.16666667*9.8+0.04465820*9.0+0.16666667*9.6+0.62200847*9.4),
+
+ // 9
+ -0.62200847*(0.62200847*5.0+0.16666667*5.6+0.04465820*5.4+0.16666667*5.8) +
+ -0.16666667*(0.16666667*5.0+0.62200847*5.6+0.16666667*5.4+0.04465820*5.8) +
+ -0.04465820*(0.04465820*5.0+0.16666667*5.6+0.62200847*5.4+0.16666667*5.8) +
+ -0.16666667*(0.16666667*5.0+0.04465820*5.6+0.16666667*5.4+0.62200847*5.8),
+
+ -0.62200847*(0.62200847*7.0+0.16666667*7.6+0.04465820*7.4+0.16666667*7.8) +
+ -0.16666667*(0.16666667*7.0+0.62200847*7.6+0.16666667*7.4+0.04465820*7.8) +
+ -0.04465820*(0.04465820*7.0+0.16666667*7.6+0.62200847*7.4+0.16666667*7.8) +
+ -0.16666667*(0.16666667*7.0+0.04465820*7.6+0.16666667*7.4+0.62200847*7.8),
+
+ -0.62200847*(0.62200847*9.0+0.16666667*9.6+0.04465820*9.4+0.16666667*9.8) +
+ -0.16666667*(0.16666667*9.0+0.62200847*9.6+0.16666667*9.4+0.04465820*9.8) +
+ -0.04465820*(0.04465820*9.0+0.16666667*9.6+0.62200847*9.4+0.16666667*9.8) +
+ -0.16666667*(0.16666667*9.0+0.04465820*9.6+0.16666667*9.4+0.62200847*9.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,
- +9.4,+5.4,-7.4, // 14
- +9.6,+5.6,-7.6, // 15
- +9.8,+5.8,-7.8, // 16
- +9.0,+5.0,-7.0, // 17
- 0.8+1.82575588523, -0.8+-0.55566483464, 0.8+0.07938069066, // 18 (constraint)
- 0.9+1.69682900001, -0.9+-0.56560966667, 0.9+0.14140241667, // 19 (constraint)
- 1.0+1.51709826228, -1.0+-0.54615537442, 1.0+0.18205179147, // 20 (constraint)
- 1.1+1.29378670385, -1.1+-0.49761027071, 1.1+0.19904410828, // 21 (constraint)
+
+ // 14
+ +0.62200847*(0.62200847*5.4+0.16666667*5.6+0.04465820*5.0+0.16666667*5.8) +
+ +0.16666667*(0.16666667*5.4+0.62200847*5.6+0.16666667*5.0+0.04465820*5.8) +
+ +0.04465820*(0.04465820*5.4+0.16666667*5.6+0.62200847*5.0+0.16666667*5.8) +
+ +0.16666667*(0.16666667*5.4+0.04465820*5.6+0.16666667*5.0+0.62200847*5.8),
+
+ +0.62200847*(0.62200847*7.4+0.16666667*7.6+0.04465820*7.0+0.16666667*7.8) +
+ +0.16666667*(0.16666667*7.4+0.62200847*7.6+0.16666667*7.0+0.04465820*7.8) +
+ +0.04465820*(0.04465820*7.4+0.16666667*7.6+0.62200847*7.0+0.16666667*7.8) +
+ +0.16666667*(0.16666667*7.4+0.04465820*7.6+0.16666667*7.0+0.62200847*7.8),
+
+ +0.62200847*(0.62200847*9.4+0.16666667*9.6+0.04465820*9.0+0.16666667*9.8) +
+ +0.16666667*(0.16666667*9.4+0.62200847*9.6+0.16666667*9.0+0.04465820*9.8) +
+ +0.04465820*(0.04465820*9.4+0.16666667*9.6+0.62200847*9.0+0.16666667*9.8) +
+ +0.16666667*(0.16666667*9.4+0.04465820*9.6+0.16666667*9.0+0.62200847*9.8),
+
+ // 15
+ +0.62200847*(0.62200847*5.6+0.16666667*5.4+0.04465820*5.8+0.16666667*5.0) +
+ +0.16666667*(0.16666667*5.6+0.62200847*5.4+0.16666667*5.8+0.04465820*5.0) +
+ +0.04465820*(0.04465820*5.6+0.16666667*5.4+0.62200847*5.8+0.16666667*5.0) +
+ +0.16666667*(0.16666667*5.6+0.04465820*5.4+0.16666667*5.8+0.62200847*5.0),
+
+ +0.62200847*(0.62200847*7.6+0.16666667*7.4+0.04465820*7.8+0.16666667*7.0) +
+ +0.16666667*(0.16666667*7.6+0.62200847*7.4+0.16666667*7.8+0.04465820*7.0) +
+ +0.04465820*(0.04465820*7.6+0.16666667*7.4+0.62200847*7.8+0.16666667*7.0) +
+ +0.16666667*(0.16666667*7.6+0.04465820*7.4+0.16666667*7.8+0.62200847*7.0),
+
+ +0.62200847*(0.62200847*9.6+0.16666667*9.4+0.04465820*9.8+0.16666667*9.0) +
+ +0.16666667*(0.16666667*9.6+0.62200847*9.4+0.16666667*9.8+0.04465820*9.0) +
+ +0.04465820*(0.04465820*9.6+0.16666667*9.4+0.62200847*9.8+0.16666667*9.0) +
+ +0.16666667*(0.16666667*9.6+0.04465820*9.4+0.16666667*9.8+0.62200847*9.0),
+
+ // 16
+ +0.62200847*(0.62200847*5.8+0.16666667*5.0+0.04465820*5.6+0.16666667*5.4) +
+ +0.16666667*(0.16666667*5.8+0.62200847*5.0+0.16666667*5.6+0.04465820*5.4) +
+ +0.04465820*(0.04465820*5.8+0.16666667*5.0+0.62200847*5.6+0.16666667*5.4) +
+ +0.16666667*(0.16666667*5.8+0.04465820*5.0+0.16666667*5.6+0.62200847*5.4),
+
+ +0.62200847*(0.62200847*7.8+0.16666667*7.0+0.04465820*7.6+0.16666667*7.4) +
+ +0.16666667*(0.16666667*7.8+0.62200847*7.0+0.16666667*7.6+0.04465820*7.4) +
+ +0.04465820*(0.04465820*7.8+0.16666667*7.0+0.62200847*7.6+0.16666667*7.4) +
+ +0.16666667*(0.16666667*7.8+0.04465820*7.0+0.16666667*7.6+0.62200847*7.4),
+
+ +0.62200847*(0.62200847*9.8+0.16666667*9.0+0.04465820*9.6+0.16666667*9.4) +
+ +0.16666667*(0.16666667*9.8+0.62200847*9.0+0.16666667*9.6+0.04465820*9.4) +
+ +0.04465820*(0.04465820*9.8+0.16666667*9.0+0.62200847*9.6+0.16666667*9.4) +
+ +0.16666667*(0.16666667*9.8+0.04465820*9.0+0.16666667*9.6+0.62200847*9.4),
+
+ // 17
+ +0.62200847*(0.62200847*5.0+0.16666667*5.6+0.04465820*5.4+0.16666667*5.8) +
+ +0.16666667*(0.16666667*5.0+0.62200847*5.6+0.16666667*5.4+0.04465820*5.8) +
+ +0.04465820*(0.04465820*5.0+0.16666667*5.6+0.62200847*5.4+0.16666667*5.8) +
+ +0.16666667*(0.16666667*5.0+0.04465820*5.6+0.16666667*5.4+0.62200847*5.8),
+
+ +0.62200847*(0.62200847*7.0+0.16666667*7.6+0.04465820*7.4+0.16666667*7.8) +
+ +0.16666667*(0.16666667*7.0+0.62200847*7.6+0.16666667*7.4+0.04465820*7.8) +
+ +0.04465820*(0.04465820*7.0+0.16666667*7.6+0.62200847*7.4+0.16666667*7.8) +
+ +0.16666667*(0.16666667*7.0+0.04465820*7.6+0.16666667*7.4+0.62200847*7.8),
+
+ +0.62200847*(0.62200847*9.0+0.16666667*9.6+0.04465820*9.4+0.16666667*9.8) +
+ +0.16666667*(0.16666667*9.0+0.62200847*9.6+0.16666667*9.4+0.04465820*9.8) +
+ +0.04465820*(0.04465820*9.0+0.16666667*9.6+0.62200847*9.4+0.16666667*9.8) +
+ +0.16666667*(0.16666667*9.0+0.04465820*9.6+0.16666667*9.4+0.62200847*9.8),
+
+ // 18 (constraint)
+ 0.62200847*(0.62200847*(5.3-4.5+0.07938069066)+0.16666667*(5.7-4.7+0.18205179147)+0.04465820*(5.9-4.8+0.19904410828)+0.16666667*(5.5-4.6+0.14140241667)) +
+ 0.16666667*(0.16666667*(5.3-4.5+0.07938069066)+0.62200847*(5.7-4.7+0.18205179147)+0.16666667*(5.9-4.8+0.19904410828)+0.04465820*(5.5-4.6+0.14140241667)) +
+ 0.04465820*(0.04465820*(5.3-4.5+0.07938069066)+0.16666667*(5.7-4.7+0.18205179147)+0.62200847*(5.9-4.8+0.19904410828)+0.16666667*(5.5-4.6+0.14140241667)) +
+ 0.16666667*(0.16666667*(5.3-4.5+0.07938069066)+0.04465820*(5.7-4.7+0.18205179147)+0.16666667*(5.9-4.8+0.19904410828)+0.62200847*(5.5-4.6+0.14140241667)),
+
+ 0.62200847*(0.62200847*(7.3-6.5+1.82575588523)+0.16666667*(7.7-6.7+1.51709826228)+0.04465820*(7.9-6.8+1.29378670385)+0.16666667*(7.5-6.6+1.69682900001)) +
+ 0.16666667*(0.16666667*(7.3-6.5+1.82575588523)+0.62200847*(7.7-6.7+1.51709826228)+0.16666667*(7.9-6.8+1.29378670385)+0.04465820*(7.5-6.6+1.69682900001)) +
+ 0.04465820*(0.04465820*(7.3-6.5+1.82575588523)+0.16666667*(7.7-6.7+1.51709826228)+0.62200847*(7.9-6.8+1.29378670385)+0.16666667*(7.5-6.6+1.69682900001)) +
+ 0.16666667*(0.16666667*(7.3-6.5+1.82575588523)+0.04465820*(7.7-6.7+1.51709826228)+0.16666667*(7.9-6.8+1.29378670385)+0.62200847*(7.5-6.6+1.69682900001)),
+
+ 0.62200847*(0.62200847*(9.3-8.5+0.55566483464)+0.16666667*(9.7-8.7+0.54615537442)+0.04465820*(9.9-8.8+0.49761027071)+0.16666667*(9.5-8.6+0.56560966667)) +
+ 0.16666667*(0.16666667*(9.3-8.5+0.55566483464)+0.62200847*(9.7-8.7+0.54615537442)+0.16666667*(9.9-8.8+0.49761027071)+0.04465820*(9.5-8.6+0.56560966667)) +
+ 0.04465820*(0.04465820*(9.3-8.5+0.55566483464)+0.16666667*(9.7-8.7+0.54615537442)+0.62200847*(9.9-8.8+0.49761027071)+0.16666667*(9.5-8.6+0.56560966667)) +
+ 0.16666667*(0.16666667*(9.3-8.5+0.55566483464)+0.04465820*(9.7-8.7+0.54615537442)+0.16666667*(9.9-8.8+0.49761027071)+0.62200847*(9.5-8.6+0.56560966667)),
+
+ // 19 (constraint)
+ 0.62200847*(0.62200847*(5.5-4.6+0.14140241667)+0.16666667*(5.3-4.5+0.07938069066)+0.04465820*(5.7-4.7+0.18205179147)+0.16666667*(5.9-4.8+0.19904410828)) +
+ 0.16666667*(0.16666667*(5.5-4.6+0.14140241667)+0.62200847*(5.3-4.5+0.07938069066)+0.16666667*(5.7-4.7+0.18205179147)+0.04465820*(5.9-4.8+0.19904410828)) +
+ 0.04465820*(0.04465820*(5.5-4.6+0.14140241667)+0.16666667*(5.3-4.5+0.07938069066)+0.62200847*(5.7-4.7+0.18205179147)+0.16666667*(5.9-4.8+0.19904410828)) +
+ 0.16666667*(0.16666667*(5.5-4.6+0.14140241667)+0.04465820*(5.3-4.5+0.07938069066)+0.16666667*(5.7-4.7+0.18205179147)+0.62200847*(5.9-4.8+0.19904410828)),
+
+ 0.62200847*(0.62200847*(7.5-6.6+1.69682900001)+0.16666667*(7.3-6.5+1.82575588523)+0.04465820*(7.7-6.7+1.51709826228)+0.16666667*(7.9-6.8+1.29378670385)) +
+ 0.16666667*(0.16666667*(7.5-6.6+1.69682900001)+0.62200847*(7.3-6.5+1.82575588523)+0.16666667*(7.7-6.7+1.51709826228)+0.04465820*(7.9-6.8+1.29378670385)) +
+ 0.04465820*(0.04465820*(7.5-6.6+1.69682900001)+0.16666667*(7.3-6.5+1.82575588523)+0.62200847*(7.7-6.7+1.51709826228)+0.16666667*(7.9-6.8+1.29378670385)) +
+ 0.16666667*(0.16666667*(7.5-6.6+1.69682900001)+0.04465820*(7.3-6.5+1.82575588523)+0.16666667*(7.7-6.7+1.51709826228)+0.62200847*(7.9-6.8+1.29378670385)),
+
+ 0.62200847*(0.62200847*(9.5-8.6+0.56560966667)+0.16666667*(9.3-8.5+0.55566483464)+0.04465820*(9.7-8.7+0.54615537442)+0.16666667*(9.9-8.8+0.49761027071)) +
+ 0.16666667*(0.16666667*(9.5-8.6+0.56560966667)+0.62200847*(9.3-8.5+0.55566483464)+0.16666667*(9.7-8.7+0.54615537442)+0.04465820*(9.9-8.8+0.49761027071)) +
+ 0.04465820*(0.04465820*(9.5-8.6+0.56560966667)+0.16666667*(9.3-8.5+0.55566483464)+0.62200847*(9.7-8.7+0.54615537442)+0.16666667*(9.9-8.8+0.49761027071)) +
+ 0.16666667*(0.16666667*(9.5-8.6+0.56560966667)+0.04465820*(9.3-8.5+0.55566483464)+0.16666667*(9.7-8.7+0.54615537442)+0.62200847*(9.9-8.8+0.49761027071)),
+
+ // 20 (constraint)
+ 0.62200847*(0.62200847*(5.7-4.7+0.18205179147)+0.16666667*(5.9-4.8+0.19904410828)+0.04465820*(5.5-4.6+0.14140241667)+0.16666667*(5.3-4.5+0.07938069066)) +
+ 0.16666667*(0.16666667*(5.7-4.7+0.18205179147)+0.62200847*(5.9-4.8+0.19904410828)+0.16666667*(5.5-4.6+0.14140241667)+0.04465820*(5.3-4.5+0.07938069066)) +
+ 0.04465820*(0.04465820*(5.7-4.7+0.18205179147)+0.16666667*(5.9-4.8+0.19904410828)+0.62200847*(5.5-4.6+0.14140241667)+0.16666667*(5.3-4.5+0.07938069066)) +
+ 0.16666667*(0.16666667*(5.7-4.7+0.18205179147)+0.04465820*(5.9-4.8+0.19904410828)+0.16666667*(5.5-4.6+0.14140241667)+0.62200847*(5.3-4.5+0.07938069066)),
+
+ 0.62200847*(0.62200847*(7.7-6.7+1.51709826228)+0.16666667*(7.9-6.8+1.29378670385)+0.04465820*(7.5-6.6+1.69682900001)+0.16666667*(7.3-6.5+1.82575588523)) +
+ 0.16666667*(0.16666667*(7.7-6.7+1.51709826228)+0.62200847*(7.9-6.8+1.29378670385)+0.16666667*(7.5-6.6+1.69682900001)+0.04465820*(7.3-6.5+1.82575588523)) +
+ 0.04465820*(0.04465820*(7.7-6.7+1.51709826228)+0.16666667*(7.9-6.8+1.29378670385)+0.62200847*(7.5-6.6+1.69682900001)+0.16666667*(7.3-6.5+1.82575588523)) +
+ 0.16666667*(0.16666667*(7.7-6.7+1.51709826228)+0.04465820*(7.9-6.8+1.29378670385)+0.16666667*(7.5-6.6+1.69682900001)+0.62200847*(7.3-6.5+1.82575588523)),
+
+ 0.62200847*(0.62200847*(9.7-8.7+0.54615537442)+0.16666667*(9.9-8.8+0.49761027071)+0.04465820*(9.5-8.6+0.56560966667)+0.16666667*(9.3-8.5+0.55566483464)) +
+ 0.16666667*(0.16666667*(9.7-8.7+0.54615537442)+0.62200847*(9.9-8.8+0.49761027071)+0.16666667*(9.5-8.6+0.56560966667)+0.04465820*(9.3-8.5+0.55566483464)) +
+ 0.04465820*(0.04465820*(9.7-8.7+0.54615537442)+0.16666667*(9.9-8.8+0.49761027071)+0.62200847*(9.5-8.6+0.56560966667)+0.16666667*(9.3-8.5+0.55566483464)) +
+ 0.16666667*(0.16666667*(9.7-8.7+0.54615537442)+0.04465820*(9.9-8.8+0.49761027071)+0.16666667*(9.5-8.6+0.56560966667)+0.62200847*(9.3-8.5+0.55566483464)),
+
+ // 21 (constraint)
+ 0.62200847*(0.62200847*(5.9-4.8+0.19904410828)+0.16666667*(5.5-4.6+0.14140241667)+0.04465820*(5.3-4.5+0.07938069066)+0.16666667*(5.7-4.7+0.18205179147)) +
+ 0.16666667*(0.16666667*(5.9-4.8+0.19904410828)+0.62200847*(5.5-4.6+0.14140241667)+0.16666667*(5.3-4.5+0.07938069066)+0.04465820*(5.7-4.7+0.18205179147)) +
+ 0.04465820*(0.04465820*(5.9-4.8+0.19904410828)+0.16666667*(5.5-4.6+0.14140241667)+0.62200847*(5.3-4.5+0.07938069066)+0.16666667*(5.7-4.7+0.18205179147)) +
+ 0.16666667*(0.16666667*(5.9-4.8+0.19904410828)+0.04465820*(5.5-4.6+0.14140241667)+0.16666667*(5.3-4.5+0.07938069066)+0.62200847*(5.7-4.7+0.18205179147)),
+
+ 0.62200847*(0.62200847*(7.9-6.8+1.29378670385)+0.16666667*(7.5-6.6+1.69682900001)+0.04465820*(7.3-6.5+1.82575588523)+0.16666667*(7.7-6.7+1.51709826228)) +
+ 0.16666667*(0.16666667*(7.9-6.8+1.29378670385)+0.62200847*(7.5-6.6+1.69682900001)+0.16666667*(7.3-6.5+1.82575588523)+0.04465820*(7.7-6.7+1.51709826228)) +
+ 0.04465820*(0.04465820*(7.9-6.8+1.29378670385)+0.16666667*(7.5-6.6+1.69682900001)+0.62200847*(7.3-6.5+1.82575588523)+0.16666667*(7.7-6.7+1.51709826228)) +
+ 0.16666667*(0.16666667*(7.9-6.8+1.29378670385)+0.04465820*(7.5-6.6+1.69682900001)+0.16666667*(7.3-6.5+1.82575588523)+0.62200847*(7.7-6.7+1.51709826228)),
+
+ 0.62200847*(0.62200847*(9.9-8.8+0.49761027071)+0.16666667*(9.5-8.6+0.56560966667)+0.04465820*(9.3-8.5+0.55566483464)+0.16666667*(9.7-8.7+0.54615537442)) +
+ 0.16666667*(0.16666667*(9.9-8.8+0.49761027071)+0.62200847*(9.5-8.6+0.56560966667)+0.16666667*(9.3-8.5+0.55566483464)+0.04465820*(9.7-8.7+0.54615537442)) +
+ 0.04465820*(0.04465820*(9.9-8.8+0.49761027071)+0.16666667*(9.5-8.6+0.56560966667)+0.62200847*(9.3-8.5+0.55566483464)+0.16666667*(9.7-8.7+0.54615537442)) +
+ 0.16666667*(0.16666667*(9.9-8.8+0.49761027071)+0.04465820*(9.5-8.6+0.56560966667)+0.16666667*(9.3-8.5+0.55566483464)+0.62200847*(9.7-8.7+0.54615537442)),
};
const double pylith::faults::CohesiveKinDataHex8::_jacobian[] = {
@@ -513,7 +876,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0, // 18
+ +1.0, 0.0, 0.0, // 18
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -533,7 +896,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +1.0, 0.0, 0.0, // 18
+ 0.0,+1.0, 0.0, // 18
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -553,7 +916,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0, 0.0, // 18
+ 0.0, 0.0,+1.0, // 18
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -574,7 +937,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0, // 19
+ +1.0, 0.0, 0.0, // 19
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 7y
@@ -594,7 +957,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +1.0, 0.0, 0.0, // 19
+ 0.0,+1.0, 0.0, // 19
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 7z
@@ -614,7 +977,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0, 0.0, // 19
+ 0.0, 0.0,+1.0, // 19
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 8x
@@ -635,7 +998,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0, // 20
+ +1.0, 0.0, 0.0, // 20
0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 8y
0.0, 0.0, 0.0,
@@ -655,7 +1018,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +1.0, 0.0, 0.0, // 20
+ 0.0,+1.0, 0.0, // 20
0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 8z
0.0, 0.0, 0.0,
@@ -675,7 +1038,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0, 0.0, // 20
+ 0.0, 0.0,+1.0, // 20
0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 9x
0.0, 0.0, 0.0,
@@ -696,7 +1059,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0, // 21
+ +1.0, 0.0, 0.0, // 21
0.0, 0.0, 0.0, // 9y
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -716,7 +1079,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +1.0, 0.0, 0.0, // 21
+ 0.0,+1.0, 0.0, // 21
0.0, 0.0, 0.0, // 9z
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -736,7 +1099,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0, 0.0, // 21
+ 0.0, 0.0,+1.0, // 21
0.0, 0.0, 0.0, // 10x
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -993,7 +1356,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0, // 18
+ -1.0, 0.0, 0.0, // 18
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1013,7 +1376,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0, 0.0, 0.0, // 18
+ 0.0,-1.0, 0.0, // 18
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1033,7 +1396,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0, 0.0, // 18
+ 0.0, 0.0,-1.0, // 18
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1054,7 +1417,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0, // 19
+ -1.0, 0.0, 0.0, // 19
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 15y
@@ -1074,7 +1437,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0, 0.0, 0.0, // 19
+ 0.0,-1.0, 0.0, // 19
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 15z
@@ -1094,7 +1457,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0, 0.0, // 19
+ 0.0, 0.0,-1.0, // 19
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 16x
@@ -1115,7 +1478,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0, // 20
+ -1.0, 0.0, 0.0, // 20
0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 16y
0.0, 0.0, 0.0,
@@ -1135,7 +1498,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0, 0.0, 0.0, // 20
+ 0.0,-1.0, 0.0, // 20
0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 16z
0.0, 0.0, 0.0,
@@ -1155,7 +1518,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0, 0.0, // 20
+ 0.0, 0.0,-1.0, // 20
0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 17x
0.0, 0.0, 0.0,
@@ -1176,7 +1539,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0, // 21
+ -1.0, 0.0, 0.0, // 21
0.0, 0.0, 0.0, // 17y
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1196,7 +1559,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0, 0.0, 0.0, // 21
+ 0.0,-1.0, 0.0, // 21
0.0, 0.0, 0.0, // 17z
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1216,12 +1579,12 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0, 0.0, // 21
+ 0.0, 0.0,-1.0, // 21
0.0, 0.0, 0.0, // 18x
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0, 0.0, // 6
+ +1.0, 0.0, 0.0, // 6
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1229,7 +1592,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0, 0.0, // 14
+ -1.0, 0.0, 0.0, // 14
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1241,7 +1604,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0, // 6
+ 0.0,+1.0, 0.0, // 6
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1249,7 +1612,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0, // 14
+ 0.0,-1.0, 0.0, // 14
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1261,7 +1624,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +1.0, 0.0, 0.0, // 6
+ 0.0, 0.0,+1.0, // 6
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1269,7 +1632,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0, 0.0, 0.0, // 14
+ 0.0, 0.0,-1.0, // 14
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1282,7 +1645,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0, 0.0, // 7
+ +1.0, 0.0, 0.0, // 7
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1290,7 +1653,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0, 0.0, // 15
+ -1.0, 0.0, 0.0, // 15
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1302,7 +1665,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0, // 7
+ 0.0,+1.0, 0.0, // 7
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1310,7 +1673,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0, // 15
+ 0.0,-1.0, 0.0, // 15
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1322,7 +1685,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +1.0, 0.0, 0.0, // 7
+ 0.0, 0.0,+1.0, // 7
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1330,7 +1693,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0, 0.0, 0.0, // 15
+ 0.0, 0.0,-1.0, // 15
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1343,7 +1706,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0, 0.0, // 8
+ +1.0, 0.0, 0.0, // 8
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1351,7 +1714,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0, 0.0, // 16
+ -1.0, 0.0, 0.0, // 16
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1363,7 +1726,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0, // 8
+ 0.0,+1.0, 0.0, // 8
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1371,7 +1734,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0, // 16
+ 0.0,-1.0, 0.0, // 16
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1383,7 +1746,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +1.0, 0.0, 0.0, // 8
+ 0.0, 0.0,+1.0, // 8
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1391,7 +1754,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0, 0.0, 0.0, // 16
+ 0.0, 0.0,-1.0, // 16
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1404,7 +1767,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0, 0.0, // 9
+ +1.0, 0.0, 0.0, // 9
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1412,7 +1775,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0, 0.0, // 17
+ -1.0, 0.0, 0.0, // 17
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1424,7 +1787,7 @@
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.0,+1.0, 0.0, // 9
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1432,7 +1795,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0, // 17
+ 0.0,-1.0, 0.0, // 17
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1444,7 +1807,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +1.0, 0.0, 0.0, // 9
+ 0.0, 0.0,+1.0, // 9
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
@@ -1452,7 +1815,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0, 0.0, 0.0, // 17
+ 0.0, 0.0,-1.0, // 17
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
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-09-22 22:14:21 UTC (rev 18963)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTet4e.cc 2011-09-23 00:11:07 UTC (rev 18964)
@@ -141,36 +141,126 @@
const double pylith::faults::CohesiveKinDataTet4e::_residual[] = {
0.0, 0.0, 0.0,
- 7.8, 3.8, 5.8, // 5
- 7.0, 3.0, 5.0, // 6
- 8.2, 4.2, 6.2, // 7
- 8.4, 4.4, 6.4, // 8
+ -(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
0.0, 0.0, 0.0,
- -7.8, -3.8, -5.8, // 10
- -7.0, -3.0, -5.0, // 11
- -8.2, -4.2, -6.2, // 12
- -8.4, -4.4, -6.4, // 13
- -0.5+1.82575588523, -0.5+-0.55566483464, -0.5+0.07938069066, // 14
- -0.6+1.69682900001, -0.6+-0.56560966667, -0.6+0.14140241667, // 15
- -0.7+1.51709826228, -0.7-0.54615537442, -0.7+0.18205179147, // 16
- -0.8+1.29378670385, -0.8-0.49761027071, -0.8+0.19904410828, // 17
+ +(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
+
+ (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
+
+ (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
+
+ (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
+
+ (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
+
};
const double pylith::faults::CohesiveKinDataTet4e::_residualIncr[] = {
0.0, 0.0, 0.0,
- 7.8, 3.8, 5.8, // 5
- 7.0, 3.0, 5.0, // 6
- 8.2, 4.2, 6.2, // 7
- 8.4, 4.4, 6.4, // 8
+ -(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
0.0, 0.0, 0.0,
- -7.8, -3.8, -5.8, // 10
- -7.0, -3.0, -5.0, // 11
- -8.2, -4.2, -6.2, // 12
- -8.4, -4.4, -6.4, // 13
- -0.5+1.82575588523, -0.5+-0.55566483464, -0.5+0.07938069066, // 14
- -0.6+1.69682900001, -0.6+-0.56560966667, -0.6+0.14140241667, // 15
- -0.7+1.51709826228, -0.7-0.54615537442, -0.7+0.18205179147, // 16
- -0.8+1.29378670385, -0.8-0.49761027071, -0.8+0.19904410828, // 17
+ +(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
+
+ (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
+
+ (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
+
+ (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
+
+ (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
+
};
const double pylith::faults::CohesiveKinDataTet4e::_jacobian[] = {
@@ -226,7 +316,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0, // 14
+ +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,
@@ -240,7 +330,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0, 0.0, 0.0, // 14
+ 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,
@@ -254,7 +344,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0, 0.0, // 14
+ 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,
@@ -269,7 +359,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0, // 15
+ +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
@@ -283,7 +373,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0, 0.0, 0.0, // 15
+ 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
@@ -297,7 +387,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0, 0.0, // 15
+ 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
@@ -312,7 +402,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0, // 16
+ +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,
@@ -326,7 +416,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0, 0.0, 0.0, // 16
+ 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,
@@ -340,7 +430,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0, 0.0, // 16
+ 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,
@@ -355,7 +445,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0, // 17
+ +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,
@@ -369,7 +459,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0, 0.0, 0.0, // 17
+ 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,
@@ -383,7 +473,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0, 0.0, // 17
+ 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,
@@ -436,7 +526,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0, // 14
+ -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,
@@ -450,7 +540,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +1.0, 0.0, 0.0, // 14
+ 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,
@@ -464,7 +554,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0, 0.0, // 14
+ 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,
@@ -479,7 +569,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0, // 15
+ -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
@@ -493,7 +583,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +1.0, 0.0, 0.0, // 15
+ 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
@@ -507,7 +597,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0, 0.0, // 15
+ 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
@@ -522,7 +612,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0, // 16
+ -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,
@@ -536,7 +626,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +1.0, 0.0, 0.0, // 16
+ 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,
@@ -550,7 +640,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0, 0.0, // 16
+ 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,
@@ -565,7 +655,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0, // 17
+ -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,
@@ -579,7 +669,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +1.0, 0.0, 0.0, // 17
+ 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,
@@ -593,14 +683,14 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0, 0.0, // 17
+ 0.0, 0.0,-1.0/3.0, // 17
0.0, 0.0, 0.0, // 14x
- 0.0,-1.0, 0.0, // 5
+ +2.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,
0.0, 0.0, 0.0,
- 0.0,+1.0, 0.0, // 10
+ -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,
@@ -609,12 +699,12 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 14y
- 0.0, 0.0,-1.0, // 5
+ 0.0,+2.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, 0.0, 0.0,
- 0.0, 0.0,+1.0, // 10
+ 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,
@@ -623,12 +713,12 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 14z
- -1.0, 0.0, 0.0, // 5
+ 0.0, 0.0,+2.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, 0.0,
- +1.0, 0.0, 0.0, // 10
+ 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,
@@ -638,12 +728,12 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 15x
0.0, 0.0, 0.0,
- 0.0,-1.0, 0.0, // 6
+ +1.0/3.0, 0.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.0, 0.0,
- 0.0,+1.0, 0.0, // 11
+ -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,
@@ -652,12 +742,12 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 15y
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0, // 6
+ 0.0,+1.0/3.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.0, 0.0,
- 0.0, 0.0,+1.0, // 11
+ 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,
@@ -666,12 +756,12 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 15z
0.0, 0.0, 0.0,
- -1.0, 0.0, 0.0, // 6
+ 0.0, 0.0,+1.0/3.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.0,
- +1.0, 0.0, 0.0, // 11
+ 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,
@@ -681,12 +771,12 @@
0.0, 0.0, 0.0, // 16x
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0, 0.0, // 7
+ +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,
- 0.0,+1.0, 0.0, // 12
+ -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,
@@ -695,12 +785,12 @@
0.0, 0.0, 0.0, // 16y
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0, // 7
+ 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, 0.0,+1.0, // 12
+ 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,
@@ -709,12 +799,12 @@
0.0, 0.0, 0.0, // 16z
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0, 0.0, 0.0, // 7
+ 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,
- +1.0, 0.0, 0.0, // 12
+ 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,
@@ -724,12 +814,12 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0, 0.0, // 8
+ +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,+1.0, 0.0, // 13
+ -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,
@@ -738,12 +828,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, // 8
+ 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,+1.0, // 13
+ 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,
@@ -752,12 +842,12 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0, 0.0, 0.0, // 8
+ 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,
- +1.0, 0.0, 0.0, // 13
+ 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,
Modified: short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTet4f.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTet4f.cc 2011-09-22 22:14:21 UTC (rev 18963)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTet4f.cc 2011-09-23 00:11:07 UTC (rev 18964)
@@ -135,30 +135,60 @@
const double pylith::faults::CohesiveKinDataTet4f::_residual[] = {
0.0, 0.0, 0.0,
- -9.7, -7.7, +8.7, // 3
- -9.9, -7.9, +8.9, // 4
- -9.1, -7.1, +8.1, // 5
+ -(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
0.0, 0.0, 0.0,
- +9.7, +7.7, -8.7, // 7
- +9.9, +7.9, -8.9, // 8
- +9.1, +7.1, -8.1, // 9
- 0.4+1.82575588523, -0.4+-0.55566483464, 0.4+0.07938069066, // 10
- 0.5+1.69682900001, -0.5+-0.56560966667, 0.5+0.14140241667, // 11
- -0.4+1.51709826228, +0.4+-0.54615537442, -0.4+0.18205179147, // 12
+ +(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
};
const double pylith::faults::CohesiveKinDataTet4f::_residualIncr[] = {
0.0, 0.0, 0.0,
- -9.7, -7.7, +8.7, // 3
- -9.9, -7.9, +8.9, // 4
- -9.1, -7.1, +8.1, // 5
+ -(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
0.0, 0.0, 0.0,
- +9.7, +7.7, -8.7, // 7
- +9.9, +7.9, -8.9, // 8
- +9.1, +7.1, -8.1, // 9
- 0.4+1.82575588523, -0.4+-0.55566483464, 0.4+0.07938069066, // 10
- 0.5+1.69682900001, -0.5+-0.56560966667, 0.5+0.14140241667, // 11
- -0.4+1.51709826228, +0.4+-0.54615537442, -0.4+0.18205179147, // 12
+ +(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
};
const double pylith::faults::CohesiveKinDataTet4f::_jacobian[] = {
@@ -203,7 +233,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0, // 10
+ +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
@@ -214,7 +244,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +1.0, 0.0, 0.0, // 10
+ 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
@@ -225,7 +255,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0, 0.0, // 10
+ 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
@@ -237,7 +267,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0, // 11
+ +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,
@@ -248,7 +278,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +1.0, 0.0, 0.0, // 11
+ 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,
@@ -259,7 +289,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0, 0.0, // 11
+ 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,
@@ -271,7 +301,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0, // 12
+ +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,
@@ -282,7 +312,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- +1.0, 0.0, 0.0, // 12
+ 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,
@@ -293,7 +323,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0, 0.0, // 12
+ 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,
@@ -335,7 +365,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0, // 10
+ -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
@@ -346,7 +376,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0, 0.0, 0.0, // 10
+ 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
@@ -357,7 +387,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0, 0.0, // 10
+ 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
@@ -369,7 +399,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0, // 11
+ -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,
@@ -380,7 +410,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0, 0.0, 0.0, // 11
+ 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,
@@ -391,7 +421,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0, 0.0, // 11
+ 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,
@@ -403,7 +433,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0, // 12
+ -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,
@@ -414,7 +444,7 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0, 0.0, 0.0, // 12
+ 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,
@@ -425,35 +455,35 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0, 0.0, // 12
+ 0.0, 0.0,-1.0/3.0, // 12
0.0, 0.0, 0.0, // 10x
- 0.0,+1.0, 0.0, // 3
+ +1.0/3.0, 0.0, 0.0, // 3
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,-1.0, 0.0, // 7
+ -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, 0.0,-1.0, // 3
+ 0.0,+1.0/3.0, 0.0, // 3
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,+1.0, // 7
+ 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
- +1.0, 0.0, 0.0, // 3
+ 0.0, 0.0,+1.0/3.0, // 3
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- -1.0, 0.0, 0.0, // 7
+ 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,
@@ -461,33 +491,33 @@
0.0, 0.0, 0.0,
0.0, 0.0, 0.0, // 11x
0.0, 0.0, 0.0,
- 0.0,+1.0, 0.0, // 4
+ +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,
- 0.0,-1.0, 0.0, // 8
+ -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, 0.0, 0.0,
- 0.0, 0.0,-1.0, // 4
+ 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, 0.0,+1.0, // 8
+ 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, 0.0,
- +1.0, 0.0, 0.0, // 4
+ 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,
- -1.0, 0.0, 0.0, // 8
+ 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,
@@ -495,33 +525,33 @@
0.0, 0.0, 0.0, // 12x
0.0, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0,+1.0, 0.0, // 5
+ +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,
- 0.0,-1.0, 0.0, // 9
+ -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, 0.0, 0.0,
0.0, 0.0, 0.0,
- 0.0, 0.0,-1.0, // 5
+ 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, 0.0,+1.0, // 9
+ 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, 0.0,
0.0, 0.0, 0.0,
- +1.0, 0.0, 0.0, // 5
+ 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,
- -1.0, 0.0, 0.0, // 9
+ 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,
More information about the CIG-COMMITS
mailing list