[cig-commits] r16616 - in short/3D/PyLith/trunk/playpen/friction: twohex8 twoquad4 twotet4 twotri3
brad at geodynamics.org
brad at geodynamics.org
Mon May 3 21:01:41 PDT 2010
Author: brad
Date: 2010-05-03 21:01:41 -0700 (Mon, 03 May 2010)
New Revision: 16616
Modified:
short/3D/PyLith/trunk/playpen/friction/twohex8/compression.cfg
short/3D/PyLith/trunk/playpen/friction/twohex8/opening.cfg
short/3D/PyLith/trunk/playpen/friction/twohex8/shear-sliding.cfg
short/3D/PyLith/trunk/playpen/friction/twohex8/shear-stick.cfg
short/3D/PyLith/trunk/playpen/friction/twoquad4/compression.cfg
short/3D/PyLith/trunk/playpen/friction/twoquad4/opening.cfg
short/3D/PyLith/trunk/playpen/friction/twoquad4/shear-sliding.cfg
short/3D/PyLith/trunk/playpen/friction/twoquad4/shear-stick.cfg
short/3D/PyLith/trunk/playpen/friction/twotet4/compression.cfg
short/3D/PyLith/trunk/playpen/friction/twotet4/opening.cfg
short/3D/PyLith/trunk/playpen/friction/twotet4/shear-sliding.cfg
short/3D/PyLith/trunk/playpen/friction/twotet4/shear-stick.cfg
short/3D/PyLith/trunk/playpen/friction/twotri3/compression.cfg
short/3D/PyLith/trunk/playpen/friction/twotri3/opening.cfg
short/3D/PyLith/trunk/playpen/friction/twotri3/shear-sliding.cfg
short/3D/PyLith/trunk/playpen/friction/twotri3/shear-stick.cfg
Log:
Updated playpen/friction parameters.
Modified: short/3D/PyLith/trunk/playpen/friction/twohex8/compression.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twohex8/compression.cfg 2010-05-04 03:28:42 UTC (rev 16615)
+++ short/3D/PyLith/trunk/playpen/friction/twohex8/compression.cfg 2010-05-04 04:01:41 UTC (rev 16616)
@@ -153,7 +153,7 @@
pc_type = asm
# Change the preconditioner settings.
-sub_pc_factor_shift_type = none
+sub_pc_factor_shift_type = nonzero
ksp_rtol = 1.0e-8
ksp_max_it = 100
@@ -171,6 +171,15 @@
#log_summary = true
#start_in_debugger = true
+# Friction sensitivity solve
+friction_pc_type = asm
+friction_sub_pc_factor_shift_type = nonzero
+friction_ksp_max_it = 25
+friction_ksp_gmres_restart = 30
+#friction_ksp_monitor = true
+#friction_ksp_view = true
+friction_ksp_converged_reason = true
+
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/playpen/friction/twohex8/opening.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twohex8/opening.cfg 2010-05-04 03:28:42 UTC (rev 16615)
+++ short/3D/PyLith/trunk/playpen/friction/twohex8/opening.cfg 2010-05-04 04:01:41 UTC (rev 16616)
@@ -150,6 +150,15 @@
#log_summary = true
#start_in_debugger = true
+# Friction sensitivity solve
+friction_pc_type = asm
+friction_sub_pc_factor_shift_type = nonzero
+friction_ksp_max_it = 25
+friction_ksp_gmres_restart = 30
+#friction_ksp_monitor = true
+#friction_ksp_view = true
+friction_ksp_converged_reason = true
+
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/playpen/friction/twohex8/shear-sliding.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twohex8/shear-sliding.cfg 2010-05-04 03:28:42 UTC (rev 16615)
+++ short/3D/PyLith/trunk/playpen/friction/twohex8/shear-sliding.cfg 2010-05-04 04:01:41 UTC (rev 16616)
@@ -232,6 +232,15 @@
#log_summary = true
#start_in_debugger = true
+# Friction sensitivity solve
+friction_pc_type = asm
+friction_sub_pc_factor_shift_type = nonzero
+friction_ksp_max_it = 25
+friction_ksp_gmres_restart = 30
+#friction_ksp_monitor = true
+#friction_ksp_view = true
+friction_ksp_converged_reason = true
+
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/playpen/friction/twohex8/shear-stick.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twohex8/shear-stick.cfg 2010-05-04 03:28:42 UTC (rev 16615)
+++ short/3D/PyLith/trunk/playpen/friction/twohex8/shear-stick.cfg 2010-05-04 04:01:41 UTC (rev 16616)
@@ -230,7 +230,7 @@
pc_type = asm
# Change the preconditioner settings.
-sub_pc_factor_shift_type = none
+sub_pc_factor_shift_type = nonzero
ksp_rtol = 1.0e-8
ksp_max_it = 100
@@ -248,6 +248,15 @@
#log_summary = true
#start_in_debugger = true
+# Friction sensitivity solve
+friction_pc_type = asm
+friction_sub_pc_factor_shift_type = nonzero
+friction_ksp_max_it = 25
+friction_ksp_gmres_restart = 30
+#friction_ksp_monitor = true
+#friction_ksp_view = true
+friction_ksp_converged_reason = true
+
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/playpen/friction/twoquad4/compression.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twoquad4/compression.cfg 2010-05-04 03:28:42 UTC (rev 16615)
+++ short/3D/PyLith/trunk/playpen/friction/twoquad4/compression.cfg 2010-05-04 04:01:41 UTC (rev 16616)
@@ -169,6 +169,15 @@
#log_summary = true
#start_in_debugger = true
+# Friction sensitivity solve
+friction_pc_type = asm
+friction_sub_pc_factor_shift_type = nonzero
+friction_ksp_max_it = 25
+friction_ksp_gmres_restart = 30
+#friction_ksp_monitor = true
+#friction_ksp_view = true
+friction_ksp_converged_reason = true
+
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/playpen/friction/twoquad4/opening.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twoquad4/opening.cfg 2010-05-04 03:28:42 UTC (rev 16615)
+++ short/3D/PyLith/trunk/playpen/friction/twoquad4/opening.cfg 2010-05-04 04:01:41 UTC (rev 16616)
@@ -192,6 +192,15 @@
#log_summary = true
#start_in_debugger = true
+# Friction sensitivity solve
+friction_pc_type = asm
+friction_sub_pc_factor_shift_type = nonzero
+friction_ksp_max_it = 25
+friction_ksp_gmres_restart = 30
+#friction_ksp_monitor = true
+#friction_ksp_view = true
+friction_ksp_converged_reason = true
+
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/playpen/friction/twoquad4/shear-sliding.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twoquad4/shear-sliding.cfg 2010-05-04 03:28:42 UTC (rev 16615)
+++ short/3D/PyLith/trunk/playpen/friction/twoquad4/shear-sliding.cfg 2010-05-04 04:01:41 UTC (rev 16616)
@@ -228,6 +228,15 @@
#log_summary = true
#start_in_debugger = true
+# Friction sensitivity solve
+friction_pc_type = asm
+friction_sub_pc_factor_shift_type = nonzero
+friction_ksp_max_it = 25
+friction_ksp_gmres_restart = 30
+#friction_ksp_monitor = true
+#friction_ksp_view = true
+friction_ksp_converged_reason = true
+
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/playpen/friction/twoquad4/shear-stick.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twoquad4/shear-stick.cfg 2010-05-04 03:28:42 UTC (rev 16615)
+++ short/3D/PyLith/trunk/playpen/friction/twoquad4/shear-stick.cfg 2010-05-04 04:01:41 UTC (rev 16616)
@@ -242,6 +242,15 @@
#log_summary = true
#start_in_debugger = true
+# Friction sensitivity solve
+friction_pc_type = asm
+friction_sub_pc_factor_shift_type = nonzero
+friction_ksp_max_it = 25
+friction_ksp_gmres_restart = 30
+#friction_ksp_monitor = true
+#friction_ksp_view = true
+friction_ksp_converged_reason = true
+
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/playpen/friction/twotet4/compression.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twotet4/compression.cfg 2010-05-04 03:28:42 UTC (rev 16615)
+++ short/3D/PyLith/trunk/playpen/friction/twotet4/compression.cfg 2010-05-04 04:01:41 UTC (rev 16616)
@@ -153,7 +153,7 @@
pc_type = asm
# Change the preconditioner settings.
-sub_pc_factor_shift_type = none
+sub_pc_factor_shift_type = nonzero
ksp_rtol = 1.0e-8
ksp_max_it = 100
@@ -171,6 +171,15 @@
#log_summary = true
#start_in_debugger = true
+# Friction sensitivity solve
+friction_pc_type = asm
+friction_sub_pc_factor_shift_type = nonzero
+friction_ksp_max_it = 25
+friction_ksp_gmres_restart = 30
+#friction_ksp_monitor = true
+#friction_ksp_view = true
+friction_ksp_converged_reason = true
+
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/playpen/friction/twotet4/opening.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twotet4/opening.cfg 2010-05-04 03:28:42 UTC (rev 16615)
+++ short/3D/PyLith/trunk/playpen/friction/twotet4/opening.cfg 2010-05-04 04:01:41 UTC (rev 16616)
@@ -131,7 +131,7 @@
pc_type = asm
# Change the preconditioner settings.
-sub_pc_factor_shift_type = none
+sub_pc_factor_shift_type = nonzero
ksp_rtol = 1.0e-8
ksp_max_it = 100
@@ -150,6 +150,15 @@
#log_summary = true
#start_in_debugger = true
+# Friction sensitivity solve
+friction_pc_type = asm
+friction_sub_pc_factor_shift_type = nonzero
+friction_ksp_max_it = 25
+friction_ksp_gmres_restart = 30
+#friction_ksp_monitor = true
+#friction_ksp_view = true
+friction_ksp_converged_reason = true
+
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/playpen/friction/twotet4/shear-sliding.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twotet4/shear-sliding.cfg 2010-05-04 03:28:42 UTC (rev 16615)
+++ short/3D/PyLith/trunk/playpen/friction/twotet4/shear-sliding.cfg 2010-05-04 04:01:41 UTC (rev 16616)
@@ -215,7 +215,7 @@
pc_type = asm
# Change the preconditioner settings.
-sub_pc_factor_shift_type = none
+sub_pc_factor_shift_type = nonzero
ksp_rtol = 1.0e-8
ksp_max_it = 100
@@ -233,6 +233,15 @@
#log_summary = true
#start_in_debugger = true
+# Friction sensitivity solve
+friction_pc_type = asm
+friction_sub_pc_factor_shift_type = nonzero
+friction_ksp_max_it = 25
+friction_ksp_gmres_restart = 30
+#friction_ksp_monitor = true
+#friction_ksp_view = true
+friction_ksp_converged_reason = true
+
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/playpen/friction/twotet4/shear-stick.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twotet4/shear-stick.cfg 2010-05-04 03:28:42 UTC (rev 16615)
+++ short/3D/PyLith/trunk/playpen/friction/twotet4/shear-stick.cfg 2010-05-04 04:01:41 UTC (rev 16616)
@@ -230,7 +230,7 @@
pc_type = asm
# Change the preconditioner settings.
-sub_pc_factor_shift_type = none
+sub_pc_factor_shift_type = nonzero
ksp_rtol = 1.0e-8
ksp_max_it = 100
@@ -248,6 +248,15 @@
#log_summary = true
#start_in_debugger = true
+# Friction sensitivity solve
+friction_pc_type = asm
+friction_sub_pc_factor_shift_type = nonzero
+friction_ksp_max_it = 25
+friction_ksp_gmres_restart = 30
+#friction_ksp_monitor = true
+#friction_ksp_view = true
+friction_ksp_converged_reason = true
+
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/playpen/friction/twotri3/compression.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twotri3/compression.cfg 2010-05-04 03:28:42 UTC (rev 16615)
+++ short/3D/PyLith/trunk/playpen/friction/twotri3/compression.cfg 2010-05-04 04:01:41 UTC (rev 16616)
@@ -35,11 +35,9 @@
formulation = pylith.problems.Implicit
formulation.solver = pylith.problems.SolverNonlinear
-# Set bc to an array with 3 boundary conditions: 'x_neg', 'y_neg' and 'x_pos'.
bc = [x_neg,y_neg,x_pos]
bc.x_pos = pylith.bc.Neumann
-# Set interfaces to an array with 1 fault: 'fault'.
interfaces = [fault]
@@ -51,58 +49,31 @@
# ----------------------------------------------------------------------
# materials
# ----------------------------------------------------------------------
-# Specify the material information for the problem.
-# The material type is isotropic elastic formulated for plane strain.
[pylithapp.timedependent.materials]
material = pylith.materials.ElasticPlaneStrain
[pylithapp.timedependent.materials.material]
-
-# We give a label of 'elastic material' to this material.
label = elastic material
-
-# The cells associated with this material are given a material ID of 0
-# in the mesh file.
id = 0
-
-# The properties for this material are given in the spatial database file
-# 'matelastic2D.spatialdb'.
db_properties.iohandler.filename = matelastic2D.spatialdb
-# We are doing 2D quadrature for a triangle.
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = triangle
# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
-# Provide information on the boundary conditions.
-
-# Boundary conditions to be applied to the negative x-side of the mesh.
[pylithapp.timedependent.bc.x_neg]
-# We are fixing the 0 (x) degree of freedom.
bc_dof = [0]
-
-# The nodes associated with this boundary condition have the name
-# 'x_neg' in the mesh file.
label = x_neg
-# Boundary conditions to be applied to the negative y-side of the mesh.
[pylithapp.timedependent.bc.y_neg]
-
-# We are fixing the 1 (y) degree of freedom.
bc_dof = [1]
-# The nodes associated with this boundary condition have the name
-# 'x_pos' in the mesh file.
label = y_neg
-# Boundary conditions to be applied to the positive x-side of the mesh.
[pylithapp.timedependent.bc.x_pos]
-
-# The nodes associated with this boundary condition have the name
-# 'x_pos' in the mesh file.
label = x_pos
db_initial = spatialdata.spatialdb.UniformDB
@@ -110,36 +81,18 @@
db_initial.values = [traction-shear,traction-normal]
db_initial.data = [0.0*MPa,-12000*MPa]
-# 1-D simplex cell in 2-D space
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = line
# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------
-# Provide information on the fault (interface).
[pylithapp.timedependent.interfaces]
-
fault = pylith.faults.FaultCohesiveDyn
-# Define fault properties.
[pylithapp.timedependent.interfaces.fault]
-
-# The nodes associated with this fault have the name 'fault' in the mesh file.
label = fault
-# NOTE: It is possible to assign an ID number to a fault (e.g.,
-# 'id = 10'). Care must be taken when doing this, however, because the
-# assigned ID will become the material ID for the cohesive element.
-# This ID must not conflict with any of the material ID numbers for
-# volume elements. The default ID for a fault is 100. If you have a
-# fault in your mesh you must:
-# 1. If you create your own fault ID, make sure it does not conflict
-# with any of you material ID's.
-# 2. If you use the default fault ID, make sure that none of your
-# material ID's are equal to 100.
-
-# The quadrature for a 2D fault is 1D with a linear shape.
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = line
@@ -160,7 +113,7 @@
pc_type = asm
# Change the preconditioner settings.
-sub_pc_factor_shift_type = none
+sub_pc_factor_shift_type = nonzero
ksp_rtol = 1.0e-8
ksp_max_it = 100
@@ -178,6 +131,15 @@
#log_summary = true
#start_in_debugger = true
+# Friction sensitivity solve
+friction_pc_type = asm
+friction_sub_pc_factor_shift_type = nonzero
+friction_ksp_max_it = 25
+friction_ksp_gmres_restart = 30
+#friction_ksp_monitor = true
+#friction_ksp_view = true
+friction_ksp_converged_reason = true
+
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/playpen/friction/twotri3/opening.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twotri3/opening.cfg 2010-05-04 03:28:42 UTC (rev 16615)
+++ short/3D/PyLith/trunk/playpen/friction/twotri3/opening.cfg 2010-05-04 04:01:41 UTC (rev 16616)
@@ -1,22 +1,9 @@
# -*- Python -*-
-
-# The settings in this file (pylithapp.cfg) will be read automatically
-# by pylith, as long as the file is placed in the run directory.
-
-# The settings in this file will override any settings in:
-# PREFIX/etc/pylithapp.cfg
-# $HOME/.pyre/pylithapp/pylithapp.cfg
-
-# The settings in this file will be overridden by any .cfg file given
-# on the command line or by any command line settings.
-
[pylithapp]
# ----------------------------------------------------------------------
# journal
# ----------------------------------------------------------------------
-# The settings below turn on journal info for the specified components.
-# If you want less output to stdout, you can turn these off.
[pylithapp.journal.info]
timedependent = 1
implicit = 1
@@ -32,14 +19,9 @@
# ----------------------------------------------------------------------
# mesh_generator
# ----------------------------------------------------------------------
-# The settings below control the mesh generation (importing mesh info).
-# Turn on debugging output for mesh generation.
[pylithapp.mesh_generator]
-debug = 0
+#debug = 1
-# This component specification means we are using PyLith ASCII format,
-# and we then specify the filename and number of space dimensions for
-# the mesh.
[pylithapp.mesh_generator.reader]
filename = twotri3.mesh
coordsys.space_dim = 2
@@ -47,24 +29,15 @@
# ----------------------------------------------------------------------
# problem
# ----------------------------------------------------------------------
-# Specify the problem settings.
-# This is a time-dependent problem, so we select this as our problem type.
-# We select a total time of 0 sec, and a time step size of 1 sec, so we
-# are performing a single time step.
-# The spatial dimension for this problem is 2.
-# For an implicit formulation (using implicit.cfg), we will perform 1
-# implicit time step from t = -1.0 to t = 0.0 (elastic solution step).
[pylithapp.timedependent]
dimension = 2
normalizer.length_scale = 1.0*m
formulation = pylith.problems.Implicit
formulation.solver = pylith.problems.SolverNonlinear
-# Set bc to an array with 2 boundary conditions: 'x_neg', 'x_pos'
bc = [x_neg,x_pos]
bc.x_pos=pylith.bc.DirichletBoundary
-# Set interfaces to an array with 1 fault: 'fault'.
interfaces = [fault]
@@ -76,51 +49,26 @@
# ----------------------------------------------------------------------
# materials
# ----------------------------------------------------------------------
-# Specify the material information for the problem.
-# The material type is isotropic elastic formulated for plane strain.
[pylithapp.timedependent.materials]
material = pylith.materials.ElasticPlaneStrain
[pylithapp.timedependent.materials.material]
-
-# We give a label of 'elastic material' to this material.
label = elastic material
-
-# The cells associated with this material are given a material ID of 0
-# in the mesh file.
id = 0
-
-# The properties for this material are given in the spatial database file
-# 'matelastic2D.spatialdb'.
db_properties.iohandler.filename = matelastic2D.spatialdb
-# We are doing 2D quadrature for a triangle.
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = triangle
# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
-# Provide information on the boundary conditions.
-
-# Boundary conditions to be applied to the negative x-side of the mesh.
[pylithapp.timedependent.bc.x_neg]
-
-# We are fixing the 0 (x) and 1 (y) degree of freedom.
bc_dof = [0,1]
-
-# The nodes associated with this boundary condition have the name
-# 'x_neg' in the mesh file.
label = x_neg
-# Boundary conditions to be applied to the positive x-side of the mesh.
[pylithapp.timedependent.bc.x_pos]
-
-# We are fixing the 1 (y) and prescribing 0 (x) degree of freedom.
bc_dof = [0,1]
-
-# The nodes associated with this boundary condition have the name
-# 'x_pos' in the mesh file.
label = x_pos
db_initial = spatialdata.spatialdb.UniformDB
@@ -131,29 +79,13 @@
# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------
-# Provide information on the fault (interface).
[pylithapp.timedependent.interfaces]
fault = pylith.faults.FaultCohesiveDyn
-# Define fault properties.
[pylithapp.timedependent.interfaces.fault]
-
-# The nodes associated with this fault have the name 'fault' in the mesh file.
label = fault
-# NOTE: It is possible to assign an ID number to a fault (e.g.,
-# 'id = 10'). Care must be taken when doing this, however, because the
-# assigned ID will become the material ID for the cohesive element.
-# This ID must not conflict with any of the material ID numbers for
-# volume elements. The default ID for a fault is 100. If you have a
-# fault in your mesh you must:
-# 1. If you create your own fault ID, make sure it does not conflict
-# with any of you material ID's.
-# 2. If you use the default fault ID, make sure that none of your
-# material ID's are equal to 100.
-
-# The quadrature for a 2D fault is 1D with a linear shape.
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = line
@@ -174,7 +106,7 @@
pc_type = asm
# Change the preconditioner settings.
-sub_pc_factor_shift_type = none
+sub_pc_factor_shift_type = nonzero
ksp_rtol = 1.0e-8
ksp_max_it = 100
@@ -192,6 +124,15 @@
#log_summary = true
#start_in_debugger = true
+# Friction sensitivity solve
+friction_pc_type = asm
+friction_sub_pc_factor_shift_type = nonzero
+friction_ksp_max_it = 25
+friction_ksp_gmres_restart = 30
+#friction_ksp_monitor = true
+#friction_ksp_view = true
+friction_ksp_converged_reason = true
+
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/playpen/friction/twotri3/shear-sliding.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twotri3/shear-sliding.cfg 2010-05-04 03:28:42 UTC (rev 16615)
+++ short/3D/PyLith/trunk/playpen/friction/twotri3/shear-sliding.cfg 2010-05-04 04:01:41 UTC (rev 16616)
@@ -1,22 +1,9 @@
# -*- Python -*-
-
-# The settings in this file (pylithapp.cfg) will be read automatically
-# by pylith, as long as the file is placed in the run directory.
-
-# The settings in this file will override any settings in:
-# PREFIX/etc/pylithapp.cfg
-# $HOME/.pyre/pylithapp/pylithapp.cfg
-
-# The settings in this file will be overridden by any .cfg file given
-# on the command line or by any command line settings.
-
[pylithapp]
# ----------------------------------------------------------------------
# journal
# ----------------------------------------------------------------------
-# The settings below turn on journal info for the specified components.
-# If you want less output to stdout, you can turn these off.
[pylithapp.journal.info]
timedependent = 1
implicit = 1
@@ -32,14 +19,9 @@
# ----------------------------------------------------------------------
# mesh_generator
# ----------------------------------------------------------------------
-# The settings below control the mesh generation (importing mesh info).
-# Turn on debugging output for mesh generation.
[pylithapp.mesh_generator]
-debug = 0
+#debug = 1
-# This component specification means we are using PyLith ASCII format,
-# and we then specify the filename and number of space dimensions for
-# the mesh.
[pylithapp.mesh_generator.reader]
filename = twotri3.mesh
coordsys.space_dim = 2
@@ -47,29 +29,19 @@
# ----------------------------------------------------------------------
# problem
# ----------------------------------------------------------------------
-# Specify the problem settings.
-# This is a time-dependent problem, so we select this as our problem type.
-# We select a total time of 0 sec, and a time step size of 1 sec, so we
-# are performing a single time step.
-# The spatial dimension for this problem is 2.
-# For an implicit formulation (using implicit.cfg), we will perform 1
-# implicit time step from t = -1.0 to t = 0.0 (elastic solution step).
[pylithapp.timedependent]
dimension = 2
normalizer.length_scale = 1.0*m
formulation = pylith.problems.Implicit
formulation.solver = pylith.problems.SolverNonlinear
-# Set bc to an array with 2 boundary conditions: 'x_neg', 'x_pos'
bc = [x_neg,x_pos,y_pos_fault,y_neg_fault]
bc.x_pos = pylith.bc.DirichletBoundary
bc.y_pos_fault = pylith.bc.Neumann
bc.y_neg_fault = pylith.bc.Neumann
-# Set interfaces to an array with 1 fault: 'fault'.
interfaces = [fault]
-
[pylithapp.timedependent.formulation.time_step]
total_time = 0.0*s
dt = 1.0*s
@@ -78,51 +50,26 @@
# ----------------------------------------------------------------------
# materials
# ----------------------------------------------------------------------
-# Specify the material information for the problem.
-# The material type is isotropic elastic formulated for plane strain.
[pylithapp.timedependent.materials]
material = pylith.materials.ElasticPlaneStrain
[pylithapp.timedependent.materials.material]
-
-# We give a label of 'elastic material' to this material.
label = elastic material
-
-# The cells associated with this material are given a material ID of 0
-# in the mesh file.
id = 0
-
-# The properties for this material are given in the spatial database file
-# 'matelastic2D.spatialdb'.
db_properties.iohandler.filename = matelastic2D.spatialdb
-# We are doing 2D quadrature for a triangle.
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = triangle
# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
-# Provide information on the boundary conditions.
-
-# Boundary conditions to be applied to the negative x-side of the mesh.
[pylithapp.timedependent.bc.x_neg]
-
-# We are fixing the 0 (x) and 1 (y) degrees of freedom.
bc_dof = [0,1]
-
-# The nodes associated with this boundary condition have the name
-# 'x_neg' in the mesh file.
label = x_neg
-# Boundary conditions to be applied to the positive x-side of the mesh.
[pylithapp.timedependent.bc.x_pos]
-
-# We are fixing the 0 (x) and prescribing 1 (y) degree of freedom.
bc_dof = [0,1]
-
-# The nodes associated with this boundary condition have the name
-# 'x_pos' in the mesh file.
label = x_pos
db_initial = spatialdata.spatialdb.UniformDB
@@ -131,12 +78,7 @@
db_initial.data = [0.0*m,1.0*m]
-
-# Boundary conditions to be applied to the positive y-side of the mesh.
[pylithapp.timedependent.bc.y_pos_fault]
-
-# The nodes associated with this boundary condition have the name
-# 'y_pos_fault' in the mesh file.
label = y_pos_fault
db_initial = spatialdata.spatialdb.UniformDB
@@ -144,15 +86,10 @@
db_initial.values = [traction-shear,traction-normal]
db_initial.data = [-6000*MPa, 0.0*MPa]
-# 1-D simplex cell in 2-D space
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = line
-# Boundary conditions to be applied to the negative y-side of the mesh.
[pylithapp.timedependent.bc.y_neg_fault]
-
-# The nodes associated with this boundary condition have the name
-# 'y_neg_fault' in the mesh file.
label = y_neg_fault
db_initial = spatialdata.spatialdb.UniformDB
@@ -160,36 +97,18 @@
db_initial.values = [traction-shear,traction-normal]
db_initial.data = [-6000*MPa, 0.0*MPa]
-# 1-D simplex cell in 2-D space
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = line
# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------
-# Provide information on the fault (interface).
[pylithapp.timedependent.interfaces]
fault = pylith.faults.FaultCohesiveDyn
-# Define fault properties.
[pylithapp.timedependent.interfaces.fault]
-
-# The nodes associated with this fault have the name 'fault' in the mesh file.
label = fault
-
-# NOTE: It is possible to assign an ID number to a fault (e.g.,
-# 'id = 10'). Care must be taken when doing this, however, because the
-# assigned ID will become the material ID for the cohesive element.
-# This ID must not conflict with any of the material ID numbers for
-# volume elements. The default ID for a fault is 100. If you have a
-# fault in your mesh you must:
-# 1. If you create your own fault ID, make sure it does not conflict
-# with any of you material ID's.
-# 2. If you use the default fault ID, make sure that none of your
-# material ID's are equal to 100.
-
-# The quadrature for a 2D fault is 1D with a linear shape.
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = line
@@ -210,7 +129,7 @@
pc_type = asm
# Change the preconditioner settings.
-sub_pc_factor_shift_type = none
+sub_pc_factor_shift_type = nonzero
ksp_rtol = 1.0e-8
ksp_max_it = 100
@@ -228,6 +147,15 @@
#log_summary = true
#start_in_debugger = true
+# Friction sensitivity solve
+friction_pc_type = asm
+friction_sub_pc_factor_shift_type = nonzero
+friction_ksp_max_it = 25
+friction_ksp_gmres_restart = 30
+#friction_ksp_monitor = true
+#friction_ksp_view = true
+friction_ksp_converged_reason = true
+
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/playpen/friction/twotri3/shear-stick.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twotri3/shear-stick.cfg 2010-05-04 03:28:42 UTC (rev 16615)
+++ short/3D/PyLith/trunk/playpen/friction/twotri3/shear-stick.cfg 2010-05-04 04:01:41 UTC (rev 16616)
@@ -1,22 +1,9 @@
# -*- Python -*-
-
-# The settings in this file (pylithapp.cfg) will be read automatically
-# by pylith, as long as the file is placed in the run directory.
-
-# The settings in this file will override any settings in:
-# PREFIX/etc/pylithapp.cfg
-# $HOME/.pyre/pylithapp/pylithapp.cfg
-
-# The settings in this file will be overridden by any .cfg file given
-# on the command line or by any command line settings.
-
[pylithapp]
# ----------------------------------------------------------------------
# journal
# ----------------------------------------------------------------------
-# The settings below turn on journal info for the specified components.
-# If you want less output to stdout, you can turn these off.
[pylithapp.journal.info]
timedependent = 1
implicit = 1
@@ -32,14 +19,9 @@
# ----------------------------------------------------------------------
# mesh_generator
# ----------------------------------------------------------------------
-# The settings below control the mesh generation (importing mesh info).
-# Turn on debugging output for mesh generation.
[pylithapp.mesh_generator]
-debug = 1
+#debug = 1
-# This component specification means we are using PyLith ASCII format,
-# and we then specify the filename and number of space dimensions for
-# the mesh.
[pylithapp.mesh_generator.reader]
filename = twotri3.mesh
coordsys.space_dim = 2
@@ -47,13 +29,6 @@
# ----------------------------------------------------------------------
# problem
# ----------------------------------------------------------------------
-# Specify the problem settings.
-# This is a time-dependent problem, so we select this as our problem type.
-# We select a total time of 0 sec, and a time step size of 1 sec, so we
-# are performing a single time step.
-# The spatial dimension for this problem is 2.
-# For an implicit formulation (using implicit.cfg), we will perform 1
-# implicit time step from t = -1.0 to t = 0.0 (elastic solution step).
[pylithapp.timedependent]
dimension = 2
normalizer.length_scale = 1.0*m
@@ -78,58 +53,28 @@
# ----------------------------------------------------------------------
# materials
# ----------------------------------------------------------------------
-# Specify the material information for the problem.
-# The material type is isotropic elastic formulated for plane strain.
[pylithapp.timedependent.materials]
material = pylith.materials.ElasticPlaneStrain
[pylithapp.timedependent.materials.material]
-
-# We give a label of 'elastic material' to this material.
label = elastic material
-
-# The cells associated with this material are given a material ID of 0
-# in the mesh file.
id = 0
-
-# The properties for this material are given in the spatial database file
-# 'matelastic2D.spatialdb'.
db_properties.iohandler.filename = matelastic2D.spatialdb
-
-# We are doing 2D quadrature for a triangle.
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = triangle
# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
-# Provide information on the boundary conditions.
-
-# Boundary conditions to be applied to the negative x-side of the mesh.
[pylithapp.timedependent.bc.x_neg]
-
-# We are fixing the 0 (x) and 1 (y) degrees of freedom.
bc_dof = [0,1]
-
-# The nodes associated with this boundary condition have the name
-# 'x_neg' in the mesh file.
label = x_neg
-# Boundary conditions to be applied to the positive x-side of the mesh.
[pylithapp.timedependent.bc.x_pos_disp]
-
-# We are fixing the 0 (x) degree of freedom.
bc_dof = [0]
-
-# The nodes associated with this boundary condition have the name
-# 'x_pos' in the mesh file.
label = x_pos
-# Boundary conditions to be applied to the positive x-side of the mesh.
[pylithapp.timedependent.bc.x_pos_tract]
-
-# The nodes associated with this boundary condition have the name
-# 'x_pos' in the mesh file.
label = x_pos
db_initial = spatialdata.spatialdb.UniformDB
@@ -142,13 +87,7 @@
quadrature.cell.shape = line
-
-
-# Boundary conditions to be applied to the positive y-side of the mesh.
[pylithapp.timedependent.bc.y_pos_fault]
-
-# The nodes associated with this boundary condition have the name
-# 'y_pos_fault' in the mesh file.
label = y_pos_fault
db_initial = spatialdata.spatialdb.UniformDB
@@ -156,15 +95,10 @@
db_initial.values = [traction-shear,traction-normal]
db_initial.data = [-1000.0*MPa, 0.0*MPa]
-# 1-D simplex cell in 2-D space
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = line
-# Boundary conditions to be applied to the negative y-side of the mesh.
[pylithapp.timedependent.bc.y_neg_fault]
-
-# The nodes associated with this boundary condition have the name
-# 'y_neg_fault' in the mesh file.
label = y_neg_fault
db_initial = spatialdata.spatialdb.UniformDB
@@ -172,7 +106,6 @@
db_initial.values = [traction-shear,traction-normal]
db_initial.data = [-1000.0*MPa, 0.0*MPa]
-# 1-D simplex cell in 2-D space
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = line
@@ -186,24 +119,8 @@
fault = pylith.faults.FaultCohesiveDyn
-# Define fault properties.
[pylithapp.timedependent.interfaces.fault]
-
-# The nodes associated with this fault have the name 'fault' in the mesh file.
label = fault
-
-# NOTE: It is possible to assign an ID number to a fault (e.g.,
-# 'id = 10'). Care must be taken when doing this, however, because the
-# assigned ID will become the material ID for the cohesive element.
-# This ID must not conflict with any of the material ID numbers for
-# volume elements. The default ID for a fault is 100. If you have a
-# fault in your mesh you must:
-# 1. If you create your own fault ID, make sure it does not conflict
-# with any of you material ID's.
-# 2. If you use the default fault ID, make sure that none of your
-# material ID's are equal to 100.
-
-# The quadrature for a 2D fault is 1D with a linear shape.
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = line
@@ -224,7 +141,7 @@
pc_type = asm
# Change the preconditioner settings.
-sub_pc_factor_shift_type = none
+sub_pc_factor_shift_type = nonzero
ksp_rtol = 1.0e-8
ksp_max_it = 100
@@ -242,6 +159,15 @@
#log_summary = true
#start_in_debugger = true
+# Friction sensitivity solve
+friction_pc_type = asm
+friction_sub_pc_factor_shift_type = nonzero
+friction_ksp_max_it = 25
+friction_ksp_gmres_restart = 30
+#friction_ksp_monitor = true
+#friction_ksp_view = true
+friction_ksp_converged_reason = true
+
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
More information about the CIG-COMMITS
mailing list