[cig-commits] r15354 - short/3D/PyLith/trunk/tests/2d/quad4
brad at geodynamics.org
brad at geodynamics.org
Fri Jun 19 12:05:03 PDT 2009
Author: brad
Date: 2009-06-19 12:05:03 -0700 (Fri, 19 Jun 2009)
New Revision: 15354
Added:
short/3D/PyLith/trunk/tests/2d/quad4/fourcells.mesh
short/3D/PyLith/trunk/tests/2d/quad4/fourcells_twofaults.cfg
Log:
Added test case for 2 faults (1 fault each processor).
Added: short/3D/PyLith/trunk/tests/2d/quad4/fourcells.mesh
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/fourcells.mesh (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/quad4/fourcells.mesh 2009-06-19 19:05:03 UTC (rev 15354)
@@ -0,0 +1,91 @@
+// Mesh (file)
+//
+// 1 --- 3 --- 5 --- 7 --- 9
+// | 0 | 1 | 2 | 3 |
+// | | | | |
+// 0 --- 2 --- 4 --- 6 --- 8
+//
+//
+//
+// Sieve mesh (1 fault)
+//
+// 15
+// 3 --- 5-14 --- 7 --- 9 --- 11
+// | 0 | | 1 | 2 | 3 |
+// | | | | |
+// 2 --- 4-12 --- 6 --- 8 --- 10
+// 13
+//
+// 13 and 15 are Lagrange vertices
+
+mesh = {
+ dimension = 2
+ vertices = {
+ dimension = 2
+ count = 10
+ coordinates = {
+ 0 -4.0 -1.0
+ 1 -4.0 +1.0
+ 2 -2.0 -1.0
+ 3 -2.0 +1.0
+ 4 0.0 -1.0
+ 5 0.0 +1.0
+ 6 2.0 -1.0
+ 7 2.0 +1.0
+ 8 4.0 -1.0
+ 9 4.0 +1.0
+ }
+ }
+ cells = {
+ num-corners = 4
+ count = 4
+ simplices = {
+ 0 0 2 3 1
+ 1 2 4 5 3
+ 2 4 6 7 5
+ 3 6 8 9 7
+ }
+ material-ids = {
+ 0 1
+ 1 1
+ 2 1
+ 3 1
+ }
+ }
+ group = {
+ name = edge-x
+ type = vertices
+ count = 2
+ indices = {
+ 0
+ 1
+ }
+ }
+ group = {
+ name = edge+x
+ type = vertices
+ count = 2
+ indices = {
+ 8
+ 9
+ }
+ }
+ group = {
+ name = fault-x
+ type = vertices
+ count = 2
+ indices = {
+ 2
+ 3
+ }
+ }
+ group = {
+ name = fault+x
+ type = vertices
+ count = 2
+ indices = {
+ 6
+ 7
+ }
+ }
+}
Added: short/3D/PyLith/trunk/tests/2d/quad4/fourcells_twofaults.cfg
===================================================================
--- short/3D/PyLith/trunk/tests/2d/quad4/fourcells_twofaults.cfg (rev 0)
+++ short/3D/PyLith/trunk/tests/2d/quad4/fourcells_twofaults.cfg 2009-06-19 19:05:03 UTC (rev 15354)
@@ -0,0 +1,143 @@
+# -*- Python -*-
+[pylithapp]
+
+[pylithapp.launcher] # WARNING: THIS IS NOT PORTABLE
+command = mpirun -np ${nodes}
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+#timedependent = 1
+#implicit = 1
+#petsc = 1
+#solverlinear = 1
+#meshiocubit = 1
+#implicitelasticity = 1
+#quadrature2d = 1
+#fiatlagrange = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+#debug = 1
+
+[pylithapp.mesh_generator.reader]
+filename = fourcells.mesh
+coordsys.space_dim = 2
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+dimension = 2
+bc = [x_neg,x_pos]
+interfaces = [fault1,fault2]
+
+[pylithapp.timedependent.formulation.time_step]
+total_time = 0.0*s
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+materials = [elastic]
+materials.elastic = pylith.materials.ElasticPlaneStrain
+
+[pylithapp.timedependent.materials.elastic]
+label = Elastic material
+id = 1
+db_properties.iohandler.filename = matprops.spatialdb
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 2
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+bc_dof = [0,1]
+label = edge-x
+db_initial.label = Dirichlet BC -x edge
+
+[pylithapp.timedependent.bc.x_neg]
+bc_dof = [0]
+label = edge+x
+db_initial.label = Dirichlet BC +x edge
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.interfaces]
+fault1 = pylith.faults.FaultCohesiveKin
+fault2 = pylith.faults.FaultCohesiveKin
+
+[pylithapp.timedependent.interfaces.fault1]
+id = 2
+label = fault-x
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+
+[pylithapp.timedependent.interfaces.fault1.eq_srcs.rupture.slip_function]
+slip = spatialdata.spatialdb.UniformDB
+slip.label = Final slip
+slip.values = [left-lateral-slip,fault-opening]
+slip.data = [+2.0*m,0.0*m]
+
+slip_time = spatialdata.spatialdb.UniformDB
+slip_time.label = Slip start time
+slip_time.values = [slip-time]
+slip_time.data = [0.0*s]
+
+[pylithapp.timedependent.interfaces.fault2]
+id = 3
+label = fault+x
+quadrature.cell = pylith.feassemble.FIATLagrange
+quadrature.cell.dimension = 1
+
+[pylithapp.timedependent.interfaces.fault2.eq_srcs.rupture.slip_function]
+slip = spatialdata.spatialdb.UniformDB
+slip.label = Final slip
+slip.values = [left-lateral-slip,fault-opening]
+slip.data = [+2.0*m,0.0*m]
+
+slip_time = spatialdata.spatialdb.UniformDB
+slip_time.label = Slip start time
+slip_time.values = [slip-time]
+slip_time.data = [0.0*s]
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_rtol = 1.0e-8
+pc_type = asm
+# Change the preconditioner settings (must turn off
+# shift_positive_definite and turn on shift_nonzero).
+sub_pc_factor_shift_positive_definite = 0
+sub_pc_factor_shift_nonzero =
+
+ksp_max_it = 100
+ksp_gmres_restart = 50
+#ksp_monitor = true
+#ksp_view = true
+#log_summary = true
+# start_in_debugger = true
+
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.output.writer]
+filename = fourcells_onefault.vtk
+
+[pylithapp.timedependent.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvgMesh
+writer.filename = fourcells_onefault-elastic.vtk
+
+[pylithapp.timedependent.interfaces.fault1.output.writer]
+filename = fourcells_onefault-fault1.vtk
+
+[pylithapp.timedependent.interfaces.fault2.output.writer]
+filename = fourcells_onefault-fault2.vtk
More information about the CIG-COMMITS
mailing list