[cig-commits] r14001 - short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology
brad at geodynamics.org
brad at geodynamics.org
Sun Feb 1 10:07:34 PST 2009
Author: brad
Date: 2009-02-01 10:07:33 -0800 (Sun, 01 Feb 2009)
New Revision: 14001
Added:
short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology/TestFieldSubMesh.cc
short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology/TestFieldSubMesh.hh
Modified:
short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology/Makefile.am
short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology/TestField.hh
short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology/TestSubMesh.hh
Log:
Added unit tests for FieldSubMesh.
Modified: short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology/Makefile.am
===================================================================
--- short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology/Makefile.am 2009-02-01 00:15:18 UTC (rev 14000)
+++ short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology/Makefile.am 2009-02-01 18:07:33 UTC (rev 14001)
@@ -21,22 +21,24 @@
# Primary source files
testtopology_SOURCES = \
- TestFieldBase.cc \
- TestField.cc \
TestMesh.cc \
TestSubMesh.cc \
TestMeshOps.cc \
+ TestFieldBase.cc \
+ TestField.cc \
+ TestFieldSubMesh.cc \
test_topology.cc
# TestFieldsManager.cc
noinst_HEADERS = \
+ TestMesh.hh \
+ TestSubMesh.hh \
+ TestMeshOps.hh \
TestFieldBase.hh \
TestField.hh \
- TestFieldsManager.hh \
- TestMesh.hh \
- TestSubMesh.hh \
- TestMeshOps.hh
+ TestFieldSubMesh.hh \
+ TestFieldsManager.hh
# Source files associated with testing data
testtopology_SOURCES +=
Modified: short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology/TestField.hh
===================================================================
--- short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology/TestField.hh 2009-02-01 00:15:18 UTC (rev 14000)
+++ short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology/TestField.hh 2009-02-01 18:07:33 UTC (rev 14001)
@@ -118,6 +118,7 @@
*
* @param mesh Finite-element mesh.
*/
+ static
void _buildMesh(Mesh* mesh);
}; // class TestField
Added: short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology/TestFieldSubMesh.cc
===================================================================
--- short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology/TestFieldSubMesh.cc (rev 0)
+++ short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology/TestFieldSubMesh.cc 2009-02-01 18:07:33 UTC (rev 14001)
@@ -0,0 +1,700 @@
+// -*- C++ -*-
+//
+// ----------------------------------------------------------------------
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ----------------------------------------------------------------------
+//
+
+#include <portinfo>
+
+#include "TestFieldSubMesh.hh" // Implementation of class methods
+#include "pylith/topology/FieldSubMesh.hh" // USES FieldSubMesh
+
+#include "pylith/topology/Mesh.hh" // USES Mesh
+#include "pylith/topology/SubMesh.hh" // USES SubMesh
+#include "pylith/utils/array.hh" // USES double_array
+
+#include "spatialdata/geocoords/CSCart.hh" // USES CSCart
+
+// ----------------------------------------------------------------------
+CPPUNIT_TEST_SUITE_REGISTRATION( pylith::topology::TestFieldSubMesh );
+
+// ----------------------------------------------------------------------
+namespace pylith {
+ namespace topology {
+ namespace _TestFieldSubMesh {
+ const int cellDim = 2;
+ const int nvertices = 4;
+ const int ncells = 2;
+ const int ncorners = 3;
+ const int cells[] = {
+ 0, 1, 3,
+ 0, 3, 2,
+ };
+ const double coordinates[] = {
+ 0.0, 0.0,
+ 1.0, 0.0,
+ 0.0, 1.0,
+ 1.0, 1.0,
+ };
+ const char* label = "bc";
+ const int groupSize = 3;
+ const int groupVertices[] = {
+ 1, 2, 3
+ };
+ const int submeshVertices[] = {
+ 3, 4, 5
+ };
+ } // _TestFieldSubMesh
+ } // topology
+} // pylith
+
+// ----------------------------------------------------------------------
+// Test constructor.
+void
+pylith::topology::TestFieldSubMesh::testConstructor(void)
+{ // testConstructor
+ Mesh mesh;
+ SubMesh submesh;
+ _buildMesh(&mesh, &submesh);
+ FieldSubMesh field(submesh);
+} // testConstructor
+
+// ----------------------------------------------------------------------
+// Test newSection().
+void
+pylith::topology::TestFieldSubMesh::testSection(void)
+{ // testSection
+ Mesh mesh;
+ SubMesh submesh;
+ _buildMesh(&mesh, &submesh);
+ FieldSubMesh field(submesh);
+
+ mesh.createSieveMesh();
+ const ALE::Obj<SieveMesh::real_section_type>& section = field.section();
+ CPPUNIT_ASSERT(section.isNull());
+} // testSection
+
+// ----------------------------------------------------------------------
+// Test spaceDim().
+void
+pylith::topology::TestFieldSubMesh::testSpaceDim(void)
+{ // testSpaceDim
+ Mesh mesh;
+ SubMesh submesh;
+ _buildMesh(&mesh, &submesh);
+ FieldSubMesh field(submesh);
+
+ CPPUNIT_ASSERT_EQUAL(_TestFieldSubMesh::cellDim, field.spaceDim());
+} // testSpaceDim
+
+// ----------------------------------------------------------------------
+// Test newSection().
+void
+pylith::topology::TestFieldSubMesh::testNewSection(void)
+{ // testNewSection
+ Mesh mesh;
+ SubMesh submesh;
+ _buildMesh(&mesh, &submesh);
+ FieldSubMesh field(submesh);
+
+ field.newSection();
+ const ALE::Obj<SieveMesh::real_section_type>& section = field.section();
+ CPPUNIT_ASSERT(!section.isNull());
+} // testNewSection
+
+// ----------------------------------------------------------------------
+// Test newSection(points).
+void
+pylith::topology::TestFieldSubMesh::testNewSectionPoints(void)
+{ // testNewSectionPoints
+ const int fiberDim = 2;
+
+ Mesh mesh;
+ SubMesh submesh;
+ _buildMesh(&mesh, &submesh);
+ const ALE::Obj<SieveMesh>& sieveMesh = submesh.sieveMesh();
+
+ FieldSubMesh field(submesh);
+ const ALE::Obj<SieveMesh::label_sequence>& vertices =
+ sieveMesh->depthStratum(0);
+ field.newSection(vertices, fiberDim);
+ const ALE::Obj<SieveMesh::real_section_type>& section = field.section();
+ CPPUNIT_ASSERT(!section.isNull());
+
+ CPPUNIT_ASSERT(!vertices.isNull());
+ for (SieveMesh::label_sequence::iterator v_iter=vertices->begin();
+ v_iter != vertices->end();
+ ++v_iter)
+ CPPUNIT_ASSERT_EQUAL(fiberDim, section->getFiberDimension(*v_iter));
+} // testNewSectionPoints
+
+// ----------------------------------------------------------------------
+// Test newSection(domain).
+void
+pylith::topology::TestFieldSubMesh::testNewSectionDomain(void)
+{ // testNewSectionDomain
+ const int fiberDim = 2;
+
+ Mesh mesh;
+ SubMesh submesh;
+ _buildMesh(&mesh, &submesh);
+ const ALE::Obj<SieveMesh>& sieveMesh = submesh.sieveMesh();
+
+ FieldSubMesh field(submesh);
+ field.newSection(FieldSubMesh::VERTICES_FIELD, fiberDim);
+
+ const ALE::Obj<SieveMesh::real_section_type>& section = field.section();
+ CPPUNIT_ASSERT(!section.isNull());
+ const ALE::Obj<SieveMesh::label_sequence>& vertices =
+ sieveMesh->depthStratum(0);
+ CPPUNIT_ASSERT(!vertices.isNull());
+ for (SieveMesh::label_sequence::iterator v_iter=vertices->begin();
+ v_iter != vertices->end();
+ ++v_iter)
+ CPPUNIT_ASSERT_EQUAL(fiberDim, section->getFiberDimension(*v_iter));
+} // testNewSectionDomain
+
+// ----------------------------------------------------------------------
+// Test newSection(chart).
+void
+pylith::topology::TestFieldSubMesh::testNewSectionChart(void)
+{ // testNewSectionChart
+ const int fiberDim = 3;
+ const int nconstraints[] = { 0, 2, 1, 3 };
+ const int constraints[] = {
+ // 0
+ 0, 3, // 1
+ 2, // 2
+ 0, 1, 2, // 3
+ };
+
+ Mesh mesh;
+ SubMesh submesh;
+ _buildMesh(&mesh, &submesh);
+ const ALE::Obj<SieveMesh>& sieveMesh = submesh.sieveMesh();
+
+ // Create field with atlas to use to create new field
+ FieldSubMesh fieldSrc(submesh);
+ fieldSrc.newSection(FieldSubMesh::VERTICES_FIELD, fiberDim);
+ const ALE::Obj<MeshRealSection>& sectionSrc = fieldSrc.section();
+ CPPUNIT_ASSERT(!sectionSrc.isNull());
+ const MeshRealSection::chart_type& chart = sectionSrc->getChart();
+
+ FieldSubMesh field(submesh);
+ field.newSection(chart, fiberDim);
+ const ALE::Obj<MeshRealSection>& section = field.section();
+ const ALE::Obj<SieveMesh::label_sequence>& vertices =
+ sieveMesh->depthStratum(0);
+ CPPUNIT_ASSERT(!vertices.isNull());
+ for (SieveMesh::label_sequence::iterator v_iter=vertices->begin();
+ v_iter != vertices->end();
+ ++v_iter)
+ CPPUNIT_ASSERT_EQUAL(fiberDim, section->getFiberDimension(*v_iter));
+} // testNewSectionChart
+
+// ----------------------------------------------------------------------
+// Test newSection(field).
+void
+pylith::topology::TestFieldSubMesh::testNewSectionField(void)
+{ // testNewSectionField
+ const int fiberDim = 3;
+ const int nconstraints[] = { 0, 2, 1, 3 };
+ const int constraints[] = {
+ // 0
+ 0, 3, // 1
+ 2, // 2
+ 0, 1, 2, // 3
+ };
+
+ Mesh mesh;
+ SubMesh submesh;
+ _buildMesh(&mesh, &submesh);
+ const ALE::Obj<SieveMesh>& sieveMesh = submesh.sieveMesh();
+
+ const ALE::Obj<SieveMesh::label_sequence>& vertices =
+ sieveMesh->depthStratum(0);
+ CPPUNIT_ASSERT(!vertices.isNull());
+
+ // Create field with atlas to use to create new field
+ FieldSubMesh fieldSrc(submesh);
+ { // Setup source field
+ fieldSrc.newSection(FieldSubMesh::VERTICES_FIELD, fiberDim);
+ const ALE::Obj<MeshRealSection>& section = fieldSrc.section();
+ int iV=0;
+
+ CPPUNIT_ASSERT(!vertices.isNull());
+ for (SieveMesh::label_sequence::iterator v_iter=vertices->begin();
+ v_iter != vertices->end();
+ ++v_iter)
+ section->addConstraintDimension(*v_iter, nconstraints[iV++]);
+ fieldSrc.allocate();
+ } // Setup source field
+
+
+ FieldSubMesh field(submesh);
+ field.newSection(fieldSrc);
+ const ALE::Obj<MeshRealSection>& section = field.section();
+ int iV = 0;
+ for (SieveMesh::label_sequence::iterator v_iter=vertices->begin();
+ v_iter != vertices->end();
+ ++v_iter) {
+ CPPUNIT_ASSERT_EQUAL(fiberDim, section->getFiberDimension(*v_iter));
+ CPPUNIT_ASSERT_EQUAL(nconstraints[iV++],
+ section->getConstraintDimension(*v_iter));
+ } // for
+} // testNewSectionField
+
+// ----------------------------------------------------------------------
+// Test clear().
+void
+pylith::topology::TestFieldSubMesh::testClear(void)
+{ // testClear
+ Mesh mesh;
+ SubMesh submesh;
+ _buildMesh(&mesh, &submesh);
+ FieldSubMesh field(submesh);
+
+ field.scale(2.0);
+ field.vectorFieldType(FieldSubMesh::TENSOR);
+ field.addDimensionOkay(true);
+
+ field.clear();
+
+ CPPUNIT_ASSERT_EQUAL(1.0, field._scale);
+ CPPUNIT_ASSERT_EQUAL(FieldSubMesh::OTHER, field._vecFieldType);
+ CPPUNIT_ASSERT_EQUAL(false, field._dimensionsOkay);
+} // testClear
+
+// ----------------------------------------------------------------------
+// Test allocate().
+void
+pylith::topology::TestFieldSubMesh::testAllocate(void)
+{ // testAllocate
+ const int fiberDim = 3;
+ const double scale = 2.0;
+ const double valuesNondim[] = {
+ 1.1, 2.2, 3.3,
+ 1.2, 2.3, 3.4,
+ 1.3, 2.4, 3.5,
+ };
+
+ Mesh mesh;
+ SubMesh submesh;
+ _buildMesh(&mesh, &submesh);
+ const ALE::Obj<SieveMesh>& sieveMesh = submesh.sieveMesh();
+ const ALE::Obj<SieveMesh::label_sequence>& vertices =
+ sieveMesh->depthStratum(0);
+
+ FieldSubMesh field(submesh);
+ field.newSection(FieldSubMesh::VERTICES_FIELD, fiberDim);
+ field.allocate();
+ const ALE::Obj<MeshRealSection>& section = field.section();
+
+ double_array values(fiberDim);
+ int i = 0;
+ for (SieveMesh::label_sequence::iterator v_iter=vertices->begin();
+ v_iter != vertices->end();
+ ++v_iter) {
+ for (int iDim=0; iDim < fiberDim; ++iDim)
+ values[iDim] = valuesNondim[i++];
+ section->updatePoint(*v_iter, &values[0]);
+ } // for
+
+ const double tolerance = 1.0e-6;
+ i = 0;
+ for (SieveMesh::label_sequence::iterator v_iter=vertices->begin();
+ v_iter != vertices->end();
+ ++v_iter) {
+ section->restrictPoint(*v_iter, &values[0], values.size());
+ for (int iDim=0; iDim < fiberDim; ++iDim) {
+ CPPUNIT_ASSERT_DOUBLES_EQUAL(valuesNondim[i++], values[iDim], tolerance);
+ } // for
+ } // for
+} // testAllocate
+
+// ----------------------------------------------------------------------
+// Test zero().
+void
+pylith::topology::TestFieldSubMesh::testZero(void)
+{ // testZero
+ const int fiberDim = 3;
+ const double scale = 2.0;
+ const double valuesNondim[] = {
+ 1.1, 2.2, 3.3,
+ 1.2, 2.3, 3.4,
+ 1.3, 2.4, 3.5,
+ };
+
+ Mesh mesh;
+ SubMesh submesh;
+ _buildMesh(&mesh, &submesh);
+ const ALE::Obj<SieveMesh>& sieveMesh = submesh.sieveMesh();
+ const ALE::Obj<SieveMesh::label_sequence>& vertices =
+ sieveMesh->depthStratum(0);
+
+ FieldSubMesh field(submesh);
+ field.newSection(FieldSubMesh::VERTICES_FIELD, fiberDim);
+ field.allocate();
+ const ALE::Obj<MeshRealSection>& section = field.section();
+
+ double_array values(fiberDim);
+ int i = 0;
+ for (SieveMesh::label_sequence::iterator v_iter=vertices->begin();
+ v_iter != vertices->end();
+ ++v_iter) {
+ for (int iDim=0; iDim < fiberDim; ++iDim)
+ values[iDim] = valuesNondim[i++];
+ section->updatePoint(*v_iter, &values[0]);
+ } // for
+
+ field.zero();
+
+ const double tolerance = 1.0e-6;
+ for (SieveMesh::label_sequence::iterator v_iter=vertices->begin();
+ v_iter != vertices->end();
+ ++v_iter) {
+ section->restrictPoint(*v_iter, &values[0], values.size());
+ for (int iDim=0; iDim < fiberDim; ++iDim) {
+ CPPUNIT_ASSERT_DOUBLES_EQUAL(0.0, values[iDim], tolerance);
+ } // for
+ } // for
+} // testZero
+
+// ----------------------------------------------------------------------
+// Test complete().
+void
+pylith::topology::TestFieldSubMesh::testComplete(void)
+{ // testComplete
+ const int fiberDim = 3;
+ const double scale = 2.0;
+ const double valuesNondim[] = {
+ 1.1, 2.2, 3.3,
+ 1.2, 2.3, 3.4,
+ 1.3, 2.4, 3.5,
+ };
+
+ Mesh mesh;
+ SubMesh submesh;
+ _buildMesh(&mesh, &submesh);
+ const ALE::Obj<SieveMesh>& sieveMesh = submesh.sieveMesh();
+ const ALE::Obj<SieveMesh::label_sequence>& vertices =
+ sieveMesh->depthStratum(0);
+
+ FieldSubMesh field(submesh);
+ field.newSection(FieldSubMesh::VERTICES_FIELD, fiberDim);
+ field.allocate();
+ const ALE::Obj<MeshRealSection>& section = field.section();
+
+ double_array values(fiberDim);
+ int i = 0;
+ for (SieveMesh::label_sequence::iterator v_iter=vertices->begin();
+ v_iter != vertices->end();
+ ++v_iter) {
+ for (int iDim=0; iDim < fiberDim; ++iDim)
+ values[iDim] = valuesNondim[i++];
+ section->updatePoint(*v_iter, &values[0]);
+ } // for
+
+ field.complete();
+
+ // Expect no change for this serial test
+ i = 0;
+ const double tolerance = 1.0e-6;
+ for (SieveMesh::label_sequence::iterator v_iter=vertices->begin();
+ v_iter != vertices->end();
+ ++v_iter) {
+ section->restrictPoint(*v_iter, &values[0], values.size());
+ for (int iDim=0; iDim < fiberDim; ++iDim) {
+ CPPUNIT_ASSERT_DOUBLES_EQUAL(valuesNondim[i++], values[iDim], tolerance);
+ } // for
+ } // for
+} // testComplete
+
+// ----------------------------------------------------------------------
+// Test copy().
+void
+pylith::topology::TestFieldSubMesh::testCopy(void)
+{ // testCopy
+ const int fiberDim = 3;
+ const double scale = 2.0;
+ const double valuesNondim[] = {
+ 1.1, 2.2, 3.3,
+ 1.2, 2.3, 3.4,
+ 1.3, 2.4, 3.5,
+ };
+
+ Mesh mesh;
+ SubMesh submesh;
+ _buildMesh(&mesh, &submesh);
+ const ALE::Obj<SieveMesh>& sieveMesh = submesh.sieveMesh();
+ const ALE::Obj<SieveMesh::label_sequence>& vertices =
+ sieveMesh->depthStratum(0);
+
+ FieldSubMesh fieldSrc(submesh);
+ { // Setup source field
+ fieldSrc.newSection(FieldSubMesh::VERTICES_FIELD, fiberDim);
+ fieldSrc.allocate();
+ const ALE::Obj<MeshRealSection>& section = fieldSrc.section();
+
+ double_array values(fiberDim);
+ int i = 0;
+ for (SieveMesh::label_sequence::iterator v_iter=vertices->begin();
+ v_iter != vertices->end();
+ ++v_iter) {
+ for (int iDim=0; iDim < fiberDim; ++iDim)
+ values[iDim] = valuesNondim[i++];
+ section->updatePoint(*v_iter, &values[0]);
+ } // for
+ } // Setup source field
+
+ FieldSubMesh field(submesh);
+ field.newSection(FieldSubMesh::VERTICES_FIELD, fiberDim);
+ field.allocate();
+ const ALE::Obj<MeshRealSection>& section = field.section();
+
+ field.copy(fieldSrc);
+
+ int i = 0;
+ double_array values(fiberDim);
+ const double tolerance = 1.0e-6;
+ for (SieveMesh::label_sequence::iterator v_iter=vertices->begin();
+ v_iter != vertices->end();
+ ++v_iter) {
+ section->restrictPoint(*v_iter, &values[0], values.size());
+ for (int iDim=0; iDim < fiberDim; ++iDim) {
+ CPPUNIT_ASSERT_DOUBLES_EQUAL(valuesNondim[i++], values[iDim], tolerance);
+ } // for
+ } // for
+} // testCopy
+
+// ----------------------------------------------------------------------
+// Test operator+=().
+void
+pylith::topology::TestFieldSubMesh::testOperatorAdd(void)
+{ // testOperateAdd
+ const int fiberDim = 3;
+ const double scale = 2.0;
+ const double valuesA[] = {
+ 1.1, 2.2, 3.3,
+ 1.2, 2.3, 3.4,
+ 1.3, 2.4, 3.5,
+ };
+ const double valuesB[] = {
+ 10.1, 20.2, 30.3,
+ 10.2, 20.3, 30.4,
+ 10.3, 20.4, 30.5,
+ };
+
+ Mesh mesh;
+ SubMesh submesh;
+ _buildMesh(&mesh, &submesh);
+ const ALE::Obj<SieveMesh>& sieveMesh = submesh.sieveMesh();
+ const ALE::Obj<SieveMesh::label_sequence>& vertices =
+ sieveMesh->depthStratum(0);
+
+ FieldSubMesh fieldSrc(submesh);
+ { // Setup source field
+ fieldSrc.newSection(FieldSubMesh::VERTICES_FIELD, fiberDim);
+ fieldSrc.allocate();
+ const ALE::Obj<MeshRealSection>& section = fieldSrc.section();
+
+ double_array values(fiberDim);
+ int i = 0;
+ for (SieveMesh::label_sequence::iterator v_iter=vertices->begin();
+ v_iter != vertices->end();
+ ++v_iter) {
+ for (int iDim=0; iDim < fiberDim; ++iDim)
+ values[iDim] = valuesA[i++];
+ section->updatePoint(*v_iter, &values[0]);
+ } // for
+ } // Setup source field
+
+ FieldSubMesh field(submesh);
+ field.newSection(FieldSubMesh::VERTICES_FIELD, fiberDim);
+ field.allocate();
+ const ALE::Obj<MeshRealSection>& section = field.section();
+ { // Setup destination field
+
+ double_array values(fiberDim);
+ int i = 0;
+ for (SieveMesh::label_sequence::iterator v_iter=vertices->begin();
+ v_iter != vertices->end();
+ ++v_iter) {
+ for (int iDim=0; iDim < fiberDim; ++iDim)
+ values[iDim] = valuesB[i++];
+ section->updatePoint(*v_iter, &values[0]);
+ } // for
+ } // Setup destination field
+
+ field += fieldSrc;
+
+ int i = 0;
+ double_array values(fiberDim);
+ const double tolerance = 1.0e-6;
+ for (SieveMesh::label_sequence::iterator v_iter=vertices->begin();
+ v_iter != vertices->end();
+ ++v_iter) {
+ section->restrictPoint(*v_iter, &values[0], values.size());
+ for (int iDim=0; iDim < fiberDim; ++iDim) {
+ const double valueE = valuesA[i] + valuesB[i];
+ ++i;
+ CPPUNIT_ASSERT_DOUBLES_EQUAL(valueE, values[iDim], tolerance);
+ } // for
+ } // for
+} // testOperateAdd
+
+// ----------------------------------------------------------------------
+// Test dimensionalize().
+void
+pylith::topology::TestFieldSubMesh::testDimensionalize(void)
+{ // testDimensionalize
+ const int fiberDim = 3;
+ const double scale = 2.0;
+ const double valuesNondim[] = {
+ 1.1, 2.2, 3.3,
+ 1.2, 2.3, 3.4,
+ 1.3, 2.4, 3.5,
+ };
+
+ Mesh mesh;
+ SubMesh submesh;
+ _buildMesh(&mesh, &submesh);
+ const ALE::Obj<SieveMesh>& sieveMesh = submesh.sieveMesh();
+
+ FieldSubMesh field(submesh);
+ field.newSection(FieldSubMesh::VERTICES_FIELD, fiberDim);
+ field.allocate();
+ const ALE::Obj<MeshRealSection>& section = field.section();
+ const ALE::Obj<SieveMesh::label_sequence>& vertices =
+ sieveMesh->depthStratum(0);
+
+ double_array values(fiberDim);
+ int i = 0;
+ for (SieveMesh::label_sequence::iterator v_iter=vertices->begin();
+ v_iter != vertices->end();
+ ++v_iter) {
+ for (int iDim=0; iDim < fiberDim; ++iDim)
+ values[iDim] = valuesNondim[i++];
+ section->updatePoint(*v_iter, &values[0]);
+ } // for
+
+ field.scale(scale);
+ field.addDimensionOkay(true);
+ field.dimensionalize();
+
+ i = 0;
+ const double tolerance = 1.0e-6;
+ for (SieveMesh::label_sequence::iterator v_iter=vertices->begin();
+ v_iter != vertices->end();
+ ++v_iter) {
+ section->restrictPoint(*v_iter, &values[0], values.size());
+ for (int iDim=0; iDim < fiberDim; ++iDim) {
+ const double valueE = valuesNondim[i++]*scale;
+ CPPUNIT_ASSERT_DOUBLES_EQUAL(valueE, values[iDim], tolerance);
+ } // for
+ } // for
+
+} // testDimensionalize
+
+// ----------------------------------------------------------------------
+// Test view().
+void
+pylith::topology::TestFieldSubMesh::testView(void)
+{ // testView
+ const int fiberDim = 3;
+ const double scale = 2.0;
+ const double valuesNondim[] = {
+ 1.1, 2.2, 3.3,
+ 1.2, 2.3, 3.4,
+ 1.3, 2.4, 3.5,
+ };
+
+ Mesh mesh;
+ SubMesh submesh;
+ _buildMesh(&mesh, &submesh);
+ const ALE::Obj<SieveMesh>& sieveMesh = submesh.sieveMesh();
+ FieldSubMesh field(submesh);
+ field.newSection(FieldSubMesh::VERTICES_FIELD, fiberDim);
+ field.allocate();
+ const ALE::Obj<MeshRealSection>& section = field.section();
+ const ALE::Obj<SieveMesh::label_sequence>& vertices =
+ sieveMesh->depthStratum(0);
+
+ double_array values(fiberDim);
+ int i = 0;
+ for (SieveMesh::label_sequence::iterator v_iter=vertices->begin();
+ v_iter != vertices->end();
+ ++v_iter) {
+ for (int iDim=0; iDim < fiberDim; ++iDim)
+ values[iDim] = valuesNondim[i++];
+ section->updatePoint(*v_iter, &values[0]);
+ } // for
+
+ field.view("Testing view");
+} // testView
+
+// ----------------------------------------------------------------------
+void
+pylith::topology::TestFieldSubMesh::_buildMesh(Mesh* mesh,
+ SubMesh* submesh)
+{ // _buildMesh
+ assert(0 != mesh);
+ assert(0 != submesh);
+
+ mesh->createSieveMesh(_TestFieldSubMesh::cellDim);
+ const ALE::Obj<SieveMesh>& sieveMesh = mesh->sieveMesh();
+
+ ALE::Obj<SieveMesh::sieve_type> sieve =
+ new SieveMesh::sieve_type(sieveMesh->comm());
+ CPPUNIT_ASSERT(!sieve.isNull());
+
+ ALE::Obj<ALE::Mesh::sieve_type> s =
+ new ALE::Mesh::sieve_type(sieve->comm(), sieve->debug());
+
+ const int cellDim = _TestFieldSubMesh::cellDim;
+ const int ncells = _TestFieldSubMesh::ncells;
+ const int* cells = _TestFieldSubMesh::cells;
+ const int nvertices = _TestFieldSubMesh::nvertices;
+ const int ncorners = _TestFieldSubMesh::ncorners;
+ const int spaceDim = _TestFieldSubMesh::cellDim;
+ const double* coordinates = _TestFieldSubMesh::coordinates;
+ const bool interpolate = false;
+ ALE::SieveBuilder<ALE::Mesh>::buildTopology(s, cellDim, ncells, (int*) cells,
+ nvertices, interpolate,
+ ncorners);
+ std::map<SieveMesh::point_type,SieveMesh::point_type> renumbering;
+ ALE::ISieveConverter::convertSieve(*s, *sieve, renumbering);
+ sieveMesh->setSieve(sieve);
+ sieveMesh->stratify();
+ ALE::SieveBuilder<SieveMesh>::buildCoordinates(sieveMesh, spaceDim,
+ coordinates);
+
+ typedef SieveMesh::int_section_type::chart_type chart_type;
+ const ALE::Obj<SieveMesh::int_section_type>& groupField =
+ sieveMesh->getIntSection(_TestFieldSubMesh::label);
+ assert(!groupField.isNull());
+
+ const int numPoints = _TestFieldSubMesh::groupSize;
+ const int numVertices = sieveMesh->depthStratum(0)->size();
+ const int numCells = sieveMesh->heightStratum(0)->size();
+ groupField->setChart(chart_type(numCells, numCells+numVertices));
+ for(int i=0; i < numPoints; ++i)
+ groupField->setFiberDimension(numCells+_TestFieldSubMesh::groupVertices[i],
+ 1);
+ sieveMesh->allocate(groupField);
+
+ spatialdata::geocoords::CSCart cs;
+ cs.setSpaceDim(spaceDim);
+ cs.initialize();
+ mesh->coordsys(&cs);
+
+ submesh->createSubMesh(*mesh, _TestFieldSubMesh::label);
+} // _buildMesh
+
+// End of file
Added: short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology/TestFieldSubMesh.hh
===================================================================
--- short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology/TestFieldSubMesh.hh (rev 0)
+++ short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology/TestFieldSubMesh.hh 2009-02-01 18:07:33 UTC (rev 14001)
@@ -0,0 +1,131 @@
+// -*- C++ -*-
+//
+// ----------------------------------------------------------------------
+//
+// Brad T. Aagaard
+// U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ----------------------------------------------------------------------
+//
+
+/**
+ * @file unittests/libtests/topology/TestfieldsubmeshSubMesh.hh
+ *
+ * @brief C++ unit testing for Field.
+ */
+
+#if !defined(pylith_topology_testfieldsubmesh_hh)
+#define pylith_topology_testfieldsubmesh_hh
+
+// Include directives ---------------------------------------------------
+#include <cppunit/extensions/HelperMacros.h>
+
+// Forward declarations -------------------------------------------------
+/// Namespace for pylith package
+namespace pylith {
+ namespace topology {
+ class TestFieldSubMesh;
+
+ class Mesh;
+ class SubMesh;
+ } // topology
+} // pylith
+
+// TestFieldSubMesh -----------------------------------------------------
+/// C++ unit testing for Field.
+class pylith::topology::TestFieldSubMesh : public CppUnit::TestFixture
+{ // class TestFieldSubMesh
+
+ // CPPUNIT TEST SUITE /////////////////////////////////////////////////
+ CPPUNIT_TEST_SUITE( TestFieldSubMesh );
+
+ CPPUNIT_TEST( testConstructor );
+ CPPUNIT_TEST( testSection );
+ CPPUNIT_TEST( testSpaceDim );
+ CPPUNIT_TEST( testNewSection );
+ CPPUNIT_TEST( testNewSectionPoints );
+ CPPUNIT_TEST( testNewSectionDomain );
+ CPPUNIT_TEST( testNewSectionChart );
+ CPPUNIT_TEST( testNewSectionField );
+ CPPUNIT_TEST( testClear );
+ CPPUNIT_TEST( testAllocate );
+ CPPUNIT_TEST( testZero );
+ CPPUNIT_TEST( testComplete );
+ CPPUNIT_TEST( testCopy );
+ CPPUNIT_TEST( testOperatorAdd );
+ CPPUNIT_TEST( testDimensionalize );
+ CPPUNIT_TEST( testView );
+
+ CPPUNIT_TEST_SUITE_END();
+
+ // PUBLIC METHODS /////////////////////////////////////////////////////
+public :
+
+ /// Test constructor.
+ void testConstructor(void);
+
+ /// Test section().
+ void testSection(void);
+
+ /// Test spaceDim().
+ void testSpaceDim(void);
+
+ /// Test newSection().
+ void testNewSection(void);
+
+ /// Test newSection(points).
+ void testNewSectionPoints(void);
+
+ /// Test newSection(domain).
+ void testNewSectionDomain(void);
+
+ /// Test newSection(chart).
+ void testNewSectionChart(void);
+
+ /// Test newSection(field).
+ void testNewSectionField(void);
+
+ /// Test clear().
+ void testClear(void);
+
+ /// Test allocate().
+ void testAllocate(void);
+
+ /// Test zero().
+ void testZero(void);
+
+ /// Test complete().
+ void testComplete(void);
+
+ /// Test copy().
+ void testCopy(void);
+
+ /// Test operator+=().
+ void testOperatorAdd(void);
+
+ /// Test dimensionalize().
+ void testDimensionalize(void);
+
+ /// Test view().
+ void testView(void);
+
+// PRIVATE METHODS /////////////////////////////////////////////////////
+private :
+
+ /** Build mesh.
+ *
+ * @param mesh Finite-element mesh.
+ * @param submesh Boundary mesh.
+ */
+ static
+ void _buildMesh(Mesh* mesh,
+ SubMesh* submesh);
+
+}; // class TestFieldSubMesh
+
+#endif // pylith_topology_testfieldsubmesh_hh
+
+
+// End of file
Modified: short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology/TestSubMesh.hh
===================================================================
--- short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology/TestSubMesh.hh 2009-02-01 00:15:18 UTC (rev 14000)
+++ short/3D/PyLith/branches/pylith-swig/unittests/libtests/topology/TestSubMesh.hh 2009-02-01 18:07:33 UTC (rev 14001)
@@ -89,6 +89,7 @@
*
* @param mesh Finite-element mesh.
*/
+ static
void _buildMesh(Mesh* mesh);
}; // class TestSubMesh
More information about the CIG-COMMITS
mailing list