[cig-commits] r16455 - in short/3D/PyLith/trunk: . examples/bar_shearwave/hex8 examples/bar_shearwave/quad4 examples/bar_shearwave/tet4 examples/bar_shearwave/tri3 playpen/friction/bar_shearwave/quad4 playpen/friction/twohex8
brad at geodynamics.org
brad at geodynamics.org
Wed Mar 24 17:15:44 PDT 2010
Author: brad
Date: 2010-03-24 17:15:43 -0700 (Wed, 24 Mar 2010)
New Revision: 16455
Modified:
short/3D/PyLith/trunk/TODO
short/3D/PyLith/trunk/examples/bar_shearwave/hex8/pylithapp.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/tet4/pylithapp.cfg
short/3D/PyLith/trunk/examples/bar_shearwave/tri3/pylithapp.cfg
short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4/shearwave.cfg
short/3D/PyLith/trunk/playpen/friction/twohex8/opening.cfg
short/3D/PyLith/trunk/playpen/friction/twohex8/shear-sliding.cfg
Log:
Improved solver settings.
Modified: short/3D/PyLith/trunk/TODO
===================================================================
--- short/3D/PyLith/trunk/TODO 2010-03-25 00:15:03 UTC (rev 16454)
+++ short/3D/PyLith/trunk/TODO 2010-03-25 00:15:43 UTC (rev 16455)
@@ -12,7 +12,6 @@
Need to pass hint when creating Jacobian.
* Friction
- + Fix bug in explicit time-stepping
+ Check output (diagnostic + simulation data)
+ Add stuff to manual
@@ -34,6 +33,10 @@
* Uniform global refinement for tets with faults
+* Move time-stepping dependent code to Explicit.cc and Implicit.cc.
+ + Explicit: compute vel(t), acc(t) at beginning of reformResidual()
+ + Implicit: compute vel(t) at begining of reformResidial()
+
* Cleanup
+ memory model
+ full-scale testing
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/hex8/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/hex8/pylithapp.cfg 2010-03-25 00:15:03 UTC (rev 16454)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/hex8/pylithapp.cfg 2010-03-25 00:15:43 UTC (rev 16455)
@@ -40,8 +40,8 @@
dimension = 3
# Change to an explicit time stepping formulation
-#formulation = pylith.problems.Explicit
-formulation = pylith.problems.ExplicitLumped
+formulation = pylith.problems.Explicit
+#formulation = pylith.problems.ExplicitLumped
# Nondimensionalize problem using wave propagation parameters.
normalizer = spatialdata.units.NondimElasticDynamic
@@ -177,12 +177,12 @@
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_atol = 1.0e-10
ksp_max_it = 50
-ksp_gmres_restart = 10
+ksp_gmres_restart = 75
#ksp_monitor = true
#ksp_view = true
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg 2010-03-25 00:15:03 UTC (rev 16454)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/quad4/pylithapp.cfg 2010-03-25 00:15:43 UTC (rev 16455)
@@ -181,16 +181,16 @@
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_atol = 1.0e-10
+ksp_atol = 1.0e-12
ksp_max_it = 50
-ksp_gmres_restart = 10
+ksp_gmres_restart = 75
-#ksp_monitor = true
-#ksp_view = true
-log_summary = true
+ksp_monitor = true
+ksp_view = true
+#log_summary = true
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/tet4/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/tet4/pylithapp.cfg 2010-03-25 00:15:03 UTC (rev 16454)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/tet4/pylithapp.cfg 2010-03-25 00:15:43 UTC (rev 16455)
@@ -41,9 +41,9 @@
dimension = 3
# Change to an explicit time stepping formulation
-#formulation = pylith.problems.Explicit
+formulation = pylith.problems.Explicit
#formulation = pylith.problems.ExplicitLumped
-formulation = pylith.problems.ExplicitLumpedTet4
+#formulation = pylith.problems.ExplicitLumpedTet4
# Nondimensionalize problem using wave propagation parameters.
normalizer = spatialdata.units.NondimElasticDynamic
@@ -83,7 +83,7 @@
# 3-D simplex cell with 2nd order quadrature
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = tetrahedron
-quadrature.cell.quad_order = 1
+quadrature.cell.quad_order = 2
[pylithapp.timedependent.materials.neg]
@@ -100,7 +100,7 @@
# 3-D simplex cell with 2nd order quadrature
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = tetrahedron
-quadrature.cell.quad_order = 1
+quadrature.cell.quad_order = 2
# ----------------------------------------------------------------------
# boundary conditions
@@ -119,7 +119,7 @@
# 2-D simplex cell in 3-D space with 2nd order quadrature
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = triangle
-quadrature.cell.quad_order = 1
+quadrature.cell.quad_order = 2
[pylithapp.timedependent.bc.x_neg]
# Absorbing boundary condition on +x face of bar
@@ -135,7 +135,7 @@
# 2-D simplex cell in 3-D space with 2nd order quadrature
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = triangle
-quadrature.cell.quad_order = 1
+quadrature.cell.quad_order = 2
[pylithapp.timedependent.bc.y_pos]
# Dirichlet boundary condition on all vertices except fault vertices
@@ -170,7 +170,7 @@
# 2-D simplex cell in 3-D space with 2nd order quadrature
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = triangle
-quadrature.cell.quad_order = 1
+quadrature.cell.quad_order = 2
# Switch to Brune slip time function
eq_srcs.rupture.slip_function = pylith.faults.BruneSlipFn
@@ -194,16 +194,16 @@
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_atol = 1.0e-10
ksp_max_it = 50
-ksp_gmres_restart = 10
+ksp_gmres_restart = 75
-#ksp_monitor = true
+ksp_monitor = true
#ksp_view = true
-log_summary = true
+#log_summary = true
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/examples/bar_shearwave/tri3/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/bar_shearwave/tri3/pylithapp.cfg 2010-03-25 00:15:03 UTC (rev 16454)
+++ short/3D/PyLith/trunk/examples/bar_shearwave/tri3/pylithapp.cfg 2010-03-25 00:15:43 UTC (rev 16455)
@@ -41,8 +41,8 @@
dimension = 2
# Change to an explicit time stepping formulation
-#formulation = pylith.problems.Explicit
-formulation = pylith.problems.ExplicitLumped
+formulation = pylith.problems.Explicit
+#formulation = pylith.problems.ExplicitLumped
# Nondimensionalize problem using wave propagation parameters.
normalizer = spatialdata.units.NondimElasticDynamic
@@ -181,7 +181,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_atol = 1.0e-10
Modified: short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4/shearwave.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4/shearwave.cfg 2010-03-25 00:15:03 UTC (rev 16454)
+++ short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4/shearwave.cfg 2010-03-25 00:15:43 UTC (rev 16455)
@@ -24,7 +24,7 @@
bc.x_neg = pylith.bc.AbsorbingDampers
[pylithapp.timedependent.formulation.time_step]
-total_time = 0.25*s
+total_time = 12.0*s
dt = 0.05*s
# ----------------------------------------------------------------------
Modified: short/3D/PyLith/trunk/playpen/friction/twohex8/opening.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twohex8/opening.cfg 2010-03-25 00:15:03 UTC (rev 16454)
+++ short/3D/PyLith/trunk/playpen/friction/twohex8/opening.cfg 2010-03-25 00:15:43 UTC (rev 16455)
@@ -131,9 +131,10 @@
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_atol = 1.0e-12
ksp_max_it = 100
ksp_gmres_restart = 50
snes_max_it = 10
@@ -154,14 +155,14 @@
# ----------------------------------------------------------------------
# Give basename for VTK domain output of solution over domain.
[pylithapp.problem.formulation.output.output.writer]
-filename = twohex8-opening.vtk
+filename = opening.vtk
# Give basename for VTK fault output.
[pylithapp.timedependent.interfaces.fault.output]
-writer.filename = twohex8-opening-fault.vtk
+writer.filename = opening-fault.vtk
vertex_info_fields = [strike_dir,normal_dir,dip_dir,initial_traction]
# Give basename for VTK output of state variables.
[pylithapp.timedependent.materials.material.output]
cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = twohex8-opening-statevars.vtk
+writer.filename = opening-statevars.vtk
Modified: short/3D/PyLith/trunk/playpen/friction/twohex8/shear-sliding.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/twohex8/shear-sliding.cfg 2010-03-25 00:15:03 UTC (rev 16454)
+++ short/3D/PyLith/trunk/playpen/friction/twohex8/shear-sliding.cfg 2010-03-25 00:15:43 UTC (rev 16455)
@@ -213,9 +213,10 @@
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_atol = 1.0e-15
ksp_max_it = 100
ksp_gmres_restart = 50
snes_max_it = 200
More information about the CIG-COMMITS
mailing list