[cig-commits] [commit] baagaard/add-examples-debugging, baagaard/add-release-2.0.3, baagaard/add-release-2.1.0, baagaard/dynrup-new-lagrange, baagaard/feature-output-station-names, baagaard/feature-progress-monitor, baagaard/fix-custom-faultpc, baagaard/fix-error-messages, baagaard/fix-faults-intersect, baagaard/fix-friction-initial-state, baagaard/update-autoconf, knepley/feature-petsc-fe, knepley/upgrade-petsc-3.5, knepley/upgrade-petsc-master, maint, master, next, williams/add-examples-grav2d, willic3/fix-plasticity: Small improvements to 2-D gravity example. (c5a9dbe)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Nov 5 15:41:13 PST 2014
- Previous message: [cig-commits] [commit] baagaard/add-examples-debugging, baagaard/add-missing-files-dist, baagaard/add-release-2.0.2, baagaard/add-release-2.0.3, baagaard/add-release-2.1.0, baagaard/dynrup-new-lagrange, baagaard/feature-output-station-names, baagaard/feature-progress-monitor, baagaard/fix-custom-faultpc, baagaard/fix-error-messages, baagaard/fix-faults-intersect, baagaard/fix-friction-initial-state, baagaard/fix-meshing-examples-trelis, baagaard/update-autoconf, knepley/feature-petsc-fe, knepley/fix-dm-composition, knepley/upgrade-petsc-3.5, knepley/upgrade-petsc-master, maint, master, next, willic3/fix-plasticity: Updated meshing README. (2c5c061)
- Next message: [cig-commits] [commit] baagaard/add-examples-debugging, baagaard/add-missing-files-dist, baagaard/add-release-2.0.2, baagaard/add-release-2.0.3, baagaard/add-release-2.1.0, baagaard/dynrup-new-lagrange, baagaard/feature-output-station-names, baagaard/feature-progress-monitor, baagaard/fix-custom-faultpc, baagaard/fix-error-messages, baagaard/fix-faults-intersect, baagaard/fix-friction-initial-state, baagaard/fix-meshing-examples-trelis, baagaard/update-autoconf, knepley/feature-petsc-fe, knepley/fix-dm-composition, knepley/upgrade-petsc-3.5, knepley/upgrade-petsc-master, maint, master, next, willic3/fix-plasticity: Updated geometry journal file based on updated surfaces (change in ids). (256700b)
- Messages sorted by:
[ date ]
[ thread ]
[ subject ]
[ author ]
Repository : https://github.com/geodynamics/pylith
On branches: baagaard/add-examples-debugging,baagaard/add-release-2.0.3,baagaard/add-release-2.1.0,baagaard/dynrup-new-lagrange,baagaard/feature-output-station-names,baagaard/feature-progress-monitor,baagaard/fix-custom-faultpc,baagaard/fix-error-messages,baagaard/fix-faults-intersect,baagaard/fix-friction-initial-state,baagaard/update-autoconf,knepley/feature-petsc-fe,knepley/upgrade-petsc-3.5,knepley/upgrade-petsc-master,maint,master,next,williams/add-examples-grav2d,willic3/fix-plasticity
Link : https://github.com/geodynamics/pylith/compare/f33c75b19fd60eedb2a3405db76a1fee333bb1d7...5b6d812b1612809fea3bd331c4e5af98c25a536a
>---------------------------------------------------------------
commit c5a9dbe5776e87a83efa0040169d00fe3571af03
Author: Brad Aagaard <baagaard at usgs.gov>
Date: Fri Jun 20 12:19:44 2014 -0700
Small improvements to 2-D gravity example.
Fixed location of inner block. Adjusted solver tolerances (KSP smaller
than SNES).
>---------------------------------------------------------------
c5a9dbe5776e87a83efa0040169d00fe3571af03
examples/2d/gravity/README | 11 ++++++-----
examples/2d/gravity/geometry.jou | 7 ++++---
examples/2d/gravity/{mesh_tri3.jou => mesh.jou} | 2 +-
examples/2d/gravity/mesh_tri3.exo | Bin 866892 -> 0 bytes
examples/2d/gravity/pylithapp.cfg | 12 +++++++-----
5 files changed, 18 insertions(+), 14 deletions(-)
diff --git a/examples/2d/gravity/README b/examples/2d/gravity/README
index 333bdf9..5f09daa 100644
--- a/examples/2d/gravity/README
+++ b/examples/2d/gravity/README
@@ -9,7 +9,7 @@ different simulations, which should generally be performed in order:
1. pylith grav_stress.cfg
This will run a single elastic time step to compute the stresses due to
'turning on' gravity. These stresses will be used as initial stresses
- in one of the other simulations. To conver the output PyLith stresses
+ in one of the other simulations. To convert the stresses PyLith outputs
to a spatial database of initial stresses, run the stress2spatialdb.py
script:
./stress2spatialdb.py
@@ -36,7 +36,8 @@ different simulations, which should generally be performed in order:
This simulation uses a finite strain formulation, and should thus
produce less vertical deformation that the infinitesimal strain version.
-The results of all the different simulations may be viewed using ParaView,
-and a ParaView state file is contained in the output directory. By loading
-this state file into ParaView and changing the filename appropriately, the
-various simulation results may be viewed.
+The results of all the different simulations may be viewed using
+ParaView. We provide a ParaView state file (output/deform_view.pvsm),
+which you can load in ParaView to view the solution over the
+domain. The visualization is intended for steps 3 and 4 in which there
+is deformation.
diff --git a/examples/2d/gravity/geometry.jou b/examples/2d/gravity/geometry.jou
index 993ab75..efc9038 100644
--- a/examples/2d/gravity/geometry.jou
+++ b/examples/2d/gravity/geometry.jou
@@ -38,10 +38,10 @@ undo off
# Create outer box (50 km x 40 km x 10 km)..
# ----------------------------------------------------------------------
#{outerBlockX=40.0*km}
-#{outerBlockY=50.0*km}
+#{outerBlockY=30.0*km}
#{outerBlockZ=10.0*km}
brick x {outerBlockX} y {outerBlockY} z {outerBlockZ}
-volume 1 move x 0.0 y {-25.0*km} z 0.0
+volume 1 move x 0.0 y {-0.5*outerBlockY} z 0.0
# ----------------------------------------------------------------------
# Create inner box (2.5 km x 5 km x 10 km)..
@@ -50,7 +50,7 @@ volume 1 move x 0.0 y {-25.0*km} z 0.0
#{innerBlockY=5.0*km}
#{innerBlockZ=10.0*km}
brick x {innerBlockX} y {innerBlockY} z {innerBlockZ}
-volume 2 move x {-40.0*km} y {-2.5*km} z 0.0
+volume 2 move x 0.0 y {-0.5*innerBlockY} z 0.0
# ----------------------------------------------------------------------
# Webcut outer block with inner block and then form a midsurface
@@ -65,3 +65,4 @@ create midsurface volume 4 surface 20 22
delete volume 3 4
imprint all
merge all
+
diff --git a/examples/2d/gravity/mesh_tri3.jou b/examples/2d/gravity/mesh.jou
similarity index 98%
rename from examples/2d/gravity/mesh_tri3.jou
rename to examples/2d/gravity/mesh.jou
index 5451801..8a0b49e 100644
--- a/examples/2d/gravity/mesh_tri3.jou
+++ b/examples/2d/gravity/mesh.jou
@@ -72,7 +72,7 @@ nodeset 14 name "face_ypos"
# ----------------------------------------------------------------------
# Export exodus file
# ----------------------------------------------------------------------
-export mesh "mesh_tri3.exo" dimension 2 overwrite
+export mesh "mesh.exo" dimension 2 overwrite
# End of file
diff --git a/examples/2d/gravity/mesh_tri3.exo b/examples/2d/gravity/mesh_tri3.exo
deleted file mode 100644
index 0899088..0000000
Binary files a/examples/2d/gravity/mesh_tri3.exo and /dev/null differ
diff --git a/examples/2d/gravity/pylithapp.cfg b/examples/2d/gravity/pylithapp.cfg
index 892e636..0d9b2b4 100644
--- a/examples/2d/gravity/pylithapp.cfg
+++ b/examples/2d/gravity/pylithapp.cfg
@@ -28,7 +28,7 @@ materials = 1
reader = pylith.meshio.MeshIOCubit
[pylithapp.mesh_generator.reader]
-filename = mesh_tri3.exo
+filename = mesh.exo
coordsys.space_dim = 2
# ----------------------------------------------------------------------
@@ -174,23 +174,25 @@ sub_pc_factor_shift_type = nonzero
# ksp_type = preonly
# Convergence parameters.
-ksp_rtol = 1.0e-8
+ksp_rtol = 1.0e-10
ksp_atol = 1.0e-12
ksp_max_it = 300
ksp_gmres_restart = 100
# Linear solver monitoring options.
ksp_monitor = true
-ksp_view = true
+#ksp_view = true
ksp_converged_reason = true
+ksp_error_if_not_converged = true
# Nonlinear solver monitoring options.
snes_rtol = 1.0e-10
-snes_atol = 1.0e-12
+snes_atol = 1.0e-11
snes_max_it = 100
snes_monitor = true
-snes_view = true
+#snes_view = true
snes_converged_reason = true
+snes_error_if_not_converged = true
# PETSc summary -- useful for performance information.
log_summary = true
- Previous message: [cig-commits] [commit] baagaard/add-examples-debugging, baagaard/add-missing-files-dist, baagaard/add-release-2.0.2, baagaard/add-release-2.0.3, baagaard/add-release-2.1.0, baagaard/dynrup-new-lagrange, baagaard/feature-output-station-names, baagaard/feature-progress-monitor, baagaard/fix-custom-faultpc, baagaard/fix-error-messages, baagaard/fix-faults-intersect, baagaard/fix-friction-initial-state, baagaard/fix-meshing-examples-trelis, baagaard/update-autoconf, knepley/feature-petsc-fe, knepley/fix-dm-composition, knepley/upgrade-petsc-3.5, knepley/upgrade-petsc-master, maint, master, next, willic3/fix-plasticity: Updated meshing README. (2c5c061)
- Next message: [cig-commits] [commit] baagaard/add-examples-debugging, baagaard/add-missing-files-dist, baagaard/add-release-2.0.2, baagaard/add-release-2.0.3, baagaard/add-release-2.1.0, baagaard/dynrup-new-lagrange, baagaard/feature-output-station-names, baagaard/feature-progress-monitor, baagaard/fix-custom-faultpc, baagaard/fix-error-messages, baagaard/fix-faults-intersect, baagaard/fix-friction-initial-state, baagaard/fix-meshing-examples-trelis, baagaard/update-autoconf, knepley/feature-petsc-fe, knepley/fix-dm-composition, knepley/upgrade-petsc-3.5, knepley/upgrade-petsc-master, maint, master, next, willic3/fix-plasticity: Updated geometry journal file based on updated surfaces (change in ids). (256700b)
- Messages sorted by:
[ date ]
[ thread ]
[ subject ]
[ author ]
More information about the CIG-COMMITS
mailing list