[cig-commits] r9264 - in short/3D/PyLith/trunk: libsrc/meshio pylith/faults tests/2d tests/2d/twofaults_quad4

brad at geodynamics.org brad at geodynamics.org
Wed Feb 6 17:59:58 PST 2008


Author: brad
Date: 2008-02-06 17:59:58 -0800 (Wed, 06 Feb 2008)
New Revision: 9264

Added:
   short/3D/PyLith/trunk/tests/2d/twofaults_quad4/
   short/3D/PyLith/trunk/tests/2d/twofaults_quad4/dislocation_para.cfg
   short/3D/PyLith/trunk/tests/2d/twofaults_quad4/dislocation_perp.cfg
   short/3D/PyLith/trunk/tests/2d/twofaults_quad4/quad4_para.mesh
   short/3D/PyLith/trunk/tests/2d/twofaults_quad4/quad4_perp.mesh
Removed:
   short/3D/PyLith/trunk/tests/2d/quad4b/
   short/3D/PyLith/trunk/tests/2d/twofaults_quad4/dislocation.cfg
   short/3D/PyLith/trunk/tests/2d/twofaults_quad4/quad4h.mesh
Modified:
   short/3D/PyLith/trunk/libsrc/meshio/DataWriterVTK.cc
   short/3D/PyLith/trunk/pylith/faults/BruneSlipFn.py
   short/3D/PyLith/trunk/pylith/faults/Fault.py
Log:
Worked on another two fault example.

Modified: short/3D/PyLith/trunk/libsrc/meshio/DataWriterVTK.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/meshio/DataWriterVTK.cc	2008-02-06 23:47:36 UTC (rev 9263)
+++ short/3D/PyLith/trunk/libsrc/meshio/DataWriterVTK.cc	2008-02-07 01:59:58 UTC (rev 9264)
@@ -151,7 +151,8 @@
     } else {
       labelName = "depth";
     }
-    const ALE::Obj<Mesh::numbering_type>& numbering = mesh->getFactory()->getNumbering(mesh, labelName, 0);
+    const ALE::Obj<Mesh::numbering_type>& numbering =
+      mesh->getFactory()->getNumbering(mesh, labelName, 0);
 
     const int fiberDim = 
       field->getFiberDimension(*mesh->getLabelStratum(labelName, 0)->begin());

Modified: short/3D/PyLith/trunk/pylith/faults/BruneSlipFn.py
===================================================================
--- short/3D/PyLith/trunk/pylith/faults/BruneSlipFn.py	2008-02-06 23:47:36 UTC (rev 9263)
+++ short/3D/PyLith/trunk/pylith/faults/BruneSlipFn.py	2008-02-07 01:59:58 UTC (rev 9264)
@@ -51,17 +51,15 @@
     from spatialdata.spatialdb.SimpleDB import SimpleDB
 
     slip = pyre.inventory.facility("slip", family="spatial_database",
-                                   factory=SimpleDB, args=["slip"])
+                                   factory=SimpleDB)
     slip.meta['tip'] = "Spatial database of slip."
 
     slipTime = pyre.inventory.facility("slip_time", family="spatial_database",
-                                       factory=SimpleDB,
-                                       args=["slip time"])
+                                       factory=SimpleDB)
     slipTime.meta['tip'] = "Spatial database of slip initiation time."
 
     slipRate = pyre.inventory.facility("slip_rate", family="spatial_database",
-                                       factory=SimpleDB,
-                                       args=["slip rate"])
+                                       factory=SimpleDB)
     slipRate.meta['tip'] = "Spatial database of peak slip rate."
 
 

Modified: short/3D/PyLith/trunk/pylith/faults/Fault.py
===================================================================
--- short/3D/PyLith/trunk/pylith/faults/Fault.py	2008-02-06 23:47:36 UTC (rev 9263)
+++ short/3D/PyLith/trunk/pylith/faults/Fault.py	2008-02-07 01:59:58 UTC (rev 9264)
@@ -103,7 +103,7 @@
 
     from spatialdata.spatialdb.SimpleDB import SimpleDB
     matDB = pyre.inventory.facility("mat_db", family="spatial_database",
-                                   factory=SimpleDB, args=["bulk materials"])
+                                   factory=SimpleDB)
     matDB.meta['tip'] = "Spatial database for bulk material properties " \
                         "(used in improving conditioning of Jacobian matrix)."
 

Copied: short/3D/PyLith/trunk/tests/2d/twofaults_quad4 (from rev 9263, short/3D/PyLith/trunk/tests/2d/quad4b)

Deleted: short/3D/PyLith/trunk/tests/2d/twofaults_quad4/dislocation.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4b/dislocation.cfg	2008-02-06 23:47:36 UTC (rev 9263)
+++ short/3D/PyLith/trunk/tests/2d/twofaults_quad4/dislocation.cfg	2008-02-07 01:59:58 UTC (rev 9264)
@@ -1,123 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# journal
-# ----------------------------------------------------------------------
-[pylithapp.journal.info]
-timedependent = 1
-explicit = 1
-implicit = 1
-petsc = 1
-solverlinear = 1
-meshioascii = 1
-homogeneous = 1
-implicitelasticity = 1
-quadrature1d = 1
-fiatsimplex = 1
-faultcohesivekin = 1
-
-# ----------------------------------------------------------------------
-# mesh_generator
-# ----------------------------------------------------------------------
-[pylithapp.mesh_generator]
-debug = 1
-
-[pylithapp.mesh_generator.importer]
-filename = quad4h.mesh
-coordsys.space_dim = 2
-
-# ----------------------------------------------------------------------
-# problem
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-total_time = 0.0*s
-default_dt = 1.0*s
-dimension = 2
-formulation = pylith.problems.Implicit
-bc = pylith.bc.BCFourSides
-interfaces = faults
-#interfaces = pylith.faults.SingleFault
-materials = materials
-
-# ----------------------------------------------------------------------
-# materials
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.materials.one]
-label = elastic material A
-id = 10
-db.iohandler.filename = matprops.spatialdb
-quadrature = pylith.feassemble.quadrature.Quadrature2D
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-
-[pylithapp.timedependent.materials.two]
-label = elastic material B
-id = 11
-db.iohandler.filename = matprops.spatialdb
-quadrature = pylith.feassemble.quadrature.Quadrature2D
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-
-[pylithapp.timedependent.materials.three]
-label = elastic material C
-id = 12
-db.iohandler.filename = matprops.spatialdb
-quadrature = pylith.feassemble.quadrature.Quadrature2D
-quadrature.cell = pylith.feassemble.FIATLagrange
-quadrature.cell.dimension = 2
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.bc.x_neg]
-fixed_dof = [0, 1]
-label = corners
-db.label = Dirichlet BC on corners
-db.iohandler.filename = dislocation_disp.spatialdb
-db.query_type = nearest
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.interfaces.one]
-id = 1
-label = faultA
-quadrature = pylith.feassemble.quadrature.Quadrature1Din2D
-quadrature.cell.shape = line
-mat_db.iohandler.filename = matprops.spatialdb
-
-[pylithapp.timedependent.interfaces.one.eq_src.slip_function]
-slip.query_type = nearest
-slip.iohandler.filename = dislocation_slipA.spatialdb
-slip.label = slip for fault A
-slip_rate.iohandler.filename = dislocation_sliprate.spatialdb
-slip_time.iohandler.filename = dislocation_sliptime.spatialdb
-
-[pylithapp.timedependent.interfaces.two]
-id = 2
-label = faultB
-quadrature = pylith.feassemble.quadrature.Quadrature1Din2D
-quadrature.cell.shape = line
-mat_db.iohandler.filename = matprops.spatialdb
-
-[pylithapp.timedependent.interfaces.two.eq_src.slip_function]
-slip.label = slip for fault B
-slip.iohandler.filename = dislocation_slipB.spatialdb
-slip_rate.iohandler.filename = dislocation_sliprate.spatialdb
-slip_time.iohandler.filename = dislocation_sliptime.spatialdb
-
-
-# ----------------------------------------------------------------------
-# PETSc
-# ----------------------------------------------------------------------
-[pylithapp.petsc]
-pc_type = asm
-ksp_rtol = 1.0e-8
-ksp_gmres_restart = 50
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-[pylithapp.problem.formulation.output.output.writer]
-filename = dislocation.vtk

Added: short/3D/PyLith/trunk/tests/2d/twofaults_quad4/dislocation_para.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/twofaults_quad4/dislocation_para.cfg	                        (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/twofaults_quad4/dislocation_para.cfg	2008-02-07 01:59:58 UTC (rev 9264)
@@ -0,0 +1,159 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+timedependent = 1
+implicit = 1
+petsc = 1
+solverlinear = 1
+meshioascii = 1
+homogeneous = 1
+implicitelasticity = 1
+quadrature1d = 1
+fiatlagrange = 1
+faultcohesivekin = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+#debug = 1
+
+[pylithapp.mesh_generator.importer]
+filename = quad4_para.mesh
+coordsys.space_dim = 2
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+total_time = 0.0*s
+default_dt = 1.0*s
+dimension = 2
+formulation = pylith.problems.Implicit
+bc = pylith.bc.BCTwoSides
+interfaces = faults
+materials = materials
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials.one]
+label = elastic material A
+id = 10
+quadrature = pylith.feassemble.quadrature.Quadrature2D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+db = spatialdata.spatialdb.UniformDB
+db.values = [vp,vs,density]
+db.data = [5291.502622129181, 3000.0, 2500.0]
+
+[pylithapp.timedependent.materials.two]
+label = elastic material B
+id = 11
+quadrature = pylith.feassemble.quadrature.Quadrature2D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+db = spatialdata.spatialdb.UniformDB
+db.values = [vp,vs,density]
+db.data = [5291.502622129181, 3000.0, 2500.0]
+
+[pylithapp.timedependent.materials.three]
+label = elastic material C
+id = 12
+quadrature = pylith.feassemble.quadrature.Quadrature2D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+db = spatialdata.spatialdb.UniformDB
+db.values = [vp,vs,density]
+db.data = [5291.502622129181, 3000.0, 2500.0]
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.neg]
+fixed_dof = [0, 1]
+label = edge -x
+db.label = Dirichlet BC on corners
+
+[pylithapp.timedependent.bc.pos]
+fixed_dof = [0, 1]
+label = edge +x
+db.label = Dirichlet BC on corners
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.one]
+id = 100
+label = faultL
+quadrature = pylith.feassemble.quadrature.Quadrature1Din2D
+quadrature.cell.shape = line
+mat_db = spatialdata.spatialdb.UniformDB
+mat_db.values = [vp,vs,density]
+mat_db.data = [5291.502622129181, 3000.0, 2500.0]
+
+[pylithapp.timedependent.interfaces.one.eq_src.slip_function]
+slip = spatialdata.spatialdb.UniformDB
+slip.values = [left-lateral-slip,fault-opening]
+slip.data = [0.2, 0.0]
+slip.label = Final slip for faultL
+
+slip_rate = spatialdata.spatialdb.UniformDB
+slip_rate.values = [slip-rate]
+slip_rate.data = [1.0e+6]
+slip_rate.label = Peak slip rate for faultL
+
+slip_time = spatialdata.spatialdb.UniformDB
+slip_time.values = [slip-time]
+slip_time.data = [-1.0]
+slip_time.label = Slip time for faultL
+
+[pylithapp.timedependent.interfaces.two]
+id = 101
+label = faultR
+quadrature = pylith.feassemble.quadrature.Quadrature1Din2D
+quadrature.cell.shape = line
+mat_db = spatialdata.spatialdb.UniformDB
+mat_db.values = [vp,vs,density]
+mat_db.data = [5291.502622129181, 3000.0, 2500.0]
+
+[pylithapp.timedependent.interfaces.two.eq_src.slip_function]
+slip = spatialdata.spatialdb.UniformDB
+slip.values = [left-lateral-slip,fault-opening]
+slip.data = [-0.2, 0.0]
+slip.label = Final slip for faultL
+
+slip_rate = spatialdata.spatialdb.UniformDB
+slip_rate.values = [slip-rate]
+slip_rate.data = [1.0e+6]
+slip_rate.label = Peak slip rate for faultL
+
+slip_time = spatialdata.spatialdb.UniformDB
+slip_time.values = [slip-time]
+slip_time.data = [-1.0]
+slip_time.label = Slip time for faultL
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+pc_type = asm
+ksp_rtol = 1.0e-8
+ksp_gmres_restart = 50
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output.writer]
+filename = dislocation-para.vtk
+
+[pylithapp.problem.interfaces.one.output.writer]
+filename = dislocation-para-faultL.vtk
+
+[pylithapp.problem.interfaces.two.output.writer]
+filename = dislocation-para-faultR.vtk

Added: short/3D/PyLith/trunk/tests/2d/twofaults_quad4/dislocation_perp.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/twofaults_quad4/dislocation_perp.cfg	                        (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/twofaults_quad4/dislocation_perp.cfg	2008-02-07 01:59:58 UTC (rev 9264)
@@ -0,0 +1,123 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+timedependent = 1
+explicit = 1
+implicit = 1
+petsc = 1
+solverlinear = 1
+meshioascii = 1
+homogeneous = 1
+implicitelasticity = 1
+quadrature1d = 1
+fiatsimplex = 1
+faultcohesivekin = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+debug = 1
+
+[pylithapp.mesh_generator.importer]
+filename = quad4h.mesh
+coordsys.space_dim = 2
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+total_time = 0.0*s
+default_dt = 1.0*s
+dimension = 2
+formulation = pylith.problems.Implicit
+bc = pylith.bc.BCFourSides
+interfaces = faults
+#interfaces = pylith.faults.SingleFault
+materials = materials
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials.one]
+label = elastic material A
+id = 10
+db.iohandler.filename = matprops.spatialdb
+quadrature = pylith.feassemble.quadrature.Quadrature2D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+
+[pylithapp.timedependent.materials.two]
+label = elastic material B
+id = 11
+db.iohandler.filename = matprops.spatialdb
+quadrature = pylith.feassemble.quadrature.Quadrature2D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+
+[pylithapp.timedependent.materials.three]
+label = elastic material C
+id = 12
+db.iohandler.filename = matprops.spatialdb
+quadrature = pylith.feassemble.quadrature.Quadrature2D
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_neg]
+fixed_dof = [0, 1]
+label = corners
+db.label = Dirichlet BC on corners
+db.iohandler.filename = dislocation_disp.spatialdb
+db.query_type = nearest
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces.one]
+id = 1
+label = faultA
+quadrature = pylith.feassemble.quadrature.Quadrature1Din2D
+quadrature.cell.shape = line
+mat_db.iohandler.filename = matprops.spatialdb
+
+[pylithapp.timedependent.interfaces.one.eq_src.slip_function]
+slip.query_type = nearest
+slip.iohandler.filename = dislocation_slipA.spatialdb
+slip.label = slip for fault A
+slip_rate.iohandler.filename = dislocation_sliprate.spatialdb
+slip_time.iohandler.filename = dislocation_sliptime.spatialdb
+
+[pylithapp.timedependent.interfaces.two]
+id = 2
+label = faultB
+quadrature = pylith.feassemble.quadrature.Quadrature1Din2D
+quadrature.cell.shape = line
+mat_db.iohandler.filename = matprops.spatialdb
+
+[pylithapp.timedependent.interfaces.two.eq_src.slip_function]
+slip.label = slip for fault B
+slip.iohandler.filename = dislocation_slipB.spatialdb
+slip_rate.iohandler.filename = dislocation_sliprate.spatialdb
+slip_time.iohandler.filename = dislocation_sliptime.spatialdb
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+pc_type = asm
+ksp_rtol = 1.0e-8
+ksp_gmres_restart = 50
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output.writer]
+filename = dislocation.vtk

Added: short/3D/PyLith/trunk/tests/2d/twofaults_quad4/quad4_para.mesh
===================================================================
--- short/3D/PyLith/trunk/tests/2d/twofaults_quad4/quad4_para.mesh	                        (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/twofaults_quad4/quad4_para.mesh	2008-02-07 01:59:58 UTC (rev 9264)
@@ -0,0 +1,96 @@
+mesh = {
+  dimension = 2
+  use-index-zero = true
+  vertices = {
+    dimension = 2
+    count = 16
+    coordinates = {
+             0     -3.0   3.0
+             1     -1.0   3.0
+             2      1.0   3.0
+             3      3.0   3.0
+             4     -3.0   1.0
+             5     -1.0   1.0
+             6      1.0   1.0
+             7      3.0   1.0
+             8     -3.0  -1.0
+             9     -1.0  -1.0
+            10      1.0  -1.0
+            11      3.0  -1.0
+            12     -3.0  -3.0
+            13     -1.0  -3.0
+            14      1.0  -3.0
+            15      3.0  -3.0
+    }
+  }
+  cells = {
+    count = 9
+    num-corners = 4
+    simplices = {
+             0       0  4  5  1
+             1       1  5  6  2
+             2       2  6  7  3
+             3       4  8  9  5
+             4       5  9 10  6
+             5       6 10 11  7
+             6       8 12 13  9
+             7       9 13 14 10
+             8      10 14 15 11
+    }
+    material-ids = {
+             0   10
+             1   11
+             2   12
+             3   10
+             4   11
+             5   12
+             6   10
+             7   11
+             8   12
+    }
+  }
+  group = {
+    name = faultL
+    type = vertices
+    count = 4
+    indices = {
+      1
+      5
+      9
+     13
+    }
+  }
+  group = {
+    name = faultR
+    type = vertices
+    count = 4
+    indices = {
+      2
+      6
+     10
+     14
+    }
+  }
+  group = {
+    name = edge -x
+    type = vertices
+    count = 4
+    indices = {
+      0
+      4
+      8
+     12
+    }
+  }
+  group = {
+    name = edge +x
+    type = vertices
+    count = 4
+    indices = {
+      3
+      7
+     11
+     15
+    }
+  }
+}

Added: short/3D/PyLith/trunk/tests/2d/twofaults_quad4/quad4_perp.mesh
===================================================================
--- short/3D/PyLith/trunk/tests/2d/twofaults_quad4/quad4_perp.mesh	                        (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/twofaults_quad4/quad4_perp.mesh	2008-02-07 01:59:58 UTC (rev 9264)
@@ -0,0 +1,82 @@
+mesh = {
+  dimension = 2
+  use-index-zero = true
+  vertices = {
+    dimension = 2
+    count = 16
+    coordinates = {
+             0     -3.0   3.0
+             1     -1.0   3.0
+             2      1.0   3.0
+             3      3.0   3.0
+             4     -3.0   1.0
+             5     -1.0   1.0
+             6      1.0   1.0
+             7      3.0   1.0
+             8     -3.0  -1.0
+             9     -1.0  -1.0
+            10      1.0  -1.0
+            11      3.0  -1.0
+            12     -3.0  -3.0
+            13     -1.0  -3.0
+            14      1.0  -3.0
+            15      3.0  -3.0
+    }
+  }
+  cells = {
+    count = 9
+    num-corners = 4
+    simplices = {
+             0       0  4  5  1
+             1       1  5  6  2
+             2       2  6  7  3
+             3       4  8  9  5
+             4       5  9 10  6
+             5       6 10 11  7
+             6       8 12 13  9
+             7       9 13 14 10
+             8      10 14 15 11
+    }
+    material-ids = {
+             0   10
+             1   10
+             2   11
+             3   10
+             4   10
+             5   11
+             6   12
+             7   12
+             8   11
+    }
+  }
+  group = {
+    name = faultA
+    type = vertices
+    count = 3
+    indices = {
+      2
+      6
+      10
+    }
+  }
+  group = {
+    name = faultB
+    type = vertices
+    count = 2
+    indices = {
+      8
+      9
+    }
+  }
+  group = {
+    name = corners
+    type = vertices
+    count = 4
+    indices = {
+      0
+      3
+     12
+     15
+    }
+  }
+}

Deleted: short/3D/PyLith/trunk/tests/2d/twofaults_quad4/quad4h.mesh
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4b/quad4h.mesh	2008-02-06 23:47:36 UTC (rev 9263)
+++ short/3D/PyLith/trunk/tests/2d/twofaults_quad4/quad4h.mesh	2008-02-07 01:59:58 UTC (rev 9264)
@@ -1,82 +0,0 @@
-mesh = {
-  dimension = 2
-  use-index-zero = true
-  vertices = {
-    dimension = 2
-    count = 16
-    coordinates = {
-             0     -3.0   3.0
-             1     -1.0   3.0
-             2      1.0   3.0
-             3      3.0   3.0
-             4     -3.0   1.0
-             5     -1.0   1.0
-             6      1.0   1.0
-             7      3.0   1.0
-             8     -3.0  -1.0
-             9     -1.0  -1.0
-            10      1.0  -1.0
-            11      3.0  -1.0
-            12     -3.0  -3.0
-            13     -1.0  -3.0
-            14      1.0  -3.0
-            15      3.0  -3.0
-    }
-  }
-  cells = {
-    count = 9
-    num-corners = 4
-    simplices = {
-             0       0  4  5  1
-             1       1  5  6  2
-             2       2  6  7  3
-             3       4  8  9  5
-             4       5  9 10  6
-             5       6 10 11  7
-             6       8 12 13  9
-             7       9 13 14 10
-             8      10 14 15 11
-    }
-    material-ids = {
-             0   10
-             1   10
-             2   11
-             3   10
-             4   10
-             5   11
-             6   12
-             7   12
-             8   11
-    }
-  }
-  group = {
-    name = faultA
-    type = vertices
-    count = 3
-    indices = {
-      2
-      6
-      10
-    }
-  }
-  group = {
-    name = faultB
-    type = vertices
-    count = 2
-    indices = {
-      8
-      9
-    }
-  }
-  group = {
-    name = corners
-    type = vertices
-    count = 4
-    indices = {
-      0
-      3
-     12
-     15
-    }
-  }
-}



More information about the cig-commits mailing list