[cig-commits] r16215 - in short/3D/PyLith/trunk: libsrc/faults unittests/libtests/faults unittests/libtests/faults/data
brad at geodynamics.org
brad at geodynamics.org
Tue Feb 2 15:40:57 PST 2010
Author: brad
Date: 2010-02-02 15:40:56 -0800 (Tue, 02 Feb 2010)
New Revision: 16215
Modified:
short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveKin.cc
short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveKin.hh
short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.cc
short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.hh
short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinHex8.hh
short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinQuad4.hh
short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinTet4.hh
short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinTri3.hh
short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveKinDataHex8.cc
short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveKinDataQuad4.cc
short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveKinDataTet4.cc
short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveKinDataTri3.cc
Log:
Updated FaultCohesiveKin unit tests for integrating residual.
Modified: short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveKin.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveKin.cc 2010-02-02 23:10:22 UTC (rev 16214)
+++ short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveKin.cc 2010-02-02 23:40:56 UTC (rev 16215)
@@ -169,15 +169,6 @@
} // splitFields
// ----------------------------------------------------------------------
-// Integrate contribution of cohesive cells to residual term that
-// require assembly across processors.
-void pylith::faults::FaultCohesiveKin::integrateResidual(const topology::Field<
- topology::Mesh>& residual,
- const double t,
- topology::SolutionFields* const fields) { // integrateResidual
-} // integrateResidual
-
-// ----------------------------------------------------------------------
// Integrate contribution of cohesive cells to residual term that do
// not require assembly across cells, vertices, or processors.
void pylith::faults::FaultCohesiveKin::integrateResidualAssembled(const topology::Field<
Modified: short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveKin.hh
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveKin.hh 2010-02-02 23:10:22 UTC (rev 16214)
+++ short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveKin.hh 2010-02-02 23:40:56 UTC (rev 16215)
@@ -139,17 +139,6 @@
void splitField(topology::Field<topology::Mesh>* field);
/** Integrate contributions to residual term (r) for operator that
- * require assembly across processors.
- *
- * @param residual Field containing values for residual
- * @param t Current time
- * @param fields Solution fields
- */
- void integrateResidual(const topology::Field<topology::Mesh>& residual,
- const double t,
- topology::SolutionFields* const fields);
-
- /** Integrate contributions to residual term (r) for operator that
* do not require assembly across cells, vertices, or processors.
*
* @param residual Field containing values for residual
Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.cc 2010-02-02 23:10:22 UTC (rev 16214)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.cc 2010-02-02 23:40:56 UTC (rev 16215)
@@ -211,7 +211,6 @@
void
pylith::faults::TestFaultCohesiveKin::testIntegrateResidual(void)
{ // testIntegrateResidual
-#if 0
topology::Mesh mesh;
FaultCohesiveKin fault;
topology::SolutionFields fields(mesh);
@@ -249,7 +248,6 @@
//residual.view("RESIDUAL"); // DEBUGGING
// Check values
- const double* valsE = _data->residual;
iVertex = 0;
const int fiberDimE = spaceDim;
const double tolerance = 1.0e-06;
@@ -263,11 +261,8 @@
for (int i=0; i < fiberDimE; ++i) {
const int index = iVertex*spaceDim+i;
- const double valE = valsE[index];
- if (fabs(valE) > tolerance)
- CPPUNIT_ASSERT_DOUBLES_EQUAL(1.0, vals[i]/valE, tolerance);
- else
- CPPUNIT_ASSERT_DOUBLES_EQUAL(valE, vals[i], tolerance);
+ const double valE = 0.0;
+ CPPUNIT_ASSERT_DOUBLES_EQUAL(valE, vals[i], tolerance);
} // for
} // for
} // Integrate residual with disp (as opposed to disp increment).
@@ -280,7 +275,6 @@
//residual.view("RESIDUAL"); // DEBUGGING
// Check values
- const double* valsE = _data->residualIncr;
iVertex = 0;
const int fiberDimE = spaceDim;
const double tolerance = 1.0e-06;
@@ -294,98 +288,14 @@
for (int i=0; i < fiberDimE; ++i) {
const int index = iVertex*spaceDim+i;
- const double valE = valsE[index];
- if (fabs(valE) > tolerance)
- CPPUNIT_ASSERT_DOUBLES_EQUAL(1.0, vals[i]/valE, tolerance);
- else
- CPPUNIT_ASSERT_DOUBLES_EQUAL(valE, vals[i], tolerance);
+ const double valE = 0.0;
+ CPPUNIT_ASSERT_DOUBLES_EQUAL(valE, vals[i], tolerance);
} // for
} // for
} // Integrate residual with disp increment.
-#endif
} // testIntegrateResidual
// ----------------------------------------------------------------------
-// Test integrateJacobian().
-void
-pylith::faults::TestFaultCohesiveKin::testIntegrateJacobian(void)
-{ // testIntegrateJacobian
- topology::Mesh mesh;
- FaultCohesiveKin fault;
- topology::SolutionFields fields(mesh);
- _initialize(&mesh, &fault, &fields);
-
- const ALE::Obj<RealSection>& dispSection = fields.get("disp(t)").section();
- CPPUNIT_ASSERT(!dispSection.isNull());
-
- const int spaceDim = _data->spaceDim;
- const ALE::Obj<SieveMesh>& sieveMesh = mesh.sieveMesh();
- CPPUNIT_ASSERT(!sieveMesh.isNull());
- const ALE::Obj<SieveMesh::label_sequence>& vertices = sieveMesh->depthStratum(0);
- CPPUNIT_ASSERT(!vertices.isNull());
- const SieveMesh::label_sequence::iterator verticesBegin = vertices->begin();
- const SieveMesh::label_sequence::iterator verticesEnd = vertices->end();
- int iVertex = 0;
- for (SieveMesh::label_sequence::iterator v_iter=verticesBegin;
- v_iter != verticesEnd;
- ++v_iter, ++iVertex) {
- dispSection->updatePoint(*v_iter, &_data->fieldT[iVertex*spaceDim]);
- } // for
-
- topology::Jacobian jacobian(fields);
-
- const double t = 2.134;
- fault.integrateJacobian(&jacobian, t, &fields);
- CPPUNIT_ASSERT_EQUAL(false, fault.needNewJacobian());
-
- jacobian.assemble("final_assembly");
-
- //MatView(jacobian, PETSC_VIEWER_STDOUT_WORLD); // DEBUGGING
-
- const double* valsE = _data->jacobian;
- const int nrowsE = dispSection->sizeWithBC();
- const int ncolsE = nrowsE;
-
- int nrows = 0;
- int ncols = 0;
- PetscMat jacobianMat = jacobian.matrix();
- MatGetSize(jacobianMat, &nrows, &ncols);
- CPPUNIT_ASSERT_EQUAL(nrowsE, nrows);
- CPPUNIT_ASSERT_EQUAL(ncolsE, ncols);
-
- PetscMat jDense;
- PetscMat jSparseAIJ;
- MatConvert(jacobianMat, MATSEQAIJ, MAT_INITIAL_MATRIX, &jSparseAIJ);
- MatConvert(jSparseAIJ, MATSEQDENSE, MAT_INITIAL_MATRIX, &jDense);
-
- double_array vals(nrows*ncols);
- int_array rows(nrows);
- int_array cols(ncols);
- for (int iRow=0; iRow < nrows; ++iRow)
- rows[iRow] = iRow;
- for (int iCol=0; iCol < ncols; ++iCol)
- cols[iCol] = iCol;
- MatGetValues(jDense, nrows, &rows[0], ncols, &cols[0], &vals[0]);
- const double tolerance = 1.0e-06;
- for (int iRow=0; iRow < nrows; ++iRow)
- for (int iCol=0; iCol < ncols; ++iCol) {
- const int index = ncols*iRow+iCol;
- const double valE = 0.0;
-#if 0 // DEBUGGING
- if (fabs(valE-vals[index]) > tolerance)
- std::cout << "ERROR: iRow: " << iRow << ", iCol: " << iCol
- << "valE: " << valE
- << ", val: " << vals[index]
- << std::endl;
-#endif // DEBUGGING
- CPPUNIT_ASSERT_DOUBLES_EQUAL(valE, vals[index], tolerance);
- } // for
- MatDestroy(jDense);
- MatDestroy(jSparseAIJ);
- CPPUNIT_ASSERT_EQUAL(false, fault.needNewJacobian());
-} // testIntegrateJacobian
-
-// ----------------------------------------------------------------------
// Test integrateResidualAssembled().
void
pylith::faults::TestFaultCohesiveKin::testIntegrateResidualAssembled(void)
@@ -422,7 +332,7 @@
fault.useSolnIncr(false);
fault.integrateResidualAssembled(residual, t, &fields);
- //residual.view("RESIDUAL"); // DEBUGGING
+ residual.view("RESIDUAL"); // DEBUGGING
// Check values
const double* valsE = _data->residual;
@@ -481,6 +391,86 @@
} // testIntegrateResidualAssembled
// ----------------------------------------------------------------------
+// Test integrateJacobian().
+void
+pylith::faults::TestFaultCohesiveKin::testIntegrateJacobian(void)
+{ // testIntegrateJacobian
+ topology::Mesh mesh;
+ FaultCohesiveKin fault;
+ topology::SolutionFields fields(mesh);
+ _initialize(&mesh, &fault, &fields);
+
+ const ALE::Obj<RealSection>& dispSection = fields.get("disp(t)").section();
+ CPPUNIT_ASSERT(!dispSection.isNull());
+
+ const int spaceDim = _data->spaceDim;
+ const ALE::Obj<SieveMesh>& sieveMesh = mesh.sieveMesh();
+ CPPUNIT_ASSERT(!sieveMesh.isNull());
+ const ALE::Obj<SieveMesh::label_sequence>& vertices = sieveMesh->depthStratum(0);
+ CPPUNIT_ASSERT(!vertices.isNull());
+ const SieveMesh::label_sequence::iterator verticesBegin = vertices->begin();
+ const SieveMesh::label_sequence::iterator verticesEnd = vertices->end();
+ int iVertex = 0;
+ for (SieveMesh::label_sequence::iterator v_iter=verticesBegin;
+ v_iter != verticesEnd;
+ ++v_iter, ++iVertex) {
+ dispSection->updatePoint(*v_iter, &_data->fieldT[iVertex*spaceDim]);
+ } // for
+
+ topology::Jacobian jacobian(fields);
+
+ const double t = 2.134;
+ fault.integrateJacobian(&jacobian, t, &fields);
+ CPPUNIT_ASSERT_EQUAL(false, fault.needNewJacobian());
+
+ jacobian.assemble("final_assembly");
+
+ //MatView(jacobian, PETSC_VIEWER_STDOUT_WORLD); // DEBUGGING
+
+ const double* valsE = _data->jacobian;
+ const int nrowsE = dispSection->sizeWithBC();
+ const int ncolsE = nrowsE;
+
+ int nrows = 0;
+ int ncols = 0;
+ PetscMat jacobianMat = jacobian.matrix();
+ MatGetSize(jacobianMat, &nrows, &ncols);
+ CPPUNIT_ASSERT_EQUAL(nrowsE, nrows);
+ CPPUNIT_ASSERT_EQUAL(ncolsE, ncols);
+
+ PetscMat jDense;
+ PetscMat jSparseAIJ;
+ MatConvert(jacobianMat, MATSEQAIJ, MAT_INITIAL_MATRIX, &jSparseAIJ);
+ MatConvert(jSparseAIJ, MATSEQDENSE, MAT_INITIAL_MATRIX, &jDense);
+
+ double_array vals(nrows*ncols);
+ int_array rows(nrows);
+ int_array cols(ncols);
+ for (int iRow=0; iRow < nrows; ++iRow)
+ rows[iRow] = iRow;
+ for (int iCol=0; iCol < ncols; ++iCol)
+ cols[iCol] = iCol;
+ MatGetValues(jDense, nrows, &rows[0], ncols, &cols[0], &vals[0]);
+ const double tolerance = 1.0e-06;
+ for (int iRow=0; iRow < nrows; ++iRow)
+ for (int iCol=0; iCol < ncols; ++iCol) {
+ const int index = ncols*iRow+iCol;
+ const double valE = 0.0;
+#if 0 // DEBUGGING
+ if (fabs(valE-vals[index]) > tolerance)
+ std::cout << "ERROR: iRow: " << iRow << ", iCol: " << iCol
+ << "valE: " << valE
+ << ", val: " << vals[index]
+ << std::endl;
+#endif // DEBUGGING
+ CPPUNIT_ASSERT_DOUBLES_EQUAL(valE, vals[index], tolerance);
+ } // for
+ MatDestroy(jDense);
+ MatDestroy(jSparseAIJ);
+ CPPUNIT_ASSERT_EQUAL(false, fault.needNewJacobian());
+} // testIntegrateJacobian
+
+// ----------------------------------------------------------------------
// Test integrateJacobianAssembled().
void
pylith::faults::TestFaultCohesiveKin::testIntegrateJacobianAssembled(void)
Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.hh
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.hh 2010-02-02 23:10:22 UTC (rev 16214)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.hh 2010-02-02 23:40:56 UTC (rev 16215)
@@ -92,12 +92,12 @@
/// Test integrateResidual().
void testIntegrateResidual(void);
+ /// Test integrateResidualAssembled().
+ void testIntegrateResidualAssembled(void);
+
/// Test integrateJacobian().
void testIntegrateJacobian(void);
- /// Test integrateResidualAssembled().
- void testIntegrateResidualAssembled(void);
-
/// Test integrateJacobianAssembled().
void testIntegrateJacobianAssembled(void);
Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinHex8.hh
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinHex8.hh 2010-02-02 23:10:22 UTC (rev 16214)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinHex8.hh 2010-02-02 23:40:56 UTC (rev 16215)
@@ -39,6 +39,7 @@
CPPUNIT_TEST( testInitialize );
CPPUNIT_TEST( testIntegrateResidual );
+ CPPUNIT_TEST( testIntegrateResidualAssembled );
CPPUNIT_TEST( testIntegrateJacobian );
CPPUNIT_TEST( testIntegrateJacobianAssembled );
CPPUNIT_TEST( testIntegrateJacobianAssembledLumped );
Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinQuad4.hh
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinQuad4.hh 2010-02-02 23:10:22 UTC (rev 16214)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinQuad4.hh 2010-02-02 23:40:56 UTC (rev 16215)
@@ -39,6 +39,7 @@
CPPUNIT_TEST( testInitialize );
CPPUNIT_TEST( testIntegrateResidual );
+ CPPUNIT_TEST( testIntegrateResidualAssembled );
CPPUNIT_TEST( testIntegrateJacobian );
CPPUNIT_TEST( testIntegrateJacobianAssembled );
CPPUNIT_TEST( testIntegrateJacobianAssembledLumped );
Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinTet4.hh
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinTet4.hh 2010-02-02 23:10:22 UTC (rev 16214)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinTet4.hh 2010-02-02 23:40:56 UTC (rev 16215)
@@ -39,6 +39,7 @@
CPPUNIT_TEST( testInitialize );
CPPUNIT_TEST( testIntegrateResidual );
+ CPPUNIT_TEST( testIntegrateResidualAssembled );
CPPUNIT_TEST( testIntegrateJacobian );
CPPUNIT_TEST( testIntegrateJacobianAssembled );
CPPUNIT_TEST( testIntegrateJacobianAssembledLumped );
Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinTri3.hh
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinTri3.hh 2010-02-02 23:10:22 UTC (rev 16214)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKinTri3.hh 2010-02-02 23:40:56 UTC (rev 16215)
@@ -39,6 +39,7 @@
CPPUNIT_TEST( testInitialize );
CPPUNIT_TEST( testIntegrateResidual );
+ CPPUNIT_TEST( testIntegrateResidualAssembled );
CPPUNIT_TEST( testIntegrateJacobian );
CPPUNIT_TEST( testIntegrateJacobianAssembled );
CPPUNIT_TEST( testIntegrateJacobianAssembledLumped );
Modified: short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveKinDataHex8.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveKinDataHex8.cc 2010-02-02 23:10:22 UTC (rev 16214)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveKinDataHex8.cc 2010-02-02 23:40:56 UTC (rev 16215)
@@ -206,10 +206,10 @@
+9.6,+5.6,+7.6, // 15
+9.8,+5.8,+7.8, // 16
+9.0,+5.0,+7.0, // 17
- 0.8, 0.8, 0.8, // 18 (constraint)
- 0.9, 0.9, 0.9, // 19 (constraint)
- 1.0, 1.0, 1.0, // 20 (constraint)
- 1.1, 1.1, 1.1, // 21 (constraint)
+ 0.8+1.82575588523, 0.8+-0.55566483464, 0.8+0.07938069066, // 18 (constraint)
+ 0.9+1.69682900001, 0.9+-0.56560966667, 0.9+0.14140241667, // 19 (constraint)
+ 1.0+1.51709826228, 1.0+-0.54615537442, 1.0+0.18205179147, // 20 (constraint)
+ 1.1+1.29378670385, 1.1+-0.49761027071, 1.1+0.19904410828, // 21 (constraint)
};
const double pylith::faults::CohesiveKinDataHex8::_residualIncr[] = {
@@ -229,10 +229,10 @@
+9.6,+5.6,+7.6, // 15
+9.8,+5.8,+7.8, // 16
+9.0,+5.0,+7.0, // 17
- 0.8, 0.8, 0.8, // 18 (constraint)
- 0.9, 0.9, 0.9, // 19 (constraint)
- 1.0, 1.0, 1.0, // 20 (constraint)
- 1.1, 1.1, 1.1, // 21 (constraint)
+ 0.8+1.82575588523, 0.8+-0.55566483464, 0.8+0.07938069066, // 18 (constraint)
+ 0.9+1.69682900001, 0.9+-0.56560966667, 0.9+0.14140241667, // 19 (constraint)
+ 1.0+1.51709826228, 1.0+-0.54615537442, 1.0+0.18205179147, // 20 (constraint)
+ 1.1+1.29378670385, 1.1+-0.49761027071, 1.1+0.19904410828, // 21 (constraint)
};
const double pylith::faults::CohesiveKinDataHex8::_jacobian[] = {
Modified: short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveKinDataQuad4.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveKinDataQuad4.cc 2010-02-02 23:10:22 UTC (rev 16214)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveKinDataQuad4.cc 2010-02-02 23:40:56 UTC (rev 16215)
@@ -156,8 +156,8 @@
0.0, 0.0,
-9.8, -8.8, // 8
-9.0, -8.0, // 9
- -0.4, -0.4, // 10
- -0.5, -0.5, // 11
+ -0.4+1.77538035254, -0.4+0.14794836271, // 10
+ -0.5+1.89546413727, -0.5+0.08241148423, // 11
};
const double pylith::faults::CohesiveKinDataQuad4::_residualIncr[] = {
@@ -169,8 +169,8 @@
0.0, 0.0,
-9.8, -8.8, // 8
-9.0, -8.0, // 9
- -0.4, -0.4, // 10
- -0.5, -0.5, // 11
+ -0.4+1.77538035254, -0.4+0.14794836271, // 10
+ -0.5+1.89546413727, -0.5+0.08241148423, // 11
};
const double pylith::faults::CohesiveKinDataQuad4::_jacobian[] = {
Modified: short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveKinDataTet4.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveKinDataTet4.cc 2010-02-02 23:10:22 UTC (rev 16214)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveKinDataTet4.cc 2010-02-02 23:40:56 UTC (rev 16215)
@@ -148,9 +148,9 @@
-9.7, -7.7, -8.7, // 7
-9.9, -7.9, -8.9, // 4
-9.1, -7.1, -8.1, // 5
- -0.4, -0.4, -0.4, // 10
- -0.5, -0.5, -0.5, // 11
- 0.4, 0.4, 0.4, // 12
+ -0.4+1.82575588523, -0.4+-0.55566483464, -0.4+0.07938069066, // 10
+ -0.5+1.69682900001, -0.5+-0.56560966667, -0.5+0.14140241667, // 11
+ 0.4+1.51709826228, 0.4+-0.54615537442, 0.4+0.18205179147, // 12
};
const double pylith::faults::CohesiveKinDataTet4::_residualIncr[] = {
@@ -162,9 +162,9 @@
-9.7, -7.7, -8.7, // 7
-9.9, -7.9, -8.9, // 4
-9.1, -7.1, -8.1, // 5
- -0.4, -0.4, -0.4, // 10
- -0.5, -0.5, -0.5, // 11
- 0.4, 0.4, 0.4, // 12
+ -0.4+1.82575588523, -0.4+-0.55566483464, -0.4+0.07938069066, // 10
+ -0.5+1.69682900001, -0.5+-0.56560966667, -0.5+0.14140241667, // 11
+ 0.4+1.51709826228, 0.4+-0.54615537442, 0.4+0.18205179147, // 12
};
const double pylith::faults::CohesiveKinDataTet4::_jacobian[] = {
Modified: short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveKinDataTri3.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveKinDataTri3.cc 2010-02-02 23:10:22 UTC (rev 16214)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveKinDataTri3.cc 2010-02-02 23:40:56 UTC (rev 16215)
@@ -151,8 +151,8 @@
0.0, 0.0,
+9.6, +8.6, // 6
+9.8, +8.8, // 7
- 0.3, 0.3, // 8
- 0.4, 0.4, // 9
+ 0.3+1.89546413727, 0.3+0.08241148423, // 8
+ 0.4+1.77538035254, 0.4+0.14794836271, // 9
};
const double pylith::faults::CohesiveKinDataTri3::_residualIncr[] = {
@@ -162,8 +162,8 @@
0.0, 0.0,
+9.6, +8.6, // 6
+9.8, +8.8, // 7
- 0.3, 0.3, // 8
- 0.4, 0.4, // 9
+ 0.3+1.89546413727, 0.3+0.08241148423, // 8
+ 0.4+1.77538035254, 0.4+0.14794836271, // 9
};
const double pylith::faults::CohesiveKinDataTri3::_jacobian[] = {
More information about the CIG-COMMITS
mailing list