[cig-commits] [commit] knepley/upgrade-petsc-interface: Updated README for initial stress case. (a5ef048)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/pylith

On branch  : knepley/upgrade-petsc-interface
Link       : https://github.com/geodynamics/pylith/compare/c55f238bb80c197e0216732910b531602b2eb7ac...1e3ea1f187463f5f843746ecb11e350472d667ca

>---------------------------------------------------------------

commit a5ef0483ed153864f53fe041dd3181308a3ffb67
Author: Charles Williams <C.Williams at gns.cri.nz>
Date:   Thu Mar 13 14:17:19 2014 +1300

    Updated README for initial stress case.
    
    Directed output to output directory and added script to run tests.


>---------------------------------------------------------------

a5ef0483ed153864f53fe041dd3181308a3ffb67
 tests/2d/powerlaw/README                       | 8 +++++++-
 tests/2d/powerlaw/axialdisp2d.cfg              | 4 ++--
 tests/2d/powerlaw/axialdisp3d.cfg              | 4 ++--
 tests/2d/powerlaw/axialdisp_init_stress_2d.cfg | 4 ++--
 tests/2d/powerlaw/axialdisp_init_stress_3d.cfg | 4 ++--
 tests/2d/powerlaw/sheardisp2d.cfg              | 4 ++--
 tests/2d/powerlaw/sheardisp3d.cfg              | 4 ++--
 7 files changed, 19 insertions(+), 13 deletions(-)

diff --git a/tests/2d/powerlaw/README b/tests/2d/powerlaw/README
index ae66da9..a287f17 100644
--- a/tests/2d/powerlaw/README
+++ b/tests/2d/powerlaw/README
@@ -1,6 +1,6 @@
 This is a simple test of the PowerLawPlaneStrain viscoelastic model. The
 results may be compared against the results from a 3D model with the
-z-displacements constrained to zero.  There are four models:
+z-displacements constrained to zero.  There are six models:
 
 Plane strain axial displacement problem:
 pylith planestrain2d.cfg axialdisp2d.cfg
@@ -8,11 +8,17 @@ pylith planestrain2d.cfg axialdisp2d.cfg
 Plane strain shear displacement problem:
 pylith planestrain2d.cfg sheardisp2d.cfg
 
+Plane strain axial displacement problem with initial stresses:
+pylith planestrain2d.cfg axialdisp_init_stress_2d.cfg
+
 3D simulated plane strain axial displacement problem:
 pylith planestrain3d.cfg axialdisp3d.cfg
 
 3D simulated plane strain shear displacement problem:
 pylith planestrain3d.cfg sheardisp3d.cfg
 
+3D simulated plane strain axial displacement problem with initial stresses:
+pylith planestrain3d.cfg axialdisp_init_stress_3d.cfg
+
 The 2D and 3D problems should yield the same displacements, strains,
 stresses, and viscous strains when the appropriate components are compared.
diff --git a/tests/2d/powerlaw/axialdisp2d.cfg b/tests/2d/powerlaw/axialdisp2d.cfg
index a067551..79446c1 100644
--- a/tests/2d/powerlaw/axialdisp2d.cfg
+++ b/tests/2d/powerlaw/axialdisp2d.cfg
@@ -91,13 +91,13 @@ db_initial.label = Dirichlet BC -y edge
 # ----------------------------------------------------------------------
 # Give basename for VTK output of solution over domain.
 [pylithapp.problem.formulation.output.output.writer]
-filename = axialdisp2d.vtk
+filename = output/axialdisp2d.vtk
 time_format = %04.0f
 time_constant = 1.0*year
 
 # Give basename for VTK output of state variables.
 [pylithapp.timedependent.materials.material.output]
 cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = axialdisp2d-statevars.vtk
+writer.filename = output/axialdisp2d-statevars.vtk
 writer.time_format = %04.0f
 writer.time_constant = 1.0*year
diff --git a/tests/2d/powerlaw/axialdisp3d.cfg b/tests/2d/powerlaw/axialdisp3d.cfg
index 611353a..b30928b 100644
--- a/tests/2d/powerlaw/axialdisp3d.cfg
+++ b/tests/2d/powerlaw/axialdisp3d.cfg
@@ -120,13 +120,13 @@ db_initial.label = Dirichlet BC -z edge
 # ----------------------------------------------------------------------
 # Give basename for VTK output of solution over domain.
 [pylithapp.problem.formulation.output.output.writer]
-filename = axialdisp3d.vtk
+filename = output/axialdisp3d.vtk
 time_format = %04.0f
 time_constant = 1.0*year
 
 # Give basename for VTK output of state variables.
 [pylithapp.timedependent.materials.material.output]
 cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = axialdisp3d-statevars.vtk
+writer.filename = output/axialdisp3d-statevars.vtk
 writer.time_format = %04.0f
 writer.time_constant = 1.0*year
diff --git a/tests/2d/powerlaw/axialdisp_init_stress_2d.cfg b/tests/2d/powerlaw/axialdisp_init_stress_2d.cfg
index e5683f8..0f24001 100644
--- a/tests/2d/powerlaw/axialdisp_init_stress_2d.cfg
+++ b/tests/2d/powerlaw/axialdisp_init_stress_2d.cfg
@@ -104,13 +104,13 @@ db_initial_state.iohandler.filename = init_state_2d.spatialdb
 # ----------------------------------------------------------------------
 # Give basename for VTK output of solution over domain.
 [pylithapp.problem.formulation.output.output.writer]
-filename = axialdisp_init_stress_2d.vtk
+filename = output/axialdisp_init_stress_2d.vtk
 time_format = %04.0f
 time_constant = 1.0*year
 
 # Give basename for VTK output of state variables.
 [pylithapp.timedependent.materials.material.output]
 cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = axialdisp_init_stress_2d-statevars.vtk
+writer.filename = output/axialdisp_init_stress_2d-statevars.vtk
 writer.time_format = %04.0f
 writer.time_constant = 1.0*year
diff --git a/tests/2d/powerlaw/axialdisp_init_stress_3d.cfg b/tests/2d/powerlaw/axialdisp_init_stress_3d.cfg
index 9a852e5..d2c180a 100644
--- a/tests/2d/powerlaw/axialdisp_init_stress_3d.cfg
+++ b/tests/2d/powerlaw/axialdisp_init_stress_3d.cfg
@@ -130,13 +130,13 @@ db_initial_stress.iohandler.filename = init_stress_3d.spatialdb
 # ----------------------------------------------------------------------
 # Give basename for VTK output of solution over domain.
 [pylithapp.problem.formulation.output.output.writer]
-filename = axialdisp_init_stress_3d.vtk
+filename = output/axialdisp_init_stress_3d.vtk
 time_format = %04.0f
 time_constant = 1.0*year
 
 # Give basename for VTK output of state variables.
 [pylithapp.timedependent.materials.material.output]
 cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = axialdisp_init_stress_3d-statevars.vtk
+writer.filename = output/axialdisp_init_stress_3d-statevars.vtk
 writer.time_format = %04.0f
 writer.time_constant = 1.0*year
diff --git a/tests/2d/powerlaw/sheardisp2d.cfg b/tests/2d/powerlaw/sheardisp2d.cfg
index 73213cf..eb5cb39 100644
--- a/tests/2d/powerlaw/sheardisp2d.cfg
+++ b/tests/2d/powerlaw/sheardisp2d.cfg
@@ -67,13 +67,13 @@ db_initial.data = [0.0*cm,1.0*cm]
 # ----------------------------------------------------------------------
 # Give basename for VTK output of solution over domain.
 [pylithapp.problem.formulation.output.output.writer]
-filename = sheardisp2d.vtk
+filename = output/sheardisp2d.vtk
 time_format = %04.0f
 time_constant = 1.0*year
 
 # Give basename for VTK output of state variables.
 [pylithapp.timedependent.materials.material.output]
 cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = sheardisp2d-statevars.vtk
+writer.filename = output/sheardisp2d-statevars.vtk
 writer.time_format = %04.0f
 writer.time_constant = 1.0*year
diff --git a/tests/2d/powerlaw/sheardisp3d.cfg b/tests/2d/powerlaw/sheardisp3d.cfg
index cfc6f06..ef81b57 100644
--- a/tests/2d/powerlaw/sheardisp3d.cfg
+++ b/tests/2d/powerlaw/sheardisp3d.cfg
@@ -94,13 +94,13 @@ db_initial.label = Dirichlet BC -z edge
 # ----------------------------------------------------------------------
 # Give basename for VTK output of solution over domain.
 [pylithapp.problem.formulation.output.output.writer]
-filename = sheardisp3d.vtk
+filename = output/sheardisp3d.vtk
 time_format = %04.0f
 time_constant = 1.0*year
 
 # Give basename for VTK output of state variables.
 [pylithapp.timedependent.materials.material.output]
 cell_filter = pylith.meshio.CellFilterAvg
-writer.filename = sheardisp3d-statevars.vtk
+writer.filename = output/sheardisp3d-statevars.vtk
 writer.time_format = %04.0f
 writer.time_constant = 1.0*year