[cig-commits] r19892 - short/3D/PyLith/benchmarks/trunk/static/performance/solvertest
brad at geodynamics.org
brad at geodynamics.org
Wed Mar 28 16:12:37 PDT 2012
Author: brad
Date: 2012-03-28 16:12:37 -0700 (Wed, 28 Mar 2012)
New Revision: 19892
Added:
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8.jou
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np001.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np002.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np004.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np008.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np016.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np032.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np064.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np128.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_scaling.jou
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4.jou
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np001.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np002.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np004.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np008.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np016.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np032.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np064.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np128.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_scaling.jou
Removed:
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/mesh_hex8.jou
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/mesh_tet4.jou
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/output/
Modified:
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/geometry.jou
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8.exo
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/setbc.jou
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4.exo
Log:
Started setting up meshes and parameter files for scaling benchmark.
Modified: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/geometry.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/geometry.jou 2012-03-28 23:09:46 UTC (rev 19891)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/geometry.jou 2012-03-28 23:12:37 UTC (rev 19892)
@@ -1,4 +1,23 @@
+# -*- Python -*- (syntax highlighting)
# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+#
+# ----------------------------------------------------------------------
+#
+# CUBIT journal file with geometry for solvertest benchmark.
+#
+# ----------------------------------------------------------------------
+# Set units to SI.
+# ----------------------------------------------------------------------
+#{Units('si')}
+#
+# ----------------------------------------------------------------------
+# Reset geometry.
+# ----------------------------------------------------------------------
+reset
+
+# ----------------------------------------------------------------------
# Create block
# ----------------------------------------------------------------------
@@ -6,35 +25,43 @@
# -36 km <= x <= 36 km
# -36 km <= y <= 36 km
# -36 km <= z <= 0 km
-reset
-brick x 72000 y 72000 z 36000
+#
+#{blockWidth=72.0*km}
+#{blockLength=72.0*km}
+#{blockHeight=36.0*km}
+brick x {blockWidth} y {blockLength} z {blockHeight}
+
# Translate block so the top is at z=0
-volume 1 move x 0 y 0 z -18000
+volume 1 move x 0 y 0 z {-0.5*blockHeight}
# ----------------------------------------------------------------------
# Create interface surfaces
# ----------------------------------------------------------------------
# Center fault
create planar surface with plane xplane offset 0
-surface 7 name "fault_middle"
+#{s=Id("surface")}
+surface {s} name "fault_middle"
rotate surface fault_middle about z angle -15
# Fault on +x side of middle fault
create planar surface with plane yplane offset 0
-surface 8 name "fault_xpos"
+#{s=Id("surface")}
+surface {s} name "fault_xpos"
rotate surface fault_xpos about z angle -45
-move surface fault_xpos x 12000 y 12000 z 0
+move surface fault_xpos x {12*km} y {12*km} z 0
# Fault on -x side of middle fault
create planar surface with plane yplane offset 0
-surface 9 name "fault_xneg"
+#{s=Id("surface")}
+surface {s} name "fault_xneg"
rotate surface fault_xneg about z angle -45
-move surface fault_xneg x -12000 y -12000 z 0
+move surface fault_xneg x {-12*km} y {-12*km} z 0
# Moho
-create planar surface with plane zplane offset -12000
-surface 10 name "material_interface"
+create planar surface with plane zplane offset {-12*km}
+#{s=Id("surface")}
+surface {s} name "material_interface"
# ----------------------------------------------------------------------
# Divide volumes using interface surfaces
@@ -61,3 +88,4 @@
imprint all with volume all
merge all
+# End of file
Modified: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8.cfg 2012-03-28 23:09:46 UTC (rev 19891)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8.cfg 2012-03-28 23:12:37 UTC (rev 19892)
@@ -1,4 +1,3 @@
-# -*- Python -*-
[pylithapp]
# ----------------------------------------------------------------------
@@ -52,10 +51,10 @@
writer.filename = output/hex8-fault-xneg.h5
[pylithapp.timedependent.materials.elastic.output]
-cell_filter = pylith.meshio.CellFilterAvgMesh
writer.filename = output/hex8-elastic.h5
[pylithapp.timedependent.materials.viscoelastic.output]
-cell_filter = pylith.meshio.CellFilterAvgMesh
writer.filename = output/hex8-viscoelastic.h5
+
+# End of file
Modified: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8.exo
===================================================================
(Binary files differ)
Copied: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8.jou (from rev 19889, short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/mesh_hex8.jou)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8.jou 2012-03-28 23:12:37 UTC (rev 19892)
@@ -0,0 +1,46 @@
+# -*- Python -*- (syntax highlighting)
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+#
+# ----------------------------------------------------------------------
+#
+# CUBIT journal file driver for creating hex8 mesh for solvertest
+# benchmark.
+#
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set scheme
+# ----------------------------------------------------------------------
+surface 23 25 41 scheme pave
+surface 33 35 57 scheme pave
+surface 27 29 67 scheme pave
+surface 39 37 77 scheme pave
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+volume all size {3.0*km}
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Generate blocks and nodesets
+# ----------------------------------------------------------------------
+playback 'setbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "hex8.exo" dimension 3 overwrite
+
+# End of file
+
+
Added: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np001.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np001.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np001.cfg 2012-03-28 23:12:37 UTC (rev 19892)
@@ -0,0 +1,42 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = hex8_np001.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/hex8_np001.h5
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/hex8_np001-groundsurf.h5
+
+[pylithapp.timedependent.interfaces.fault_middle.output]
+writer.filename = output/hex8_np001-fault-middle.h5
+
+[pylithapp.timedependent.interfaces.fault_xpos.output]
+writer.filename = output/hex8_np001-fault-xpos.h5
+
+[pylithapp.timedependent.interfaces.fault_xneg.output]
+writer.filename = output/hex8_np001-fault-xneg.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/hex8_np001-elastic.h5
+
+[pylithapp.timedependent.materials.viscoelastic.output]
+writer.filename = output/hex8_np001-viscoelastic.h5
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_monitor_true_residual = False
+log_summary_python = hex8_np001.py
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np002.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np002.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np002.cfg 2012-03-28 23:12:37 UTC (rev 19892)
@@ -0,0 +1,42 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = hex8_np002.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/hex8_np002.h5
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/hex8_np002-groundsurf.h5
+
+[pylithapp.timedependent.interfaces.fault_middle.output]
+writer.filename = output/hex8_np002-fault-middle.h5
+
+[pylithapp.timedependent.interfaces.fault_xpos.output]
+writer.filename = output/hex8_np002-fault-xpos.h5
+
+[pylithapp.timedependent.interfaces.fault_xneg.output]
+writer.filename = output/hex8_np002-fault-xneg.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/hex8_np002-elastic.h5
+
+[pylithapp.timedependent.materials.viscoelastic.output]
+writer.filename = output/hex8_np002-viscoelastic.h5
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_monitor_true_residual = False
+log_summary_python = hex8_np002.py
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np004.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np004.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np004.cfg 2012-03-28 23:12:37 UTC (rev 19892)
@@ -0,0 +1,42 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = hex8_np001.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/hex8_np004.h5
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/hex8_np004-groundsurf.h5
+
+[pylithapp.timedependent.interfaces.fault_middle.output]
+writer.filename = output/hex8_np004-fault-middle.h5
+
+[pylithapp.timedependent.interfaces.fault_xpos.output]
+writer.filename = output/hex8_np004-fault-xpos.h5
+
+[pylithapp.timedependent.interfaces.fault_xneg.output]
+writer.filename = output/hex8_np004-fault-xneg.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/hex8_np004-elastic.h5
+
+[pylithapp.timedependent.materials.viscoelastic.output]
+writer.filename = output/hex8_np004-viscoelastic.h5
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_monitor_true_residual = False
+log_summary_python = hex8_np004.py
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np008.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np008.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np008.cfg 2012-03-28 23:12:37 UTC (rev 19892)
@@ -0,0 +1,42 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = hex8_np008.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/hex8_np008.h5
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/hex8_np008-groundsurf.h5
+
+[pylithapp.timedependent.interfaces.fault_middle.output]
+writer.filename = output/hex8_np008-fault-middle.h5
+
+[pylithapp.timedependent.interfaces.fault_xpos.output]
+writer.filename = output/hex8_np008-fault-xpos.h5
+
+[pylithapp.timedependent.interfaces.fault_xneg.output]
+writer.filename = output/hex8_np008-fault-xneg.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/hex8_np008-elastic.h5
+
+[pylithapp.timedependent.materials.viscoelastic.output]
+writer.filename = output/hex8_np008-viscoelastic.h5
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_monitor_true_residual = False
+log_summary_python = hex8_np008.py
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np016.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np016.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np016.cfg 2012-03-28 23:12:37 UTC (rev 19892)
@@ -0,0 +1,42 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = hex8_np016.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/hex8_np016.h5
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/hex8_np016-groundsurf.h5
+
+[pylithapp.timedependent.interfaces.fault_middle.output]
+writer.filename = output/hex8_np016-fault-middle.h5
+
+[pylithapp.timedependent.interfaces.fault_xpos.output]
+writer.filename = output/hex8_np016-fault-xpos.h5
+
+[pylithapp.timedependent.interfaces.fault_xneg.output]
+writer.filename = output/hex8_np016-fault-xneg.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/hex8_np016-elastic.h5
+
+[pylithapp.timedependent.materials.viscoelastic.output]
+writer.filename = output/hex8_np016-viscoelastic.h5
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_monitor_true_residual = False
+log_summary_python = hex8_np016.py
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np032.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np032.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np032.cfg 2012-03-28 23:12:37 UTC (rev 19892)
@@ -0,0 +1,42 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = hex8_np032.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/hex8_np032.h5
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/hex8_np032-groundsurf.h5
+
+[pylithapp.timedependent.interfaces.fault_middle.output]
+writer.filename = output/hex8_np032-fault-middle.h5
+
+[pylithapp.timedependent.interfaces.fault_xpos.output]
+writer.filename = output/hex8_np032-fault-xpos.h5
+
+[pylithapp.timedependent.interfaces.fault_xneg.output]
+writer.filename = output/hex8_np032-fault-xneg.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/hex8_np032-elastic.h5
+
+[pylithapp.timedependent.materials.viscoelastic.output]
+writer.filename = output/hex8_np032-viscoelastic.h5
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_monitor_true_residual = False
+log_summary_python = hex8_np032.py
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np064.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np064.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np064.cfg 2012-03-28 23:12:37 UTC (rev 19892)
@@ -0,0 +1,42 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = hex8_np064.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/hex8_np064.h5
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/hex8_np064-groundsurf.h5
+
+[pylithapp.timedependent.interfaces.fault_middle.output]
+writer.filename = output/hex8_np064-fault-middle.h5
+
+[pylithapp.timedependent.interfaces.fault_xpos.output]
+writer.filename = output/hex8_np064-fault-xpos.h5
+
+[pylithapp.timedependent.interfaces.fault_xneg.output]
+writer.filename = output/hex8_np064-fault-xneg.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/hex8_np064-elastic.h5
+
+[pylithapp.timedependent.materials.viscoelastic.output]
+writer.filename = output/hex8_np064-viscoelastic.h5
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_monitor_true_residual = False
+log_summary_python = hex8_np064.py
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np128.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np128.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_np128.cfg 2012-03-28 23:12:37 UTC (rev 19892)
@@ -0,0 +1,42 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = hex8_np128.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/hex8_np128.h5
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/hex8_np128-groundsurf.h5
+
+[pylithapp.timedependent.interfaces.fault_middle.output]
+writer.filename = output/hex8_np128-fault-middle.h5
+
+[pylithapp.timedependent.interfaces.fault_xpos.output]
+writer.filename = output/hex8_np128-fault-xpos.h5
+
+[pylithapp.timedependent.interfaces.fault_xneg.output]
+writer.filename = output/hex8_np128-fault-xneg.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/hex8_np128-elastic.h5
+
+[pylithapp.timedependent.materials.viscoelastic.output]
+writer.filename = output/hex8_np128-viscoelastic.h5
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_monitor_true_residual = False
+log_summary_python = hex8_np128.py
+
+
+# End of file
Copied: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_scaling.jou (from rev 19889, short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/mesh_hex8.jou)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_scaling.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8_scaling.jou 2012-03-28 23:12:37 UTC (rev 19892)
@@ -0,0 +1,59 @@
+# -*- Python -*- (syntax highlighting)
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+#
+# ----------------------------------------------------------------------
+#
+# CUBIT journal file driver for creating hex8 meshes for solvertest
+# scaling benchmark.
+#
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set scheme
+# ----------------------------------------------------------------------
+surface 23 25 41 scheme pave
+surface 33 35 57 scheme pave
+surface 27 29 67 scheme pave
+surface 39 37 77 scheme pave
+
+# ----------------------------------------------------------------------
+# Loop to create meshes for scaling benchmark.
+#{dx0=1500.0*m}
+#{iloop=0}
+#{Loop(8)}
+ #{if(iloop==0)}
+ #{dx=1520.0*m}
+ #{else}
+ #{dx=dx0/(2.0**iloop)**(1.0/3.0)}
+ #{endif}
+
+ #{if(iloop==0)} {meshout="hex8_np001.exo"} {endif}
+ #{if(iloop==1)} {meshout="hex8_np002.exo"} {endif}
+ #{if(iloop==2)} {meshout="hex8_np004.exo"} {endif}
+ #{if(iloop==3)} {meshout="hex8_np008.exo"} {endif}
+ #{if(iloop==4)} {meshout="hex8_np016.exo"} {endif}
+ #{if(iloop==5)} {meshout="hex8_np032.exo"} {endif}
+ #{if(iloop==6)} {meshout="hex8_np064.exo"} {endif}
+ #{if(iloop==7)} {meshout="hex8_np128.exo"} {endif}
+
+ delete mesh volume all propagate
+
+ volume all size {dx}
+ mesh volume all
+ playback 'setbc.jou'
+ export mesh {Quote(meshout)} dimension 3 overwrite
+
+ #{tt=iloop+1}
+ #{iloop=tt}
+#{endloop}
+
+
+# End of file
+
+
+
Deleted: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/mesh_hex8.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/mesh_hex8.jou 2012-03-28 23:09:46 UTC (rev 19891)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/mesh_hex8.jou 2012-03-28 23:12:37 UTC (rev 19892)
@@ -1,32 +0,0 @@
-# ----------------------------------------------------------------------
-# Generate geometry
-# ----------------------------------------------------------------------
-playback 'geometry.jou'
-
-# ----------------------------------------------------------------------
-# Set scheme
-# ----------------------------------------------------------------------
-surface 23 25 41 scheme pave
-surface 33 35 57 scheme pave
-surface 27 29 67 scheme pave
-surface 39 37 77 scheme pave
-
-# ----------------------------------------------------------------------
-# Set discretization size
-# ----------------------------------------------------------------------
-volume all size 3000
-
-# ----------------------------------------------------------------------
-# Generate the mesh
-# ----------------------------------------------------------------------
-mesh volume all
-
-# ----------------------------------------------------------------------
-# Generate blocks and nodesets
-# ----------------------------------------------------------------------
-playback 'setbc.jou'
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "hex8.exo" dimension 3 overwrite
Deleted: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/mesh_tet4.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/mesh_tet4.jou 2012-03-28 23:09:46 UTC (rev 19891)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/mesh_tet4.jou 2012-03-28 23:12:37 UTC (rev 19892)
@@ -1,30 +0,0 @@
-# ----------------------------------------------------------------------
-# Generate geometry
-# ----------------------------------------------------------------------
-playback 'geometry.jou'
-
-# ----------------------------------------------------------------------
-# Set scheme
-# ----------------------------------------------------------------------
-volume all scheme tetmesh
-
-# ----------------------------------------------------------------------
-# Set discretization size
-# ----------------------------------------------------------------------
-volume all size 3000
-
-# ----------------------------------------------------------------------
-# Generate the mesh
-# ----------------------------------------------------------------------
-mesh volume all
-
-# ----------------------------------------------------------------------
-# Generate blocks and nodesets
-# ----------------------------------------------------------------------
-playback 'setbc.jou'
-
-# ----------------------------------------------------------------------
-# Export exodus file
-# ----------------------------------------------------------------------
-export mesh "tet4.exo" dimension 3 overwrite
-
Modified: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/pylithapp.cfg 2012-03-28 23:09:46 UTC (rev 19891)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/pylithapp.cfg 2012-03-28 23:12:37 UTC (rev 19892)
@@ -1,4 +1,3 @@
-# -*- Python -*-
[pylithapp]
# ----------------------------------------------------------------------
@@ -23,9 +22,9 @@
# mesh_generator
# ----------------------------------------------------------------------
[pylithapp.mesh_generator]
-#debug = 1
distributor.write_partition = True
-distributor.data_writer.filename = output/partition.vtk
+distributor.data_writer = pylith.meshio.DataWriterHDF5ExtMesh
+distributor.data_writer.filename = output/partition.h5
reader = pylith.meshio.MeshIOCubit
reorder_mesh = True
@@ -163,38 +162,40 @@
[pylithapp.problem.formulation.output.domain]
output_freq = time_step
time_step = 0.999*year
-writer = pylith.meshio.DataWriterHDF5Mesh
+writer = pylith.meshio.DataWriterHDF5ExtMesh
[pylithapp.problem.formulation.output.subdomain]
label = face_zpos
output_freq = time_step
time_step = 0.999*year
-writer = pylith.meshio.DataWriterHDF5SubMesh
+writer = pylith.meshio.DataWriterHDF5ExtSubMesh
[pylithapp.problem.interfaces.fault_middle.output]
output_freq = time_step
time_step = 0.999*year
-writer = pylith.meshio.DataWriterHDF5SubSubMesh
+writer = pylith.meshio.DataWriterHDF5ExtSubSubMesh
[pylithapp.problem.interfaces.fault_xpos.output]
output_freq = time_step
time_step = 0.999*year
-writer = pylith.meshio.DataWriterHDF5SubSubMesh
+writer = pylith.meshio.DataWriterHDF5ExtSubSubMesh
[pylithapp.problem.interfaces.fault_xneg.output]
output_freq = time_step
time_step = 0.999*year
-writer = pylith.meshio.DataWriterHDF5SubSubMesh
+writer = pylith.meshio.DataWriterHDF5ExtSubSubMesh
[pylithapp.problem.materials.elastic.output]
+cell_filter = pylith.meshio.CellFilterAvgMesh
output_freq = time_step
time_step = 0.999*year
-writer = pylith.meshio.DataWriterHDF5Mesh
+writer = pylith.meshio.DataWriterHDF5ExtMesh
[pylithapp.problem.materials.viscoelastic.output]
+cell_filter = pylith.meshio.CellFilterAvgMesh
output_freq = time_step
time_step = 0.999*year
-writer = pylith.meshio.DataWriterHDF5Mesh
+writer = pylith.meshio.DataWriterHDF5ExtMesh
# ----------------------------------------------------------------------
# PETSc
@@ -211,3 +212,5 @@
ksp_monitor_true_residual = true
log_summary = true
+
+# End of file
Modified: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/setbc.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/setbc.jou 2012-03-28 23:09:46 UTC (rev 19891)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/setbc.jou 2012-03-28 23:12:37 UTC (rev 19892)
@@ -1,4 +1,13 @@
+# -*- Python -*- (syntax highlighting)
# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+#
+# ----------------------------------------------------------------------
+#
+# CUBIT journal file with boundary condition info for solvertest benchmark.
+#
+# ----------------------------------------------------------------------
# Create blocks for materials
# ----------------------------------------------------------------------
block 1 volume elastic_xpos_yneg elastic_xpos_ypos elastic_xneg_yneg elastic_xneg_ypos
@@ -9,60 +18,42 @@
# ----------------------------------------------------------------------
# Create nodeset for +x face
# ----------------------------------------------------------------------
-group "face_xpos" add node in surface 45
-group "face_xpos" add node in surface 47
-group "face_xpos" add node in surface 69
-group "face_xpos" add node in surface 76
+group "face_xpos" add node in surface 45 47 69 76
nodeset 10 group face_xpos
nodeset 10 name "face xpos"
# ----------------------------------------------------------------------
# Create nodeset for -x face
# ----------------------------------------------------------------------
-group "face_xneg" add node in surface 59
-group "face_xneg" add node in surface 65
-group "face_xneg" add node in surface 81
-group "face_xneg" add node in surface 83
+group "face_xneg" add node in surface 59 65 81 83
nodeset 11 group face_xneg
nodeset 11 name "face xneg"
# ----------------------------------------------------------------------
# Create nodeset for +y face
# ----------------------------------------------------------------------
-group "face_ypos" add node in surface 42
-group "face_ypos" add node in surface 50
-group "face_ypos" add node in surface 58
-group "face_ypos" add node in surface 66
+group "face_ypos" add node in surface 42 50 58 66
nodeset 12 group face_ypos
nodeset 12 name "face ypos"
# ----------------------------------------------------------------------
# Create nodeset for -y face
# ----------------------------------------------------------------------
-group "face_yneg" add node in surface 68
-group "face_yneg" add node in surface 78
-group "face_yneg" add node in surface 75
-group "face_yneg" add node in surface 86
+group "face_yneg" add node in surface 68 78 75 86
nodeset 13 group face_yneg
nodeset 13 name "face yneg"
# ----------------------------------------------------------------------
# Create nodeset for -z face
# ----------------------------------------------------------------------
-group "face_zneg" add node in surface 23
-group "face_zneg" add node in surface 29
-group "face_zneg" add node in surface 35
-group "face_zneg" add node in surface 37
+group "face_zneg" add node in surface 23 29 35 37
nodeset 14 group face_zneg
nodeset 14 name "face zneg"
# ----------------------------------------------------------------------
# Create nodeset for +z face
# ----------------------------------------------------------------------
-group "face_zpos" add node in surface 25
-group "face_zpos" add node in surface 27
-group "face_zpos" add node in surface 33
-group "face_zpos" add node in surface 39
+group "face_zpos" add node in surface 25 27 33 39
nodeset 15 group face_zpos
nodeset 15 name "face zpos"
@@ -92,3 +83,4 @@
nodeset 102 name "fault_xneg"
+# End of file
\ No newline at end of file
Modified: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4.cfg 2012-03-28 23:09:46 UTC (rev 19891)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4.cfg 2012-03-28 23:12:37 UTC (rev 19892)
@@ -1,4 +1,3 @@
-# -*- Python -*-
[pylithapp]
# ----------------------------------------------------------------------
@@ -13,10 +12,12 @@
[pylithapp.timedependent.materials.elastic]
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = tetrahedron
+quadrature.cell.quad_order = 1
[pylithapp.timedependent.materials.viscoelastic]
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = tetrahedron
+quadrature.cell.quad_order = 1
# ----------------------------------------------------------------------
# faults
@@ -24,14 +25,17 @@
[pylithapp.timedependent.interfaces.fault_middle]
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 1
[pylithapp.timedependent.interfaces.fault_xpos]
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 1
[pylithapp.timedependent.interfaces.fault_xneg]
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 1
# ----------------------------------------------------------------------
# output
@@ -52,10 +56,10 @@
writer.filename = output/tet4-fault-xneg.h5
[pylithapp.timedependent.materials.elastic.output]
-cell_filter = pylith.meshio.CellFilterAvgMesh
writer.filename = output/tet4-elastic.h5
[pylithapp.timedependent.materials.viscoelastic.output]
-cell_filter = pylith.meshio.CellFilterAvgMesh
writer.filename = output/tet4-viscoelastic.h5
+
+# End of file
Modified: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4.exo
===================================================================
(Binary files differ)
Copied: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4.jou (from rev 19889, short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/mesh_tet4.jou)
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4.jou 2012-03-28 23:12:37 UTC (rev 19892)
@@ -0,0 +1,53 @@
+# -*- Python -*- (syntax highlighting)
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+#
+# ----------------------------------------------------------------------
+#
+# CUBIT journal file driver for creating tet4 mesh for solvertest
+# benchmark.
+#
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set scheme
+# ----------------------------------------------------------------------
+volume all scheme tetmesh
+
+# ----------------------------------------------------------------------
+# Set discretization size
+# ----------------------------------------------------------------------
+volume all size {3.0*km}
+
+# ----------------------------------------------------------------------
+# Generate the mesh
+# ----------------------------------------------------------------------
+mesh volume all
+
+# ----------------------------------------------------------------------
+# Smooth mesh to improve quality
+# ----------------------------------------------------------------------
+#{condnum=2.0}
+#{loop(4)}
+cleanup volume all
+volume all smooth scheme condition number beta {condnum} cpu 10
+smooth volume all
+#{condnum=condnum-0.1}
+#{endloop}
+
+# ----------------------------------------------------------------------
+# Generate blocks and nodesets
+# ----------------------------------------------------------------------
+playback 'setbc.jou'
+
+# ----------------------------------------------------------------------
+# Export exodus file
+# ----------------------------------------------------------------------
+export mesh "tet4.exo" dimension 3 overwrite
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np001.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np001.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np001.cfg 2012-03-28 23:12:37 UTC (rev 19892)
@@ -0,0 +1,42 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tet4_np001.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/tet4_np001.h5
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/tet4_np001-groundsurf.h5
+
+[pylithapp.timedependent.interfaces.fault_middle.output]
+writer.filename = output/tet4_np001-fault-middle.h5
+
+[pylithapp.timedependent.interfaces.fault_xpos.output]
+writer.filename = output/tet4_np001-fault-xpos.h5
+
+[pylithapp.timedependent.interfaces.fault_xneg.output]
+writer.filename = output/tet4_np001-fault-xneg.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tet4_np001-elastic.h5
+
+[pylithapp.timedependent.materials.viscoelastic.output]
+writer.filename = output/tet4_np001-viscoelastic.h5
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_monitor_true_residual = False
+log_summary_python = tet4_np001.py
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np002.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np002.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np002.cfg 2012-03-28 23:12:37 UTC (rev 19892)
@@ -0,0 +1,42 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tet4_np002.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/tet4_np002.h5
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/tet4_np002-groundsurf.h5
+
+[pylithapp.timedependent.interfaces.fault_middle.output]
+writer.filename = output/tet4_np002-fault-middle.h5
+
+[pylithapp.timedependent.interfaces.fault_xpos.output]
+writer.filename = output/tet4_np002-fault-xpos.h5
+
+[pylithapp.timedependent.interfaces.fault_xneg.output]
+writer.filename = output/tet4_np002-fault-xneg.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tet4_np002-elastic.h5
+
+[pylithapp.timedependent.materials.viscoelastic.output]
+writer.filename = output/tet4_np002-viscoelastic.h5
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_monitor_true_residual = False
+log_summary_python = tet4_np002.py
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np004.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np004.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np004.cfg 2012-03-28 23:12:37 UTC (rev 19892)
@@ -0,0 +1,42 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tet4_np004.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/tet4_np004.h5
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/tet4_np004-groundsurf.h5
+
+[pylithapp.timedependent.interfaces.fault_middle.output]
+writer.filename = output/tet4_np004-fault-middle.h5
+
+[pylithapp.timedependent.interfaces.fault_xpos.output]
+writer.filename = output/tet4_np004-fault-xpos.h5
+
+[pylithapp.timedependent.interfaces.fault_xneg.output]
+writer.filename = output/tet4_np004-fault-xneg.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tet4_np004-elastic.h5
+
+[pylithapp.timedependent.materials.viscoelastic.output]
+writer.filename = output/tet4_np004-viscoelastic.h5
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_monitor_true_residual = False
+log_summary_python = tet4_np004.py
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np008.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np008.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np008.cfg 2012-03-28 23:12:37 UTC (rev 19892)
@@ -0,0 +1,42 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tet4_np008.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/tet4_np008.h5
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/tet4_np008-groundsurf.h5
+
+[pylithapp.timedependent.interfaces.fault_middle.output]
+writer.filename = output/tet4_np008-fault-middle.h5
+
+[pylithapp.timedependent.interfaces.fault_xpos.output]
+writer.filename = output/tet4_np008-fault-xpos.h5
+
+[pylithapp.timedependent.interfaces.fault_xneg.output]
+writer.filename = output/tet4_np008-fault-xneg.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tet4_np008-elastic.h5
+
+[pylithapp.timedependent.materials.viscoelastic.output]
+writer.filename = output/tet4_np008-viscoelastic.h5
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_monitor_true_residual = False
+log_summary_python = tet4_np008.py
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np016.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np016.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np016.cfg 2012-03-28 23:12:37 UTC (rev 19892)
@@ -0,0 +1,42 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tet4_np016.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/tet4_np016.h5
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/tet4_np016-groundsurf.h5
+
+[pylithapp.timedependent.interfaces.fault_middle.output]
+writer.filename = output/tet4_np016-fault-middle.h5
+
+[pylithapp.timedependent.interfaces.fault_xpos.output]
+writer.filename = output/tet4_np016-fault-xpos.h5
+
+[pylithapp.timedependent.interfaces.fault_xneg.output]
+writer.filename = output/tet4_np016-fault-xneg.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tet4_np016-elastic.h5
+
+[pylithapp.timedependent.materials.viscoelastic.output]
+writer.filename = output/tet4_np016-viscoelastic.h5
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_monitor_true_residual = False
+log_summary_python = tet4_np016.py
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np032.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np032.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np032.cfg 2012-03-28 23:12:37 UTC (rev 19892)
@@ -0,0 +1,42 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tet4_np032.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/tet4_np032.h5
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/tet4_np032-groundsurf.h5
+
+[pylithapp.timedependent.interfaces.fault_middle.output]
+writer.filename = output/tet4_np032-fault-middle.h5
+
+[pylithapp.timedependent.interfaces.fault_xpos.output]
+writer.filename = output/tet4_np032-fault-xpos.h5
+
+[pylithapp.timedependent.interfaces.fault_xneg.output]
+writer.filename = output/tet4_np032-fault-xneg.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tet4_np032-elastic.h5
+
+[pylithapp.timedependent.materials.viscoelastic.output]
+writer.filename = output/tet4_np032-viscoelastic.h5
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_monitor_true_residual = False
+log_summary_python = tet4_np032.py
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np064.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np064.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np064.cfg 2012-03-28 23:12:37 UTC (rev 19892)
@@ -0,0 +1,42 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tet4_np064.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/tet4_np064.h5
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/tet4_np064-groundsurf.h5
+
+[pylithapp.timedependent.interfaces.fault_middle.output]
+writer.filename = output/tet4_np064-fault-middle.h5
+
+[pylithapp.timedependent.interfaces.fault_xpos.output]
+writer.filename = output/tet4_np064-fault-xpos.h5
+
+[pylithapp.timedependent.interfaces.fault_xneg.output]
+writer.filename = output/tet4_np064-fault-xneg.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tet4_np064-elastic.h5
+
+[pylithapp.timedependent.materials.viscoelastic.output]
+writer.filename = output/tet4_np064-viscoelastic.h5
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_monitor_true_residual = False
+log_summary_python = tet4_np064.py
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np128.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np128.cfg (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_np128.cfg 2012-03-28 23:12:37 UTC (rev 19892)
@@ -0,0 +1,42 @@
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator.reader]
+filename = tet4_np128.exo
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/tet4_np128.h5
+
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/tet4_np128-groundsurf.h5
+
+[pylithapp.timedependent.interfaces.fault_middle.output]
+writer.filename = output/tet4_np128-fault-middle.h5
+
+[pylithapp.timedependent.interfaces.fault_xpos.output]
+writer.filename = output/tet4_np128-fault-xpos.h5
+
+[pylithapp.timedependent.interfaces.fault_xneg.output]
+writer.filename = output/tet4_np128-fault-xneg.h5
+
+[pylithapp.timedependent.materials.elastic.output]
+writer.filename = output/tet4_np128-elastic.h5
+
+[pylithapp.timedependent.materials.viscoelastic.output]
+writer.filename = output/tet4_np128-viscoelastic.h5
+
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_monitor_true_residual = False
+log_summary_python = tet4_np128.py
+
+
+# End of file
Added: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_scaling.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_scaling.jou (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4_scaling.jou 2012-03-28 23:12:37 UTC (rev 19892)
@@ -0,0 +1,68 @@
+# -*- Python -*- (syntax highlighting)
+# ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+#
+# ----------------------------------------------------------------------
+#
+# CUBIT journal file driver for creating hex8 meshes for solvertest
+# scaling benchmark.
+#
+# ----------------------------------------------------------------------
+# Generate geometry
+# ----------------------------------------------------------------------
+playback 'geometry.jou'
+
+# ----------------------------------------------------------------------
+# Set scheme
+# ----------------------------------------------------------------------
+volume all scheme tetmesh
+
+# ----------------------------------------------------------------------
+# Loop to create meshes for scaling benchmark.
+#{dx0=1500.0*m}
+#{iloop=0}
+#{Loop(6)}
+ #{if(iloop==0)}
+ #{dx=1520.0*m}
+ #{else}
+ #{dx=dx0/(2.0**iloop)**(1.0/3.0)}
+ #{endif}
+
+ #{if(iloop==0)} {meshout="tet4_np001.exo"} {endif}
+ #{if(iloop==1)} {meshout="tet4_np002.exo"} {endif}
+ #{if(iloop==2)} {meshout="tet4_np004.exo"} {endif}
+ #{if(iloop==3)} {meshout="tet4_np008.exo"} {endif}
+ #{if(iloop==4)} {meshout="tet4_np016.exo"} {endif}
+ #{if(iloop==5)} {meshout="tet4_np032.exo"} {endif}
+ #{if(iloop==6)} {meshout="tet4_np064.exo"} {endif}
+ #{if(iloop==7)} {meshout="tet4_np128.exo"} {endif}
+
+ delete mesh volume all propagate
+
+ volume all size {dx}
+ mesh volume all
+
+# ----------------------------------------------------------------------
+# Smooth mesh to improve quality
+# ----------------------------------------------------------------------
+#{condnum=2.1}
+#{loop(4)}
+cleanup volume all
+volume all smooth scheme condition number beta {condnum} cpu 10
+smooth volume all
+#{condnum=condnum-0.1}
+#{endloop}
+
+ playback 'setbc.jou'
+ export mesh {Quote(meshout)} dimension 3 overwrite
+
+ #{tt=iloop+1}
+ #{iloop=tt}
+#{endloop}
+
+
+# End of file
+
+
+
More information about the CIG-COMMITS
mailing list