[cig-commits] r6611 - in short/3D/PyLith/trunk: libsrc/faults unittests/libtests/faults unittests/libtests/faults/data

brad at geodynamics.org brad at geodynamics.org
Thu Apr 19 14:48:52 PDT 2007


Author: brad
Date: 2007-04-19 14:48:51 -0700 (Thu, 19 Apr 2007)
New Revision: 6611

Added:
   short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataHex8.cc
   short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataHex8.hh
   short/3D/PyLith/trunk/unittests/libtests/faults/data/meshHex8A.txt
Modified:
   short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc
   short/3D/PyLith/trunk/unittests/libtests/faults/Makefile.am
   short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesive.cc
   short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesive.hh
   short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataLine2.cc
   short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataQuad4.cc
   short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataTet4.cc
   short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataTri3.cc
   short/3D/PyLith/trunk/unittests/libtests/faults/data/Makefile.am
Log:
Added test for creating cohesive elements for hex8 element. Fixed unit tests to use Sieve numbering (original cells then vertices then cohesive cells).

Modified: short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc	2007-04-19 21:26:49 UTC (rev 6610)
+++ short/3D/PyLith/trunk/libsrc/faults/CohesiveTopology.cc	2007-04-19 21:48:51 UTC (rev 6611)
@@ -200,7 +200,8 @@
     }
   } // for
   mesh->stratify();
-  mesh->view("Mesh with Cohesive Elements");
+  if (debug)
+    mesh->view("Mesh with Cohesive Elements");
 
   // Fix coordinates
   const ALE::Obj<real_section_type>& coordinates = 

Modified: short/3D/PyLith/trunk/unittests/libtests/faults/Makefile.am
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/Makefile.am	2007-04-19 21:26:49 UTC (rev 6610)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/Makefile.am	2007-04-19 21:48:51 UTC (rev 6611)
@@ -35,6 +35,7 @@
 	data/CohesiveDataLine2.cc \
 	data/CohesiveDataTri3.cc \
 	data/CohesiveDataQuad4.cc \
+	data/CohesiveDataHex8.cc \
 	data/CohesiveDataTet4.cc
 
 noinst_HEADERS += \
@@ -42,6 +43,7 @@
 	data/CohesiveDataLine2.hh \
 	data/CohesiveDataTri3.hh \
 	data/CohesiveDataQuad4.hh \
+	data/CohesiveDataHex8.hh \
 	data/CohesiveDataTet4.hh
 
 testfaults_LDFLAGS = $(PETSC_LIB)

Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesive.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesive.cc	2007-04-19 21:26:49 UTC (rev 6610)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesive.cc	2007-04-19 21:48:51 UTC (rev 6611)
@@ -24,8 +24,8 @@
 #include "data/CohesiveDataTri3.hh" // USES CohesiveDataTri3
 #include "data/CohesiveDataQuad4.hh" // USES CohesiveDataQuad4
 #include "data/CohesiveDataTet4.hh" // USES CohesiveDataTet4
+#include "data/CohesiveDataHex8.hh" // USES CohesiveDataHex8
 
-
 // ----------------------------------------------------------------------
 CPPUNIT_TEST_SUITE_REGISTRATION( pylith::faults::TestFaultCohesive );
 
@@ -66,6 +66,15 @@
 } // testAdjustTopologyTet4
 
 // ----------------------------------------------------------------------
+// Test adjustTopology() with 3-D hexahedral element.
+void
+pylith::faults::TestFaultCohesive::testAdjustTopologyHex8(void)
+{ // testAdjustTopologyHex8
+  CohesiveDataHex8 data;
+  _testAdjustTopology(data);
+} // testAdjustTopologyHex8
+
+// ----------------------------------------------------------------------
 // Test adjustTopology().
 void
 pylith::faults::TestFaultCohesive::_testAdjustTopology(const CohesiveData& data)
@@ -73,7 +82,7 @@
   ALE::Obj<ALE::Mesh> mesh;
   meshio::MeshIOAscii iohandler;
   iohandler.filename(data.filename);
-  iohandler.debug(true);
+  iohandler.debug(false);
   iohandler.interpolate(false);
   iohandler.read(&mesh);
 
@@ -82,8 +91,6 @@
   fault.label("fault");
   fault.adjustTopology(&mesh);
 
-  mesh->view("Mesh");
-
   CPPUNIT_ASSERT_EQUAL(data.cellDim, mesh->getDimension());
 
   // Check vertices
@@ -119,19 +126,19 @@
   const int numCells = cells->size();
   CPPUNIT_ASSERT_EQUAL(data.numCells, numCells);
 
+  int iCell = 0;
   i = 0;
-  const int offset = numCells;
   for(Mesh::label_sequence::iterator c_iter = cells->begin();
       c_iter != cells->end();
       ++c_iter) {
     const int numCorners = sieve->nCone(*c_iter, mesh->depth())->size();
-    CPPUNIT_ASSERT_EQUAL(data.numCorners[i], numCorners);
+    CPPUNIT_ASSERT_EQUAL(data.numCorners[iCell++], numCorners);
     const ALE::Obj<sieve_type::traits::coneSequence>& cone = 
       sieve->cone(*c_iter);
     for(sieve_type::traits::coneSequence::iterator v_iter = cone->begin();
 	v_iter != cone->end();
 	++v_iter)
-      CPPUNIT_ASSERT_EQUAL(data.cells[i++], *v_iter - offset);
+      CPPUNIT_ASSERT_EQUAL(data.cells[i++], *v_iter);
   } // for
 
   // check materials
@@ -166,11 +173,10 @@
     const int numPoints = chart.size();
     int_array points(numPoints);
     int i = 0;
-    const int offset = ("vertex" == groupType) ? numCells : 0;
     for(int_section_type::chart_type::iterator c_iter = chart.begin();
 	c_iter != chart.end();
 	++c_iter)
-      points[i++] = *c_iter - offset;
+      points[i++] = *c_iter;
 
     CPPUNIT_ASSERT_EQUAL(std::string(data.groupNames[iGroup]), *name);
     CPPUNIT_ASSERT_EQUAL(std::string(data.groupTypes[iGroup]), groupType);

Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesive.hh
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesive.hh	2007-04-19 21:26:49 UTC (rev 6610)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesive.hh	2007-04-19 21:48:51 UTC (rev 6611)
@@ -44,6 +44,7 @@
   CPPUNIT_TEST( testAdjustTopologyTri3 );
   CPPUNIT_TEST( testAdjustTopologyQuad4 );
   CPPUNIT_TEST( testAdjustTopologyTet4 );
+  CPPUNIT_TEST( testAdjustTopologyHex8 );
   CPPUNIT_TEST_SUITE_END();
 
   // PUBLIC METHODS /////////////////////////////////////////////////////
@@ -61,6 +62,9 @@
   /// Test adjustTopology() with 3-D tetrahedral element.
   void testAdjustTopologyTet4(void);
 
+  /// Test adjustTopology() with 3-D hexahedral element.
+  void testAdjustTopologyHex8(void);
+
   // PROTECTED METHODS //////////////////////////////////////////////////
 public :
 

Added: short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataHex8.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataHex8.cc	2007-04-19 21:26:49 UTC (rev 6610)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataHex8.cc	2007-04-19 21:48:51 UTC (rev 6611)
@@ -0,0 +1,119 @@
+// -*- C++ -*-
+//
+// ======================================================================
+//
+//                           Brad T. Aagaard
+//                        U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ======================================================================
+//
+
+/* Original mesh
+ *
+ * Cells are 0-1 and vertices are 2-13.
+ *
+ *       2,3,4,5 -------- 6,7,8,9 -------- 10,11,12,13
+ *
+ *                        ^^^^^^^ Vertices forming fault
+ *
+ * After adding cohesive elements
+ *
+ * Cells are 0-1,16 and vertices are 4-15.
+ *
+ *       2,3,4,5 -------- 6,7,8,9 -- 14,15,16,17 -------- 10,11,12,13
+ *
+ *                        ^^^^^^^^^^^^^^^^^^^^^^ Cohesive element
+ *
+ */
+
+#include "CohesiveDataHex8.hh"
+
+const int pylith::faults::CohesiveDataHex8::_numVertices = 16;
+
+const int pylith::faults::CohesiveDataHex8::_spaceDim = 3;
+
+const int pylith::faults::CohesiveDataHex8::_numCells = 3;
+
+const int pylith::faults::CohesiveDataHex8::_cellDim = 3;
+
+const double pylith::faults::CohesiveDataHex8::_vertices[] = {
+  -2.0, -1.0, -1.0,
+  -2.0,  1.0, -1.0,
+  -2.0, -1.0,  1.0,
+  -2.0,  1.0,  1.0,
+   0.0, -1.0, -1.0,
+   0.0,  1.0, -1.0,
+   0.0, -1.0,  1.0,
+   0.0,  1.0,  1.0,
+   2.0, -1.0, -1.0,
+   2.0,  1.0, -1.0,
+   2.0, -1.0,  1.0,
+   2.0,  1.0,  1.0,
+   0.0, -1.0, -1.0,
+   0.0,  1.0, -1.0,
+   0.0, -1.0,  1.0,
+   0.0,  1.0,  1.0,
+};
+
+const int pylith::faults::CohesiveDataHex8::_numCorners[] = {
+  8,
+  8,
+  8
+};
+
+const int pylith::faults::CohesiveDataHex8::_cells[] = {
+  2,  4,  5,  3,  6,  8,  9,  7,
+ 14, 16, 17, 15, 10, 12, 13, 11,
+  6,  8,  9,  7, 14, 16, 17, 15,
+};
+
+const int pylith::faults::CohesiveDataHex8::_materialIds[] = {
+  0,  0,
+  1
+};
+
+const int pylith::faults::CohesiveDataHex8::_numGroups = 2;
+
+const int pylith::faults::CohesiveDataHex8::_groupSizes[] = 
+  { 8, 8 };
+
+const int pylith::faults::CohesiveDataHex8::_groups[] = {
+  6, 7, 8, 9, 14, 15, 16, 17,
+  4, 5, 8, 9, 12, 13, 16, 17
+};
+
+const char* pylith::faults::CohesiveDataHex8::_groupNames[] = {
+  "fault", "output"
+};
+
+const char* pylith::faults::CohesiveDataHex8::_groupTypes[] = {
+  "vertex", "vertex"
+};
+
+const char* pylith::faults::CohesiveDataHex8::_filename = "data/meshHex8A.txt";
+
+pylith::faults::CohesiveDataHex8::CohesiveDataHex8(void)
+{ // constructor
+  numVertices = _numVertices;
+  spaceDim = _spaceDim;
+  numCells = _numCells;
+  cellDim = _cellDim;
+  vertices = const_cast<double*>(_vertices);
+  numCorners = const_cast<int*>(_numCorners);
+  cells = const_cast<int*>(_cells);
+  materialIds = const_cast<int*>(_materialIds);
+  groups = const_cast<int*>(_groups);
+  groupSizes = const_cast<int*>(_groupSizes);
+  groupNames = const_cast<char**>(_groupNames);
+  groupTypes = const_cast<char**>(_groupTypes);
+  numGroups = _numGroups;
+  filename = const_cast<char*>(_filename);
+} // constructor
+
+pylith::faults::CohesiveDataHex8::~CohesiveDataHex8(void)
+{}
+
+
+// End of file

Added: short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataHex8.hh
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataHex8.hh	2007-04-19 21:26:49 UTC (rev 6610)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataHex8.hh	2007-04-19 21:48:51 UTC (rev 6611)
@@ -0,0 +1,60 @@
+// -*- C++ -*-
+//
+// ======================================================================
+//
+//                           Brad T. Aagaard
+//                        U.S. Geological Survey
+//
+// {LicenseText}
+//
+// ======================================================================
+//
+
+#if !defined(pylith_faults_cohesivedatahex8_hh)
+#define pylith_faults_cohesivedatahex8_hh
+
+#include "CohesiveData.hh"
+
+namespace pylith {
+  namespace faults {
+     class CohesiveDataHex8;
+  } // pylith
+} // faults
+
+class pylith::faults::CohesiveDataHex8 : public CohesiveData
+{
+
+// PUBLIC METHODS ///////////////////////////////////////////////////////
+public: 
+
+  /// Constructor
+  CohesiveDataHex8(void);
+
+  /// Destructor
+  ~CohesiveDataHex8(void);
+
+// PRIVATE MEMBERS //////////////////////////////////////////////////////
+private:
+
+  static const int _numVertices; ///< Number of vertices
+  static const int _spaceDim; ///< Number of dimensions in vertex coordinates
+  static const int _numCells; ///< Number of cells
+  static const int _cellDim; ///< Number of dimensions associated with cell
+
+  static const double _vertices[]; ///< Pointer to coordinates of vertices
+  static const int _numCorners[]; ///< Number of vertices in cell
+  static const int _cells[]; ///< Pointer to indices of vertices in cells
+  static const int _materialIds[]; ///< Pointer to cell material identifiers
+
+  static const int _groups[]; ///< Groups of points
+  static const int _groupSizes[]; ///< Sizes of groups
+  static const char* _groupNames[]; ///< Array of group names
+  static const char* _groupTypes[]; ///< Array of group types
+  static const int _numGroups; ///< Number of groups
+
+  static const char* _filename; ///< Filename of input mesh
+};
+
+#endif // pylith_faults_cohesivedatahex8_hh
+
+// End of file

Modified: short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataLine2.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataLine2.cc	2007-04-19 21:26:49 UTC (rev 6610)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataLine2.cc	2007-04-19 21:48:51 UTC (rev 6611)
@@ -12,10 +12,12 @@
 
 /* Original mesh
  *
- *   0 -------- 1 -------- 2
+ *   Cells are 0-1, vertices are 2-4.
  *
+ *   2 -------- 3 -------- 4
+ *
  * After adding cohesive elements
- *   0 -------- 1 -- 3 -------- 2
+ *   2 -------- 3 -- 5 -------- 4
  */
 
 #include "CohesiveDataLine2.hh"
@@ -42,9 +44,9 @@
 };
 
 const int pylith::faults::CohesiveDataLine2::_cells[] = {
-       0,  1,
-       3,  2,
-       1,  3,
+       2,  3,
+       5,  4,
+       3,  5,
 };
 
 const int pylith::faults::CohesiveDataLine2::_materialIds[] = {
@@ -55,11 +57,11 @@
 const int pylith::faults::CohesiveDataLine2::_numGroups = 2;
 
 const int pylith::faults::CohesiveDataLine2::_groupSizes[] = 
-  { 1, 2 };
+  { 2, 3 };
 
 const int pylith::faults::CohesiveDataLine2::_groups[] = {
-  1,
-  0, 1
+  3, 5,
+  2, 3, 5
 };
 
 const char* pylith::faults::CohesiveDataLine2::_groupNames[] = {

Modified: short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataQuad4.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataQuad4.cc	2007-04-19 21:26:49 UTC (rev 6610)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataQuad4.cc	2007-04-19 21:48:51 UTC (rev 6611)
@@ -12,7 +12,9 @@
 
 /* Original mesh
  *
- *       1 -------- 3 -------- 5
+ * Cells are 0-1, vertices are 2-7.
+ *
+ *       3 -------- 5 -------- 7
  *       |          |          |
  *       |          |          |
  *       |          |          |
@@ -21,11 +23,13 @@
  *       |          |          |
  *       |          |          |
  *       |          |          |
- *       0 -------- 2 -------- 4
+ *       2 -------- 4 -------- 6
  *
  * After adding cohesive elements
  *
- *       1 -------- 3 -- 7 -------- 5
+ * Cells are 0-1,10 vertices are 2-9.
+ *
+ *       3 -------- 5 -- 9 -------- 7
  *       |          |    |          |
  *       |          |    |          |
  *       |          |    |          |
@@ -34,7 +38,7 @@
  *       |          |    |          |
  *       |          |    |          |
  *       |          |    |          |
- *       0 -------- 2 -- 6 -------- 4
+ *       2 -------- 4 -- 8 -------- 6
  */
 
 #include "CohesiveDataQuad4.hh"
@@ -65,9 +69,9 @@
 };
 
 const int pylith::faults::CohesiveDataQuad4::_cells[] = {
-  0,  2,  3,  1,
-  4,  5,  7,  6,
-  2,  3,  6,  7,
+  2,  4,  5,  3,
+  6,  7,  9,  8,
+  4,  5,  8,  9,
 };
 
 const int pylith::faults::CohesiveDataQuad4::_materialIds[] = {
@@ -78,11 +82,11 @@
 const int pylith::faults::CohesiveDataQuad4::_numGroups = 2;
 
 const int pylith::faults::CohesiveDataQuad4::_groupSizes[] = 
-  { 2, 3 };
+  { 4, 4 };
 
 const int pylith::faults::CohesiveDataQuad4::_groups[] = {
-  2, 3,
-  1, 3, 5
+  4, 5, 8, 9,
+  3, 5, 7, 9
 };
 
 const char* pylith::faults::CohesiveDataQuad4::_groupNames[] = {

Modified: short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataTet4.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataTet4.cc	2007-04-19 21:26:49 UTC (rev 6610)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataTet4.cc	2007-04-19 21:48:51 UTC (rev 6611)
@@ -12,14 +12,18 @@
 
 /* Original mesh
  *
- * 0   1,2,3  4
+ * Cells are 0-1, vertices are 2-6.
  *
+ * 2   3,4,5  6
+ *
  *     ^^^^^ Face in x-y plane
  *
  * After adding cohesive elements
  *
- * 0   1,2,3  5,6,7   4
+ * Cells are 0-1,10, vertices are 2-9.
  *
+ * 2   3,4,5  7,8,9   6
+ *
  *     ^^^^^^^^^^^^ Cohesive element in x-y plane.
  */
 
@@ -51,9 +55,9 @@
 };
 
 const int pylith::faults::CohesiveDataTet4::_cells[] = {
-  1,  2,  3,  0,
-  5,  7,  6,  4,
-  1,  3,  2,  5,  7,  6
+  3,  4,  5,  2,
+  7,  9,  8,  6,
+  3,  4,  5,  7,  8,  9
 };
 
 const int pylith::faults::CohesiveDataTet4::_materialIds[] = {
@@ -64,11 +68,11 @@
 const int pylith::faults::CohesiveDataTet4::_numGroups = 2;
 
 const int pylith::faults::CohesiveDataTet4::_groupSizes[] = 
-  { 3, 3 };
+  { 6, 5 };
 
 const int pylith::faults::CohesiveDataTet4::_groups[] = {
-  1, 2, 3,
-  0, 2, 3
+  3, 4, 5, 7, 8, 9,
+  2, 4, 5, 8, 9
 };
 
 const char* pylith::faults::CohesiveDataTet4::_groupNames[] = {

Modified: short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataTri3.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataTri3.cc	2007-04-19 21:26:49 UTC (rev 6610)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/data/CohesiveDataTri3.cc	2007-04-19 21:48:51 UTC (rev 6611)
@@ -12,32 +12,36 @@
 
 /* Original mesh
  *
- *              1
+ * Cells are 0-1, vertices are 2-5.
+ *
+ *              3
  *             /|\
  *            / | \
  *           /  |  \
  *          /   |   \
- *         0    |    3
+ *         2    |    5
  *          \   |   /
  *           \  |  /
  *            \ | /
  *             \|/
- *              2
+ *              4
  *
  *
  * After adding cohesive elements
  *
- *              1 -- 4
+ * Cells are 0-1, 8, vertices are 2-7.
+ *
+ *              3 -- 6
  *             /|    |\
  *            / |    | \
  *           /  |    |  \
  *          /   |    |   \
- *         0    |    |    3
+ *         2    |    |    5
  *          \   |    |   /
  *           \  |    |  /
  *            \ |    | /
  *             \|    |/
- *              2 -- 5
+ *              4 -- 7
  */
 
 #include "CohesiveDataTri3.hh"
@@ -66,9 +70,9 @@
 };
 
 const int pylith::faults::CohesiveDataTri3::_cells[] = {
-  0,  2,  1,
-  4,  5,  3,
-  1,  2,  4, 5
+  2,  4,  3,
+  6,  7,  5,
+  3,  4,  6, 7
 };
 
 const int pylith::faults::CohesiveDataTri3::_materialIds[] = {
@@ -79,11 +83,11 @@
 const int pylith::faults::CohesiveDataTri3::_numGroups = 2;
 
 const int pylith::faults::CohesiveDataTri3::_groupSizes[] = 
-  { 2, 3 };
+  { 4, 5 };
 
 const int pylith::faults::CohesiveDataTri3::_groups[] = {
-  1, 2,
-  1, 2, 3
+  3, 4, 6, 7,
+  3, 4, 5, 6, 7
 };
 
 const char* pylith::faults::CohesiveDataTri3::_groupNames[] = {

Modified: short/3D/PyLith/trunk/unittests/libtests/faults/data/Makefile.am
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/data/Makefile.am	2007-04-19 21:26:49 UTC (rev 6610)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/data/Makefile.am	2007-04-19 21:48:51 UTC (rev 6611)
@@ -14,7 +14,8 @@
 	meshLine.txt \
 	meshTri3A.txt \
 	meshQuad4A.txt \
-	meshTet4A.txt
+	meshTet4A.txt \
+	meshHex8A.txt
 
 noinst_TMP =
 

Added: short/3D/PyLith/trunk/unittests/libtests/faults/data/meshHex8A.txt
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/data/meshHex8A.txt	2007-04-19 21:26:49 UTC (rev 6610)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/data/meshHex8A.txt	2007-04-19 21:48:51 UTC (rev 6611)
@@ -0,0 +1,58 @@
+mesh = {
+  dimension = 3
+  use-index-zero = true
+  vertices = {
+    dimension = 3
+    count = 12
+    coordinates = {
+             0     -2.0 -1.0 -1.0
+             1     -2.0  1.0 -1.0
+             2     -2.0 -1.0  1.0
+             3     -2.0  1.0  1.0
+             4      0.0 -1.0 -1.0
+             5      0.0  1.0 -1.0
+             6      0.0 -1.0  1.0
+             7      0.0  1.0  1.0
+             8      2.0 -1.0 -1.0
+             9      2.0  1.0 -1.0
+            10      2.0 -1.0  1.0
+            11      2.0  1.0  1.0
+    }
+  }
+  cells = {
+    count = 2
+    num-corners = 8
+    simplices = {
+             0       0  2  3  1  4  6  7  5
+             1       4  6  7  5  8 10 11  9
+    }
+    material-ids = {
+             0   0
+             1   0
+    }
+  }
+  group = {
+    name = fault
+    type = vertices
+    count = 4
+    indices = {
+      4
+      5
+      6
+      7
+    }
+  }
+  group = {
+    name = output
+    type = vertices
+    count = 6
+    indices = {
+      2
+      3
+      6
+      7
+     10
+     11
+    }
+  }
+}



More information about the cig-commits mailing list