[cig-commits] r18975 - in short/3D/PyLith/branches/v1.6-revisedfault: libsrc/pylith/faults playpen/faultpc unittests/libtests/faults unittests/libtests/faults/data

brad at geodynamics.org brad at geodynamics.org
Tue Sep 27 10:54:56 PDT 2011


Author: brad
Date: 2011-09-27 10:54:55 -0700 (Tue, 27 Sep 2011)
New Revision: 18975

Modified:
   short/3D/PyLith/branches/v1.6-revisedfault/libsrc/pylith/faults/FaultCohesiveLagrange.cc
   short/3D/PyLith/branches/v1.6-revisedfault/playpen/faultpc/checkfaultpc.py
   short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/TestFaultCohesiveKin.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/CohesiveKinDataLine2.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/CohesiveKinDataTet4f.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
   short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinSrcsDataHex8.cc
   short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinSrcsDataQuad4.cc
   short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinSrcsDataTet4.cc
   short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinSrcsDataTri3.cc
   short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/adjustsoln.py
Log:
Fixed bug in integrating Jacobian for fault (jB indices depend on jBasis) and adjusting solution space. Updated test data accordingly.

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-26 21:38:49 UTC (rev 18974)
+++ short/3D/PyLith/branches/v1.6-revisedfault/libsrc/pylith/faults/FaultCohesiveLagrange.cc	2011-09-27 17:54:55 UTC (rev 18975)
@@ -373,7 +373,7 @@
 	      (dispTpdtCell[jBP + iDim] - dispTpdtCell[jBN + iDim] -
 	       slipCell[jBasis*spaceDim+iDim]);
 
-#if 1
+#if 0
 	    std::cout << "iBasis: " << iBasis
 		      << ", jBasis: " << jBasis
 		      << ", iDim: " << iDim
@@ -532,13 +532,19 @@
       const double wt = quadWts[iQuad] * jacobianDet[iQuad];
       for (int iBasis=0, iQ=iQuad*numBasis; iBasis < numBasis; ++iBasis) {
         const double valI = wt*basis[iQ+iBasis];
+
+	// First index for positive, negative, and Lagrange vertices
+	const int iBN = 0*numBasis*spaceDim + iBasis*spaceDim;
+	const int iBP = 1*numBasis*spaceDim + iBasis*spaceDim;
+	const int iBL = 2*numBasis*spaceDim + iBasis*spaceDim;
+
         for (int jBasis=0; jBasis < numBasis; ++jBasis) {
           const double valIJ = valI * basis[iQ+jBasis];
 
 	  // First index for positive, negative, and Lagrange vertices
-	  const int iBN = 0*numBasis*spaceDim + iBasis*spaceDim;
-	  const int iBP = 1*numBasis*spaceDim + iBasis*spaceDim;
-	  const int iBL = 2*numBasis*spaceDim + iBasis*spaceDim;
+	  const int jBN = 0*numBasis*spaceDim + jBasis*spaceDim;
+	  const int jBP = 1*numBasis*spaceDim + jBasis*spaceDim;
+	  const int jBL = 2*numBasis*spaceDim + jBasis*spaceDim;
 
 	  for (int iDim=0; iDim < spaceDim; ++iDim) {
 	    // Add entries to Jacobian at (i,j) where
@@ -549,15 +555,15 @@
 
 	    // Indices for negative vertex
             const int iN = (iBN + iDim) * rowSize; // row
-            const int jN = (iBN + iDim); // col
+            const int jN = (jBN + iDim); // col
 
 	    // Indices for positive vertex
             const int iP = (iBP + iDim) * rowSize; // row
-            const int jP = (iBP + iDim); // col
+            const int jP = (jBP + iDim); // col
 
 	    // Indices for Lagrange vertex
             const int iL = (iBL + iDim) * rowSize; // row
-            const int jL = (iBL + iDim); // col
+            const int jL = (jBL + iDim); // col
 
             jacobianCell[iN + jL] += valIJ;
             jacobianCell[iL + jN] += valIJ;
@@ -1413,9 +1419,9 @@
 
     for (int iDim=0; iDim < spaceDim; ++iDim) {
       assert(jacobianVertexL[iDim] > 0.0);
-      const double Sinv = (jacobianVertexP[iDim] + jacobianVertexN[iDim]) /
-	(jacobianVertexL[iDim]*jacobianVertexL[iDim]);
-      dispTIncrVertexL[iDim] = Sinv * 
+      const double S = (1.0/jacobianVertexP[iDim] + 1.0/jacobianVertexN[iDim]) *
+	jacobianVertexL[iDim]*jacobianVertexL[iDim];
+      dispTIncrVertexL[iDim] = 1.0/S * 
 	(-residualVertexL[iDim] +
 	 jacobianVertexL[iDim] * 
 	 (dispTIncrVertexP[iDim] - dispTIncrVertexN[iDim]));
@@ -2014,8 +2020,6 @@
     _fields->get("slip").section();
   assert(!slipSection.isNull());
 
-  slipSection->view("SLIP (FAULT)");
-
   const ALE::Obj<RealSection>& orientationSection =
     _fields->get("orientation").section();
   assert(!orientationSection.isNull());
@@ -2045,7 +2049,7 @@
   
   PetscLogFlops(numVertices * (2*spaceDim*spaceDim) );
 
-#if 1 // DEBUGGING
+#if 0 // DEBUGGING
   slipSection->view("SLIP (GLOBAL)");
 #endif
 } // _slipFaultToGlobal

Modified: short/3D/PyLith/branches/v1.6-revisedfault/playpen/faultpc/checkfaultpc.py
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/playpen/faultpc/checkfaultpc.py	2011-09-26 21:38:49 UTC (rev 18974)
+++ short/3D/PyLith/branches/v1.6-revisedfault/playpen/faultpc/checkfaultpc.py	2011-09-27 17:54:55 UTC (rev 18975)
@@ -3,7 +3,7 @@
 import numpy
 import numpy.linalg as linalg
 
-A = numpy.array([[4/3.0, -0.5, 1/6.0, 0,   0,0,0,0],
+K = numpy.array([[4/3.0, -0.5, 1/6.0, 0,   0,0,0,0],
                  [-0.5, 4/3.0, 0, -5/6.0,  0,0,0,0],
                  [1/6.0, 0, 4/3.0, 0.5,    0,0,0,0],
                  [0, -5/6.0, 0.5, 4/3.0,   0,0,0,0],
@@ -12,47 +12,47 @@
                  [0,0,0,0, 1/6.0, 0, 4/3.0, -0.5],
                  [0,0,0,0, 0, -5/6.0, -0.5, 4/3.0]],
                 dtype=numpy.float64)
-Ai = linalg.inv(A)
+Ki = linalg.inv(K)
 
-C = numpy.array([[0,-1,  0,0,  0,1, 0,0],
-                 [-1,0,  0,0,  1,0, 0,0],
-                 [0,0,  0,-1,  0,0, 0,1],
-                 [0,0, -1,0,   0,0, 1,0]],
+L = numpy.array([[1,0, 0,0, -1,0, 0,0],
+                 [0,1, 0,0, 0,-1, 0,0],
+                 [0,0, 1,0,  0,0, -1,0],
+                 [0,0, 0,1,  0,0, 0,-1]],
                 dtype=numpy.float64)
 Z = numpy.zeros( (4,4), dtype=numpy.float64)
 
-J = numpy.vstack( (numpy.hstack( (A, C.transpose()) ),
-                   numpy.hstack( (C, Z) ) ) )
-Jinv = linalg.inv(J)
+A = numpy.vstack( (numpy.hstack( (K, L.transpose()) ),
+                   numpy.hstack( (L, Z) ) ) )
+Ainv = linalg.inv(A)
 
-# Compute [C] [A]^(-1) [C]^T and its inverse.
-CAC = numpy.dot(numpy.dot(C, Ai), C.transpose())
-CACi = numpy.linalg.inv(CAC)
+# Compute [L] [K]^(-1) [L]^T and its inverse.
+LKiL = numpy.dot(numpy.dot(L, Ki), L.transpose())
+LKiLi = numpy.linalg.inv(LKiL)
 
-# Compute diagonal approximation of CAC and its inverse
-Aid = 1.0 / A.diagonal() * numpy.identity(A.shape[0])
-CACd = numpy.dot(numpy.dot(C, Aid), C.transpose())
-CACdi = 1.0 / CACd.diagonal() * numpy.identity(CACd.shape[0])
+# Compute diagonal approximation of LKiL and its inverse
+Kid = 1.0 / K.diagonal() * numpy.identity(K.shape[0])
+LKiLd = numpy.dot(numpy.dot(L, Kid), L.transpose())
+LKiLdi = 1.0 / LKiLd.diagonal() * numpy.identity(LKiLd.shape[0])
 
 # Compute preconditioner using full matrices (no approximations)
-P = J
+P = A
 Pi = numpy.linalg.inv(P)
 
 # Compute condition number
-evals, evecs = numpy.linalg.eig(numpy.dot(J, Pi))
+evals, evecs = numpy.linalg.eig(numpy.dot(A, Pi))
 print numpy.abs(evals)
 print numpy.max(numpy.abs(evals))/numpy.min(numpy.abs(evals))
 
 # Compute preconditioner using diagonal approximations (but full A)
-Pd = numpy.zeros(J.shape)
-Pd[0:8,0:8] = A
+Pd = numpy.zeros(A.shape)
+Pd[0:8,0:8] = K
 Pd[0:8,8:12] = 0.0
 Pd[8:12,0:8] = 0.0
-Pd[8:12,8:12] = -CACd
+Pd[8:12,8:12] = -LKiLd
 
 Pdi = numpy.linalg.inv(Pd)
 
 # Compute condition number for diagonal approximations
-evals, evecs = numpy.linalg.eig(numpy.dot(Pdi, J))
+evals, evecs = numpy.linalg.eig(numpy.dot(Pdi, A))
 print numpy.abs(evals)
 print numpy.max(numpy.abs(evals))/numpy.min(numpy.abs(evals))

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-09-26 21:38:49 UTC (rev 18974)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/TestFaultCohesiveKin.cc	2011-09-27 17:54:55 UTC (rev 18975)
@@ -492,7 +492,7 @@
     
     for (int iDim=0; iDim < spaceDim; ++iDim) {
       const double valE = _data->jacobianLumped[iVertex*spaceDim+iDim];
-#if 1 // debugging
+#if 0 // debugging
       std::cout << "vertex: " << *v_iter << ", iDim: " << iDim
 		<< ", valE: " << valE
 		<< ", val: " << vals[iDim]
@@ -601,7 +601,6 @@
     const double* solutionVertex = solutionSection->restrictPoint(*v_iter);
     CPPUNIT_ASSERT(0 != solutionVertex);
     for (int iDim=0; iDim < spaceDim; ++iDim, ++i) {
-      std::cout << "valE: " << solutionE[i] << ", val: " << solutionVertex[iDim] << std::endl;
       if (0.0 != solutionE[i])
         CPPUNIT_ASSERT_DOUBLES_EQUAL(1.0, solutionVertex[iDim]/solutionE[i],
           tolerance);

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-26 21:38:49 UTC (rev 18974)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataHex8.cc	2011-09-27 17:54:55 UTC (rev 18975)
@@ -876,10 +876,10 @@
   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, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+ +4.0/9.0, 0.0, 0.0, // 18
+ +2.0/9.0, 0.0, 0.0, // 19
+ +2.0/9.0, 0.0, 0.0, // 20
+ +1.0/9.0, 0.0, 0.0, // 21
   0.0, 0.0, 0.0, // 6y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -896,10 +896,10 @@
   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, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+  0.0,+4.0/9.0, 0.0, // 18
+  0.0,+2.0/9.0, 0.0, // 19
+  0.0,+2.0/9.0, 0.0, // 20
+  0.0,+1.0/9.0, 0.0, // 21
   0.0, 0.0, 0.0, // 6z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -916,10 +916,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, // 18
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+  0.0, 0.0,+4.0/9.0, // 18
+  0.0, 0.0,+2.0/9.0, // 19
+  0.0, 0.0,+2.0/9.0, // 20
+  0.0, 0.0,+1.0/9.0, // 21
   0.0, 0.0, 0.0, // 7x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -936,10 +936,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,
- +1.0, 0.0, 0.0, // 19
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+ +2.0/9.0, 0.0, 0.0, // 18
+ +4.0/9.0, 0.0, 0.0, // 19
+ +1.0/9.0, 0.0, 0.0, // 20
+ +2.0/9.0, 0.0, 0.0, // 21
   0.0, 0.0, 0.0, // 7y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -956,10 +956,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,+1.0, 0.0, // 19
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+  0.0,+2.0/9.0, 0.0, // 18
+  0.0,+4.0/9.0, 0.0, // 19
+  0.0,+1.0/9.0, 0.0, // 20
+  0.0,+2.0/9.0, 0.0, // 21
   0.0, 0.0, 0.0, // 7z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -976,10 +976,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,+1.0, // 19
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+  0.0, 0.0,+2.0/9.0, // 18
+  0.0, 0.0,+4.0/9.0, // 19
+  0.0, 0.0,+1.0/9.0, // 20
+  0.0, 0.0,+2.0/9.0, // 21
   0.0, 0.0, 0.0, // 8x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -996,10 +996,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,
- +1.0, 0.0, 0.0, // 20
-  0.0, 0.0, 0.0,
+ +2.0/9.0, 0.0, 0.0, // 18
+ +1.0/9.0, 0.0, 0.0, // 19
+ +4.0/9.0, 0.0, 0.0, // 20
+ +2.0/9.0, 0.0, 0.0, // 21
   0.0, 0.0, 0.0, // 8y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1016,10 +1016,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,+1.0, 0.0, // 20
-  0.0, 0.0, 0.0,
+  0.0,+2.0/9.0, 0.0, // 18
+  0.0,+1.0/9.0, 0.0, // 19
+  0.0,+4.0/9.0, 0.0, // 20
+  0.0,+2.0/9.0, 0.0, // 21
   0.0, 0.0, 0.0, // 8z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1036,10 +1036,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,+1.0, // 20
-  0.0, 0.0, 0.0,
+  0.0, 0.0,+2.0/9.0, // 18
+  0.0, 0.0,+1.0/9.0, // 19
+  0.0, 0.0,+4.0/9.0, // 20
+  0.0, 0.0,+2.0/9.0, // 21
   0.0, 0.0, 0.0, // 9x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1056,10 +1056,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,
- +1.0, 0.0, 0.0, // 21
+ +1.0/9.0, 0.0, 0.0, // 18
+ +2.0/9.0, 0.0, 0.0, // 19
+ +2.0/9.0, 0.0, 0.0, // 20
+ +4.0/9.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,
@@ -1076,10 +1076,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,+1.0, 0.0, // 21
+  0.0,+1.0/9.0, 0.0, // 18
+  0.0,+2.0/9.0, 0.0, // 19
+  0.0,+2.0/9.0, 0.0, // 20
+  0.0,+4.0/9.0, 0.0, // 21
   0.0, 0.0, 0.0, // 9z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1096,10 +1096,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,+1.0, // 21
+  0.0, 0.0,+1.0/9.0, // 18
+  0.0, 0.0,+2.0/9.0, // 19
+  0.0, 0.0,+2.0/9.0, // 20
+  0.0, 0.0,+4.0/9.0, // 21
   0.0, 0.0, 0.0, // 10x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1356,10 +1356,10 @@
   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, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+ -4.0/9.0, 0.0, 0.0, // 18
+ -2.0/9.0, 0.0, 0.0, // 19
+ -2.0/9.0, 0.0, 0.0, // 20
+ -1.0/9.0, 0.0, 0.0, // 21
   0.0, 0.0, 0.0, // 14y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1376,10 +1376,10 @@
   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, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+  0.0,-4.0/9.0, 0.0, // 18
+  0.0,-2.0/9.0, 0.0, // 19
+  0.0,-2.0/9.0, 0.0, // 20
+  0.0,-1.0/9.0, 0.0, // 21
   0.0, 0.0, 0.0, // 14z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1396,10 +1396,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, // 18
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+  0.0, 0.0,-4.0/9.0, // 18
+  0.0, 0.0,-2.0/9.0, // 19
+  0.0, 0.0,-2.0/9.0, // 20
+  0.0, 0.0,-1.0/9.0, // 21
   0.0, 0.0, 0.0, // 15x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1416,10 +1416,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,
- -1.0, 0.0, 0.0, // 19
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+ -2.0/9.0, 0.0, 0.0, // 18
+ -4.0/9.0, 0.0, 0.0, // 19
+ -1.0/9.0, 0.0, 0.0, // 20
+ -2.0/9.0, 0.0, 0.0, // 21
   0.0, 0.0, 0.0, // 15y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1436,10 +1436,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,-1.0, 0.0, // 19
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+  0.0,-2.0/9.0, 0.0, // 18
+  0.0,-4.0/9.0, 0.0, // 19
+  0.0,-1.0/9.0, 0.0, // 20
+  0.0,-2.0/9.0, 0.0, // 21
   0.0, 0.0, 0.0, // 15z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1456,10 +1456,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,-1.0, // 19
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+  0.0, 0.0,-2.0/9.0, // 18
+  0.0, 0.0,-4.0/9.0, // 19
+  0.0, 0.0,-1.0/9.0, // 20
+  0.0, 0.0,-2.0/9.0, // 21
   0.0, 0.0, 0.0, // 16x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1476,10 +1476,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,
- -1.0, 0.0, 0.0, // 20
-  0.0, 0.0, 0.0,
+ -2.0/9.0, 0.0, 0.0, // 18
+ -1.0/9.0, 0.0, 0.0, // 19
+ -4.0/9.0, 0.0, 0.0, // 20
+ -2.0/9.0, 0.0, 0.0, // 21
   0.0, 0.0, 0.0, // 16y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1496,10 +1496,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,-1.0, 0.0, // 20
-  0.0, 0.0, 0.0,
+  0.0,-2.0/9.0, 0.0, // 18
+  0.0,-1.0/9.0, 0.0, // 19
+  0.0,-4.0/9.0, 0.0, // 20
+  0.0,-2.0/9.0, 0.0, // 21
   0.0, 0.0, 0.0, // 16z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1516,10 +1516,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,-1.0, // 20
-  0.0, 0.0, 0.0,
+  0.0, 0.0,-2.0/9.0, // 18
+  0.0, 0.0,-1.0/9.0, // 19
+  0.0, 0.0,-4.0/9.0, // 20
+  0.0, 0.0,-2.0/9.0, // 21
   0.0, 0.0, 0.0, // 17x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1536,10 +1536,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,
- -1.0, 0.0, 0.0, // 21
+ -1.0/9.0, 0.0, 0.0, // 18
+ -2.0/9.0, 0.0, 0.0, // 19
+ -2.0/9.0, 0.0, 0.0, // 20
+ -4.0/9.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,
@@ -1556,10 +1556,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,-1.0, 0.0, // 21
+  0.0,-1.0/9.0, 0.0, // 18
+  0.0,-2.0/9.0, 0.0, // 19
+  0.0,-2.0/9.0, 0.0, // 20
+  0.0,-4.0/9.0, 0.0, // 21
   0.0, 0.0, 0.0, // 17z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1576,250 +1576,250 @@
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
   0.0, 0.0, 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, // 21
+  0.0, 0.0,-1.0/9.0, // 18
+  0.0, 0.0,-2.0/9.0, // 19
+  0.0, 0.0,-2.0/9.0, // 20
+  0.0, 0.0,-4.0/9.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,
- +1.0, 0.0, 0.0, // 6
+ +4.0/9.0, 0.0, 0.0, // 6
+ +2.0/9.0, 0.0, 0.0, // 7
+ +2.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,
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
+ -4.0/9.0, 0.0, 0.0, // 14
+ -2.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,
   0.0, 0.0, 0.0,
- -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,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 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.0,
   0.0, 0.0, 0.0,
-  0.0,+1.0, 0.0, // 6
+  0.0,+4.0/9.0, 0.0, // 6
+  0.0,+2.0/9.0, 0.0, // 7
+  0.0,+2.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, 0.0, 0.0,
   0.0, 0.0, 0.0,
+  0.0,-4.0/9.0, 0.0, // 14
+  0.0,-2.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, 0.0, 0.0,
-  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,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
   0.0, 0.0, 0.0, // 18z
   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, 0.0,+4.0/9.0, // 6
+  0.0, 0.0,+2.0/9.0, // 7
+  0.0, 0.0,+2.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, 0.0,
   0.0, 0.0, 0.0,
+  0.0, 0.0,-4.0/9.0, // 14
+  0.0, 0.0,-2.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, 0.0,
-  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,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
   0.0, 0.0, 0.0, // 19x
   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, // 6
+ +4.0/9.0, 0.0, 0.0, // 7
+ +1.0/9.0, 0.0, 0.0, // 8
+ +2.0/9.0, 0.0, 0.0, // 9
   0.0, 0.0, 0.0,
- +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,
+ -2.0/9.0, 0.0, 0.0, // 14
+ -4.0/9.0, 0.0, 0.0, // 15
+ -1.0/9.0, 0.0, 0.0, // 16
+ -2.0/9.0, 0.0, 0.0, // 17
   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, // 15
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
   0.0, 0.0, 0.0, // 19y
   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, // 6
+  0.0,+4.0/9.0, 0.0, // 7
+  0.0,+1.0/9.0, 0.0, // 8
+  0.0,+2.0/9.0, 0.0, // 9
   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,
   0.0, 0.0, 0.0,
+  0.0,-2.0/9.0, 0.0, // 14
+  0.0,-4.0/9.0, 0.0, // 15
+  0.0,-1.0/9.0, 0.0, // 16
+  0.0,-2.0/9.0, 0.0, // 17
   0.0, 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, // 15
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
   0.0, 0.0, 0.0, // 19z
   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, // 6
+  0.0, 0.0,+4.0/9.0, // 7
+  0.0, 0.0,+1.0/9.0, // 8
+  0.0, 0.0,+2.0/9.0, // 9
   0.0, 0.0, 0.0,
-  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,
+  0.0, 0.0,-2.0/9.0, // 14
+  0.0, 0.0,-4.0/9.0, // 15
+  0.0, 0.0,-1.0/9.0, // 16
+  0.0, 0.0,-2.0/9.0, // 17
   0.0, 0.0, 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, // 15
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
   0.0, 0.0, 0.0, // 20x
   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, // 6
+ +1.0/9.0, 0.0, 0.0, // 7
+ +4.0/9.0, 0.0, 0.0, // 8
+ +2.0/9.0, 0.0, 0.0, // 9
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
- +1.0, 0.0, 0.0, // 8
   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
+ -4.0/9.0, 0.0, 0.0, // 16
+ -2.0/9.0, 0.0, 0.0, // 17
   0.0, 0.0, 0.0,
   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, // 16
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
   0.0, 0.0, 0.0, // 20y
   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, // 6
+  0.0,+1.0/9.0, 0.0, // 7
+  0.0,+4.0/9.0, 0.0, // 8
+  0.0,+2.0/9.0, 0.0, // 9
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
-  0.0,+1.0, 0.0, // 8
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
+  0.0,-2.0/9.0, 0.0, // 6
+  0.0,-1.0/9.0, 0.0, // 7
+  0.0,-4.0/9.0, 0.0, // 8
+  0.0,-2.0/9.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,
-  0.0, 0.0, 0.0,
-  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,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
   0.0, 0.0, 0.0, // 20z
   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, // 6
+  0.0, 0.0,+1.0/9.0, // 7
+  0.0, 0.0,+4.0/9.0, // 8
+  0.0, 0.0,+2.0/9.0, // 9
   0.0, 0.0, 0.0,
   0.0, 0.0, 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,-2.0/9.0, // 14
+  0.0, 0.0,-1.0/9.0, // 15
+  0.0, 0.0,-4.0/9.0, // 16
+  0.0, 0.0,-2.0/9.0, // 17
   0.0, 0.0, 0.0,
   0.0, 0.0, 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, // 16
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
   0.0, 0.0, 0.0, // 21x
   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, // 6
+ +2.0/9.0, 0.0, 0.0, // 7
+ +2.0/9.0, 0.0, 0.0, // 8
+ +4.0/9.0, 0.0, 0.0, // 9
   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, 0.0,
+ -1.0/9.0, 0.0, 0.0, // 14
+ -2.0/9.0, 0.0, 0.0, // 15
+ -2.0/9.0, 0.0, 0.0, // 16
+ -4.0/9.0, 0.0, 0.0, // 17
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
   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, // 17
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
   0.0, 0.0, 0.0, // 21y
   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, // 6
+  0.0,+2.0/9.0, 0.0, // 7
+  0.0,+2.0/9.0, 0.0, // 8
+  0.0,+4.0/9.0, 0.0, // 9
   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
   0.0, 0.0, 0.0,
+  0.0,-1.0/9.0, 0.0, // 14
+  0.0,-2.0/9.0, 0.0, // 15
+  0.0,-2.0/9.0, 0.0, // 16
+  0.0,-4.0/9.0, 0.0, // 17
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
   0.0, 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, // 17
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
   0.0, 0.0, 0.0, // 21z
   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, // 6
+  0.0, 0.0,+2.0/9.0, // 7
+  0.0, 0.0,+2.0/9.0, // 8
+  0.0, 0.0,+4.0/9.0, // 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.0, 0.0, 0.0,
+  0.0, 0.0,-1.0/9.0, // 14
+  0.0, 0.0,-2.0/9.0, // 15
+  0.0, 0.0,-2.0/9.0, // 16
+  0.0, 0.0,-4.0/9.0, // 17
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
   0.0, 0.0, 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, // 17
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
 };
 
 const double pylith::faults::CohesiveKinDataHex8::_fieldIncrAdjusted[] = {
@@ -1827,22 +1827,22 @@
   3.2, 4.2, 5.2,
   3.3, 4.3, 5.3,
   3.4, 4.4, 5.4,
-  1.205353233, -0.669560195189, 2.4212859344, // 6
-  1.31708438297, -0.532975672131, 2.53979527703, // 7
-  1.44967569358, -0.34540170722, 2.68738185699, // 8
-  1.6022763433, -0.117061939422, 2.86063862983, // 9
+  2.92926515362, 3.21420377798, 4.73424841481, // 6
+  3.02986498638, 3.31809140342, 4.8357449183, // 7
+  3.1374189234, 3.4386495732, 4.94684546425, // 8
+  3.25097042324, 3.57163244464, 5.06569642907, // 9
   3.9, 4.9, 5.9,
   3.0, 4.0, 5.0,
   3.1, 4.1, 5.1,
   3.2, 4.2, 5.2,
-  5.94766934654, 10.2648771483, 8.85236238339, // 14
-  5.9351099915, 9.97517405027, 8.76421837117, // 15
-  5.95032430642, 9.74540170722, 8.71261814301, // 16
-  5.98205399056, 9.55826920577, 8.68465814016, // 17
-  -3.4419701505, -7.75434029278, -4.6180710984, // 18
-  -3.65266498725, -8.21276107541, -4.89632755675, // 19
-  -3.82555132091, -8.57718290227, -5.12145084311, // 20
-  -3.95590258206, -8.85071149096, -5.29085046631, // 21
+  3.95854020737, 5.78361102541, 6.18355952138, // 14
+  4.10814401453, 5.86736916968, 6.31520542048, // 15
+  4.2625810766, 5.9613504268, 6.45315453575, // 16
+  4.42013328325, 6.06371663139, 6.59565601456, // 17
+  -0.856102269575, -1.92869433303, -1.14862737779, // 18
+  -0.912216021789, -2.05105375452, -1.22280813072, // 19
+  -0.956387830228, -2.14429572557, -1.28036271078, // 20
+  -0.988253238176, -2.21106159964, -1.32174642767, // 21
 };
 
 pylith::faults::CohesiveKinDataHex8::CohesiveKinDataHex8(void)

Modified: short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataLine2.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataLine2.cc	2011-09-26 21:38:49 UTC (rev 18974)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataLine2.cc	2011-09-27 17:54:55 UTC (rev 18975)
@@ -154,10 +154,10 @@
 
 const double pylith::faults::CohesiveKinDataLine2::_fieldIncrAdjusted[] = {
   1.1,
-  5.16324319611, // 3
+  2.18893781075, // 3
   1.3,
-  -2.23297292977, // 5
-  8.719135031441999, // 6
+  0.49347367348, // 5
+  2.17566318365, // 6
 };
 
 pylith::faults::CohesiveKinDataLine2::CohesiveKinDataLine2(void)

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-09-26 21:38:49 UTC (rev 18974)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataQuad4.cc	2011-09-27 17:54:55 UTC (rev 18975)
@@ -247,8 +247,8 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
- +1.0, 0.0, // 10
-  0.0, 0.0,
+ +0.5, 0.0, // 10
+ +0.5, 0.0, // 11
   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,+1.0, //  10
-  0.0, 0.0,
+  0.0,+0.5, // 10
+  0.0,+0.5, // 11
   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.0, 0.0,
- +1.0, 0.0, //  11
+ +0.5, 0.0, // 10
+ +0.5, 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.0,
-  0.0,+1.0, //  11
+  0.0,+0.5, // 10
+  0.0,+0.5, // 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,
- -1.0, 0.0, //  10
-  0.0, 0.0,
+ -0.5, 0.0, // 10
+ -0.5, 0.0, // 11
   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,-1.0, // 10
-  0.0, 0.0,
+  0.0,-0.5, // 10
+  0.0,-0.5, // 11
   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.0, 0.0,
- -1.0, 0.0, // 11
+ -0.5, 0.0, // 10
+ -0.5, 0.0, // 11
   0.0, 0.0, // 9y
   0.0, 0.0,
   0.0, 0.0,
@@ -357,61 +357,61 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
-  0.0, 0.0,
-  0.0,-1.0, // 11
+  0.0,-0.5, // 10
+  0.0,-0.5, // 11
   0.0, 0.0, // 10x
   0.0, 0.0,
- +1.0, 0.0, //  4
+ +0.5, 0.0, // 4
+ +0.5, 0.0, // 5
   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,+1.0, //  4
+  0.0,+0.5, // 4
+  0.0,+0.5, // 5
   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[] = {
   3.1, 4.1,
   3.2, 4.2,
-  3.45041520801, 8.42020518055, // 4
-  3.78935267675, 8.84420957728, // 5
+  3.33693528997, 5.31173927211, // 4
+  3.4951035923, 5.48554614101, // 5
   3.5, 4.5,
   3.6, 4.6,
-  3.58497660564, 1.54925486193, // 8
-  3.61310855397, 1.62531925885, // 9
-  0.19553977041, 5.35626673472, // 10
-  0.545093747451, 6.22189340819, // 11
+  3.6717553665, 3.92631702721, // 8
+  3.82992366883, 4.1001238961, // 9
+  0.0480158769562, 1.31526105375, // 10
+  0.133145029221, 1.51976459741, // 11
 };
 
 pylith::faults::CohesiveKinDataQuad4::CohesiveKinDataQuad4(void)

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-09-26 21:38:49 UTC (rev 18974)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataQuad4e.cc	2011-09-27 17:54:55 UTC (rev 18975)
@@ -327,9 +327,9 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
- +1.0, 0.0, // 16
+ +0.5, 0.0, // 16
+ +0.5, 0.0, // 17
   0.0, 0.0,
-  0.0, 0.0,
   0.0, 0.0, // 6y
   0.0, 0.0,
   0.0, 0.0,
@@ -342,9 +342,9 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
-  0.0,+1.0, // 16
+  0.0,+0.5, // 16
+  0.0,+0.5, // 17
   0.0, 0.0,
-  0.0, 0.0,
   0.0, 0.0, // 7x
   0.0, 0.0,
   0.0, 0.0,
@@ -357,9 +357,9 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
-  0.0, 0.0,
- +2.0, 0.0, // 17
-  0.0, 0.0,
+ +0.5, 0.0, // 16
+ +1.0, 0.0, // 17
+ +0.5, 0.0, // 18
   0.0, 0.0, // 7y
   0.0, 0.0,
   0.0, 0.0,
@@ -372,9 +372,9 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
-  0.0, 0.0,
-  0.0,+2.0, // 17
-  0.0, 0.0,
+  0.0,+0.5, // 16
+  0.0,+1.0, // 17
+  0.0,+0.5, // 18
   0.0, 0.0, // 8x
   0.0, 0.0,
   0.0, 0.0,
@@ -478,8 +478,8 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
-  0.0, 0.0,
- +1.0, 0.0, // 18
+ +0.5, 0.0, // 17
+ +0.5, 0.0, // 18
   0.0, 0.0, // 11y
   0.0, 0.0,
   0.0, 0.0,
@@ -493,8 +493,8 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
-  0.0, 0.0,
-  0.0,+1.0, // 18
+  0.0,+0.5, // 17
+  0.0,+0.5, // 18
   0.0, 0.0, // 12x
   0.0, 0.0,
   0.0, 0.0,
@@ -537,9 +537,9 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
- -1.0, 0.0, // 16
+ -0.5, 0.0, // 16
+ -0.5, 0.0, // 17
   0.0, 0.0,
-  0.0, 0.0,
   0.0, 0.0, // 13y
   0.0, 0.0,
   0.0, 0.0,
@@ -552,9 +552,9 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
-  0.0,-1.0, // 16
+  0.0,-0.5, // 16
+  0.0,-0.5, // 17
   0.0, 0.0,
-  0.0, 0.0,
   0.0, 0.0, // 14x
   0.0, 0.0,
   0.0, 0.0,
@@ -567,9 +567,9 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
-  0.0, 0.0,
- -2.0, 0.0, // 17
-  0.0, 0.0,
+ -0.5, 0.0, // 16
+ -1.0, 0.0, // 17
+ -0.5, 0.0, // 18
   0.0, 0.0, // 14y
   0.0, 0.0,
   0.0, 0.0,
@@ -582,9 +582,9 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
-  0.0, 0.0,
-  0.0,-2.0, // 17
-  0.0, 0.0,
+  0.0,-0.5, // 16
+  0.0,-1.0, // 17
+  0.0,-0.5, // 18
   0.0, 0.0, // 15x
   0.0, 0.0,
   0.0, 0.0,
@@ -598,8 +598,8 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
-  0.0, 0.0,
- -1.0, 0.0, // 18
+ -0.5, 0.0, // 17
+ -0.5, 0.0, // 18
   0.0, 0.0, // 15y
   0.0, 0.0,
   0.0, 0.0,
@@ -613,98 +613,98 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
-  0.0, 0.0,
-  0.0,-1.0, // 18
+  0.0,-0.5, // 17
+  0.0,-0.5, // 18
   0.0, 0.0, // 16x
   0.0, 0.0,
- +1.0, 0.0, // 6
+ +0.5, 0.0, // 6
+ +0.5, 0.0, // 7
   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,+1.0, // 6
+  0.0,+0.5, // 6
+  0.0,+0.5, // 7
   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-09-26 21:38:49 UTC (rev 18974)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTet4.cc	2011-09-27 17:54:55 UTC (rev 18975)
@@ -261,9 +261,9 @@
   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, // 10
-  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, // 3y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -272,9 +272,9 @@
   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, // 10
-  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, // 3z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -283,9 +283,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/3.0, // 10
-  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, // 4x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -294,9 +294,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,
- +1.0/3.0, 0.0, 0.0, // 11
-  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, // 4y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -305,9 +305,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,
-  0.0,+1.0/3.0, 0.0, // 11
-  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, // 4z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -316,9 +316,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,
-  0.0, 0.0,+1.0/3.0, // 11
-  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, // 5x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -327,9 +327,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,
-  0.0, 0.0, 0.0,
- +1.0/3.0, 0.0, 0.0, // 12
+ +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, // 5y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -338,9 +338,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,
-  0.0, 0.0, 0.0,
-  0.0,+1.0/3.0, 0.0, // 12
+  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, // 5z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -349,9 +349,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,
-  0.0, 0.0, 0.0,
-  0.0, 0.0,+1.0/3.0, // 12
+  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, // 6x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -393,9 +393,9 @@
   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, // 10
-  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, // 7y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -404,9 +404,9 @@
   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, // 10
-  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, // 7z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -415,9 +415,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/3.0, // 10
-  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, // 8x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -426,9 +426,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,
- -1.0/3.0, 0.0, 0.0, // 11
-  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, // 8y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -437,9 +437,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,
-  0.0,-1.0/3.0, 0.0, // 11
-  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, // 8z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -448,9 +448,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,
-  0.0, 0.0,-1.0/3.0, // 11
-  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, // 9x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -459,9 +459,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,
-  0.0, 0.0, 0.0,
- -1.0/3.0, 0.0, 0.0, // 12
+ -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, // 9y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -470,9 +470,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,
-  0.0, 0.0, 0.0,
-  0.0,-1.0/3.0, 0.0, // 12
+  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, // 9z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -481,122 +481,122 @@
   0.0, 0.0, 0.0,
   0.0, 0.0, 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, // 12
+  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, // 10x
- +1.0/3.0, 0.0, 0.0, // 3
+ +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,
   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/3.0, 0.0, // 3
+  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, 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/3.0, // 3
+  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, 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/3.0, 0.0, 0.0, // 4
+ -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,
   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/3.0, 0.0, // 4
+  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, 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/3.0, // 4
+  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, 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[] = {
   3.1, 4.1, 5.1,
-  4.05776047684, 8.66419800363, 4.4475545411, // 3
-  4.41507389392, 9.10077378393, 4.76947881673, // 4
-  2.65876279931, 6.30839037001, 3.47575435673, // 5
+  3.41006379025, 5.29327298048, 5.01572764272, // 3
+  3.57151643203, 5.46897093178, 5.17082002405, // 4
+  3.21983818039, 4.8638448816, 4.93230140615, // 5
   3.5, 4.5, 5.5,
-  2.95667964237, 1.25185149728, 6.16433409418, // 7
-  2.99466885439, 1.33277448939, 6.18315418792, // 8
-  4.03773208096, 1.32825348198, 7.69394390058, // 9
-  3.08793771664, 16.0711128131, -2.70880365204, // 10
-  4.34878818628, 18.7230177573, -2.06903261476, // 11
-  -3.11319624288, 8.01523955405, -8.08183170175, // 12
+  3.44245215731, 3.78004526464, 5.73820426796, // 7
+  3.60390479909, 3.95574321594, 5.8932966493, // 8
+  3.25222654746, 3.35061716576, 5.65477803139, // 9
+  0.756229644891, 3.93578272973, -0.663380486215, // 10
+  1.0589140849, 4.55898663393, -0.503801906196, // 11
+  -0.756679642367, 1.94814850272, -1.96433409418, // 12
 };
 
 pylith::faults::CohesiveKinDataTet4::CohesiveKinDataTet4(void)

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-26 21:38:49 UTC (rev 18974)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTet4e.cc	2011-09-27 17:54:55 UTC (rev 18975)
@@ -350,10 +350,10 @@
   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, // 14
-  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, // 5y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -364,10 +364,10 @@
   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, // 14
-  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, // 5z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -378,10 +378,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/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,+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, // 6x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -392,10 +392,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,
@@ -406,10 +406,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,
@@ -420,10 +420,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,
@@ -434,10 +434,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,
-  +2.0/3.0, 0.0, 0.0, // 16
-  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, // 7y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -448,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, 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,+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, // 7z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -462,10 +462,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,+2.0/3.0, // 16
-  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, // 8x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -476,10 +476,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,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
- +1.0/3.0, 0.0, 0.0, // 17
+ +1.0/9.0, 0.0, 0.0, // 16
+ +1.0/9.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,
@@ -490,10 +490,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, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0,+1.0/3.0, 0.0, // 17
+  0.0,+1.0/9.0, 0.0, // 16
+  0.0,+1.0/9.0, 0.0, // 17
   0.0, 0.0, 0.0, // 8z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -504,10 +504,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, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0,+1.0/3.0, // 17
+  0.0, 0.0,+1.0/9.0, // 16
+  0.0, 0.0,+1.0/9.0, // 17
   0.0, 0.0, 0.0, // 9x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -560,10 +560,10 @@
   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, // 14
-  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, // 10y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -574,10 +574,10 @@
   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, // 14
-  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, // 10z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -588,10 +588,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/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,-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, // 11x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -602,10 +602,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,
@@ -616,10 +616,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,
@@ -630,10 +630,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,
@@ -644,10 +644,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,
- -2.0/3.0, 0.0, 0.0, // 16
-  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, // 12y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -658,10 +658,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,-2.0/3.0, 0.0, // 16
-  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, // 12z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -672,10 +672,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,-2.0/3.0, // 16
-  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, // 13x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -686,10 +686,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,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
- -1.0/3.0, 0.0, 0.0, // 17
+ -1.0/9.0, 0.0, 0.0, // 16
+ -1.0/9.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,
@@ -700,10 +700,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, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0,-1.0/3.0, 0.0, // 17
+  0.0,-1.0/9.0, 0.0, // 16
+  0.0,-1.0/9.0, 0.0, // 17
   0.0, 0.0, 0.0, // 13z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -714,178 +714,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, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0,-1.0/3.0, // 17
+  0.0, 0.0,-1.0/9.0, // 16
+  0.0, 0.0,-1.0/9.0, // 17
   0.0, 0.0, 0.0, // 14x
- +2.0/3.0, 0.0, 0.0, // 5
+ +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,
   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/3.0, 0.0, // 5
+  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, 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/3.0, // 5
+  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, 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,// 17x
+ +1.0/9.0, 0.0, 0.0, // 5
   0.0, 0.0, 0.0,
-  0.0, 0.0,-2.0/3.0, // 12
+ +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
   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.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,
- -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/3.0, 0.0, // 8
+  0.0,-1.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, 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/3.0, // 8
+  0.0, 0.0,-1.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, 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/CohesiveKinDataTet4f.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTet4f.cc	2011-09-26 21:38:49 UTC (rev 18974)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTet4f.cc	2011-09-27 17:54:55 UTC (rev 18975)
@@ -261,9 +261,9 @@
   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, // 10
-  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, // 3y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -272,9 +272,9 @@
   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, // 10
-  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, // 3z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -283,9 +283,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/3.0, // 10
-  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, // 4x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -294,9 +294,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,
- +1.0/3.0, 0.0, 0.0, // 11
-  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, // 4y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -305,9 +305,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,
-  0.0,+1.0/3.0, 0.0, // 11
-  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, // 4z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -316,9 +316,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,
-  0.0, 0.0,+1.0/3.0, // 11
-  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, // 5x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -327,9 +327,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,
-  0.0, 0.0, 0.0,
- +1.0/3.0, 0.0, 0.0, // 12
+ +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, // 5y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -338,9 +338,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,
-  0.0, 0.0, 0.0,
-  0.0,+1.0/3.0, 0.0, // 12
+  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, // 5z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -349,9 +349,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,
-  0.0, 0.0, 0.0,
-  0.0, 0.0,+1.0/3.0, // 12
+  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, // 6x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -393,9 +393,9 @@
   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, // 10
-  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, // 7y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -404,9 +404,9 @@
   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, // 10
-  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, // 7z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -415,9 +415,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/3.0, // 10
-  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, // 8x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -426,9 +426,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,
- -1.0/3.0, 0.0, 0.0, // 11
-  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, // 8y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -437,9 +437,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,
-  0.0,-1.0/3.0, 0.0, // 11
-  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, // 8z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -448,9 +448,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,
-  0.0, 0.0,-1.0/3.0, // 11
-  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, // 9x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -459,9 +459,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,
-  0.0, 0.0, 0.0,
- -1.0/3.0, 0.0, 0.0, // 12
+ -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, // 9y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -470,9 +470,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,
-  0.0, 0.0, 0.0,
-  0.0,-1.0/3.0, 0.0, // 12
+  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, // 9z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -481,108 +481,108 @@
   0.0, 0.0, 0.0,
   0.0, 0.0, 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, // 12
+  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, // 10x
- +1.0/3.0, 0.0, 0.0, // 3
+ +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,
   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/3.0, 0.0, // 3
+  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, 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/3.0, // 3
+  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, 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/3.0, 0.0, 0.0, // 4
+ -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,
   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/3.0, 0.0, // 4
+  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, 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/3.0, // 4
+  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, 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,
 };
 
 pylith::faults::CohesiveKinDataTet4f::CohesiveKinDataTet4f(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-09-26 21:38:49 UTC (rev 18974)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTri3.cc	2011-09-27 17:54:55 UTC (rev 18975)
@@ -213,32 +213,32 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
- +1.0, 0.0, // 8
-  0.0, 0.0,
+ +0.5, 0.0, // 8
+ +0.5, 0.0, // 9
   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,+1.0, //  8
-  0.0, 0.0,
+  0.0,+0.5, // 8
+  0.0,+0.5, // 9
   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.0, 0.0,
- +1.0, 0.0, //  9
+ +0.5, 0.0, // 8
+ +0.5, 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.0,
-  0.0,+1.0, //  9
+  0.0,+0.5, // 8
+  0.0,+0.5, // 9
   0.0, 0.0, // 5x
   0.0, 0.0,
   0.0, 0.0,
@@ -261,77 +261,77 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
- -1.0, 0.0, //  8
-  0.0, 0.0,
+ -0.5, 0.0, // 8
+ -0.5, 0.0, // 9
   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,-1.0, //  8
-  0.0, 0.0,
+  0.0,-0.5, // 8
+  0.0,-0.5, // 9
   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.0, 0.0,
- -1.0, 0.0, //  9
+ -0.5, 0.0, // 8
+ -0.5, 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.0,
-  0.0,-1.0, //  9
+  0.0,-0.5, // 8
+  0.0,-0.5, // 9
 
   0.0, 0.0, // 8x
- +1.0, 0.0, //  3
+ +0.5, 0.0, // 3
+ +0.5, 0.0, // 4
   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,+1.0, //  3
+  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, 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,
- +1.0, 0.0, //  4
+ -0.5, 0.0, // 6
+ -0.5, 0.0, // 7
   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,+1.0, //  4
+  0.0,-0.5, // 6
+  0.0,-0.5, // 7
   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[] = {
   3.1, 4.1,
-  2.82834517219, -0.0422000510363, // 3
-  3.14958479199, 0.17979481945, // 4
+  3.10823337585, 3.15254319728, // 3
+  3.26306471003, 3.28826072789, // 4
   3.4, 4.4,
-  3.79732386225, 7.89376004083, // 6
-  3.81502339436, 7.85074513807, // 7
-  -0.445985793369, -5.09064006124, // 8
-  -0.19553977041, -5.35626673471, // 9
+  3.57341329932, 5.33796544218, // 6
+  3.7282446335, 5.47368297279, // 7
+  -0.11011994898, -1.25694816327, // 8
+  -0.0480158769562, -1.31526105375, // 9
 };
 
 pylith::faults::CohesiveKinDataTri3::CohesiveKinDataTri3(void)

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-09-26 21:38:49 UTC (rev 18974)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinDataTri3d.cc	2011-09-27 17:54:55 UTC (rev 18975)
@@ -280,9 +280,9 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
- +2.0, 0.0, // 13
-  0.0, 0.0,
-  0.0, 0.0,
+ +1.0, 0.0, // 13
+ +0.5, 0.0, // 14
+ +0.5, 0.0, // 15
   0.0, 0.0, // 5y
   0.0, 0.0,
   0.0, 0.0,
@@ -292,9 +292,9 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
-  0.0,+2.0, // 13
-  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, 0.0, // 6x
   0.0, 0.0,
   0.0, 0.0,
@@ -304,9 +304,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 +316,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 +352,9 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
+ +0.5, 0.0, // 13
   0.0, 0.0,
-  0.0, 0.0,
- +1.0, 0.0, // 15
+ +0.5, 0.0, // 15
   0.0, 0.0, // 8y
   0.0, 0.0,
   0.0, 0.0,
@@ -364,9 +364,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.0,
-  0.0,+1.0, // 15
+  0.0,+0.5, // 15
   0.0, 0.0, // 9x
   0.0, 0.0,
   0.0, 0.0,
@@ -400,9 +400,9 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
- -2.0, 0.0, // 13
-  0.0, 0.0,
-  0.0, 0.0,
+ -1.0, 0.0, // 13
+ -0.5, 0.0, // 14
+ -0.5, 0.0, // 15
   0.0, 0.0, // 10y
   0.0, 0.0,
   0.0, 0.0,
@@ -412,9 +412,9 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
-  0.0,-2.0, // 13
-  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, 0.0, // 11x
   0.0, 0.0,
   0.0, 0.0,
@@ -424,9 +424,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 +436,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 +448,9 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
+ -0.5, 0.0, // 13
   0.0, 0.0,
-  0.0, 0.0,
- -1.0, 0.0, // 15
+ -0.5, 0.0, // 15
   0.0, 0.0, // 12y
   0.0, 0.0,
   0.0, 0.0,
@@ -460,81 +460,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.0,
-  0.0,-1.0, // 15
+  0.0,-0.5, // 15
   0.0, 0.0, // 13x
- +2.0, 0.0, // 5
+ +1.0, 0.0, // 5
+ +0.5, 0.0, // 6
   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,+2.0, // 5
+  0.0,+1.0, // 5
+  0.0,+0.5, // 6
   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,
- +1.0, 0.0, // 8
+ -0.5, 0.0, // 10
   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,+1.0, // 8
+  0.0,-0.5, // 10
   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)

Modified: short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinSrcsDataHex8.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinSrcsDataHex8.cc	2011-09-26 21:38:49 UTC (rev 18974)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinSrcsDataHex8.cc	2011-09-27 17:54:55 UTC (rev 18975)
@@ -876,10 +876,10 @@
   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, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+ +4.0/9.0, 0.0, 0.0, // 18
+ +2.0/9.0, 0.0, 0.0, // 19
+ +2.0/9.0, 0.0, 0.0, // 20
+ +1.0/9.0, 0.0, 0.0, // 21
   0.0, 0.0, 0.0, // 6y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -896,10 +896,10 @@
   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, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+  0.0,+4.0/9.0, 0.0, // 18
+  0.0,+2.0/9.0, 0.0, // 19
+  0.0,+2.0/9.0, 0.0, // 20
+  0.0,+1.0/9.0, 0.0, // 21
   0.0, 0.0, 0.0, // 6z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -916,10 +916,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, // 18
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+  0.0, 0.0,+4.0/9.0, // 18
+  0.0, 0.0,+2.0/9.0, // 19
+  0.0, 0.0,+2.0/9.0, // 20
+  0.0, 0.0,+1.0/9.0, // 21
   0.0, 0.0, 0.0, // 7x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -936,10 +936,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,
- +1.0, 0.0, 0.0, // 19
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+ +2.0/9.0, 0.0, 0.0, // 18
+ +4.0/9.0, 0.0, 0.0, // 19
+ +1.0/9.0, 0.0, 0.0, // 20
+ +2.0/9.0, 0.0, 0.0, // 21
   0.0, 0.0, 0.0, // 7y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -956,10 +956,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,+1.0, 0.0, // 19
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+  0.0,+2.0/9.0, 0.0, // 18
+  0.0,+4.0/9.0, 0.0, // 19
+  0.0,+1.0/9.0, 0.0, // 20
+  0.0,+2.0/9.0, 0.0, // 21
   0.0, 0.0, 0.0, // 7z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -976,10 +976,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,+1.0, // 19
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+  0.0, 0.0,+2.0/9.0, // 18
+  0.0, 0.0,+4.0/9.0, // 19
+  0.0, 0.0,+1.0/9.0, // 20
+  0.0, 0.0,+2.0/9.0, // 21
   0.0, 0.0, 0.0, // 8x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -996,10 +996,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,
- +1.0, 0.0, 0.0, // 20
-  0.0, 0.0, 0.0,
+ +2.0/9.0, 0.0, 0.0, // 18
+ +1.0/9.0, 0.0, 0.0, // 19
+ +4.0/9.0, 0.0, 0.0, // 20
+ +2.0/9.0, 0.0, 0.0, // 21
   0.0, 0.0, 0.0, // 8y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1016,10 +1016,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,+1.0, 0.0, // 20
-  0.0, 0.0, 0.0,
+  0.0,+2.0/9.0, 0.0, // 18
+  0.0,+1.0/9.0, 0.0, // 19
+  0.0,+4.0/9.0, 0.0, // 20
+  0.0,+2.0/9.0, 0.0, // 21
   0.0, 0.0, 0.0, // 8z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1036,10 +1036,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,+1.0, // 20
-  0.0, 0.0, 0.0,
+  0.0, 0.0,+2.0/9.0, // 18
+  0.0, 0.0,+1.0/9.0, // 19
+  0.0, 0.0,+4.0/9.0, // 20
+  0.0, 0.0,+2.0/9.0, // 21
   0.0, 0.0, 0.0, // 9x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1056,10 +1056,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,
- +1.0, 0.0, 0.0, // 21
+ +1.0/9.0, 0.0, 0.0, // 18
+ +2.0/9.0, 0.0, 0.0, // 19
+ +2.0/9.0, 0.0, 0.0, // 20
+ +4.0/9.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,
@@ -1076,10 +1076,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,+1.0, 0.0, // 21
+  0.0,+1.0/9.0, 0.0, // 18
+  0.0,+2.0/9.0, 0.0, // 19
+  0.0,+2.0/9.0, 0.0, // 20
+  0.0,+4.0/9.0, 0.0, // 21
   0.0, 0.0, 0.0, // 9z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1096,10 +1096,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,+1.0, // 21
+  0.0, 0.0,+1.0/9.0, // 18
+  0.0, 0.0,+2.0/9.0, // 19
+  0.0, 0.0,+2.0/9.0, // 20
+  0.0, 0.0,+4.0/9.0, // 21
   0.0, 0.0, 0.0, // 10x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1356,10 +1356,10 @@
   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, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+ -4.0/9.0, 0.0, 0.0, // 18
+ -2.0/9.0, 0.0, 0.0, // 19
+ -2.0/9.0, 0.0, 0.0, // 20
+ -1.0/9.0, 0.0, 0.0, // 21
   0.0, 0.0, 0.0, // 14y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1376,10 +1376,10 @@
   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, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+  0.0,-4.0/9.0, 0.0, // 18
+  0.0,-2.0/9.0, 0.0, // 19
+  0.0,-2.0/9.0, 0.0, // 20
+  0.0,-1.0/9.0, 0.0, // 21
   0.0, 0.0, 0.0, // 14z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1396,10 +1396,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, // 18
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+  0.0, 0.0,-4.0/9.0, // 18
+  0.0, 0.0,-2.0/9.0, // 19
+  0.0, 0.0,-2.0/9.0, // 20
+  0.0, 0.0,-1.0/9.0, // 21
   0.0, 0.0, 0.0, // 15x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1416,10 +1416,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,
- -1.0, 0.0, 0.0, // 19
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+ -2.0/9.0, 0.0, 0.0, // 18
+ -4.0/9.0, 0.0, 0.0, // 19
+ -1.0/9.0, 0.0, 0.0, // 20
+ -2.0/9.0, 0.0, 0.0, // 21
   0.0, 0.0, 0.0, // 15y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1436,10 +1436,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,-1.0, 0.0, // 19
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+  0.0,-2.0/9.0, 0.0, // 18
+  0.0,-4.0/9.0, 0.0, // 19
+  0.0,-1.0/9.0, 0.0, // 20
+  0.0,-2.0/9.0, 0.0, // 21
   0.0, 0.0, 0.0, // 15z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1456,10 +1456,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,-1.0, // 19
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
+  0.0, 0.0,-2.0/9.0, // 18
+  0.0, 0.0,-4.0/9.0, // 19
+  0.0, 0.0,-1.0/9.0, // 20
+  0.0, 0.0,-2.0/9.0, // 21
   0.0, 0.0, 0.0, // 16x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1476,10 +1476,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,
- -1.0, 0.0, 0.0, // 20
-  0.0, 0.0, 0.0,
+ -2.0/9.0, 0.0, 0.0, // 18
+ -1.0/9.0, 0.0, 0.0, // 19
+ -4.0/9.0, 0.0, 0.0, // 20
+ -2.0/9.0, 0.0, 0.0, // 21
   0.0, 0.0, 0.0, // 16y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1496,10 +1496,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,-1.0, 0.0, // 20
-  0.0, 0.0, 0.0,
+  0.0,-2.0/9.0, 0.0, // 18
+  0.0,-1.0/9.0, 0.0, // 19
+  0.0,-4.0/9.0, 0.0, // 20
+  0.0,-2.0/9.0, 0.0, // 21
   0.0, 0.0, 0.0, // 16z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1516,10 +1516,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,-1.0, // 20
-  0.0, 0.0, 0.0,
+  0.0, 0.0,-2.0/9.0, // 18
+  0.0, 0.0,-1.0/9.0, // 19
+  0.0, 0.0,-4.0/9.0, // 20
+  0.0, 0.0,-2.0/9.0, // 21
   0.0, 0.0, 0.0, // 17x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1536,10 +1536,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,
- -1.0, 0.0, 0.0, // 21
+ -1.0/9.0, 0.0, 0.0, // 18
+ -2.0/9.0, 0.0, 0.0, // 19
+ -2.0/9.0, 0.0, 0.0, // 20
+ -4.0/9.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,
@@ -1556,10 +1556,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,-1.0, 0.0, // 21
+  0.0,-1.0/9.0, 0.0, // 18
+  0.0,-2.0/9.0, 0.0, // 19
+  0.0,-2.0/9.0, 0.0, // 20
+  0.0,-4.0/9.0, 0.0, // 21
   0.0, 0.0, 0.0, // 17z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -1576,250 +1576,250 @@
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
   0.0, 0.0, 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, // 21
+  0.0, 0.0,-1.0/9.0, // 18
+  0.0, 0.0,-2.0/9.0, // 19
+  0.0, 0.0,-2.0/9.0, // 20
+  0.0, 0.0,-4.0/9.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,
- +1.0, 0.0, 0.0, // 6
+ +4.0/9.0, 0.0, 0.0, // 6
+ +2.0/9.0, 0.0, 0.0, // 7
+ +2.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,
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
+ -4.0/9.0, 0.0, 0.0, // 14
+ -2.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,
   0.0, 0.0, 0.0,
- -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,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 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.0,
   0.0, 0.0, 0.0,
-  0.0,+1.0, 0.0, // 6
+  0.0,+4.0/9.0, 0.0, // 6
+  0.0,+2.0/9.0, 0.0, // 7
+  0.0,+2.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, 0.0, 0.0,
   0.0, 0.0, 0.0,
+  0.0,-4.0/9.0, 0.0, // 14
+  0.0,-2.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, 0.0, 0.0,
-  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,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
   0.0, 0.0, 0.0, // 18z
   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, 0.0,+4.0/9.0, // 6
+  0.0, 0.0,+2.0/9.0, // 7
+  0.0, 0.0,+2.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, 0.0,
   0.0, 0.0, 0.0,
+  0.0, 0.0,-4.0/9.0, // 14
+  0.0, 0.0,-2.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, 0.0,
-  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,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
   0.0, 0.0, 0.0, // 19x
   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, // 6
+ +4.0/9.0, 0.0, 0.0, // 7
+ +1.0/9.0, 0.0, 0.0, // 8
+ +2.0/9.0, 0.0, 0.0, // 9
   0.0, 0.0, 0.0,
- +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,
+ -2.0/9.0, 0.0, 0.0, // 14
+ -4.0/9.0, 0.0, 0.0, // 15
+ -1.0/9.0, 0.0, 0.0, // 16
+ -2.0/9.0, 0.0, 0.0, // 17
   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, // 15
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
   0.0, 0.0, 0.0, // 19y
   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, // 6
+  0.0,+4.0/9.0, 0.0, // 7
+  0.0,+1.0/9.0, 0.0, // 8
+  0.0,+2.0/9.0, 0.0, // 9
   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,
   0.0, 0.0, 0.0,
+  0.0,-2.0/9.0, 0.0, // 14
+  0.0,-4.0/9.0, 0.0, // 15
+  0.0,-1.0/9.0, 0.0, // 16
+  0.0,-2.0/9.0, 0.0, // 17
   0.0, 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, // 15
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
   0.0, 0.0, 0.0, // 19z
   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, // 6
+  0.0, 0.0,+4.0/9.0, // 7
+  0.0, 0.0,+1.0/9.0, // 8
+  0.0, 0.0,+2.0/9.0, // 9
   0.0, 0.0, 0.0,
-  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,
+  0.0, 0.0,-2.0/9.0, // 14
+  0.0, 0.0,-4.0/9.0, // 15
+  0.0, 0.0,-1.0/9.0, // 16
+  0.0, 0.0,-2.0/9.0, // 17
   0.0, 0.0, 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, // 15
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
   0.0, 0.0, 0.0, // 20x
   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, // 6
+ +1.0/9.0, 0.0, 0.0, // 7
+ +4.0/9.0, 0.0, 0.0, // 8
+ +2.0/9.0, 0.0, 0.0, // 9
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
- +1.0, 0.0, 0.0, // 8
   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
+ -4.0/9.0, 0.0, 0.0, // 16
+ -2.0/9.0, 0.0, 0.0, // 17
   0.0, 0.0, 0.0,
   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, // 16
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
   0.0, 0.0, 0.0, // 20y
   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, // 6
+  0.0,+1.0/9.0, 0.0, // 7
+  0.0,+4.0/9.0, 0.0, // 8
+  0.0,+2.0/9.0, 0.0, // 9
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
-  0.0,+1.0, 0.0, // 8
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
+  0.0,-2.0/9.0, 0.0, // 6
+  0.0,-1.0/9.0, 0.0, // 7
+  0.0,-4.0/9.0, 0.0, // 8
+  0.0,-2.0/9.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,
-  0.0, 0.0, 0.0,
-  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,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
   0.0, 0.0, 0.0, // 20z
   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, // 6
+  0.0, 0.0,+1.0/9.0, // 7
+  0.0, 0.0,+4.0/9.0, // 8
+  0.0, 0.0,+2.0/9.0, // 9
   0.0, 0.0, 0.0,
   0.0, 0.0, 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,-2.0/9.0, // 14
+  0.0, 0.0,-1.0/9.0, // 15
+  0.0, 0.0,-4.0/9.0, // 16
+  0.0, 0.0,-2.0/9.0, // 17
   0.0, 0.0, 0.0,
   0.0, 0.0, 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, // 16
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
   0.0, 0.0, 0.0, // 21x
   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, // 6
+ +2.0/9.0, 0.0, 0.0, // 7
+ +2.0/9.0, 0.0, 0.0, // 8
+ +4.0/9.0, 0.0, 0.0, // 9
   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, 0.0,
+ -1.0/9.0, 0.0, 0.0, // 14
+ -2.0/9.0, 0.0, 0.0, // 15
+ -2.0/9.0, 0.0, 0.0, // 16
+ -4.0/9.0, 0.0, 0.0, // 17
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
   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, // 17
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
   0.0, 0.0, 0.0, // 21y
   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, // 6
+  0.0,+2.0/9.0, 0.0, // 7
+  0.0,+2.0/9.0, 0.0, // 8
+  0.0,+4.0/9.0, 0.0, // 9
   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
   0.0, 0.0, 0.0,
+  0.0,-1.0/9.0, 0.0, // 14
+  0.0,-2.0/9.0, 0.0, // 15
+  0.0,-2.0/9.0, 0.0, // 16
+  0.0,-4.0/9.0, 0.0, // 17
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
   0.0, 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, // 17
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
-  0.0, 0.0, 0.0,
   0.0, 0.0, 0.0, // 21z
   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, // 6
+  0.0, 0.0,+2.0/9.0, // 7
+  0.0, 0.0,+2.0/9.0, // 8
+  0.0, 0.0,+4.0/9.0, // 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.0, 0.0, 0.0,
+  0.0, 0.0,-1.0/9.0, // 14
+  0.0, 0.0,-2.0/9.0, // 15
+  0.0, 0.0,-2.0/9.0, // 16
+  0.0, 0.0,-4.0/9.0, // 17
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
   0.0, 0.0, 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, // 17
-  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::CohesiveKinSrcsDataHex8::CohesiveKinSrcsDataHex8(void)

Modified: short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinSrcsDataQuad4.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinSrcsDataQuad4.cc	2011-09-26 21:38:49 UTC (rev 18974)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinSrcsDataQuad4.cc	2011-09-27 17:54:55 UTC (rev 18975)
@@ -253,8 +253,8 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
- +1.0, 0.0, // 10
-  0.0, 0.0,
+ +0.5, 0.0, // 10
+ +0.5, 0.0, // 11
   0.0, 0.0, // 4y
   0.0, 0.0,
   0.0, 0.0,
@@ -263,8 +263,8 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
-  0.0,+1.0, //  10
-  0.0, 0.0,
+  0.0,+0.5, // 10
+  0.0,+0.5, // 11
   0.0, 0.0, // 5x
   0.0, 0.0,
   0.0, 0.0,
@@ -273,8 +273,8 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
-  0.0, 0.0,
- +1.0, 0.0, //  11
+ +0.5, 0.0, // 10
+ +0.5, 0.0, // 11
   0.0, 0.0, // 5y
   0.0, 0.0,
   0.0, 0.0,
@@ -283,8 +283,8 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
-  0.0, 0.0,
-  0.0,+1.0, //  11
+  0.0,+0.5, // 10
+  0.0,+0.5, // 11
   0.0, 0.0, // 6x
   0.0, 0.0,
   0.0, 0.0,
@@ -333,8 +333,8 @@
   0.0, 0.0, 
   0.0, 0.0,
   0.0, 0.0,
- -1.0, 0.0, //  10
-  0.0, 0.0,
+ -0.5, 0.0, // 10
+ -0.5, 0.0, // 11
   0.0, 0.0, // 8y
   0.0, 0.0,
   0.0, 0.0,
@@ -343,8 +343,8 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
-  0.0,-1.0, // 10
-  0.0, 0.0,
+  0.0,-0.5, // 10
+  0.0,-0.5, // 11
   0.0, 0.0, // 9x
   0.0, 0.0,
   0.0, 0.0,
@@ -353,8 +353,8 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
-  0.0, 0.0,
- -1.0, 0.0, // 11
+ -0.5, 0.0, // 10
+ -0.5, 0.0, // 11
   0.0, 0.0, // 9y
   0.0, 0.0,
   0.0, 0.0,
@@ -363,48 +363,48 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
-  0.0, 0.0,
-  0.0,-1.0, // 11
+  0.0,-0.5, // 10
+  0.0,-0.5, // 11
   0.0, 0.0, // 10x
   0.0, 0.0,
- +1.0, 0.0, //  4
+ +0.5, 0.0, // 4
+ +0.5, 0.0, // 5
   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,+1.0, //  4
+  0.0,+0.5, // 4
+  0.0,+0.5, // 5
   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,
 };
 
 pylith::faults::CohesiveKinSrcsDataQuad4::CohesiveKinSrcsDataQuad4(void)

Modified: short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinSrcsDataTet4.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinSrcsDataTet4.cc	2011-09-26 21:38:49 UTC (rev 18974)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinSrcsDataTet4.cc	2011-09-27 17:54:55 UTC (rev 18975)
@@ -260,9 +260,9 @@
   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, // 10
-  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, // 3y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -271,9 +271,9 @@
   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, // 10
-  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, // 3z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -282,9 +282,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/3.0, // 10
-  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, // 4x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -293,9 +293,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,
- +1.0/3.0, 0.0, 0.0, // 11
-  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, // 4y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -304,9 +304,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,
-  0.0,+1.0/3.0, 0.0, // 11
-  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, // 4z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -315,9 +315,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,
-  0.0, 0.0,+1.0/3.0, // 11
-  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, // 5x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -326,9 +326,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,
-  0.0, 0.0, 0.0,
- +1.0/3.0, 0.0, 0.0, // 12
+ +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, // 5y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -337,9 +337,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,
-  0.0, 0.0, 0.0,
-  0.0,+1.0/3.0, 0.0, // 12
+  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, // 5z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -348,9 +348,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,
-  0.0, 0.0, 0.0,
-  0.0, 0.0,+1.0/3.0, // 12
+  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, // 6x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -392,9 +392,9 @@
   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, // 10
-  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, // 7y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -403,9 +403,9 @@
   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, // 10
-  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, // 7z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -414,9 +414,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/3.0, // 10
-  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, // 8x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -425,9 +425,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,
- -1.0/3.0, 0.0, 0.0, // 11
-  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, // 8y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -436,9 +436,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,
-  0.0,-1.0/3.0, 0.0, // 11
-  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, // 8z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -447,9 +447,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,
-  0.0, 0.0,-1.0/3.0, // 11
-  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, // 9x
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -458,9 +458,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,
-  0.0, 0.0, 0.0,
- -1.0/3.0, 0.0, 0.0, // 12
+ -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, // 9y
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -469,9 +469,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,
-  0.0, 0.0, 0.0,
-  0.0,-1.0/3.0, 0.0, // 12
+  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, // 9z
   0.0, 0.0, 0.0,
   0.0, 0.0, 0.0,
@@ -480,108 +480,108 @@
   0.0, 0.0, 0.0,
   0.0, 0.0, 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, // 12
+  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, // 10x
- +1.0/3.0, 0.0, 0.0, // 3
+ +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,
   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/3.0, 0.0, // 3
+  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, 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/3.0, // 3
+  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, 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/3.0, 0.0, 0.0, // 4
+ -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,
   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/3.0, 0.0, // 4
+  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, 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/3.0, // 4
+  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, 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,
 };
 
 pylith::faults::CohesiveKinSrcsDataTet4::CohesiveKinSrcsDataTet4(void)

Modified: short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinSrcsDataTri3.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinSrcsDataTri3.cc	2011-09-26 21:38:49 UTC (rev 18974)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/CohesiveKinSrcsDataTri3.cc	2011-09-27 17:54:55 UTC (rev 18975)
@@ -232,32 +232,32 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
- +1.0, 0.0, // 8
-  0.0, 0.0,
+ +0.5, 0.0, // 8
+ +0.5, 0.0, // 9
   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,+1.0, //  8
-  0.0, 0.0,
+  0.0,+0.5, // 8
+  0.0,+0.5, // 9
   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.0, 0.0,
- +1.0, 0.0, //  9
+ +0.5, 0.0, // 8
+ +0.5, 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.0,
-  0.0,+1.0, //  9
+  0.0,+0.5, // 8
+  0.0,+0.5, // 9
   0.0, 0.0, // 5x
   0.0, 0.0,
   0.0, 0.0,
@@ -280,66 +280,66 @@
   0.0, 0.0,
   0.0, 0.0,
   0.0, 0.0,
- -1.0, 0.0, //  8
-  0.0, 0.0,
+ -0.5, 0.0, // 8
+ -0.5, 0.0, // 9
   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,-1.0, //  8
-  0.0, 0.0,
+  0.0,-0.5, // 8
+  0.0,-0.5, // 9
   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.0, 0.0,
- -1.0, 0.0, //  9
+ -0.5, 0.0, // 8
+ -0.5, 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.0,
-  0.0,-1.0, //  9
+  0.0,-0.5, // 8
+  0.0,-0.5, // 9
 
   0.0, 0.0, // 8x
- +1.0, 0.0, //  3
+ +0.5, 0.0, // 3
+ +0.5, 0.0, // 4
   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,+1.0, //  3
+  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, 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,
- +1.0, 0.0, //  4
+ -0.5, 0.0, // 6
+ -0.5, 0.0, // 7
   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,+1.0, //  4
+  0.0,-0.5, // 6
+  0.0,-0.5, // 7
   0.0, 0.0,
   0.0, 0.0,
-  0.0,-1.0, //  7
-  0.0, 0.0,
-  0.0, 0.0,
 };
 
 pylith::faults::CohesiveKinSrcsDataTri3::CohesiveKinSrcsDataTri3(void)

Modified: short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/adjustsoln.py
===================================================================
--- short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/adjustsoln.py	2011-09-26 21:38:49 UTC (rev 18974)
+++ short/3D/PyLith/branches/v1.6-revisedfault/unittests/libtests/faults/data/adjustsoln.py	2011-09-27 17:54:55 UTC (rev 18975)
@@ -1,5 +1,5 @@
 test = "hex8"
-vertex = 3
+vertex = 
 
 # ----------------------------------------------------------------------
 if test == "line2":
@@ -7,7 +7,7 @@
     jN = 2.2; duN = 1.2
     jP = 2.4; duP = 1.4
 
-    Sinv = (jN + jP) /jL**2
+    Sinv = 1.0/(jL**2 * (1.0/jN + 1.0/jP))
     duL = Sinv * (-rL + jL*(duP-duN));
 
     dduN = jL / jN * duL
@@ -59,7 +59,7 @@
             jN = 1.4; duNx = 3.4; duNy = 4.4;
             jP = 1.9; duPx = 3.9; duPy = 4.9;
 
-    Sinv = (jN + jP) / (jL*jL)
+    Sinv = 1.0/(jL**2 * (1.0/jN + 1.0/jP))
     duLx = Sinv * (-rLx + jL*(duPx - duNx))
     duLy = Sinv * (-rLy + jL*(duPy - duNy))
 
@@ -188,7 +188,7 @@
             jN = 1.8; duNx = 3.8; duNy = 4.8; duNz = 5.8;
             jP = 1.9; duPx = 3.9; duPy = 4.9; duPz = 5.9;
             
-    Sinv = (jN + jP) / (jL*jL)
+    Sinv = 1.0/(jL**2 * (1.0/jN + 1.0/jP))
     duLx = Sinv * (-rLx + jL*(duPx - duNx))
     duLy = Sinv * (-rLy + jL*(duPy - duNy))
     duLz = Sinv * (-rLz + jL*(duPz - duNz))



More information about the CIG-COMMITS mailing list