[cig-commits] [commit] knepley/upgrade-petsc-interface: Modified static plasticity examples to run without errors. (0e1bd7c)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Sun Mar 9 20:16:51 PDT 2014
Repository : ssh://geoshell/pylith
On branch : knepley/upgrade-petsc-interface
Link : https://github.com/geodynamics/pylith/compare/b82388a245b770e4c791bf0acee5ffc180bcba70...0e1bd7c2889e7d4b3a442ab5150ac175e0babed7
>---------------------------------------------------------------
commit 0e1bd7c2889e7d4b3a442ab5150ac175e0babed7
Author: Charles Williams <C.Williams at gns.cri.nz>
Date: Mon Mar 10 16:16:24 2014 +1300
Modified static plasticity examples to run without errors.
It was necessary to reduce the amount of displacement for the axial example and increase the amount of displacement for the shear example.
It was also necessary to allow tensile yield and to change some PETSc settings.
>---------------------------------------------------------------
0e1bd7c2889e7d4b3a442ab5150ac175e0babed7
tests/2d/plasticity/static/axialdisp_fine_2d.cfg | 4 ++--
tests/2d/plasticity/static/axialdisp_fine_3d.cfg | 4 ++--
tests/2d/plasticity/static/axialdisp_twocells_2d.cfg | 4 ++--
tests/2d/plasticity/static/axialdisp_twocells_3d.cfg | 4 ++--
tests/2d/plasticity/static/planestrain_fine_2d.cfg | 1 +
tests/2d/plasticity/static/planestrain_fine_3d.cfg | 1 +
tests/2d/plasticity/static/planestrain_twocells_2d.cfg | 1 +
tests/2d/plasticity/static/planestrain_twocells_3d.cfg | 1 +
tests/2d/plasticity/static/pylithapp.cfg | 12 ++++++------
tests/2d/plasticity/static/sheardisp_fine_2d.cfg | 6 +++---
tests/2d/plasticity/static/sheardisp_fine_3d.cfg | 4 ++--
tests/2d/plasticity/static/sheardisp_twocells_2d.cfg | 4 ++--
tests/2d/plasticity/static/sheardisp_twocells_3d.cfg | 4 ++--
13 files changed, 27 insertions(+), 23 deletions(-)
diff --git a/tests/2d/plasticity/static/axialdisp_fine_2d.cfg b/tests/2d/plasticity/static/axialdisp_fine_2d.cfg
index 89a5b3f..8442912 100644
--- a/tests/2d/plasticity/static/axialdisp_fine_2d.cfg
+++ b/tests/2d/plasticity/static/axialdisp_fine_2d.cfg
@@ -51,7 +51,7 @@ db_initial.label = Dirichlet BC -x edge
# Assign the displacement BC values
db_initial.values = [displacement-x]
-db_initial.data = [1.00*cm]
+db_initial.data = [0.10*cm]
# Boundary conditions to be applied to the positive x-side of the mesh.
[pylithapp.timedependent.bc.x_pos]
@@ -72,7 +72,7 @@ db_initial.label = Dirichlet BC +x edge
# Assign the displacement BC values
db_initial.values = [displacement-x]
-db_initial.data = [-1.00*cm]
+db_initial.data = [-0.10*cm]
# Boundary conditions to be applied to the bottom boundary of the mesh.
[pylithapp.timedependent.bc.y_neg]
diff --git a/tests/2d/plasticity/static/axialdisp_fine_3d.cfg b/tests/2d/plasticity/static/axialdisp_fine_3d.cfg
index 29cacd2..161ecd2 100644
--- a/tests/2d/plasticity/static/axialdisp_fine_3d.cfg
+++ b/tests/2d/plasticity/static/axialdisp_fine_3d.cfg
@@ -54,7 +54,7 @@ db_initial.label = Dirichlet BC -x edge
# Assign the displacement BC values
db_initial.values = [displacement-x]
-db_initial.data = [1.00*cm]
+db_initial.data = [0.10*cm]
# Boundary conditions to be applied to the positive x-side of the mesh.
[pylithapp.timedependent.bc.x_pos]
@@ -75,7 +75,7 @@ db_initial.label = Dirichlet BC +x edge
# Assign the displacement BC values
db_initial.values = [displacement-x]
-db_initial.data = [-1.00*cm]
+db_initial.data = [-0.10*cm]
# Boundary conditions to be applied to the bottom boundary of the mesh.
[pylithapp.timedependent.bc.y_neg]
diff --git a/tests/2d/plasticity/static/axialdisp_twocells_2d.cfg b/tests/2d/plasticity/static/axialdisp_twocells_2d.cfg
index a176b0e..5d69339 100644
--- a/tests/2d/plasticity/static/axialdisp_twocells_2d.cfg
+++ b/tests/2d/plasticity/static/axialdisp_twocells_2d.cfg
@@ -51,7 +51,7 @@ db_initial.label = Dirichlet BC -x edge
# Assign the displacement BC values
db_initial.values = [displacement-x]
-db_initial.data = [1.00*cm]
+db_initial.data = [0.10*cm]
# Boundary conditions to be applied to the positive x-side of the mesh.
[pylithapp.timedependent.bc.x_pos]
@@ -72,7 +72,7 @@ db_initial.label = Dirichlet BC +x edge
# Assign the displacement BC values
db_initial.values = [displacement-x]
-db_initial.data = [-1.00*cm]
+db_initial.data = [-0.10*cm]
# Boundary conditions to be applied to the bottom boundary of the mesh.
[pylithapp.timedependent.bc.y_neg]
diff --git a/tests/2d/plasticity/static/axialdisp_twocells_3d.cfg b/tests/2d/plasticity/static/axialdisp_twocells_3d.cfg
index c8e3653..785cf06 100644
--- a/tests/2d/plasticity/static/axialdisp_twocells_3d.cfg
+++ b/tests/2d/plasticity/static/axialdisp_twocells_3d.cfg
@@ -54,7 +54,7 @@ db_initial.label = Dirichlet BC -x edge
# Assign the displacement BC values
db_initial.values = [displacement-x]
-db_initial.data = [1.00*cm]
+db_initial.data = [0.10*cm]
# Boundary conditions to be applied to the positive x-side of the mesh.
[pylithapp.timedependent.bc.x_pos]
@@ -75,7 +75,7 @@ db_initial.label = Dirichlet BC +x edge
# Assign the displacement BC values
db_initial.values = [displacement-x]
-db_initial.data = [-1.00*cm]
+db_initial.data = [-0.10*cm]
# Boundary conditions to be applied to the bottom boundary of the mesh.
[pylithapp.timedependent.bc.y_neg]
diff --git a/tests/2d/plasticity/static/planestrain_fine_2d.cfg b/tests/2d/plasticity/static/planestrain_fine_2d.cfg
index f731fa3..9913a21 100644
--- a/tests/2d/plasticity/static/planestrain_fine_2d.cfg
+++ b/tests/2d/plasticity/static/planestrain_fine_2d.cfg
@@ -29,5 +29,6 @@ dimension = 2
# Other material settings are defined in pylithapp.cfg.
[pylithapp.timedependent.materials]
material = pylith.materials.DruckerPragerPlaneStrain
+material.allow_tensile_yield = True
# The cell dimension for this material is 2.
material.quadrature.cell.dimension = 2
diff --git a/tests/2d/plasticity/static/planestrain_fine_3d.cfg b/tests/2d/plasticity/static/planestrain_fine_3d.cfg
index c43dae4..7eea926 100644
--- a/tests/2d/plasticity/static/planestrain_fine_3d.cfg
+++ b/tests/2d/plasticity/static/planestrain_fine_3d.cfg
@@ -29,5 +29,6 @@ dimension = 3
# Other material settings are defined in pylithapp.cfg.
[pylithapp.timedependent.materials]
material = pylith.materials.DruckerPrager3D
+material.allow_tensile_yield = True
# The cell dimension for this material is 3.
material.quadrature.cell.dimension = 3
diff --git a/tests/2d/plasticity/static/planestrain_twocells_2d.cfg b/tests/2d/plasticity/static/planestrain_twocells_2d.cfg
index b46755c..cba7087 100644
--- a/tests/2d/plasticity/static/planestrain_twocells_2d.cfg
+++ b/tests/2d/plasticity/static/planestrain_twocells_2d.cfg
@@ -29,5 +29,6 @@ dimension = 2
# Other material settings are defined in pylithapp.cfg.
[pylithapp.timedependent.materials]
material = pylith.materials.DruckerPragerPlaneStrain
+material.allow_tensile_yield = True
# The cell dimension for this material is 2.
material.quadrature.cell.dimension = 2
diff --git a/tests/2d/plasticity/static/planestrain_twocells_3d.cfg b/tests/2d/plasticity/static/planestrain_twocells_3d.cfg
index 558d534..bd0cdca 100644
--- a/tests/2d/plasticity/static/planestrain_twocells_3d.cfg
+++ b/tests/2d/plasticity/static/planestrain_twocells_3d.cfg
@@ -29,5 +29,6 @@ dimension = 3
# Other material settings are defined in pylithapp.cfg.
[pylithapp.timedependent.materials]
material = pylith.materials.DruckerPrager3D
+material.allow_tensile_yield = True
# The cell dimension for this material is 3.
material.quadrature.cell.dimension = 3
diff --git a/tests/2d/plasticity/static/pylithapp.cfg b/tests/2d/plasticity/static/pylithapp.cfg
index 16006c8..22235da 100644
--- a/tests/2d/plasticity/static/pylithapp.cfg
+++ b/tests/2d/plasticity/static/pylithapp.cfg
@@ -83,10 +83,10 @@ pc_type = asm
# Change the preconditioner settings.
sub_pc_factor_shift_type = nonzero
-ksp_rtol = 1.0e-10
-ksp_atol = 1.0e-20
-ksp_max_it = 500
-ksp_gmres_restart = 50
+ksp_rtol = 1.0e-15
+ksp_atol = 1.0e-18
+ksp_max_it = 5000
+ksp_gmres_restart = 500
ksp_monitor = true
ksp_view = true
@@ -94,8 +94,8 @@ ksp_converged_reason = true
ksp_error_if_not_converged = true
snes_rtol = 1.0e-10
-snes_atol = 1.0e-16
-snes_max_it = 100
+snes_atol = 1.0e-12
+snes_max_it = 400
snes_monitor = true
snes_ls_monitor = true
snes_view = true
diff --git a/tests/2d/plasticity/static/sheardisp_fine_2d.cfg b/tests/2d/plasticity/static/sheardisp_fine_2d.cfg
index b018927..f062a86 100644
--- a/tests/2d/plasticity/static/sheardisp_fine_2d.cfg
+++ b/tests/2d/plasticity/static/sheardisp_fine_2d.cfg
@@ -39,7 +39,7 @@ db_initial.label = Dirichlet BC -x edge
# Assign the displacement BC values.
db_initial.values = [displacement-x,displacement-y]
-db_initial.data = [0.0*cm,-0.0001*cm]
+db_initial.data = [0.0*cm,-0.05*cm]
# Boundary conditions to be applied to the positive x-side of the mesh.
@@ -61,7 +61,7 @@ db_initial.label = Dirichlet BC +x edge
# Assign the displacement BC values.
db_initial.values = [displacement-x,displacement-y]
-db_initial.data = [0.0*cm,0.0001*cm]
+db_initial.data = [0.0*cm,0.05*cm]
# ----------------------------------------------------------------------
# PETSc
@@ -69,7 +69,7 @@ db_initial.data = [0.0*cm,0.0001*cm]
[pylithapp.petsc]
ksp_rtol = 1.0e-18
snes_rtol = 1.0e-10
-snes_atol = 1.0e-18
+snes_atol = 1.0e-12
# ----------------------------------------------------------------------
# output
diff --git a/tests/2d/plasticity/static/sheardisp_fine_3d.cfg b/tests/2d/plasticity/static/sheardisp_fine_3d.cfg
index d9865c6..e4dd633 100644
--- a/tests/2d/plasticity/static/sheardisp_fine_3d.cfg
+++ b/tests/2d/plasticity/static/sheardisp_fine_3d.cfg
@@ -40,7 +40,7 @@ db_initial.label = Dirichlet BC -x edge
# Assign the displacement BC values.
db_initial.values = [displacement-x,displacement-y]
-db_initial.data = [0.0*cm,-0.0001*cm]
+db_initial.data = [0.0*cm,-0.05*cm]
# Boundary conditions to be applied to the positive x-side of the mesh.
@@ -62,7 +62,7 @@ db_initial.label = Dirichlet BC +x edge
# Assign the displacement BC values.
db_initial.values = [displacement-x,displacement-y]
-db_initial.data = [0.0*cm,0.0001*cm]
+db_initial.data = [0.0*cm,0.05*cm]
# Boundary conditions to be applied to the front of the mesh.
[pylithapp.timedependent.bc.z_pos]
diff --git a/tests/2d/plasticity/static/sheardisp_twocells_2d.cfg b/tests/2d/plasticity/static/sheardisp_twocells_2d.cfg
index 4b8a885..acd3b9d 100644
--- a/tests/2d/plasticity/static/sheardisp_twocells_2d.cfg
+++ b/tests/2d/plasticity/static/sheardisp_twocells_2d.cfg
@@ -39,7 +39,7 @@ db_initial.label = Dirichlet BC -x edge
# Assign the displacement BC values.
db_initial.values = [displacement-x,displacement-y]
-db_initial.data = [0.0*cm,-0.0001*cm]
+db_initial.data = [0.0*cm,-0.05*cm]
# Boundary conditions to be applied to the positive x-side of the mesh.
@@ -61,7 +61,7 @@ db_initial.label = Dirichlet BC +x edge
# Assign the displacement BC values.
db_initial.values = [displacement-x,displacement-y]
-db_initial.data = [0.0*cm,0.0001*cm]
+db_initial.data = [0.0*cm,0.05*cm]
# ----------------------------------------------------------------------
# PETSc
diff --git a/tests/2d/plasticity/static/sheardisp_twocells_3d.cfg b/tests/2d/plasticity/static/sheardisp_twocells_3d.cfg
index 1fcf06d..9839540 100644
--- a/tests/2d/plasticity/static/sheardisp_twocells_3d.cfg
+++ b/tests/2d/plasticity/static/sheardisp_twocells_3d.cfg
@@ -40,7 +40,7 @@ db_initial.label = Dirichlet BC -x edge
# Assign the displacement BC values.
db_initial.values = [displacement-x,displacement-y]
-db_initial.data = [0.0*cm,-0.0001*cm]
+db_initial.data = [0.0*cm,-0.05*cm]
# Boundary conditions to be applied to the positive x-side of the mesh.
@@ -62,7 +62,7 @@ db_initial.label = Dirichlet BC +x edge
# Assign the displacement BC values.
db_initial.values = [displacement-x,displacement-y]
-db_initial.data = [0.0*cm,0.0001*cm]
+db_initial.data = [0.0*cm,0.05*cm]
# Boundary conditions to be applied to the front of the mesh.
[pylithapp.timedependent.bc.z_pos]
More information about the CIG-COMMITS
mailing list