[cig-commits] r16068 - in short/3D/PyLith/branches/pylith-friction: libsrc/feassemble unittests/libtests/feassemble unittests/libtests/feassemble/data
brad at geodynamics.org
brad at geodynamics.org
Fri Dec 4 10:58:07 PST 2009
Author: brad
Date: 2009-12-04 10:58:03 -0800 (Fri, 04 Dec 2009)
New Revision: 16068
Added:
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav1DLinear.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav1DLinear.hh
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav1DQuadratic.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav1DQuadratic.hh
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav2DLinear.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav2DLinear.hh
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav2DQuadratic.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav2DQuadratic.hh
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav3DLinear.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav3DLinear.hh
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav3DQuadratic.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav3DQuadratic.hh
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData1DLinear.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData1DLinear.hh
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData1DQuadratic.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData1DQuadratic.hh
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData2DLinear.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData2DLinear.hh
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData2DQuadratic.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData2DQuadratic.hh
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData3DLinear.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData3DLinear.hh
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData3DQuadratic.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData3DQuadratic.hh
Modified:
short/3D/PyLith/branches/pylith-friction/libsrc/feassemble/ElasticityExplicit.cc
short/3D/PyLith/branches/pylith-friction/libsrc/feassemble/ElasticityImplicit.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/Makefile.am
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit1DLinear.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit1DQuadratic.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit2DLinear.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit2DQuadratic.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit3DLinear.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit3DQuadratic.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicit.py
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitData1DLinear.cc
short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/generate.sh
Log:
Added unit tests for gravity with explicit integration of elasticity terms.
Modified: short/3D/PyLith/branches/pylith-friction/libsrc/feassemble/ElasticityExplicit.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/libsrc/feassemble/ElasticityExplicit.cc 2009-12-04 17:44:24 UTC (rev 16067)
+++ short/3D/PyLith/branches/pylith-friction/libsrc/feassemble/ElasticityExplicit.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -247,7 +247,41 @@
const double_array& basis = _quadrature->basis();
const double_array& basisDeriv = _quadrature->basisDeriv();
const double_array& jacobianDet = _quadrature->jacobianDet();
+ const double_array& quadPtsNondim = _quadrature->quadPts();
+ // Compute body force vector if gravity is being used.
+ if (0 != _gravityField) {
+ _logger->eventBegin(computeEvent);
+ const spatialdata::geocoords::CoordSys* cs = fields->mesh().coordsys();
+ assert(0 != cs);
+
+ // Get density at quadrature points for this cell
+ const double_array& density = _material->calcDensity();
+
+ quadPtsGlobal = quadPtsNondim;
+ _normalizer->dimensionalize(&quadPtsGlobal[0], quadPtsGlobal.size(),
+ lengthScale);
+
+ // Compute action for element body forces
+ for (int iQuad = 0; iQuad < numQuadPts; ++iQuad) {
+ const int err = _gravityField->query(&gravVec[0], gravVec.size(),
+ &quadPtsGlobal[0], spaceDim, cs);
+ if (err)
+ throw std::runtime_error("Unable to get gravity vector for point.");
+ _normalizer->nondimensionalize(&gravVec[0], gravVec.size(),
+ gravityScale);
+ const double wt = quadWts[iQuad] * jacobianDet[iQuad] * density[iQuad];
+ for (int iBasis = 0, iQ = iQuad * numBasis; iBasis < numBasis; ++iBasis) {
+ const double valI = wt * basis[iQ + iBasis];
+ for (int iDim = 0; iDim < spaceDim; ++iDim) {
+ _cellVector[iBasis * spaceDim + iDim] += valI * gravVec[iDim];
+ } // for
+ } // for
+ } // for
+ PetscLogFlops(numQuadPts * (2 + numBasis * (1 + 2 * spaceDim)));
+ _logger->eventEnd(computeEvent);
+ } // if
+
// Compute action for inertial terms
_logger->eventBegin(computeEvent);
const double_array& density = _material->calcDensity();
Modified: short/3D/PyLith/branches/pylith-friction/libsrc/feassemble/ElasticityImplicit.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/libsrc/feassemble/ElasticityImplicit.cc 2009-12-04 17:44:24 UTC (rev 16067)
+++ short/3D/PyLith/branches/pylith-friction/libsrc/feassemble/ElasticityImplicit.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -258,33 +258,32 @@
_logger->eventBegin(computeEvent);
const spatialdata::geocoords::CoordSys* cs = fields->mesh().coordsys();
assert(0 != cs);
-
+
// Get density at quadrature points for this cell
const double_array& density = _material->calcDensity();
quadPtsGlobal = quadPtsNondim;
_normalizer->dimensionalize(&quadPtsGlobal[0], quadPtsGlobal.size(),
- lengthScale);
+ lengthScale);
// Compute action for element body forces
- for (int iQuad=0; iQuad < numQuadPts; ++iQuad) {
- const int err = _gravityField->query(&gravVec[0], gravVec.size(),
- &quadPtsGlobal[0], spaceDim, cs);
- if (err)
- throw std::runtime_error("Unable to get gravity vector for point.");
- _normalizer->nondimensionalize(&gravVec[0], gravVec.size(),
- gravityScale);
- const double wt = quadWts[iQuad] * jacobianDet[iQuad] * density[iQuad];
- for (int iBasis=0, iQ=iQuad*numBasis;
- iBasis < numBasis; ++iBasis) {
- const double valI = wt*basis[iQ+iBasis];
- for (int iDim=0; iDim < spaceDim; ++iDim) {
- _cellVector[iBasis*spaceDim+iDim] += valI*gravVec[iDim];
- } // for
- } // for
+ for (int iQuad = 0; iQuad < numQuadPts; ++iQuad) {
+ const int err = _gravityField->query(&gravVec[0], gravVec.size(),
+ &quadPtsGlobal[0], spaceDim, cs);
+ if (err)
+ throw std::runtime_error("Unable to get gravity vector for point.");
+ _normalizer->nondimensionalize(&gravVec[0], gravVec.size(),
+ gravityScale);
+ const double wt = quadWts[iQuad] * jacobianDet[iQuad] * density[iQuad];
+ for (int iBasis = 0, iQ = iQuad * numBasis; iBasis < numBasis; ++iBasis) {
+ const double valI = wt * basis[iQ + iBasis];
+ for (int iDim = 0; iDim < spaceDim; ++iDim) {
+ _cellVector[iBasis * spaceDim + iDim] += valI * gravVec[iDim];
+ } // for
+ } // for
} // for
- PetscLogFlops(numQuadPts*(2+numBasis*(1+2*spaceDim)));
- _logger->eventEnd(computeEvent);
+ PetscLogFlops(numQuadPts * (2 + numBasis * (1 + 2 * spaceDim)));
+ _logger->eventEnd(computeEvent);
} // if
// Compute B(transpose) * sigma, first computing strains
Modified: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/Makefile.am
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/Makefile.am 2009-12-04 17:44:24 UTC (rev 16067)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/Makefile.am 2009-12-04 18:58:03 UTC (rev 16068)
@@ -53,6 +53,12 @@
TestElasticityExplicit2DQuadratic.cc \
TestElasticityExplicit3DLinear.cc \
TestElasticityExplicit3DQuadratic.cc \
+ TestElasticityExplicitGrav1DLinear.cc \
+ TestElasticityExplicitGrav1DQuadratic.cc \
+ TestElasticityExplicitGrav2DLinear.cc \
+ TestElasticityExplicitGrav2DQuadratic.cc \
+ TestElasticityExplicitGrav3DLinear.cc \
+ TestElasticityExplicitGrav3DQuadratic.cc \
TestElasticityImplicit.cc \
TestElasticityImplicit1DLinear.cc \
TestElasticityImplicit1DQuadratic.cc \
@@ -93,6 +99,12 @@
TestElasticityExplicit2DQuadratic.hh \
TestElasticityExplicit3DLinear.hh \
TestElasticityExplicit3DQuadratic.hh \
+ TestElasticityExplicitGrav1DLinear.hh \
+ TestElasticityExplicitGrav1DQuadratic.hh \
+ TestElasticityExplicitGrav2DLinear.hh \
+ TestElasticityExplicitGrav2DQuadratic.hh \
+ TestElasticityExplicitGrav3DLinear.hh \
+ TestElasticityExplicitGrav3DQuadratic.hh \
TestElasticityImplicit.hh \
TestElasticityImplicit1DLinear.hh \
TestElasticityImplicit1DQuadratic.hh \
@@ -166,6 +178,12 @@
data/ElasticityExplicitData2DQuadratic.cc \
data/ElasticityExplicitData3DLinear.cc \
data/ElasticityExplicitData3DQuadratic.cc \
+ data/ElasticityExplicitGravData1DLinear.cc \
+ data/ElasticityExplicitGravData1DQuadratic.cc \
+ data/ElasticityExplicitGravData2DLinear.cc \
+ data/ElasticityExplicitGravData2DQuadratic.cc \
+ data/ElasticityExplicitGravData3DLinear.cc \
+ data/ElasticityExplicitGravData3DQuadratic.cc \
data/ElasticityImplicitData1DLinear.cc \
data/ElasticityImplicitData1DQuadratic.cc \
data/ElasticityImplicitData2DLinear.cc \
@@ -228,6 +246,12 @@
data/ElasticityExplicitData2DQuadratic.hh \
data/ElasticityExplicitData3DLinear.hh \
data/ElasticityExplicitData3DQuadratic.hh \
+ data/ElasticityExplicitGravData1DLinear.hh \
+ data/ElasticityExplicitGravData1DQuadratic.hh \
+ data/ElasticityExplicitGravData2DLinear.hh \
+ data/ElasticityExplicitGravData2DQuadratic.hh \
+ data/ElasticityExplicitGravData3DLinear.hh \
+ data/ElasticityExplicitGravData3DQuadratic.hh \
data/ElasticityImplicitData1DLinear.hh \
data/ElasticityImplicitData1DQuadratic.hh \
data/ElasticityImplicitData2DLinear.hh \
Modified: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit.cc 2009-12-04 17:44:24 UTC (rev 16067)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -181,6 +181,13 @@
const int size = residualSection->sizeWithBC();
CPPUNIT_ASSERT_EQUAL(sizeE, size);
+#if 0
+ residual.view("RESIDUAL");
+ std::cout << "EXPECTED RESIDUAL" << std::endl;
+ for (int i=0; i < size; ++i)
+ std::cout << " " << valsE[i] << std::endl;
+#endif
+
const double tolerance = 1.0e-06;
for (int i=0; i < size; ++i)
if (fabs(valsE[i]) > 1.0)
@@ -408,9 +415,6 @@
mesh->coordsys(&cs);
mesh->nondimensionalize(normalizer);
- // Setup gravityField
- _gravityField = 0;
-
// Setup material
spatialdata::spatialdb::SimpleIOAscii iohandler;
iohandler.filename(_data->matDBFilename);
Modified: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit1DLinear.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit1DLinear.cc 2009-12-04 17:44:24 UTC (rev 16067)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit1DLinear.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -32,6 +32,7 @@
TestElasticityExplicit::setUp();
_data = new ElasticityExplicitData1DLinear();
+ _gravityField = 0;
GeometryLine1D geometry;
CPPUNIT_ASSERT(0 != _quadrature);
_quadrature->refGeometry(&geometry);
Modified: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit1DQuadratic.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit1DQuadratic.cc 2009-12-04 17:44:24 UTC (rev 16067)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit1DQuadratic.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -32,6 +32,7 @@
TestElasticityExplicit::setUp();
_data = new ElasticityExplicitData1DQuadratic();
+ _gravityField = 0;
GeometryLine1D geometry;
CPPUNIT_ASSERT(0 != _quadrature);
_quadrature->refGeometry(&geometry);
Modified: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit2DLinear.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit2DLinear.cc 2009-12-04 17:44:24 UTC (rev 16067)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit2DLinear.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -32,6 +32,7 @@
TestElasticityExplicit::setUp();
_data = new ElasticityExplicitData2DLinear();
+ _gravityField = 0;
CPPUNIT_ASSERT(0 != _quadrature);
GeometryTri2D geometry;
_quadrature->refGeometry(&geometry);
Modified: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit2DQuadratic.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit2DQuadratic.cc 2009-12-04 17:44:24 UTC (rev 16067)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit2DQuadratic.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -32,6 +32,7 @@
TestElasticityExplicit::setUp();
_data = new ElasticityExplicitData2DQuadratic();
+ _gravityField = 0;
CPPUNIT_ASSERT(0 != _quadrature);
GeometryTri2D geometry;
_quadrature->refGeometry(&geometry);
Modified: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit3DLinear.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit3DLinear.cc 2009-12-04 17:44:24 UTC (rev 16067)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit3DLinear.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -32,6 +32,7 @@
TestElasticityExplicit::setUp();
_data = new ElasticityExplicitData3DLinear();
+ _gravityField = 0;
CPPUNIT_ASSERT(0 != _quadrature);
GeometryTet3D geometry;
_quadrature->refGeometry(&geometry);
Modified: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit3DQuadratic.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit3DQuadratic.cc 2009-12-04 17:44:24 UTC (rev 16067)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicit3DQuadratic.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -32,6 +32,7 @@
TestElasticityExplicit::setUp();
_data = new ElasticityExplicitData3DQuadratic();
+ _gravityField = 0;
CPPUNIT_ASSERT(0 != _quadrature);
GeometryTet3D geometry;
_quadrature->refGeometry(&geometry);
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav1DLinear.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav1DLinear.cc (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav1DLinear.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,56 @@
+// -*- C++ -*-
+//
+// ----------------------------------------------------------------------
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ----------------------------------------------------------------------
+//
+
+#include <portinfo>
+
+#include "TestElasticityExplicitGrav1DLinear.hh" // Implementation of class methods
+
+#include "data/ElasticityExplicitGravData1DLinear.hh"
+
+#include "pylith/topology/Mesh.hh" // USES Quadrature<Mesh>
+#include "pylith/feassemble/Quadrature.hh" // USES Quadrature
+#include "pylith/feassemble/GeometryLine1D.hh" // USES GeometryLine1D
+#include "pylith/materials/ElasticStrain1D.hh" // USES ElasticStrain1D
+#include "spatialdata/spatialdb/GravityField.hh" // USES GravityField
+
+// ----------------------------------------------------------------------
+CPPUNIT_TEST_SUITE_REGISTRATION( pylith::feassemble::TestElasticityExplicitGrav1DLinear );
+
+// ----------------------------------------------------------------------
+// Setup testing data.
+void
+pylith::feassemble::TestElasticityExplicitGrav1DLinear::setUp(void)
+{ // setUp
+ TestElasticityExplicit::setUp();
+
+ _data = new ElasticityExplicitGravData1DLinear();
+ _gravityField = new spatialdata::spatialdb::GravityField();
+ CPPUNIT_ASSERT(0 != _quadrature);
+ CPPUNIT_ASSERT(0 != _gravityField);
+ GeometryLine1D geometry;
+ _quadrature->refGeometry(&geometry);
+
+ const double g = 1.0e8;
+ const double gravityDir[] = { -1.0, 0.0, 0.0};
+ _gravityField->gravAcceleration(g);
+ _gravityField->gravityDir(gravityDir[0], gravityDir[1], gravityDir[2]);
+
+ _material = new materials::ElasticStrain1D;
+ CPPUNIT_ASSERT(0 != _material);
+
+ CPPUNIT_ASSERT(0 != _data);
+ CPPUNIT_ASSERT_EQUAL(std::string("ElasticStrain1D"),
+ std::string(_data->matType));
+} // setUp
+
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav1DLinear.cc
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav1DLinear.hh
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav1DLinear.hh (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav1DLinear.hh 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,59 @@
+// -*- C++ -*-
+//
+// ----------------------------------------------------------------------
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ----------------------------------------------------------------------
+//
+
+/**
+ * @file unittests/libtests/feassemble/TestElasticityExplicitGrav1DLinear.hh
+ *
+ * @brief C++ TestElasticityExplicit object
+ *
+ * C++ unit testing for ElasticityExplicit with 1-D linear cells and gravity.
+ */
+
+#if !defined(pylith_feassemble_testelasticityexplicitgrav1dlinear_hh)
+#define pylith_feassemble_testelasticityexplicitgrav1dlinear_hh
+
+#include "TestElasticityExplicit.hh" // ISA TestElasticityExplicit
+
+/// Namespace for pylith package
+namespace pylith {
+ namespace feassemble {
+ class TestElasticityExplicitGrav1DLinear;
+ } // feassemble
+} // pylith
+
+/// C++ unit testing for ElasticityExplicit
+class pylith::feassemble::TestElasticityExplicitGrav1DLinear :
+ public TestElasticityExplicit
+{ // class TestElasticityExplicit1DLinear
+
+ // CPPUNIT TEST SUITE /////////////////////////////////////////////////
+ CPPUNIT_TEST_SUITE( TestElasticityExplicitGrav1DLinear );
+
+ CPPUNIT_TEST( testInitialize );
+ CPPUNIT_TEST( testIntegrateResidual );
+ CPPUNIT_TEST( testIntegrateJacobian );
+ CPPUNIT_TEST( testUpdateStateVars );
+
+ CPPUNIT_TEST_SUITE_END();
+
+ // PUBLIC METHODS /////////////////////////////////////////////////////
+public :
+
+ /// Setup testing data.
+ void setUp(void);
+
+}; // class TestElasticityExplicitGrav1DLinear
+
+#endif // pylith_feassemble_testelasticityexplicitgrav1dlinear_hh
+
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav1DLinear.hh
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav1DQuadratic.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav1DQuadratic.cc (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav1DQuadratic.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,56 @@
+// -*- C++ -*-
+//
+// ----------------------------------------------------------------------
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ----------------------------------------------------------------------
+//
+
+#include <portinfo>
+
+#include "TestElasticityExplicitGrav1DQuadratic.hh" // Implementation of class methods
+
+#include "data/ElasticityExplicitGravData1DQuadratic.hh"
+
+#include "pylith/topology/Mesh.hh" // USES Quadrature<Mesh>
+#include "pylith/feassemble/Quadrature.hh" // USES Quadrature
+#include "pylith/feassemble/GeometryLine1D.hh" // USES GeometryLine1D
+#include "pylith/materials/ElasticStrain1D.hh" // USES ElasticStrain1D
+#include "spatialdata/spatialdb/GravityField.hh" // USES GravityField
+
+// ----------------------------------------------------------------------
+CPPUNIT_TEST_SUITE_REGISTRATION( pylith::feassemble::TestElasticityExplicitGrav1DQuadratic );
+
+// ----------------------------------------------------------------------
+// Setup testing data.
+void
+pylith::feassemble::TestElasticityExplicitGrav1DQuadratic::setUp(void)
+{ // setUp
+ TestElasticityExplicit::setUp();
+
+ _data = new ElasticityExplicitGravData1DQuadratic();
+ _gravityField = new spatialdata::spatialdb::GravityField();
+ CPPUNIT_ASSERT(0 != _quadrature);
+ CPPUNIT_ASSERT(0 != _gravityField);
+ GeometryLine1D geometry;
+ _quadrature->refGeometry(&geometry);
+
+ const double g = 1.0e8;
+ const double gravityDir[] = { -1.0, 0.0, 0.0};
+ _gravityField->gravAcceleration(g);
+ _gravityField->gravityDir(gravityDir[0], gravityDir[1], gravityDir[2]);
+
+ _material = new materials::ElasticStrain1D;
+ CPPUNIT_ASSERT(0 != _material);
+
+ CPPUNIT_ASSERT(0 != _data);
+ CPPUNIT_ASSERT_EQUAL(std::string("ElasticStrain1D"),
+ std::string(_data->matType));
+} // setUp
+
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav1DQuadratic.cc
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav1DQuadratic.hh
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav1DQuadratic.hh (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav1DQuadratic.hh 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,59 @@
+// -*- C++ -*-
+//
+// ----------------------------------------------------------------------
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ----------------------------------------------------------------------
+//
+
+/**
+ * @file unittests/libtests/feassemble/TestElasticityExplicitGrav1DQuadratic.hh
+ *
+ * @brief C++ TestElasticityExplicit object
+ *
+ * C++ unit testing for ElasticityExplicit with 1-D quadratic cells and gravity.
+ */
+
+#if !defined(pylith_feassemble_testelasticityexplicitgrav1dquadratic_hh)
+#define pylith_feassemble_testelasticityexplicitgrav1dquadratic_hh
+
+#include "TestElasticityExplicit.hh" // ISA TestElasticityExplicit
+
+/// Namespace for pylith package
+namespace pylith {
+ namespace feassemble {
+ class TestElasticityExplicitGrav1DQuadratic;
+ } // feassemble
+} // pylith
+
+/// C++ unit testing for ElasticityExplicit
+class pylith::feassemble::TestElasticityExplicitGrav1DQuadratic :
+ public TestElasticityExplicit
+{ // class TestElasticityExplicit1DQuadratic
+
+ // CPPUNIT TEST SUITE /////////////////////////////////////////////////
+ CPPUNIT_TEST_SUITE( TestElasticityExplicitGrav1DQuadratic );
+
+ CPPUNIT_TEST( testInitialize );
+ CPPUNIT_TEST( testIntegrateResidual );
+ CPPUNIT_TEST( testIntegrateJacobian );
+ CPPUNIT_TEST( testUpdateStateVars );
+
+ CPPUNIT_TEST_SUITE_END();
+
+ // PUBLIC METHODS /////////////////////////////////////////////////////
+public :
+
+ /// Setup testing data.
+ void setUp(void);
+
+}; // class TestElasticityExplicitGrav1DQuadratic
+
+#endif // pylith_feassemble_testelasticityexplicitgrav1dquadratic_hh
+
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav1DQuadratic.hh
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav2DLinear.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav2DLinear.cc (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav2DLinear.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,56 @@
+// -*- C++ -*-
+//
+// ----------------------------------------------------------------------
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ----------------------------------------------------------------------
+//
+
+#include <portinfo>
+
+#include "TestElasticityExplicitGrav2DLinear.hh" // Implementation of class methods
+
+#include "data/ElasticityExplicitGravData2DLinear.hh"
+
+#include "pylith/topology/Mesh.hh" // USES Quadrature<Mesh>
+#include "pylith/feassemble/Quadrature.hh" // USES Quadrature
+#include "pylith/feassemble/GeometryTri2D.hh" // USES GeometryTri2D
+#include "pylith/materials/ElasticPlaneStrain.hh" // USES ElasticPlaneStrain
+#include "spatialdata/spatialdb/GravityField.hh" // USES GravityField
+
+// ----------------------------------------------------------------------
+CPPUNIT_TEST_SUITE_REGISTRATION( pylith::feassemble::TestElasticityExplicitGrav2DLinear );
+
+// ----------------------------------------------------------------------
+// Setup testing data.
+void
+pylith::feassemble::TestElasticityExplicitGrav2DLinear::setUp(void)
+{ // setUp
+ TestElasticityExplicit::setUp();
+
+ _data = new ElasticityExplicitGravData2DLinear();
+ _gravityField = new spatialdata::spatialdb::GravityField();
+ CPPUNIT_ASSERT(0 != _quadrature);
+ CPPUNIT_ASSERT(0 != _gravityField);
+ GeometryTri2D geometry;
+ _quadrature->refGeometry(&geometry);
+
+ const double g = 1.0e8;
+ const double gravityDir[] = { 0.0, -1.0, 0.0 };
+ _gravityField->gravAcceleration(g);
+ _gravityField->gravityDir(gravityDir[0], gravityDir[1], gravityDir[2]);
+
+ _material = new materials::ElasticPlaneStrain;
+ CPPUNIT_ASSERT(0 != _material);
+
+ CPPUNIT_ASSERT(0 != _data);
+ CPPUNIT_ASSERT_EQUAL(std::string("ElasticPlaneStrain"),
+ std::string(_data->matType));
+} // setUp
+
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav2DLinear.cc
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav2DLinear.hh
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav2DLinear.hh (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav2DLinear.hh 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,59 @@
+// -*- C++ -*-
+//
+// ----------------------------------------------------------------------
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ----------------------------------------------------------------------
+//
+
+/**
+ * @file unittests/libtests/feassemble/TestElasticityExplicitGrav2DLinear.hh
+ *
+ * @brief C++ TestElasticityExplicit object
+ *
+ * C++ unit testing for ElasticityExplicit with 2-D linear cells and gravity.
+ */
+
+#if !defined(pylith_feassemble_testelasticityexplicitgrav2dlinear_hh)
+#define pylith_feassemble_testelasticityexplicitgrav2dlinear_hh
+
+#include "TestElasticityExplicit.hh" // ISA TestElasticityExplicit
+
+/// Namespace for pylith package
+namespace pylith {
+ namespace feassemble {
+ class TestElasticityExplicitGrav2DLinear;
+ } // feassemble
+} // pylith
+
+/// C++ unit testing for ElasticityExplicit
+class pylith::feassemble::TestElasticityExplicitGrav2DLinear :
+ public TestElasticityExplicit
+{ // class TestElasticityExplicit2DLinear
+
+ // CPPUNIT TEST SUITE /////////////////////////////////////////////////
+ CPPUNIT_TEST_SUITE( TestElasticityExplicitGrav2DLinear );
+
+ CPPUNIT_TEST( testInitialize );
+ CPPUNIT_TEST( testIntegrateResidual );
+ CPPUNIT_TEST( testIntegrateJacobian );
+ CPPUNIT_TEST( testUpdateStateVars );
+
+ CPPUNIT_TEST_SUITE_END();
+
+ // PUBLIC METHODS /////////////////////////////////////////////////////
+public :
+
+ /// Setup testing data.
+ void setUp(void);
+
+}; // class TestElasticityExplicitGrav2DLinear
+
+#endif // pylith_feassemble_testelasticityexplicitgrav2dlinear_hh
+
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav2DLinear.hh
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav2DQuadratic.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav2DQuadratic.cc (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav2DQuadratic.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,56 @@
+// -*- C++ -*-
+//
+// ----------------------------------------------------------------------
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ----------------------------------------------------------------------
+//
+
+#include <portinfo>
+
+#include "TestElasticityExplicitGrav2DQuadratic.hh" // Implementation of class methods
+
+#include "data/ElasticityExplicitGravData2DQuadratic.hh"
+
+#include "pylith/topology/Mesh.hh" // USES Quadrature<Mesh>
+#include "pylith/feassemble/Quadrature.hh" // USES Quadrature
+#include "pylith/feassemble/GeometryTri2D.hh" // USES GeometryTri2D
+#include "pylith/materials/ElasticPlaneStrain.hh" // USES ElasticPlaneStrain
+#include "spatialdata/spatialdb/GravityField.hh" // USES GravityField
+
+// ----------------------------------------------------------------------
+CPPUNIT_TEST_SUITE_REGISTRATION( pylith::feassemble::TestElasticityExplicitGrav2DQuadratic );
+
+// ----------------------------------------------------------------------
+// Setup testing data.
+void
+pylith::feassemble::TestElasticityExplicitGrav2DQuadratic::setUp(void)
+{ // setUp
+ TestElasticityExplicit::setUp();
+
+ _data = new ElasticityExplicitGravData2DQuadratic();
+ _gravityField = new spatialdata::spatialdb::GravityField();
+ CPPUNIT_ASSERT(0 != _quadrature);
+ CPPUNIT_ASSERT(0 != _gravityField);
+ GeometryTri2D geometry;
+ _quadrature->refGeometry(&geometry);
+
+ const double g = 1.0e8;
+ const double gravityDir[] = { 0.0, -1.0, 0.0 };
+ _gravityField->gravAcceleration(g);
+ _gravityField->gravityDir(gravityDir[0], gravityDir[1], gravityDir[2]);
+
+ _material = new materials::ElasticPlaneStrain;
+ CPPUNIT_ASSERT(0 != _material);
+
+ CPPUNIT_ASSERT(0 != _data);
+ CPPUNIT_ASSERT_EQUAL(std::string("ElasticPlaneStrain"),
+ std::string(_data->matType));
+} // setUp
+
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav2DQuadratic.cc
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav2DQuadratic.hh
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav2DQuadratic.hh (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav2DQuadratic.hh 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,59 @@
+// -*- C++ -*-
+//
+// ----------------------------------------------------------------------
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ----------------------------------------------------------------------
+//
+
+/**
+ * @file unittests/libtests/feassemble/TestElasticityExplicitGrav2DQuadratic.hh
+ *
+ * @brief C++ TestElasticityExplicit object
+ *
+ * C++ unit testing for ElasticityExplicit with 2-D quadratic cells and gravity.
+ */
+
+#if !defined(pylith_feassemble_testelasticityexplicitgrav2dquadratic_hh)
+#define pylith_feassemble_testelasticityexplicitgrav2dquadratic_hh
+
+#include "TestElasticityExplicit.hh" // ISA TestElasticityExplicit
+
+/// Namespace for pylith package
+namespace pylith {
+ namespace feassemble {
+ class TestElasticityExplicitGrav2DQuadratic;
+ } // feassemble
+} // pylith
+
+/// C++ unit testing for ElasticityExplicit
+class pylith::feassemble::TestElasticityExplicitGrav2DQuadratic :
+ public TestElasticityExplicit
+{ // class TestElasticityExplicit2DQuadratic
+
+ // CPPUNIT TEST SUITE /////////////////////////////////////////////////
+ CPPUNIT_TEST_SUITE( TestElasticityExplicitGrav2DQuadratic );
+
+ CPPUNIT_TEST( testInitialize );
+ CPPUNIT_TEST( testIntegrateResidual );
+ CPPUNIT_TEST( testIntegrateJacobian );
+ CPPUNIT_TEST( testUpdateStateVars );
+
+ CPPUNIT_TEST_SUITE_END();
+
+ // PUBLIC METHODS /////////////////////////////////////////////////////
+public :
+
+ /// Setup testing data.
+ void setUp(void);
+
+}; // class TestElasticityExplicitGrav2DQuadratic
+
+#endif // pylith_feassemble_testelasticityexplicitgrav2dquadratic_hh
+
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav2DQuadratic.hh
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav3DLinear.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav3DLinear.cc (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav3DLinear.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,55 @@
+// -*- C++ -*-
+//
+// ----------------------------------------------------------------------
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ----------------------------------------------------------------------
+//
+
+#include <portinfo>
+
+#include "TestElasticityExplicitGrav3DLinear.hh" // Implementation of class methods
+
+#include "data/ElasticityExplicitGravData3DLinear.hh"
+
+#include "pylith/topology/Mesh.hh" // USES Quadrature<Mesh>
+#include "pylith/feassemble/Quadrature.hh" // USES Quadrature
+#include "pylith/feassemble/Quadrature3D.hh" // USES Quadrature3D
+#include "pylith/feassemble/GeometryTet3D.hh" // USES GeometryTet3D
+#include "pylith/materials/ElasticIsotropic3D.hh" // USES ElasticIsotropic3D
+#include "spatialdata/spatialdb/GravityField.hh" // USES GravityField
+
+// ----------------------------------------------------------------------
+CPPUNIT_TEST_SUITE_REGISTRATION( pylith::feassemble::TestElasticityExplicitGrav3DLinear );
+
+// ----------------------------------------------------------------------
+// Setup testing data.
+void
+pylith::feassemble::TestElasticityExplicitGrav3DLinear::setUp(void)
+{ // setUp
+ TestElasticityExplicit::setUp();
+
+ _data = new ElasticityExplicitGravData3DLinear();
+ _gravityField = new spatialdata::spatialdb::GravityField();
+ CPPUNIT_ASSERT(0 != _quadrature);
+ CPPUNIT_ASSERT(0 != _gravityField);
+ GeometryTet3D geometry;
+ _quadrature->refGeometry(&geometry);
+
+ const double g = 1.0e8;
+ _gravityField->gravAcceleration(g);
+
+ _material = new materials::ElasticIsotropic3D;
+ CPPUNIT_ASSERT(0 != _material);
+
+ CPPUNIT_ASSERT(0 != _data);
+ CPPUNIT_ASSERT_EQUAL(std::string("ElasticIsotropic3D"),
+ std::string(_data->matType));
+} // setUp
+
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav3DLinear.cc
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav3DLinear.hh
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav3DLinear.hh (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav3DLinear.hh 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,59 @@
+// -*- C++ -*-
+//
+// ----------------------------------------------------------------------
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ----------------------------------------------------------------------
+//
+
+/**
+ * @file unittests/libtests/feassemble/TestElasticityExplicitGrav3DLinear.hh
+ *
+ * @brief C++ TestElasticityExplicit object
+ *
+ * C++ unit testing for ElasticityExplicit with 1-D linear cells and gravity.
+ */
+
+#if !defined(pylith_feassemble_testelasticityexplicitgrav3dlinear_hh)
+#define pylith_feassemble_testelasticityexplicitgrav3dlinear_hh
+
+#include "TestElasticityExplicit.hh" // ISA TestElasticityExplicit
+
+/// Namespace for pylith package
+namespace pylith {
+ namespace feassemble {
+ class TestElasticityExplicitGrav3DLinear;
+ } // feassemble
+} // pylith
+
+/// C++ unit testing for ElasticityExplicit
+class pylith::feassemble::TestElasticityExplicitGrav3DLinear :
+ public TestElasticityExplicit
+{ // class TestElasticityExplicit3DLinear
+
+ // CPPUNIT TEST SUITE /////////////////////////////////////////////////
+ CPPUNIT_TEST_SUITE( TestElasticityExplicitGrav3DLinear );
+
+ CPPUNIT_TEST( testInitialize );
+ CPPUNIT_TEST( testIntegrateResidual );
+ CPPUNIT_TEST( testIntegrateJacobian );
+ CPPUNIT_TEST( testUpdateStateVars );
+
+ CPPUNIT_TEST_SUITE_END();
+
+ // PUBLIC METHODS /////////////////////////////////////////////////////
+public :
+
+ /// Setup testing data.
+ void setUp(void);
+
+}; // class TestElasticityExplicitGrav3DLinear
+
+#endif // pylith_feassemble_testelasticityexplicitgrav3dlinear_hh
+
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav3DLinear.hh
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav3DQuadratic.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav3DQuadratic.cc (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav3DQuadratic.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,55 @@
+// -*- C++ -*-
+//
+// ----------------------------------------------------------------------
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ----------------------------------------------------------------------
+//
+
+#include <portinfo>
+
+#include "TestElasticityExplicitGrav3DQuadratic.hh" // Implementation of class methods
+
+#include "data/ElasticityExplicitGravData3DQuadratic.hh"
+
+#include "pylith/topology/Mesh.hh" // USES Quadrature<Mesh>
+#include "pylith/feassemble/Quadrature.hh" // USES Quadrature
+#include "pylith/feassemble/Quadrature3D.hh" // USES Quadrature3D
+#include "pylith/feassemble/GeometryTet3D.hh" // USES GeometryTet3D
+#include "pylith/materials/ElasticIsotropic3D.hh" // USES ElasticIsotropic3D
+#include "spatialdata/spatialdb/GravityField.hh" // USES GravityField
+
+// ----------------------------------------------------------------------
+CPPUNIT_TEST_SUITE_REGISTRATION( pylith::feassemble::TestElasticityExplicitGrav3DQuadratic );
+
+// ----------------------------------------------------------------------
+// Setup testing data.
+void
+pylith::feassemble::TestElasticityExplicitGrav3DQuadratic::setUp(void)
+{ // setUp
+ TestElasticityExplicit::setUp();
+
+ _data = new ElasticityExplicitGravData3DQuadratic();
+ _gravityField = new spatialdata::spatialdb::GravityField();
+ CPPUNIT_ASSERT(0 != _quadrature);
+ CPPUNIT_ASSERT(0 != _gravityField);
+ GeometryTet3D geometry;
+ _quadrature->refGeometry(&geometry);
+
+ const double g = 1.0e8;
+ _gravityField->gravAcceleration(g);
+
+ _material = new materials::ElasticIsotropic3D;
+ CPPUNIT_ASSERT(0 != _material);
+
+ CPPUNIT_ASSERT(0 != _data);
+ CPPUNIT_ASSERT_EQUAL(std::string("ElasticIsotropic3D"),
+ std::string(_data->matType));
+} // setUp
+
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav3DQuadratic.cc
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav3DQuadratic.hh
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav3DQuadratic.hh (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav3DQuadratic.hh 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,59 @@
+// -*- C++ -*-
+//
+// ----------------------------------------------------------------------
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ----------------------------------------------------------------------
+//
+
+/**
+ * @file unittests/libtests/feassemble/TestElasticityExplicitGrav3DQuadratic.hh
+ *
+ * @brief C++ TestElasticityExplicit object
+ *
+ * C++ unit testing for ElasticityExplicit with 1-D quadratic cells and gravity.
+ */
+
+#if !defined(pylith_feassemble_testelasticityexplicitgrav3dquadratic_hh)
+#define pylith_feassemble_testelasticityexplicitgrav3dquadratic_hh
+
+#include "TestElasticityExplicit.hh" // ISA TestElasticityExplicit
+
+/// Namespace for pylith package
+namespace pylith {
+ namespace feassemble {
+ class TestElasticityExplicitGrav3DQuadratic;
+ } // feassemble
+} // pylith
+
+/// C++ unit testing for ElasticityExplicit
+class pylith::feassemble::TestElasticityExplicitGrav3DQuadratic :
+ public TestElasticityExplicit
+{ // class TestElasticityExplicit3DQuadratic
+
+ // CPPUNIT TEST SUITE /////////////////////////////////////////////////
+ CPPUNIT_TEST_SUITE( TestElasticityExplicitGrav3DQuadratic );
+
+ CPPUNIT_TEST( testInitialize );
+ CPPUNIT_TEST( testIntegrateResidual );
+ CPPUNIT_TEST( testIntegrateJacobian );
+ CPPUNIT_TEST( testUpdateStateVars );
+
+ CPPUNIT_TEST_SUITE_END();
+
+ // PUBLIC METHODS /////////////////////////////////////////////////////
+public :
+
+ /// Setup testing data.
+ void setUp(void);
+
+}; // class TestElasticityExplicitGrav3DQuadratic
+
+#endif // pylith_feassemble_testelasticityexplicitgrav3dquadratic_hh
+
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/TestElasticityExplicitGrav3DQuadratic.hh
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Modified: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicit.py
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicit.py 2009-12-04 17:44:24 UTC (rev 16067)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicit.py 2009-12-04 18:58:03 UTC (rev 16068)
@@ -53,7 +53,7 @@
dispResult = integrator.fieldT - integrator.fieldTmdt
residual = 1.0/integrator.dt**2 * numpy.dot(M, dispResult) - \
numpy.dot(K, integrator.fieldT)
- return residual
+ return residual.flatten()
def calculateJacobian(self, integrator):
Modified: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitData1DLinear.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitData1DLinear.cc 2009-12-04 17:44:24 UTC (rev 16067)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitData1DLinear.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -89,7 +89,7 @@
};
const double pylith::feassemble::ElasticityExplicitData1DLinear::_valsResidual[] = {
- 1.60042188e+10,
+ 1.60042187e+10,
-1.59957812e+10,
};
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData1DLinear.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData1DLinear.cc (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData1DLinear.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,135 @@
+// -*- C++ -*-
+//
+// ======================================================================
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ======================================================================
+//
+
+// DO NOT EDIT THIS FILE
+// This file was generated from python application integratorelasticity.
+
+#include "ElasticityExplicitGravData1DLinear.hh"
+
+const int pylith::feassemble::ElasticityExplicitGravData1DLinear::_spaceDim = 1;
+
+const int pylith::feassemble::ElasticityExplicitGravData1DLinear::_cellDim = 1;
+
+const int pylith::feassemble::ElasticityExplicitGravData1DLinear::_numVertices = 2;
+
+const int pylith::feassemble::ElasticityExplicitGravData1DLinear::_numCells = 1;
+
+const int pylith::feassemble::ElasticityExplicitGravData1DLinear::_numBasis = 2;
+
+const int pylith::feassemble::ElasticityExplicitGravData1DLinear::_numQuadPts = 1;
+
+const char* pylith::feassemble::ElasticityExplicitGravData1DLinear::_matType = "ElasticStrain1D";
+
+const char* pylith::feassemble::ElasticityExplicitGravData1DLinear::_matDBFilename = "data/elasticstrain1d.spatialdb";
+
+const int pylith::feassemble::ElasticityExplicitGravData1DLinear::_matId = 0;
+
+const char* pylith::feassemble::ElasticityExplicitGravData1DLinear::_matLabel = "elastic strain 1-D";
+
+const double pylith::feassemble::ElasticityExplicitGravData1DLinear::_dt = 1.00000000e-02;
+
+const double pylith::feassemble::ElasticityExplicitGravData1DLinear::_gravityVec[] = {
+ -1.00000000e+08,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DLinear::_vertices[] = {
+ -2.50000000e-01,
+ 2.00000000e+00,
+};
+
+const int pylith::feassemble::ElasticityExplicitGravData1DLinear::_cells[] = {
+0,1,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DLinear::_verticesRef[] = {
+ -1.00000000e+00,
+ 1.00000000e+00,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DLinear::_quadPts[] = {
+ 0.00000000e+00,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DLinear::_quadWts[] = {
+ 2.00000000e+00,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DLinear::_basis[] = {
+ 5.00000000e-01,
+ 5.00000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DLinear::_basisDerivRef[] = {
+ -5.00000000e-01,
+ 5.00000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DLinear::_fieldTIncr[] = {
+ 1.20000000e+00,
+ 1.70000000e+00,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DLinear::_fieldT[] = {
+ 1.10000000e+00,
+ 1.50000000e+00,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DLinear::_fieldTmdt[] = {
+ 1.00000000e+00,
+ 1.30000000e+00,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DLinear::_valsResidual[] = {
+ -2.65245781e+11,
+ -2.97245781e+11,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DLinear::_valsJacobian[] = {
+ 1.40625000e+07,
+ 1.40625000e+07,
+ 1.40625000e+07,
+ 1.40625000e+07,
+};
+
+pylith::feassemble::ElasticityExplicitGravData1DLinear::ElasticityExplicitGravData1DLinear(void)
+{ // constructor
+ spaceDim = _spaceDim;
+ cellDim = _cellDim;
+ numVertices = _numVertices;
+ numCells = _numCells;
+ numBasis = _numBasis;
+ numQuadPts = _numQuadPts;
+ matType = const_cast<char*>(_matType);
+ matDBFilename = const_cast<char*>(_matDBFilename);
+ matId = _matId;
+ matLabel = const_cast<char*>(_matLabel);
+ dt = _dt;
+ gravityVec = const_cast<double*>(_gravityVec);
+ vertices = const_cast<double*>(_vertices);
+ cells = const_cast<int*>(_cells);
+ verticesRef = const_cast<double*>(_verticesRef);
+ quadPts = const_cast<double*>(_quadPts);
+ quadWts = const_cast<double*>(_quadWts);
+ basis = const_cast<double*>(_basis);
+ basisDerivRef = const_cast<double*>(_basisDerivRef);
+ fieldTIncr = const_cast<double*>(_fieldTIncr);
+ fieldT = const_cast<double*>(_fieldT);
+ fieldTmdt = const_cast<double*>(_fieldTmdt);
+ valsResidual = const_cast<double*>(_valsResidual);
+ valsJacobian = const_cast<double*>(_valsJacobian);
+} // constructor
+
+pylith::feassemble::ElasticityExplicitGravData1DLinear::~ElasticityExplicitGravData1DLinear(void)
+{}
+
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData1DLinear.cc
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData1DLinear.hh
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData1DLinear.hh (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData1DLinear.hh 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,92 @@
+// -*- C++ -*-
+//
+// ======================================================================
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ======================================================================
+//
+
+// DO NOT EDIT THIS FILE
+// This file was generated from python application integratorelasticity.
+
+#if !defined(pylith_feassemble_elasticityexplicitgravdata1dlinear_hh)
+#define pylith_feassemble_elasticityexplicitgravdata1dlinear_hh
+
+#include "IntegratorData.hh"
+
+namespace pylith {
+ namespace feassemble {
+ class ElasticityExplicitGravData1DLinear;
+ } // pylith
+} // feassemble
+
+class pylith::feassemble::ElasticityExplicitGravData1DLinear : public IntegratorData
+{
+
+public:
+
+ /// Constructor
+ ElasticityExplicitGravData1DLinear(void);
+
+ /// Destructor
+ ~ElasticityExplicitGravData1DLinear(void);
+
+private:
+
+ static const int _spaceDim;
+
+ static const int _cellDim;
+
+ static const int _numVertices;
+
+ static const int _numCells;
+
+ static const int _numBasis;
+
+ static const int _numQuadPts;
+
+ static const char* _matType;
+
+ static const char* _matDBFilename;
+
+ static const int _matId;
+
+ static const char* _matLabel;
+
+ static const double _dt;
+
+ static const double _gravityVec[];
+
+ static const double _vertices[];
+
+ static const int _cells[];
+
+ static const double _verticesRef[];
+
+ static const double _quadPts[];
+
+ static const double _quadWts[];
+
+ static const double _basis[];
+
+ static const double _basisDerivRef[];
+
+ static const double _fieldTIncr[];
+
+ static const double _fieldT[];
+
+ static const double _fieldTmdt[];
+
+ static const double _valsResidual[];
+
+ static const double _valsJacobian[];
+
+};
+
+#endif // pylith_feassemble_elasticityexplicitgravdata1dlinear_hh
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData1DLinear.hh
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData1DQuadratic.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData1DQuadratic.cc (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData1DQuadratic.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,155 @@
+// -*- C++ -*-
+//
+// ======================================================================
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ======================================================================
+//
+
+// DO NOT EDIT THIS FILE
+// This file was generated from python application integratorelasticity.
+
+#include "ElasticityExplicitGravData1DQuadratic.hh"
+
+const int pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_spaceDim = 1;
+
+const int pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_cellDim = 1;
+
+const int pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_numVertices = 3;
+
+const int pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_numCells = 1;
+
+const int pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_numBasis = 3;
+
+const int pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_numQuadPts = 2;
+
+const char* pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_matType = "ElasticStrain1D";
+
+const char* pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_matDBFilename = "data/elasticstrain1d.spatialdb";
+
+const int pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_matId = 0;
+
+const char* pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_matLabel = "elastic strain 1-D";
+
+const double pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_dt = 1.00000000e-02;
+
+const double pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_gravityVec[] = {
+ -1.00000000e+08,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_vertices[] = {
+ -2.50000000e-01,
+ 2.00000000e+00,
+ 8.75000000e-01,
+};
+
+const int pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_cells[] = {
+0,1,2,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_verticesRef[] = {
+ -1.00000000e+00,
+ 1.00000000e+00,
+ 0.00000000e+00,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_quadPts[] = {
+ -5.77350269e-01,
+ 5.77350269e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_quadWts[] = {
+ 1.00000000e+00, 1.00000000e+00,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_basis[] = {
+ 4.55341801e-01,
+ -1.22008468e-01,
+ 6.66666667e-01,
+ -1.22008468e-01,
+ 4.55341801e-01,
+ 6.66666667e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_basisDerivRef[] = {
+ -1.07735027e+00,
+ -7.73502692e-02,
+ 1.15470054e+00,
+ 7.73502692e-02,
+ 1.07735027e+00,
+ -1.15470054e+00,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_fieldTIncr[] = {
+ 1.20000000e+00,
+ 0.00000000e+00,
+ 1.70000000e+00,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_fieldT[] = {
+ 1.10000000e+00,
+ 1.00000000e-01,
+ 1.50000000e+00,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_fieldTmdt[] = {
+ 1.00000000e+00,
+ 1.00000000e-01,
+ 1.30000000e+00,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_valsResidual[] = {
+ -3.77481250e+10,
+ 4.22509375e+10,
+ -5.66994375e+11,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData1DQuadratic::_valsJacobian[] = {
+ 6.25000000e+06,
+ -3.12500000e+06,
+ 6.25000000e+06,
+ -3.12500000e+06,
+ 6.25000000e+06,
+ 6.25000000e+06,
+ 6.25000000e+06,
+ 6.25000000e+06,
+ 2.50000000e+07,
+};
+
+pylith::feassemble::ElasticityExplicitGravData1DQuadratic::ElasticityExplicitGravData1DQuadratic(void)
+{ // constructor
+ spaceDim = _spaceDim;
+ cellDim = _cellDim;
+ numVertices = _numVertices;
+ numCells = _numCells;
+ numBasis = _numBasis;
+ numQuadPts = _numQuadPts;
+ matType = const_cast<char*>(_matType);
+ matDBFilename = const_cast<char*>(_matDBFilename);
+ matId = _matId;
+ matLabel = const_cast<char*>(_matLabel);
+ dt = _dt;
+ gravityVec = const_cast<double*>(_gravityVec);
+ vertices = const_cast<double*>(_vertices);
+ cells = const_cast<int*>(_cells);
+ verticesRef = const_cast<double*>(_verticesRef);
+ quadPts = const_cast<double*>(_quadPts);
+ quadWts = const_cast<double*>(_quadWts);
+ basis = const_cast<double*>(_basis);
+ basisDerivRef = const_cast<double*>(_basisDerivRef);
+ fieldTIncr = const_cast<double*>(_fieldTIncr);
+ fieldT = const_cast<double*>(_fieldT);
+ fieldTmdt = const_cast<double*>(_fieldTmdt);
+ valsResidual = const_cast<double*>(_valsResidual);
+ valsJacobian = const_cast<double*>(_valsJacobian);
+} // constructor
+
+pylith::feassemble::ElasticityExplicitGravData1DQuadratic::~ElasticityExplicitGravData1DQuadratic(void)
+{}
+
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData1DQuadratic.cc
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData1DQuadratic.hh
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData1DQuadratic.hh (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData1DQuadratic.hh 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,92 @@
+// -*- C++ -*-
+//
+// ======================================================================
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ======================================================================
+//
+
+// DO NOT EDIT THIS FILE
+// This file was generated from python application integratorelasticity.
+
+#if !defined(pylith_feassemble_elasticityexplicitgravdata1dquadratic_hh)
+#define pylith_feassemble_elasticityexplicitgravdata1dquadratic_hh
+
+#include "IntegratorData.hh"
+
+namespace pylith {
+ namespace feassemble {
+ class ElasticityExplicitGravData1DQuadratic;
+ } // pylith
+} // feassemble
+
+class pylith::feassemble::ElasticityExplicitGravData1DQuadratic : public IntegratorData
+{
+
+public:
+
+ /// Constructor
+ ElasticityExplicitGravData1DQuadratic(void);
+
+ /// Destructor
+ ~ElasticityExplicitGravData1DQuadratic(void);
+
+private:
+
+ static const int _spaceDim;
+
+ static const int _cellDim;
+
+ static const int _numVertices;
+
+ static const int _numCells;
+
+ static const int _numBasis;
+
+ static const int _numQuadPts;
+
+ static const char* _matType;
+
+ static const char* _matDBFilename;
+
+ static const int _matId;
+
+ static const char* _matLabel;
+
+ static const double _dt;
+
+ static const double _gravityVec[];
+
+ static const double _vertices[];
+
+ static const int _cells[];
+
+ static const double _verticesRef[];
+
+ static const double _quadPts[];
+
+ static const double _quadWts[];
+
+ static const double _basis[];
+
+ static const double _basisDerivRef[];
+
+ static const double _fieldTIncr[];
+
+ static const double _fieldT[];
+
+ static const double _fieldTmdt[];
+
+ static const double _valsResidual[];
+
+ static const double _valsJacobian[];
+
+};
+
+#endif // pylith_feassemble_elasticityexplicitgravdata1dquadratic_hh
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData1DQuadratic.hh
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData2DLinear.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData2DLinear.cc (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData2DLinear.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,155 @@
+// -*- C++ -*-
+//
+// ======================================================================
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ======================================================================
+//
+
+// DO NOT EDIT THIS FILE
+// This file was generated from python application integratorelasticity.
+
+#include "ElasticityExplicitGravData2DLinear.hh"
+
+const int pylith::feassemble::ElasticityExplicitGravData2DLinear::_spaceDim = 2;
+
+const int pylith::feassemble::ElasticityExplicitGravData2DLinear::_cellDim = 2;
+
+const int pylith::feassemble::ElasticityExplicitGravData2DLinear::_numVertices = 3;
+
+const int pylith::feassemble::ElasticityExplicitGravData2DLinear::_numCells = 1;
+
+const int pylith::feassemble::ElasticityExplicitGravData2DLinear::_numBasis = 3;
+
+const int pylith::feassemble::ElasticityExplicitGravData2DLinear::_numQuadPts = 1;
+
+const char* pylith::feassemble::ElasticityExplicitGravData2DLinear::_matType = "ElasticPlaneStrain";
+
+const char* pylith::feassemble::ElasticityExplicitGravData2DLinear::_matDBFilename = "data/elasticplanestrain.spatialdb";
+
+const int pylith::feassemble::ElasticityExplicitGravData2DLinear::_matId = 0;
+
+const char* pylith::feassemble::ElasticityExplicitGravData2DLinear::_matLabel = "elastic strain 2-D";
+
+const double pylith::feassemble::ElasticityExplicitGravData2DLinear::_dt = 1.00000000e-02;
+
+const double pylith::feassemble::ElasticityExplicitGravData2DLinear::_gravityVec[] = {
+ 0.00000000e+00, -1.00000000e+08,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DLinear::_vertices[] = {
+ 2.00000000e-01, -4.00000000e-01,
+ 3.00000000e-01, 5.00000000e-01,
+ -1.00000000e+00, -2.00000000e-01,
+};
+
+const int pylith::feassemble::ElasticityExplicitGravData2DLinear::_cells[] = {
+0,1,2,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DLinear::_verticesRef[] = {
+ -1.00000000e+00, -1.00000000e+00,
+ 1.00000000e+00, -1.00000000e+00,
+ -1.00000000e+00, 1.00000000e+00,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DLinear::_quadPts[] = {
+ -3.33333333e-01, -3.33333333e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DLinear::_quadWts[] = {
+ 5.00000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DLinear::_basis[] = {
+ 3.33333333e-01, 3.33333333e-01,
+ 3.33333333e-01,};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DLinear::_basisDerivRef[] = {
+ -5.00000000e-01, -5.00000000e-01,
+ 5.00000000e-01, 0.00000000e+00,
+ 0.00000000e+00, 5.00000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DLinear::_fieldTIncr[] = {
+ 1.30000000e+00, -9.00000000e-01,
+ 1.40000000e+00, 1.50000000e+00,
+ 5.00000000e-01, -9.00000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DLinear::_fieldT[] = {
+ 1.60000000e+00, -8.00000000e-01,
+ 9.00000000e-01, 7.00000000e-01,
+ -2.00000000e-01, -1.10000000e+00,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DLinear::_fieldTmdt[] = {
+ 8.00000000e-01, 1.00000000e-01,
+ 5.00000000e-01, 3.00000000e-01,
+ -1.00000000e-01, -6.00000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DLinear::_valsResidual[] = {
+ -1.66787702e+10, 1.94562705e+10,
+ -2.30412532e+09, -3.86705051e+10,
+ 1.89841559e+10, -1.51619113e+10,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DLinear::_valsJacobian[] = {
+ 3.81944444e+05, 0.00000000e+00,
+ 3.81944444e+05, 0.00000000e+00,
+ 3.81944444e+05, 0.00000000e+00,
+ 0.00000000e+00, 3.81944444e+05,
+ 0.00000000e+00, 3.81944444e+05,
+ 0.00000000e+00, 3.81944444e+05,
+ 3.81944444e+05, 0.00000000e+00,
+ 3.81944444e+05, 0.00000000e+00,
+ 3.81944444e+05, 0.00000000e+00,
+ 0.00000000e+00, 3.81944444e+05,
+ 0.00000000e+00, 3.81944444e+05,
+ 0.00000000e+00, 3.81944444e+05,
+ 3.81944444e+05, 0.00000000e+00,
+ 3.81944444e+05, 0.00000000e+00,
+ 3.81944444e+05, 0.00000000e+00,
+ 0.00000000e+00, 3.81944444e+05,
+ 0.00000000e+00, 3.81944444e+05,
+ 0.00000000e+00, 3.81944444e+05,
+};
+
+pylith::feassemble::ElasticityExplicitGravData2DLinear::ElasticityExplicitGravData2DLinear(void)
+{ // constructor
+ spaceDim = _spaceDim;
+ cellDim = _cellDim;
+ numVertices = _numVertices;
+ numCells = _numCells;
+ numBasis = _numBasis;
+ numQuadPts = _numQuadPts;
+ matType = const_cast<char*>(_matType);
+ matDBFilename = const_cast<char*>(_matDBFilename);
+ matId = _matId;
+ matLabel = const_cast<char*>(_matLabel);
+ dt = _dt;
+ gravityVec = const_cast<double*>(_gravityVec);
+ vertices = const_cast<double*>(_vertices);
+ cells = const_cast<int*>(_cells);
+ verticesRef = const_cast<double*>(_verticesRef);
+ quadPts = const_cast<double*>(_quadPts);
+ quadWts = const_cast<double*>(_quadWts);
+ basis = const_cast<double*>(_basis);
+ basisDerivRef = const_cast<double*>(_basisDerivRef);
+ fieldTIncr = const_cast<double*>(_fieldTIncr);
+ fieldT = const_cast<double*>(_fieldT);
+ fieldTmdt = const_cast<double*>(_fieldTmdt);
+ valsResidual = const_cast<double*>(_valsResidual);
+ valsJacobian = const_cast<double*>(_valsJacobian);
+} // constructor
+
+pylith::feassemble::ElasticityExplicitGravData2DLinear::~ElasticityExplicitGravData2DLinear(void)
+{}
+
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData2DLinear.cc
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData2DLinear.hh
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData2DLinear.hh (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData2DLinear.hh 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,92 @@
+// -*- C++ -*-
+//
+// ======================================================================
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ======================================================================
+//
+
+// DO NOT EDIT THIS FILE
+// This file was generated from python application integratorelasticity.
+
+#if !defined(pylith_feassemble_elasticityexplicitgravdata2dlinear_hh)
+#define pylith_feassemble_elasticityexplicitgravdata2dlinear_hh
+
+#include "IntegratorData.hh"
+
+namespace pylith {
+ namespace feassemble {
+ class ElasticityExplicitGravData2DLinear;
+ } // pylith
+} // feassemble
+
+class pylith::feassemble::ElasticityExplicitGravData2DLinear : public IntegratorData
+{
+
+public:
+
+ /// Constructor
+ ElasticityExplicitGravData2DLinear(void);
+
+ /// Destructor
+ ~ElasticityExplicitGravData2DLinear(void);
+
+private:
+
+ static const int _spaceDim;
+
+ static const int _cellDim;
+
+ static const int _numVertices;
+
+ static const int _numCells;
+
+ static const int _numBasis;
+
+ static const int _numQuadPts;
+
+ static const char* _matType;
+
+ static const char* _matDBFilename;
+
+ static const int _matId;
+
+ static const char* _matLabel;
+
+ static const double _dt;
+
+ static const double _gravityVec[];
+
+ static const double _vertices[];
+
+ static const int _cells[];
+
+ static const double _verticesRef[];
+
+ static const double _quadPts[];
+
+ static const double _quadWts[];
+
+ static const double _basis[];
+
+ static const double _basisDerivRef[];
+
+ static const double _fieldTIncr[];
+
+ static const double _fieldT[];
+
+ static const double _fieldTmdt[];
+
+ static const double _valsResidual[];
+
+ static const double _valsJacobian[];
+
+};
+
+#endif // pylith_feassemble_elasticityexplicitgravdata2dlinear_hh
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData2DLinear.hh
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData2DQuadratic.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData2DQuadratic.cc (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData2DQuadratic.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,282 @@
+// -*- C++ -*-
+//
+// ======================================================================
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ======================================================================
+//
+
+// DO NOT EDIT THIS FILE
+// This file was generated from python application integratorelasticity.
+
+#include "ElasticityExplicitGravData2DQuadratic.hh"
+
+const int pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_spaceDim = 2;
+
+const int pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_cellDim = 2;
+
+const int pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_numVertices = 6;
+
+const int pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_numCells = 1;
+
+const int pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_numBasis = 6;
+
+const int pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_numQuadPts = 6;
+
+const char* pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_matType = "ElasticPlaneStrain";
+
+const char* pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_matDBFilename = "data/elasticplanestrain.spatialdb";
+
+const int pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_matId = 0;
+
+const char* pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_matLabel = "elastic strain 2-D";
+
+const double pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_dt = 1.00000000e-02;
+
+const double pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_gravityVec[] = {
+ 0.00000000e+00, -1.00000000e+08,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_vertices[] = {
+ -1.00000000e+00, -1.00000000e+00,
+ 1.00000000e+00, 2.00000000e-01,
+ -1.50000000e+00, 5.00000000e-01,
+ -2.50000000e-01, 3.50000000e-01,
+ -1.25000000e+00, -2.50000000e-01,
+ 0.00000000e+00, -4.00000000e-01,
+};
+
+const int pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_cells[] = {
+0,1,2,3,4,5,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_verticesRef[] = {
+ -1.00000000e+00, -1.00000000e+00,
+ 1.00000000e+00, -1.00000000e+00,
+ -1.00000000e+00, 1.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00,
+ -1.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, -1.00000000e+00,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_quadPts[] = {
+ -7.50000000e-01, -7.50000000e-01,
+ 7.50000000e-01, -7.50000000e-01,
+ -7.50000000e-01, 7.50000000e-01,
+ 0.00000000e+00, -7.50000000e-01,
+ -7.50000000e-01, 0.00000000e+00,
+ 2.50000000e-01, 2.50000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_quadWts[] = {
+ 3.33333333e-01, 3.33333333e-01, 3.33333333e-01, 3.33333333e-01, 3.33333333e-01, 3.33333333e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_basis[] = {
+ 3.75000000e-01, -9.37500000e-02,
+ -9.37500000e-02, 6.25000000e-02,
+ 3.75000000e-01, 3.75000000e-01,
+ 0.00000000e+00, 6.56250000e-01,
+ -9.37500000e-02, 4.37500000e-01,
+ -0.00000000e+00, -0.00000000e+00,
+ 0.00000000e+00, -9.37500000e-02,
+ 6.56250000e-01, 4.37500000e-01,
+ -0.00000000e+00, -0.00000000e+00,
+ -9.37500000e-02, 0.00000000e+00,
+ -9.37500000e-02, 2.50000000e-01,
+ 1.87500000e-01, 7.50000000e-01,
+ -9.37500000e-02, -9.37500000e-02,
+ 0.00000000e+00, 2.50000000e-01,
+ 7.50000000e-01, 1.87500000e-01,
+ 3.75000000e-01, 1.56250000e-01,
+ 1.56250000e-01, 1.56250000e+00,
+ -6.25000000e-01, -6.25000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_basisDerivRef[] = {
+ -1.00000000e+00, -1.00000000e+00,
+ -2.50000000e-01, 0.00000000e+00,
+ 0.00000000e+00, -2.50000000e-01,
+ 2.50000000e-01, 2.50000000e-01,
+ -2.50000000e-01, 1.25000000e+00,
+ 1.25000000e+00, -2.50000000e-01,
+ 5.00000000e-01, 5.00000000e-01,
+ 1.25000000e+00, 0.00000000e+00,
+ 0.00000000e+00, -2.50000000e-01,
+ 2.50000000e-01, 1.75000000e+00,
+ -2.50000000e-01, -2.50000000e-01,
+ -1.75000000e+00, -1.75000000e+00,
+ 5.00000000e-01, 5.00000000e-01,
+ -2.50000000e-01, 0.00000000e+00,
+ 0.00000000e+00, 1.25000000e+00,
+ 1.75000000e+00, 2.50000000e-01,
+ -1.75000000e+00, -1.75000000e+00,
+ -2.50000000e-01, -2.50000000e-01,
+ -2.50000000e-01, -2.50000000e-01,
+ 5.00000000e-01, 0.00000000e+00,
+ 0.00000000e+00, -2.50000000e-01,
+ 2.50000000e-01, 1.00000000e+00,
+ -2.50000000e-01, 5.00000000e-01,
+ -2.50000000e-01, -1.00000000e+00,
+ -2.50000000e-01, -2.50000000e-01,
+ -2.50000000e-01, 0.00000000e+00,
+ 0.00000000e+00, 5.00000000e-01,
+ 1.00000000e+00, 2.50000000e-01,
+ -1.00000000e+00, -2.50000000e-01,
+ 5.00000000e-01, -2.50000000e-01,
+ 1.00000000e+00, 1.00000000e+00,
+ 7.50000000e-01, 0.00000000e+00,
+ 0.00000000e+00, 7.50000000e-01,
+ 1.25000000e+00, 1.25000000e+00,
+ -1.25000000e+00, -1.75000000e+00,
+ -1.75000000e+00, -1.25000000e+00,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_fieldTIncr[] = {
+ -4.00000000e-01, -6.00000000e-01,
+ 7.00000000e-01, 8.00000000e-01,
+ 0.00000000e+00, 2.00000000e-01,
+ -5.00000000e-01, -4.00000000e-01,
+ 3.00000000e-01, 9.00000000e-01,
+ -3.00000000e-01, -9.00000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_fieldT[] = {
+ -3.00000000e-01, -4.00000000e-01,
+ 5.00000000e-01, 6.00000000e-01,
+ 0.00000000e+00, 1.00000000e-01,
+ -2.00000000e-01, -3.00000000e-01,
+ 2.00000000e-01, 3.00000000e-01,
+ -1.00000000e-01, -2.00000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_fieldTmdt[] = {
+ -2.00000000e-01, -3.00000000e-01,
+ 3.00000000e-01, 4.00000000e-01,
+ 0.00000000e+00, -1.00000000e-01,
+ -3.00000000e-01, -2.00000000e-01,
+ 1.00000000e-01, 4.00000000e-01,
+ -2.00000000e-01, -6.00000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_valsResidual[] = {
+ 1.06595078e+10, 4.79704396e+10,
+ -2.93254295e+10, -4.45162488e+10,
+ 4.90758945e+09, 2.18761966e+10,
+ 6.30081264e+10, -1.22877719e+11,
+ -6.40775573e+10, -2.36370117e+11,
+ 1.48314194e+10, -1.16082082e+11,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData2DQuadratic::_valsJacobian[] = {
+ 2.24121094e+06, 0.00000000e+00,
+ 2.41699219e+05, 0.00000000e+00,
+ 2.41699219e+05, 0.00000000e+00,
+ 4.21875000e+06, 0.00000000e+00,
+ -1.36230469e+06, 0.00000000e+00,
+ -1.36230469e+06, 0.00000000e+00,
+ 0.00000000e+00, 2.24121094e+06,
+ 0.00000000e+00, 2.41699219e+05,
+ 0.00000000e+00, 2.41699219e+05,
+ 0.00000000e+00, 4.21875000e+06,
+ 0.00000000e+00, -1.36230469e+06,
+ 0.00000000e+00, -1.36230469e+06,
+ 2.41699219e+05, 0.00000000e+00,
+ 3.61083984e+06, 0.00000000e+00,
+ -6.73828125e+05, 0.00000000e+00,
+ 3.45703125e+06, 0.00000000e+00,
+ -1.52343750e+06, 0.00000000e+00,
+ -1.12792969e+06, 0.00000000e+00,
+ 0.00000000e+00, 2.41699219e+05,
+ 0.00000000e+00, 3.61083984e+06,
+ 0.00000000e+00, -6.73828125e+05,
+ 0.00000000e+00, 3.45703125e+06,
+ 0.00000000e+00, -1.52343750e+06,
+ 0.00000000e+00, -1.12792969e+06,
+ 2.41699219e+05, 0.00000000e+00,
+ -6.73828125e+05, 0.00000000e+00,
+ 3.61083984e+06, 0.00000000e+00,
+ 3.45703125e+06, 0.00000000e+00,
+ -1.12792969e+06, 0.00000000e+00,
+ -1.52343750e+06, 0.00000000e+00,
+ 0.00000000e+00, 2.41699219e+05,
+ 0.00000000e+00, -6.73828125e+05,
+ 0.00000000e+00, 3.61083984e+06,
+ 0.00000000e+00, 3.45703125e+06,
+ 0.00000000e+00, -1.12792969e+06,
+ 0.00000000e+00, -1.52343750e+06,
+ 4.21875000e+06, 0.00000000e+00,
+ 3.45703125e+06, 0.00000000e+00,
+ 3.45703125e+06, 0.00000000e+00,
+ 2.21484375e+07, 0.00000000e+00,
+ -5.39062500e+06, 0.00000000e+00,
+ -5.39062500e+06, 0.00000000e+00,
+ 0.00000000e+00, 4.21875000e+06,
+ 0.00000000e+00, 3.45703125e+06,
+ 0.00000000e+00, 3.45703125e+06,
+ 0.00000000e+00, 2.21484375e+07,
+ 0.00000000e+00, -5.39062500e+06,
+ 0.00000000e+00, -5.39062500e+06,
+ -1.36230469e+06, 0.00000000e+00,
+ -1.52343750e+06, 0.00000000e+00,
+ -1.12792969e+06, 0.00000000e+00,
+ -5.39062500e+06, 0.00000000e+00,
+ 8.46679688e+06, 0.00000000e+00,
+ 6.09375000e+06, 0.00000000e+00,
+ 0.00000000e+00, -1.36230469e+06,
+ 0.00000000e+00, -1.52343750e+06,
+ 0.00000000e+00, -1.12792969e+06,
+ 0.00000000e+00, -5.39062500e+06,
+ 0.00000000e+00, 8.46679688e+06,
+ 0.00000000e+00, 6.09375000e+06,
+ -1.36230469e+06, 0.00000000e+00,
+ -1.12792969e+06, 0.00000000e+00,
+ -1.52343750e+06, 0.00000000e+00,
+ -5.39062500e+06, 0.00000000e+00,
+ 6.09375000e+06, 0.00000000e+00,
+ 8.46679688e+06, 0.00000000e+00,
+ 0.00000000e+00, -1.36230469e+06,
+ 0.00000000e+00, -1.12792969e+06,
+ 0.00000000e+00, -1.52343750e+06,
+ 0.00000000e+00, -5.39062500e+06,
+ 0.00000000e+00, 6.09375000e+06,
+ 0.00000000e+00, 8.46679688e+06,
+};
+
+pylith::feassemble::ElasticityExplicitGravData2DQuadratic::ElasticityExplicitGravData2DQuadratic(void)
+{ // constructor
+ spaceDim = _spaceDim;
+ cellDim = _cellDim;
+ numVertices = _numVertices;
+ numCells = _numCells;
+ numBasis = _numBasis;
+ numQuadPts = _numQuadPts;
+ matType = const_cast<char*>(_matType);
+ matDBFilename = const_cast<char*>(_matDBFilename);
+ matId = _matId;
+ matLabel = const_cast<char*>(_matLabel);
+ dt = _dt;
+ gravityVec = const_cast<double*>(_gravityVec);
+ vertices = const_cast<double*>(_vertices);
+ cells = const_cast<int*>(_cells);
+ verticesRef = const_cast<double*>(_verticesRef);
+ quadPts = const_cast<double*>(_quadPts);
+ quadWts = const_cast<double*>(_quadWts);
+ basis = const_cast<double*>(_basis);
+ basisDerivRef = const_cast<double*>(_basisDerivRef);
+ fieldTIncr = const_cast<double*>(_fieldTIncr);
+ fieldT = const_cast<double*>(_fieldT);
+ fieldTmdt = const_cast<double*>(_fieldTmdt);
+ valsResidual = const_cast<double*>(_valsResidual);
+ valsJacobian = const_cast<double*>(_valsJacobian);
+} // constructor
+
+pylith::feassemble::ElasticityExplicitGravData2DQuadratic::~ElasticityExplicitGravData2DQuadratic(void)
+{}
+
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData2DQuadratic.cc
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData2DQuadratic.hh
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData2DQuadratic.hh (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData2DQuadratic.hh 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,92 @@
+// -*- C++ -*-
+//
+// ======================================================================
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ======================================================================
+//
+
+// DO NOT EDIT THIS FILE
+// This file was generated from python application integratorelasticity.
+
+#if !defined(pylith_feassemble_elasticityexplicitgravdata2dquadratic_hh)
+#define pylith_feassemble_elasticityexplicitgravdata2dquadratic_hh
+
+#include "IntegratorData.hh"
+
+namespace pylith {
+ namespace feassemble {
+ class ElasticityExplicitGravData2DQuadratic;
+ } // pylith
+} // feassemble
+
+class pylith::feassemble::ElasticityExplicitGravData2DQuadratic : public IntegratorData
+{
+
+public:
+
+ /// Constructor
+ ElasticityExplicitGravData2DQuadratic(void);
+
+ /// Destructor
+ ~ElasticityExplicitGravData2DQuadratic(void);
+
+private:
+
+ static const int _spaceDim;
+
+ static const int _cellDim;
+
+ static const int _numVertices;
+
+ static const int _numCells;
+
+ static const int _numBasis;
+
+ static const int _numQuadPts;
+
+ static const char* _matType;
+
+ static const char* _matDBFilename;
+
+ static const int _matId;
+
+ static const char* _matLabel;
+
+ static const double _dt;
+
+ static const double _gravityVec[];
+
+ static const double _vertices[];
+
+ static const int _cells[];
+
+ static const double _verticesRef[];
+
+ static const double _quadPts[];
+
+ static const double _quadWts[];
+
+ static const double _basis[];
+
+ static const double _basisDerivRef[];
+
+ static const double _fieldTIncr[];
+
+ static const double _fieldT[];
+
+ static const double _fieldTmdt[];
+
+ static const double _valsResidual[];
+
+ static const double _valsJacobian[];
+
+};
+
+#endif // pylith_feassemble_elasticityexplicitgravdata2dquadratic_hh
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData2DQuadratic.hh
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData3DLinear.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData3DLinear.cc (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData3DLinear.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,192 @@
+// -*- C++ -*-
+//
+// ======================================================================
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ======================================================================
+//
+
+// DO NOT EDIT THIS FILE
+// This file was generated from python application integratorelasticity.
+
+#include "ElasticityExplicitGravData3DLinear.hh"
+
+const int pylith::feassemble::ElasticityExplicitGravData3DLinear::_spaceDim = 3;
+
+const int pylith::feassemble::ElasticityExplicitGravData3DLinear::_cellDim = 3;
+
+const int pylith::feassemble::ElasticityExplicitGravData3DLinear::_numVertices = 4;
+
+const int pylith::feassemble::ElasticityExplicitGravData3DLinear::_numCells = 1;
+
+const int pylith::feassemble::ElasticityExplicitGravData3DLinear::_numBasis = 4;
+
+const int pylith::feassemble::ElasticityExplicitGravData3DLinear::_numQuadPts = 1;
+
+const char* pylith::feassemble::ElasticityExplicitGravData3DLinear::_matType = "ElasticIsotropic3D";
+
+const char* pylith::feassemble::ElasticityExplicitGravData3DLinear::_matDBFilename = "data/elasticisotropic3d.spatialdb";
+
+const int pylith::feassemble::ElasticityExplicitGravData3DLinear::_matId = 0;
+
+const char* pylith::feassemble::ElasticityExplicitGravData3DLinear::_matLabel = "elastic isotropic 3-D";
+
+const double pylith::feassemble::ElasticityExplicitGravData3DLinear::_dt = 1.00000000e-02;
+
+const double pylith::feassemble::ElasticityExplicitGravData3DLinear::_gravityVec[] = {
+ 0.00000000e+00, 0.00000000e+00, -1.00000000e+08,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DLinear::_vertices[] = {
+ -5.00000000e-01, -1.00000000e+00, -5.00000000e-01,
+ 2.00000000e+00, -5.00000000e-01, -4.00000000e-01,
+ 1.00000000e+00, -1.00000000e-01, -3.00000000e-01,
+ -2.00000000e-01, 5.00000000e-01, 2.00000000e+00,
+};
+
+const int pylith::feassemble::ElasticityExplicitGravData3DLinear::_cells[] = {
+0,1,2,3,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DLinear::_verticesRef[] = {
+ -1.00000000e+00, -1.00000000e+00, -1.00000000e+00,
+ 1.00000000e+00, -1.00000000e+00, -1.00000000e+00,
+ -1.00000000e+00, 1.00000000e+00, -1.00000000e+00,
+ -1.00000000e+00, -1.00000000e+00, 1.00000000e+00,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DLinear::_quadPts[] = {
+ -5.00000000e-01, -5.00000000e-01, -5.00000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DLinear::_quadWts[] = {
+ 1.33333333e+00,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DLinear::_basis[] = {
+ 2.50000000e-01, 2.50000000e-01, 2.50000000e-01,
+ 2.50000000e-01,};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DLinear::_basisDerivRef[] = {
+ -5.00000000e-01, -5.00000000e-01, -5.00000000e-01,
+ 5.00000000e-01, 0.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, 5.00000000e-01, 0.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00, 5.00000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DLinear::_fieldTIncr[] = {
+ 3.00000000e-01, 2.00000000e-01, -5.00000000e-01,
+ -3.00000000e-01, -4.00000000e-01, -6.00000000e-01,
+ 2.00000000e-01, 6.00000000e-01, 3.00000000e-01,
+ -6.00000000e-01, -1.00000000e-01, -3.00000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DLinear::_fieldT[] = {
+ 8.00000000e-01, 1.00000000e-01, -6.00000000e-01,
+ -1.00000000e-01, -2.00000000e-01, -5.00000000e-01,
+ 1.00000000e-01, 7.00000000e-01, 2.00000000e-01,
+ -5.00000000e-01, -0.00000000e+00, -2.00000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DLinear::_fieldTmdt[] = {
+ 1.00000000e-01, 1.00000000e-01, -3.00000000e-01,
+ -2.00000000e-01, -1.00000000e-01, -5.00000000e-01,
+ 2.00000000e-01, 4.00000000e-01, 1.00000000e-01,
+ -4.00000000e-01, -1.00000000e-01, -1.00000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DLinear::_valsResidual[] = {
+ -6.07582772e+09, 3.62539077e+10, -3.04285979e+10,
+ -4.00897518e+09, 6.66855682e+10, -1.17054566e+10,
+ 6.67380053e+09, -1.05592747e+11, -6.51248804e+10,
+ 3.41301986e+09, 2.65427945e+09, -2.72420737e+10,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DLinear::_valsJacobian[] = {
+ 8.40625000e+05, 0.00000000e+00, 0.00000000e+00,
+ 8.40625000e+05, 0.00000000e+00, 0.00000000e+00,
+ 8.40625000e+05, 0.00000000e+00, 0.00000000e+00,
+ 8.40625000e+05, 0.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, 8.40625000e+05, 0.00000000e+00,
+ 0.00000000e+00, 8.40625000e+05, 0.00000000e+00,
+ 0.00000000e+00, 8.40625000e+05, 0.00000000e+00,
+ 0.00000000e+00, 8.40625000e+05, 0.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00, 8.40625000e+05,
+ 0.00000000e+00, 0.00000000e+00, 8.40625000e+05,
+ 0.00000000e+00, 0.00000000e+00, 8.40625000e+05,
+ 0.00000000e+00, 0.00000000e+00, 8.40625000e+05,
+ 8.40625000e+05, 0.00000000e+00, 0.00000000e+00,
+ 8.40625000e+05, 0.00000000e+00, 0.00000000e+00,
+ 8.40625000e+05, 0.00000000e+00, 0.00000000e+00,
+ 8.40625000e+05, 0.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, 8.40625000e+05, 0.00000000e+00,
+ 0.00000000e+00, 8.40625000e+05, 0.00000000e+00,
+ 0.00000000e+00, 8.40625000e+05, 0.00000000e+00,
+ 0.00000000e+00, 8.40625000e+05, 0.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00, 8.40625000e+05,
+ 0.00000000e+00, 0.00000000e+00, 8.40625000e+05,
+ 0.00000000e+00, 0.00000000e+00, 8.40625000e+05,
+ 0.00000000e+00, 0.00000000e+00, 8.40625000e+05,
+ 8.40625000e+05, 0.00000000e+00, 0.00000000e+00,
+ 8.40625000e+05, 0.00000000e+00, 0.00000000e+00,
+ 8.40625000e+05, 0.00000000e+00, 0.00000000e+00,
+ 8.40625000e+05, 0.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, 8.40625000e+05, 0.00000000e+00,
+ 0.00000000e+00, 8.40625000e+05, 0.00000000e+00,
+ 0.00000000e+00, 8.40625000e+05, 0.00000000e+00,
+ 0.00000000e+00, 8.40625000e+05, 0.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00, 8.40625000e+05,
+ 0.00000000e+00, 0.00000000e+00, 8.40625000e+05,
+ 0.00000000e+00, 0.00000000e+00, 8.40625000e+05,
+ 0.00000000e+00, 0.00000000e+00, 8.40625000e+05,
+ 8.40625000e+05, 0.00000000e+00, 0.00000000e+00,
+ 8.40625000e+05, 0.00000000e+00, 0.00000000e+00,
+ 8.40625000e+05, 0.00000000e+00, 0.00000000e+00,
+ 8.40625000e+05, 0.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, 8.40625000e+05, 0.00000000e+00,
+ 0.00000000e+00, 8.40625000e+05, 0.00000000e+00,
+ 0.00000000e+00, 8.40625000e+05, 0.00000000e+00,
+ 0.00000000e+00, 8.40625000e+05, 0.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00, 8.40625000e+05,
+ 0.00000000e+00, 0.00000000e+00, 8.40625000e+05,
+ 0.00000000e+00, 0.00000000e+00, 8.40625000e+05,
+ 0.00000000e+00, 0.00000000e+00, 8.40625000e+05,
+};
+
+pylith::feassemble::ElasticityExplicitGravData3DLinear::ElasticityExplicitGravData3DLinear(void)
+{ // constructor
+ spaceDim = _spaceDim;
+ cellDim = _cellDim;
+ numVertices = _numVertices;
+ numCells = _numCells;
+ numBasis = _numBasis;
+ numQuadPts = _numQuadPts;
+ matType = const_cast<char*>(_matType);
+ matDBFilename = const_cast<char*>(_matDBFilename);
+ matId = _matId;
+ matLabel = const_cast<char*>(_matLabel);
+ dt = _dt;
+ gravityVec = const_cast<double*>(_gravityVec);
+ vertices = const_cast<double*>(_vertices);
+ cells = const_cast<int*>(_cells);
+ verticesRef = const_cast<double*>(_verticesRef);
+ quadPts = const_cast<double*>(_quadPts);
+ quadWts = const_cast<double*>(_quadWts);
+ basis = const_cast<double*>(_basis);
+ basisDerivRef = const_cast<double*>(_basisDerivRef);
+ fieldTIncr = const_cast<double*>(_fieldTIncr);
+ fieldT = const_cast<double*>(_fieldT);
+ fieldTmdt = const_cast<double*>(_fieldTmdt);
+ valsResidual = const_cast<double*>(_valsResidual);
+ valsJacobian = const_cast<double*>(_valsJacobian);
+} // constructor
+
+pylith::feassemble::ElasticityExplicitGravData3DLinear::~ElasticityExplicitGravData3DLinear(void)
+{}
+
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData3DLinear.cc
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData3DLinear.hh
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData3DLinear.hh (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData3DLinear.hh 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,92 @@
+// -*- C++ -*-
+//
+// ======================================================================
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ======================================================================
+//
+
+// DO NOT EDIT THIS FILE
+// This file was generated from python application integratorelasticity.
+
+#if !defined(pylith_feassemble_elasticityexplicitgravdata3dlinear_hh)
+#define pylith_feassemble_elasticityexplicitgravdata3dlinear_hh
+
+#include "IntegratorData.hh"
+
+namespace pylith {
+ namespace feassemble {
+ class ElasticityExplicitGravData3DLinear;
+ } // pylith
+} // feassemble
+
+class pylith::feassemble::ElasticityExplicitGravData3DLinear : public IntegratorData
+{
+
+public:
+
+ /// Constructor
+ ElasticityExplicitGravData3DLinear(void);
+
+ /// Destructor
+ ~ElasticityExplicitGravData3DLinear(void);
+
+private:
+
+ static const int _spaceDim;
+
+ static const int _cellDim;
+
+ static const int _numVertices;
+
+ static const int _numCells;
+
+ static const int _numBasis;
+
+ static const int _numQuadPts;
+
+ static const char* _matType;
+
+ static const char* _matDBFilename;
+
+ static const int _matId;
+
+ static const char* _matLabel;
+
+ static const double _dt;
+
+ static const double _gravityVec[];
+
+ static const double _vertices[];
+
+ static const int _cells[];
+
+ static const double _verticesRef[];
+
+ static const double _quadPts[];
+
+ static const double _quadWts[];
+
+ static const double _basis[];
+
+ static const double _basisDerivRef[];
+
+ static const double _fieldTIncr[];
+
+ static const double _fieldT[];
+
+ static const double _fieldTmdt[];
+
+ static const double _valsResidual[];
+
+ static const double _valsJacobian[];
+
+};
+
+#endif // pylith_feassemble_elasticityexplicitgravdata3dlinear_hh
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData3DLinear.hh
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData3DQuadratic.cc
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData3DQuadratic.cc (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData3DQuadratic.cc 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,531 @@
+// -*- C++ -*-
+//
+// ======================================================================
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ======================================================================
+//
+
+// DO NOT EDIT THIS FILE
+// This file was generated from python application integratorelasticity.
+
+#include "ElasticityExplicitGravData3DQuadratic.hh"
+
+const int pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_spaceDim = 3;
+
+const int pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_cellDim = 3;
+
+const int pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_numVertices = 10;
+
+const int pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_numCells = 1;
+
+const int pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_numBasis = 10;
+
+const int pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_numQuadPts = 4;
+
+const char* pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_matType = "ElasticIsotropic3D";
+
+const char* pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_matDBFilename = "data/elasticisotropic3d.spatialdb";
+
+const int pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_matId = 0;
+
+const char* pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_matLabel = "elastic isotropic 3-D";
+
+const double pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_dt = 1.00000000e-02;
+
+const double pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_gravityVec[] = {
+ 0.00000000e+00, 0.00000000e+00, -1.00000000e+08,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_vertices[] = {
+ -5.00000000e-01, -2.00000000e+00, -1.00000000e+00,
+ 2.00000000e+00, -2.00000000e+00, -5.00000000e-01,
+ 1.00000000e+00, 1.00000000e+00, 0.00000000e+00,
+ 2.00000000e-01, 5.00000000e-01, 2.00000000e+00,
+ 1.50000000e+00, -5.00000000e-01, -2.50000000e-01,
+ 2.50000000e-01, -5.00000000e-01, -5.00000000e-01,
+ 7.50000000e-01, -2.00000000e+00, -7.50000000e-01,
+ -1.50000000e-01, -7.50000000e-01, 5.00000000e-01,
+ 1.10000000e+00, -7.50000000e-01, 7.50000000e-01,
+ 6.00000000e-01, 7.50000000e-01, 1.00000000e+00,
+};
+
+const int pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_cells[] = {
+0,1,2,3,4,5,6,7,8,9,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_verticesRef[] = {
+ -1.00000000e+00, -1.00000000e+00, -1.00000000e+00,
+ 1.00000000e+00, -1.00000000e+00, -1.00000000e+00,
+ -1.00000000e+00, 1.00000000e+00, -1.00000000e+00,
+ -1.00000000e+00, -1.00000000e+00, 1.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00, -1.00000000e+00,
+ -1.00000000e+00, 0.00000000e+00, -1.00000000e+00,
+ 0.00000000e+00, -1.00000000e+00, -1.00000000e+00,
+ -1.00000000e+00, -1.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, -1.00000000e+00, 0.00000000e+00,
+ -1.00000000e+00, 0.00000000e+00, 0.00000000e+00,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_quadPts[] = {
+ -8.00000000e-01, -8.00000000e-01, -8.00000000e-01,
+ 5.00000000e-01, -8.00000000e-01, -8.00000000e-01,
+ -8.00000000e-01, 5.00000000e-01, -8.00000000e-01,
+ -8.00000000e-01, -8.00000000e-01, 5.00000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_quadWts[] = {
+ 3.33333333e-01, 3.33333333e-01, 3.33333333e-01, 3.33333333e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_basis[] = {
+ 2.80000000e-01, -8.00000000e-02, -8.00000000e-02,
+ -8.00000000e-02, 4.00000000e-02, 2.80000000e-01,
+ 2.80000000e-01, 2.80000000e-01, 4.00000000e-02,
+ 4.00000000e-02, -4.50000000e-02, 3.75000000e-01,
+ -8.00000000e-02, -8.00000000e-02, 3.00000000e-01,
+ 2.00000000e-02, 1.50000000e-01, 2.00000000e-02,
+ 3.00000000e-01, 4.00000000e-02, -4.50000000e-02,
+ -8.00000000e-02, 3.75000000e-01, -8.00000000e-02,
+ 3.00000000e-01, 1.50000000e-01, 2.00000000e-02,
+ 2.00000000e-02, 4.00000000e-02, 3.00000000e-01,
+ -4.50000000e-02, -8.00000000e-02, -8.00000000e-02,
+ 3.75000000e-01, 4.00000000e-02, 2.00000000e-02,
+ 2.00000000e-02, 1.50000000e-01, 3.00000000e-01,
+ 3.00000000e-01,};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_basisDerivRef[] = {
+ -9.00000000e-01, -9.00000000e-01, -9.00000000e-01,
+ -3.00000000e-01, 0.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, -3.00000000e-01, 0.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00, -3.00000000e-01,
+ 2.00000000e-01, 2.00000000e-01, 0.00000000e+00,
+ -2.00000000e-01, 1.20000000e+00, -2.00000000e-01,
+ 1.20000000e+00, -2.00000000e-01, -2.00000000e-01,
+ -2.00000000e-01, -2.00000000e-01, 1.20000000e+00,
+ 2.00000000e-01, 0.00000000e+00, 2.00000000e-01,
+ 0.00000000e+00, 2.00000000e-01, 2.00000000e-01,
+ 4.00000000e-01, 4.00000000e-01, 4.00000000e-01,
+ 1.00000000e+00, 0.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, -3.00000000e-01, 0.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00, -3.00000000e-01,
+ 2.00000000e-01, 1.50000000e+00, 0.00000000e+00,
+ -2.00000000e-01, -1.00000000e-01, -2.00000000e-01,
+ -1.40000000e+00, -1.50000000e+00, -1.50000000e+00,
+ -2.00000000e-01, -2.00000000e-01, -1.00000000e-01,
+ 2.00000000e-01, 0.00000000e+00, 1.50000000e+00,
+ 0.00000000e+00, 2.00000000e-01, 2.00000000e-01,
+ 4.00000000e-01, 4.00000000e-01, 4.00000000e-01,
+ -3.00000000e-01, 0.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, 1.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00, -3.00000000e-01,
+ 1.50000000e+00, 2.00000000e-01, 0.00000000e+00,
+ -1.50000000e+00, -1.40000000e+00, -1.50000000e+00,
+ -1.00000000e-01, -2.00000000e-01, -2.00000000e-01,
+ -2.00000000e-01, -2.00000000e-01, -1.00000000e-01,
+ 2.00000000e-01, 0.00000000e+00, 2.00000000e-01,
+ 0.00000000e+00, 2.00000000e-01, 1.50000000e+00,
+ 4.00000000e-01, 4.00000000e-01, 4.00000000e-01,
+ -3.00000000e-01, 0.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, -3.00000000e-01, 0.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00, 1.00000000e+00,
+ 2.00000000e-01, 2.00000000e-01, 0.00000000e+00,
+ -2.00000000e-01, -1.00000000e-01, -2.00000000e-01,
+ -1.00000000e-01, -2.00000000e-01, -2.00000000e-01,
+ -1.50000000e+00, -1.50000000e+00, -1.40000000e+00,
+ 1.50000000e+00, 0.00000000e+00, 2.00000000e-01,
+ 0.00000000e+00, 1.50000000e+00, 2.00000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_fieldTIncr[] = {
+ 3.00000000e-01, -4.00000000e-01, -4.00000000e-01,
+ -6.00000000e-01, 8.00000000e-01, 2.00000000e-01,
+ 5.00000000e-01, 5.00000000e-01, 7.00000000e-01,
+ -7.00000000e-01, -5.00000000e-01, -7.00000000e-01,
+ -6.00000000e-01, -3.00000000e-01, 8.00000000e-01,
+ -4.00000000e-01, -8.00000000e-01, -5.00000000e-01,
+ 7.00000000e-01, 8.00000000e-01, -5.00000000e-01,
+ -5.00000000e-01, -5.00000000e-01, -7.00000000e-01,
+ -3.00000000e-01, -9.00000000e-01, 8.00000000e-01,
+ -1.00000000e-01, 5.00000000e-01, -9.00000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_fieldT[] = {
+ 1.00000000e-01, -2.00000000e-01, -6.00000000e-01,
+ -3.00000000e-01, 4.00000000e-01, 9.00000000e-01,
+ 6.00000000e-01, 8.00000000e-01, 5.00000000e-01,
+ -8.00000000e-01, -6.00000000e-01, -8.00000000e-01,
+ -0.00000000e+00, -2.00000000e-01, 6.00000000e-01,
+ -4.00000000e-01, -7.00000000e-01, -2.00000000e-01,
+ 7.00000000e-01, 6.00000000e-01, -1.00000000e-01,
+ -4.00000000e-01, -3.00000000e-01, -3.00000000e-01,
+ -7.00000000e-01, -6.00000000e-01, 1.00000000e-01,
+ -9.00000000e-01, 3.00000000e-01, -8.00000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_fieldTmdt[] = {
+ 2.00000000e-01, -3.00000000e-01, -1.00000000e-01,
+ -4.00000000e-01, 2.00000000e-01, 3.00000000e-01,
+ -5.00000000e-01, 2.00000000e-01, 2.00000000e-01,
+ -3.00000000e-01, -8.00000000e-01, -3.00000000e-01,
+ -5.00000000e-01, -9.00000000e-01, 4.00000000e-01,
+ -3.00000000e-01, -6.00000000e-01, -8.00000000e-01,
+ 9.00000000e-01, 5.00000000e-01, -2.00000000e-01,
+ -7.00000000e-01, -2.00000000e-01, -9.00000000e-01,
+ -5.00000000e-01, -8.00000000e-01, 4.00000000e-01,
+ -4.00000000e-01, 5.00000000e-01, -7.00000000e-01,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_valsResidual[] = {
+ 2.17446982e+10, -9.08226289e+09, -5.09583471e+09,
+ 6.22338118e+10, -2.47879968e+10, -1.65884560e+10,
+ -4.60076154e+10, -5.34101030e+10, -5.55279242e+10,
+ -1.00555635e+10, 4.19850265e+10, -6.27466022e+10,
+ -8.43197122e+09, 6.33093784e+10, -1.24366385e+11,
+ 6.62258554e+10, 1.09239192e+11, -1.20877067e+11,
+ -9.59444688e+10, -5.29496681e+10, -1.08716467e+11,
+ -1.25318752e+10, -6.06765547e+10, -7.32262485e+10,
+ -3.24446843e+10, 5.14096428e+10, -1.13619631e+11,
+ 5.52108169e+10, -6.50263642e+10, -3.06855968e+10,
+};
+
+const double pylith::feassemble::ElasticityExplicitGravData3DQuadratic::_valsJacobian[] = {
+ 1.50251107e+06, 0.00000000e+00, 0.00000000e+00,
+ -5.70500651e+05, 0.00000000e+00, 0.00000000e+00,
+ -5.70500651e+05, 0.00000000e+00, 0.00000000e+00,
+ -5.70500651e+05, 0.00000000e+00, 0.00000000e+00,
+ -3.13041667e+05, 0.00000000e+00, 0.00000000e+00,
+ 1.24238411e+06, 0.00000000e+00, 0.00000000e+00,
+ 1.24238411e+06, 0.00000000e+00, 0.00000000e+00,
+ 1.24238411e+06, 0.00000000e+00, 0.00000000e+00,
+ -3.13041667e+05, 0.00000000e+00, 0.00000000e+00,
+ -3.13041667e+05, 0.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, 1.50251107e+06, 0.00000000e+00,
+ 0.00000000e+00, -5.70500651e+05, 0.00000000e+00,
+ 0.00000000e+00, -5.70500651e+05, 0.00000000e+00,
+ 0.00000000e+00, -5.70500651e+05, 0.00000000e+00,
+ 0.00000000e+00, -3.13041667e+05, 0.00000000e+00,
+ 0.00000000e+00, 1.24238411e+06, 0.00000000e+00,
+ 0.00000000e+00, 1.24238411e+06, 0.00000000e+00,
+ 0.00000000e+00, 1.24238411e+06, 0.00000000e+00,
+ 0.00000000e+00, -3.13041667e+05, 0.00000000e+00,
+ 0.00000000e+00, -3.13041667e+05, 0.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00, 1.50251107e+06,
+ 0.00000000e+00, 0.00000000e+00, -5.70500651e+05,
+ 0.00000000e+00, 0.00000000e+00, -5.70500651e+05,
+ 0.00000000e+00, 0.00000000e+00, -5.70500651e+05,
+ 0.00000000e+00, 0.00000000e+00, -3.13041667e+05,
+ 0.00000000e+00, 0.00000000e+00, 1.24238411e+06,
+ 0.00000000e+00, 0.00000000e+00, 1.24238411e+06,
+ 0.00000000e+00, 0.00000000e+00, 1.24238411e+06,
+ 0.00000000e+00, 0.00000000e+00, -3.13041667e+05,
+ 0.00000000e+00, 0.00000000e+00, -3.13041667e+05,
+ -5.70500651e+05, 0.00000000e+00, 0.00000000e+00,
+ 2.84272070e+06, 0.00000000e+00, 0.00000000e+00,
+ -8.39520833e+05, 0.00000000e+00, 0.00000000e+00,
+ -8.39520833e+05, 0.00000000e+00, 0.00000000e+00,
+ 1.46026823e+06, 0.00000000e+00, 0.00000000e+00,
+ -5.06914062e+05, 0.00000000e+00, 0.00000000e+00,
+ 5.45154948e+05, 0.00000000e+00, 0.00000000e+00,
+ -5.06914062e+05, 0.00000000e+00, 0.00000000e+00,
+ 1.46026823e+06, 0.00000000e+00, 0.00000000e+00,
+ -6.43869792e+05, 0.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, -5.70500651e+05, 0.00000000e+00,
+ 0.00000000e+00, 2.84272070e+06, 0.00000000e+00,
+ 0.00000000e+00, -8.39520833e+05, 0.00000000e+00,
+ 0.00000000e+00, -8.39520833e+05, 0.00000000e+00,
+ 0.00000000e+00, 1.46026823e+06, 0.00000000e+00,
+ 0.00000000e+00, -5.06914062e+05, 0.00000000e+00,
+ 0.00000000e+00, 5.45154948e+05, 0.00000000e+00,
+ 0.00000000e+00, -5.06914062e+05, 0.00000000e+00,
+ 0.00000000e+00, 1.46026823e+06, 0.00000000e+00,
+ 0.00000000e+00, -6.43869792e+05, 0.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00, -5.70500651e+05,
+ 0.00000000e+00, 0.00000000e+00, 2.84272070e+06,
+ 0.00000000e+00, 0.00000000e+00, -8.39520833e+05,
+ 0.00000000e+00, 0.00000000e+00, -8.39520833e+05,
+ 0.00000000e+00, 0.00000000e+00, 1.46026823e+06,
+ 0.00000000e+00, 0.00000000e+00, -5.06914062e+05,
+ 0.00000000e+00, 0.00000000e+00, 5.45154948e+05,
+ 0.00000000e+00, 0.00000000e+00, -5.06914062e+05,
+ 0.00000000e+00, 0.00000000e+00, 1.46026823e+06,
+ 0.00000000e+00, 0.00000000e+00, -6.43869792e+05,
+ -5.70500651e+05, 0.00000000e+00, 0.00000000e+00,
+ -8.39520833e+05, 0.00000000e+00, 0.00000000e+00,
+ 2.84272070e+06, 0.00000000e+00, 0.00000000e+00,
+ -8.39520833e+05, 0.00000000e+00, 0.00000000e+00,
+ 1.46026823e+06, 0.00000000e+00, 0.00000000e+00,
+ 5.45154948e+05, 0.00000000e+00, 0.00000000e+00,
+ -5.06914062e+05, 0.00000000e+00, 0.00000000e+00,
+ -5.06914062e+05, 0.00000000e+00, 0.00000000e+00,
+ -6.43869792e+05, 0.00000000e+00, 0.00000000e+00,
+ 1.46026823e+06, 0.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, -5.70500651e+05, 0.00000000e+00,
+ 0.00000000e+00, -8.39520833e+05, 0.00000000e+00,
+ 0.00000000e+00, 2.84272070e+06, 0.00000000e+00,
+ 0.00000000e+00, -8.39520833e+05, 0.00000000e+00,
+ 0.00000000e+00, 1.46026823e+06, 0.00000000e+00,
+ 0.00000000e+00, 5.45154948e+05, 0.00000000e+00,
+ 0.00000000e+00, -5.06914062e+05, 0.00000000e+00,
+ 0.00000000e+00, -5.06914062e+05, 0.00000000e+00,
+ 0.00000000e+00, -6.43869792e+05, 0.00000000e+00,
+ 0.00000000e+00, 1.46026823e+06, 0.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00, -5.70500651e+05,
+ 0.00000000e+00, 0.00000000e+00, -8.39520833e+05,
+ 0.00000000e+00, 0.00000000e+00, 2.84272070e+06,
+ 0.00000000e+00, 0.00000000e+00, -8.39520833e+05,
+ 0.00000000e+00, 0.00000000e+00, 1.46026823e+06,
+ 0.00000000e+00, 0.00000000e+00, 5.45154948e+05,
+ 0.00000000e+00, 0.00000000e+00, -5.06914062e+05,
+ 0.00000000e+00, 0.00000000e+00, -5.06914062e+05,
+ 0.00000000e+00, 0.00000000e+00, -6.43869792e+05,
+ 0.00000000e+00, 0.00000000e+00, 1.46026823e+06,
+ -5.70500651e+05, 0.00000000e+00, 0.00000000e+00,
+ -8.39520833e+05, 0.00000000e+00, 0.00000000e+00,
+ -8.39520833e+05, 0.00000000e+00, 0.00000000e+00,
+ 2.84272070e+06, 0.00000000e+00, 0.00000000e+00,
+ -6.43869792e+05, 0.00000000e+00, 0.00000000e+00,
+ -5.06914062e+05, 0.00000000e+00, 0.00000000e+00,
+ -5.06914062e+05, 0.00000000e+00, 0.00000000e+00,
+ 5.45154948e+05, 0.00000000e+00, 0.00000000e+00,
+ 1.46026823e+06, 0.00000000e+00, 0.00000000e+00,
+ 1.46026823e+06, 0.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, -5.70500651e+05, 0.00000000e+00,
+ 0.00000000e+00, -8.39520833e+05, 0.00000000e+00,
+ 0.00000000e+00, -8.39520833e+05, 0.00000000e+00,
+ 0.00000000e+00, 2.84272070e+06, 0.00000000e+00,
+ 0.00000000e+00, -6.43869792e+05, 0.00000000e+00,
+ 0.00000000e+00, -5.06914062e+05, 0.00000000e+00,
+ 0.00000000e+00, -5.06914062e+05, 0.00000000e+00,
+ 0.00000000e+00, 5.45154948e+05, 0.00000000e+00,
+ 0.00000000e+00, 1.46026823e+06, 0.00000000e+00,
+ 0.00000000e+00, 1.46026823e+06, 0.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00, -5.70500651e+05,
+ 0.00000000e+00, 0.00000000e+00, -8.39520833e+05,
+ 0.00000000e+00, 0.00000000e+00, -8.39520833e+05,
+ 0.00000000e+00, 0.00000000e+00, 2.84272070e+06,
+ 0.00000000e+00, 0.00000000e+00, -6.43869792e+05,
+ 0.00000000e+00, 0.00000000e+00, -5.06914062e+05,
+ 0.00000000e+00, 0.00000000e+00, -5.06914062e+05,
+ 0.00000000e+00, 0.00000000e+00, 5.45154948e+05,
+ 0.00000000e+00, 0.00000000e+00, 1.46026823e+06,
+ 0.00000000e+00, 0.00000000e+00, 1.46026823e+06,
+ -3.13041667e+05, 0.00000000e+00, 0.00000000e+00,
+ 1.46026823e+06, 0.00000000e+00, 0.00000000e+00,
+ 1.46026823e+06, 0.00000000e+00, 0.00000000e+00,
+ -6.43869792e+05, 0.00000000e+00, 0.00000000e+00,
+ 3.25847917e+06, 0.00000000e+00, 0.00000000e+00,
+ 1.12054688e+06, 0.00000000e+00, 0.00000000e+00,
+ 1.12054688e+06, 0.00000000e+00, 0.00000000e+00,
+ 5.19364583e+05, 0.00000000e+00, 0.00000000e+00,
+ 2.05611458e+06, 0.00000000e+00, 0.00000000e+00,
+ 2.05611458e+06, 0.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, -3.13041667e+05, 0.00000000e+00,
+ 0.00000000e+00, 1.46026823e+06, 0.00000000e+00,
+ 0.00000000e+00, 1.46026823e+06, 0.00000000e+00,
+ 0.00000000e+00, -6.43869792e+05, 0.00000000e+00,
+ 0.00000000e+00, 3.25847917e+06, 0.00000000e+00,
+ 0.00000000e+00, 1.12054688e+06, 0.00000000e+00,
+ 0.00000000e+00, 1.12054688e+06, 0.00000000e+00,
+ 0.00000000e+00, 5.19364583e+05, 0.00000000e+00,
+ 0.00000000e+00, 2.05611458e+06, 0.00000000e+00,
+ 0.00000000e+00, 2.05611458e+06, 0.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00, -3.13041667e+05,
+ 0.00000000e+00, 0.00000000e+00, 1.46026823e+06,
+ 0.00000000e+00, 0.00000000e+00, 1.46026823e+06,
+ 0.00000000e+00, 0.00000000e+00, -6.43869792e+05,
+ 0.00000000e+00, 0.00000000e+00, 3.25847917e+06,
+ 0.00000000e+00, 0.00000000e+00, 1.12054688e+06,
+ 0.00000000e+00, 0.00000000e+00, 1.12054688e+06,
+ 0.00000000e+00, 0.00000000e+00, 5.19364583e+05,
+ 0.00000000e+00, 0.00000000e+00, 2.05611458e+06,
+ 0.00000000e+00, 0.00000000e+00, 2.05611458e+06,
+ 1.24238411e+06, 0.00000000e+00, 0.00000000e+00,
+ -5.06914062e+05, 0.00000000e+00, 0.00000000e+00,
+ 5.45154948e+05, 0.00000000e+00, 0.00000000e+00,
+ -5.06914062e+05, 0.00000000e+00, 0.00000000e+00,
+ 1.12054688e+06, 0.00000000e+00, 0.00000000e+00,
+ 1.80888281e+06, 0.00000000e+00, 0.00000000e+00,
+ 1.50829167e+06, 0.00000000e+00, 0.00000000e+00,
+ 1.50829167e+06, 0.00000000e+00, 0.00000000e+00,
+ 5.19364583e+05, 0.00000000e+00, 0.00000000e+00,
+ 1.12054688e+06, 0.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, 1.24238411e+06, 0.00000000e+00,
+ 0.00000000e+00, -5.06914062e+05, 0.00000000e+00,
+ 0.00000000e+00, 5.45154948e+05, 0.00000000e+00,
+ 0.00000000e+00, -5.06914062e+05, 0.00000000e+00,
+ 0.00000000e+00, 1.12054688e+06, 0.00000000e+00,
+ 0.00000000e+00, 1.80888281e+06, 0.00000000e+00,
+ 0.00000000e+00, 1.50829167e+06, 0.00000000e+00,
+ 0.00000000e+00, 1.50829167e+06, 0.00000000e+00,
+ 0.00000000e+00, 5.19364583e+05, 0.00000000e+00,
+ 0.00000000e+00, 1.12054688e+06, 0.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00, 1.24238411e+06,
+ 0.00000000e+00, 0.00000000e+00, -5.06914062e+05,
+ 0.00000000e+00, 0.00000000e+00, 5.45154948e+05,
+ 0.00000000e+00, 0.00000000e+00, -5.06914062e+05,
+ 0.00000000e+00, 0.00000000e+00, 1.12054688e+06,
+ 0.00000000e+00, 0.00000000e+00, 1.80888281e+06,
+ 0.00000000e+00, 0.00000000e+00, 1.50829167e+06,
+ 0.00000000e+00, 0.00000000e+00, 1.50829167e+06,
+ 0.00000000e+00, 0.00000000e+00, 5.19364583e+05,
+ 0.00000000e+00, 0.00000000e+00, 1.12054688e+06,
+ 1.24238411e+06, 0.00000000e+00, 0.00000000e+00,
+ 5.45154948e+05, 0.00000000e+00, 0.00000000e+00,
+ -5.06914062e+05, 0.00000000e+00, 0.00000000e+00,
+ -5.06914062e+05, 0.00000000e+00, 0.00000000e+00,
+ 1.12054688e+06, 0.00000000e+00, 0.00000000e+00,
+ 1.50829167e+06, 0.00000000e+00, 0.00000000e+00,
+ 1.80888281e+06, 0.00000000e+00, 0.00000000e+00,
+ 1.50829167e+06, 0.00000000e+00, 0.00000000e+00,
+ 1.12054688e+06, 0.00000000e+00, 0.00000000e+00,
+ 5.19364583e+05, 0.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, 1.24238411e+06, 0.00000000e+00,
+ 0.00000000e+00, 5.45154948e+05, 0.00000000e+00,
+ 0.00000000e+00, -5.06914062e+05, 0.00000000e+00,
+ 0.00000000e+00, -5.06914062e+05, 0.00000000e+00,
+ 0.00000000e+00, 1.12054688e+06, 0.00000000e+00,
+ 0.00000000e+00, 1.50829167e+06, 0.00000000e+00,
+ 0.00000000e+00, 1.80888281e+06, 0.00000000e+00,
+ 0.00000000e+00, 1.50829167e+06, 0.00000000e+00,
+ 0.00000000e+00, 1.12054688e+06, 0.00000000e+00,
+ 0.00000000e+00, 5.19364583e+05, 0.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00, 1.24238411e+06,
+ 0.00000000e+00, 0.00000000e+00, 5.45154948e+05,
+ 0.00000000e+00, 0.00000000e+00, -5.06914062e+05,
+ 0.00000000e+00, 0.00000000e+00, -5.06914062e+05,
+ 0.00000000e+00, 0.00000000e+00, 1.12054688e+06,
+ 0.00000000e+00, 0.00000000e+00, 1.50829167e+06,
+ 0.00000000e+00, 0.00000000e+00, 1.80888281e+06,
+ 0.00000000e+00, 0.00000000e+00, 1.50829167e+06,
+ 0.00000000e+00, 0.00000000e+00, 1.12054688e+06,
+ 0.00000000e+00, 0.00000000e+00, 5.19364583e+05,
+ 1.24238411e+06, 0.00000000e+00, 0.00000000e+00,
+ -5.06914062e+05, 0.00000000e+00, 0.00000000e+00,
+ -5.06914062e+05, 0.00000000e+00, 0.00000000e+00,
+ 5.45154948e+05, 0.00000000e+00, 0.00000000e+00,
+ 5.19364583e+05, 0.00000000e+00, 0.00000000e+00,
+ 1.50829167e+06, 0.00000000e+00, 0.00000000e+00,
+ 1.50829167e+06, 0.00000000e+00, 0.00000000e+00,
+ 1.80888281e+06, 0.00000000e+00, 0.00000000e+00,
+ 1.12054688e+06, 0.00000000e+00, 0.00000000e+00,
+ 1.12054688e+06, 0.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, 1.24238411e+06, 0.00000000e+00,
+ 0.00000000e+00, -5.06914062e+05, 0.00000000e+00,
+ 0.00000000e+00, -5.06914062e+05, 0.00000000e+00,
+ 0.00000000e+00, 5.45154948e+05, 0.00000000e+00,
+ 0.00000000e+00, 5.19364583e+05, 0.00000000e+00,
+ 0.00000000e+00, 1.50829167e+06, 0.00000000e+00,
+ 0.00000000e+00, 1.50829167e+06, 0.00000000e+00,
+ 0.00000000e+00, 1.80888281e+06, 0.00000000e+00,
+ 0.00000000e+00, 1.12054688e+06, 0.00000000e+00,
+ 0.00000000e+00, 1.12054688e+06, 0.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00, 1.24238411e+06,
+ 0.00000000e+00, 0.00000000e+00, -5.06914062e+05,
+ 0.00000000e+00, 0.00000000e+00, -5.06914062e+05,
+ 0.00000000e+00, 0.00000000e+00, 5.45154948e+05,
+ 0.00000000e+00, 0.00000000e+00, 5.19364583e+05,
+ 0.00000000e+00, 0.00000000e+00, 1.50829167e+06,
+ 0.00000000e+00, 0.00000000e+00, 1.50829167e+06,
+ 0.00000000e+00, 0.00000000e+00, 1.80888281e+06,
+ 0.00000000e+00, 0.00000000e+00, 1.12054688e+06,
+ 0.00000000e+00, 0.00000000e+00, 1.12054688e+06,
+ -3.13041667e+05, 0.00000000e+00, 0.00000000e+00,
+ 1.46026823e+06, 0.00000000e+00, 0.00000000e+00,
+ -6.43869792e+05, 0.00000000e+00, 0.00000000e+00,
+ 1.46026823e+06, 0.00000000e+00, 0.00000000e+00,
+ 2.05611458e+06, 0.00000000e+00, 0.00000000e+00,
+ 5.19364583e+05, 0.00000000e+00, 0.00000000e+00,
+ 1.12054688e+06, 0.00000000e+00, 0.00000000e+00,
+ 1.12054688e+06, 0.00000000e+00, 0.00000000e+00,
+ 3.25847917e+06, 0.00000000e+00, 0.00000000e+00,
+ 2.05611458e+06, 0.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, -3.13041667e+05, 0.00000000e+00,
+ 0.00000000e+00, 1.46026823e+06, 0.00000000e+00,
+ 0.00000000e+00, -6.43869792e+05, 0.00000000e+00,
+ 0.00000000e+00, 1.46026823e+06, 0.00000000e+00,
+ 0.00000000e+00, 2.05611458e+06, 0.00000000e+00,
+ 0.00000000e+00, 5.19364583e+05, 0.00000000e+00,
+ 0.00000000e+00, 1.12054688e+06, 0.00000000e+00,
+ 0.00000000e+00, 1.12054688e+06, 0.00000000e+00,
+ 0.00000000e+00, 3.25847917e+06, 0.00000000e+00,
+ 0.00000000e+00, 2.05611458e+06, 0.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00, -3.13041667e+05,
+ 0.00000000e+00, 0.00000000e+00, 1.46026823e+06,
+ 0.00000000e+00, 0.00000000e+00, -6.43869792e+05,
+ 0.00000000e+00, 0.00000000e+00, 1.46026823e+06,
+ 0.00000000e+00, 0.00000000e+00, 2.05611458e+06,
+ 0.00000000e+00, 0.00000000e+00, 5.19364583e+05,
+ 0.00000000e+00, 0.00000000e+00, 1.12054688e+06,
+ 0.00000000e+00, 0.00000000e+00, 1.12054688e+06,
+ 0.00000000e+00, 0.00000000e+00, 3.25847917e+06,
+ 0.00000000e+00, 0.00000000e+00, 2.05611458e+06,
+ -3.13041667e+05, 0.00000000e+00, 0.00000000e+00,
+ -6.43869792e+05, 0.00000000e+00, 0.00000000e+00,
+ 1.46026823e+06, 0.00000000e+00, 0.00000000e+00,
+ 1.46026823e+06, 0.00000000e+00, 0.00000000e+00,
+ 2.05611458e+06, 0.00000000e+00, 0.00000000e+00,
+ 1.12054688e+06, 0.00000000e+00, 0.00000000e+00,
+ 5.19364583e+05, 0.00000000e+00, 0.00000000e+00,
+ 1.12054688e+06, 0.00000000e+00, 0.00000000e+00,
+ 2.05611458e+06, 0.00000000e+00, 0.00000000e+00,
+ 3.25847917e+06, 0.00000000e+00, 0.00000000e+00,
+ 0.00000000e+00, -3.13041667e+05, 0.00000000e+00,
+ 0.00000000e+00, -6.43869792e+05, 0.00000000e+00,
+ 0.00000000e+00, 1.46026823e+06, 0.00000000e+00,
+ 0.00000000e+00, 1.46026823e+06, 0.00000000e+00,
+ 0.00000000e+00, 2.05611458e+06, 0.00000000e+00,
+ 0.00000000e+00, 1.12054688e+06, 0.00000000e+00,
+ 0.00000000e+00, 5.19364583e+05, 0.00000000e+00,
+ 0.00000000e+00, 1.12054688e+06, 0.00000000e+00,
+ 0.00000000e+00, 2.05611458e+06, 0.00000000e+00,
+ 0.00000000e+00, 3.25847917e+06, 0.00000000e+00,
+ 0.00000000e+00, 0.00000000e+00, -3.13041667e+05,
+ 0.00000000e+00, 0.00000000e+00, -6.43869792e+05,
+ 0.00000000e+00, 0.00000000e+00, 1.46026823e+06,
+ 0.00000000e+00, 0.00000000e+00, 1.46026823e+06,
+ 0.00000000e+00, 0.00000000e+00, 2.05611458e+06,
+ 0.00000000e+00, 0.00000000e+00, 1.12054688e+06,
+ 0.00000000e+00, 0.00000000e+00, 5.19364583e+05,
+ 0.00000000e+00, 0.00000000e+00, 1.12054688e+06,
+ 0.00000000e+00, 0.00000000e+00, 2.05611458e+06,
+ 0.00000000e+00, 0.00000000e+00, 3.25847917e+06,
+};
+
+pylith::feassemble::ElasticityExplicitGravData3DQuadratic::ElasticityExplicitGravData3DQuadratic(void)
+{ // constructor
+ spaceDim = _spaceDim;
+ cellDim = _cellDim;
+ numVertices = _numVertices;
+ numCells = _numCells;
+ numBasis = _numBasis;
+ numQuadPts = _numQuadPts;
+ matType = const_cast<char*>(_matType);
+ matDBFilename = const_cast<char*>(_matDBFilename);
+ matId = _matId;
+ matLabel = const_cast<char*>(_matLabel);
+ dt = _dt;
+ gravityVec = const_cast<double*>(_gravityVec);
+ vertices = const_cast<double*>(_vertices);
+ cells = const_cast<int*>(_cells);
+ verticesRef = const_cast<double*>(_verticesRef);
+ quadPts = const_cast<double*>(_quadPts);
+ quadWts = const_cast<double*>(_quadWts);
+ basis = const_cast<double*>(_basis);
+ basisDerivRef = const_cast<double*>(_basisDerivRef);
+ fieldTIncr = const_cast<double*>(_fieldTIncr);
+ fieldT = const_cast<double*>(_fieldT);
+ fieldTmdt = const_cast<double*>(_fieldTmdt);
+ valsResidual = const_cast<double*>(_valsResidual);
+ valsJacobian = const_cast<double*>(_valsJacobian);
+} // constructor
+
+pylith::feassemble::ElasticityExplicitGravData3DQuadratic::~ElasticityExplicitGravData3DQuadratic(void)
+{}
+
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData3DQuadratic.cc
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Added: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData3DQuadratic.hh
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData3DQuadratic.hh (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData3DQuadratic.hh 2009-12-04 18:58:03 UTC (rev 16068)
@@ -0,0 +1,92 @@
+// -*- C++ -*-
+//
+// ======================================================================
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ======================================================================
+//
+
+// DO NOT EDIT THIS FILE
+// This file was generated from python application integratorelasticity.
+
+#if !defined(pylith_feassemble_elasticityexplicitgravdata3dquadratic_hh)
+#define pylith_feassemble_elasticityexplicitgravdata3dquadratic_hh
+
+#include "IntegratorData.hh"
+
+namespace pylith {
+ namespace feassemble {
+ class ElasticityExplicitGravData3DQuadratic;
+ } // pylith
+} // feassemble
+
+class pylith::feassemble::ElasticityExplicitGravData3DQuadratic : public IntegratorData
+{
+
+public:
+
+ /// Constructor
+ ElasticityExplicitGravData3DQuadratic(void);
+
+ /// Destructor
+ ~ElasticityExplicitGravData3DQuadratic(void);
+
+private:
+
+ static const int _spaceDim;
+
+ static const int _cellDim;
+
+ static const int _numVertices;
+
+ static const int _numCells;
+
+ static const int _numBasis;
+
+ static const int _numQuadPts;
+
+ static const char* _matType;
+
+ static const char* _matDBFilename;
+
+ static const int _matId;
+
+ static const char* _matLabel;
+
+ static const double _dt;
+
+ static const double _gravityVec[];
+
+ static const double _vertices[];
+
+ static const int _cells[];
+
+ static const double _verticesRef[];
+
+ static const double _quadPts[];
+
+ static const double _quadWts[];
+
+ static const double _basis[];
+
+ static const double _basisDerivRef[];
+
+ static const double _fieldTIncr[];
+
+ static const double _fieldT[];
+
+ static const double _fieldTmdt[];
+
+ static const double _valsResidual[];
+
+ static const double _valsJacobian[];
+
+};
+
+#endif // pylith_feassemble_elasticityexplicitgravdata3dquadratic_hh
+
+// End of file
Property changes on: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/ElasticityExplicitGravData3DQuadratic.hh
___________________________________________________________________
Name: svn:mime-type
+ text/plain
Modified: short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/generate.sh
===================================================================
--- short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/generate.sh 2009-12-04 17:44:24 UTC (rev 16067)
+++ short/3D/PyLith/branches/pylith-friction/unittests/libtests/feassemble/data/generate.sh 2009-12-04 18:58:03 UTC (rev 16068)
@@ -156,6 +156,28 @@
--material=MaterialElasticStrain1D \
--solution=Solution1DQuadratic
+ python IntegratorElasticity.py \
+ --use_gravity=True \
+ --formulation=ElasticityExplicit \
+ --data.namespace=pylith,feassemble \
+ --data.object=ElasticityExplicitGravData1DLinear \
+ --data.parent=IntegratorData \
+ --mesh=Mesh1DLinear \
+ --quadrature=Quadrature1DLinear \
+ --material=MaterialElasticStrain1D \
+ --solution=Solution1DLinear
+
+ python IntegratorElasticity.py \
+ --use_gravity=True \
+ --formulation=ElasticityExplicit \
+ --data.namespace=pylith,feassemble \
+ --data.object=ElasticityExplicitGravData1DQuadratic \
+ --data.parent=IntegratorData \
+ --mesh=Mesh1DQuadratic \
+ --quadrature=Quadrature1DQuadratic \
+ --material=MaterialElasticStrain1D \
+ --solution=Solution1DQuadratic
+
# 2-D ----------------------------------------------------------------
python IntegratorElasticity.py \
@@ -178,7 +200,30 @@
--material=ElasticPlaneStrain \
--solution=Solution2DQuadratic
+ python IntegratorElasticity.py \
+ --use_gravity=True \
+ --formulation=ElasticityExplicit \
+ --data.namespace=pylith,feassemble \
+ --data.object=ElasticityExplicitGravData2DLinear \
+ --data.parent=IntegratorData \
+ --mesh=Mesh2DLinear \
+ --quadrature=Quadrature2DLinear \
+ --material=ElasticPlaneStrain \
+ --solution=Solution2DLinear
+ python IntegratorElasticity.py \
+ --use_gravity=True \
+ --formulation=ElasticityExplicit \
+ --data.namespace=pylith,feassemble \
+ --data.object=ElasticityExplicitGravData2DQuadratic \
+ --data.parent=IntegratorData \
+ --mesh=Mesh2DQuadratic \
+ --quadrature=Quadrature2DQuadratic \
+ --material=ElasticPlaneStrain \
+ --solution=Solution2DQuadratic
+
+
+
# 3-D ----------------------------------------------------------------
python IntegratorElasticity.py \
@@ -201,6 +246,27 @@
--material=ElasticIsotropic3D \
--solution=Solution3DQuadratic
+ python IntegratorElasticity.py \
+ --use_gravity=True \
+ --formulation=ElasticityExplicit \
+ --data.namespace=pylith,feassemble \
+ --data.object=ElasticityExplicitGravData3DLinear \
+ --data.parent=IntegratorData \
+ --mesh=Mesh3DLinear \
+ --quadrature=Quadrature3DLinear \
+ --material=ElasticIsotropic3D \
+ --solution=Solution3DLinear
+
+ python IntegratorElasticity.py \
+ --use_gravity=True \
+ --formulation=ElasticityExplicit \
+ --data.namespace=pylith,feassemble \
+ --data.object=ElasticityExplicitGravData3DQuadratic \
+ --data.parent=IntegratorData \
+ --mesh=Mesh3DQuadratic \
+ --quadrature=Quadrature3DQuadratic \
+ --material=ElasticIsotropic3D \
+ --solution=Solution3DQuadratic
fi
@@ -339,6 +405,7 @@
--quadrature=Quadrature3DQuadratic \
--material=ElasticIsotropic3D \
--solution=Solution3DQuadratic
+
fi
# //////////////////////////////////////////////////////////////////////
More information about the CIG-COMMITS
mailing list