[cig-commits] r19019 - short/3D/PyLith/benchmarks/trunk/static/performance/solvertest
brad at geodynamics.org
brad at geodynamics.org
Wed Oct 5 16:19:17 PDT 2011
Author: brad
Date: 2011-10-05 16:19:16 -0700 (Wed, 05 Oct 2011)
New Revision: 19019
Modified:
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/asm.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/pylithapp.cfg
short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4.cfg
Log:
Switched to HDF5 output.
Modified: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/asm.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/asm.cfg 2011-10-05 23:17:27 UTC (rev 19018)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/asm.cfg 2011-10-05 23:19:16 UTC (rev 19019)
@@ -1,6 +1,6 @@
[pylithapp.petsc]
ksp_type = gmres
-ksp_gmres_restart = 500
+ksp_gmres_restart = 100
pc_type = asm
sub_pc_factor_shift_type = nonzero
Modified: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8.cfg 2011-10-05 23:17:27 UTC (rev 19018)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/hex8.cfg 2011-10-05 23:19:16 UTC (rev 19019)
@@ -24,38 +24,44 @@
[pylithapp.timedependent.interfaces.fault_middle]
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 2
+#quadrature.cell.quad_order = 1
+#quadrature.cell.collocate_quad = True
[pylithapp.timedependent.interfaces.fault_xpos]
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 2
+#quadrature.cell.quad_order = 1
+#quadrature.cell.collocate_quad = True
[pylithapp.timedependent.interfaces.fault_xneg]
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 2
+#quadrature.cell.quad_order = 1
+#quadrature.cell.collocate_quad = True
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
[pylithapp.problem.formulation.output.domain]
-writer.filename = output/hex8.vtk
+writer.filename = output/hex8.h5
[pylithapp.problem.formulation.output.subdomain]
-writer.filename = output/hex8-groundsurf.vtk
+writer.filename = output/hex8-groundsurf.h5
[pylithapp.timedependent.interfaces.fault_middle.output]
-writer.filename = output/hex8-fault-middle.vtk
+writer.filename = output/hex8-fault-middle.h5
[pylithapp.timedependent.interfaces.fault_xpos.output]
-writer.filename = output/hex8-fault-xpos.vtk
+writer.filename = output/hex8-fault-xpos.h5
[pylithapp.timedependent.interfaces.fault_xneg.output]
-writer.filename = output/hex8-fault-xneg.vtk
+writer.filename = output/hex8-fault-xneg.h5
[pylithapp.timedependent.materials.elastic.output]
cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = output/hex8-elastic.vtk
+writer.filename = output/hex8-elastic.h5
[pylithapp.timedependent.materials.viscoelastic.output]
cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = output/hex8-viscoelastic.vtk
+writer.filename = output/hex8-viscoelastic.h5
Modified: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/pylithapp.cfg 2011-10-05 23:17:27 UTC (rev 19018)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/pylithapp.cfg 2011-10-05 23:19:16 UTC (rev 19019)
@@ -43,7 +43,7 @@
bc = [x_neg,x_pos,y_neg,y_pos,z_neg]
interfaces = [fault_middle,fault_xpos,fault_xneg]
-normalizer.length_scale = 1.0*m
+normalizer.length_scale = 3.0*km
[pylithapp.timedependent.formulation]
output = [domain,subdomain]
@@ -117,10 +117,9 @@
[pylithapp.timedependent.interfaces.fault_middle.eq_srcs.rupture.slip_function]
slip.label = Final slip on middle fault
-slip = spatialdata.spatialdb.UniformDB
-slip.values = [left-lateral-slip,reverse-slip,fault-opening]
-slip.data = [-1.0*m,0.0*m,0.0*m]
+slip = spatialdata.spatialdb.SimpleDB
slip.label = Final slip for fault middle
+slip.iohandler.filename = finalslip_middle.spatialdb
slip_time = spatialdata.spatialdb.UniformDB
slip_time.values = [slip-time]
@@ -133,10 +132,9 @@
[pylithapp.timedependent.interfaces.fault_xpos.eq_srcs.rupture.slip_function]
slip.label = Final slip on +x fault
-slip = spatialdata.spatialdb.UniformDB
-slip.values = [left-lateral-slip,reverse-slip,fault-opening]
-slip.data = [-0.5*m,0.0*m,0.0*m]
+slip = spatialdata.spatialdb.SimpleDB
slip.label = Final slip for fault xpos
+slip.iohandler.filename = finalslip_xpos.spatialdb
slip_time = spatialdata.spatialdb.UniformDB
slip_time.values = [slip-time]
@@ -150,10 +148,9 @@
[pylithapp.timedependent.interfaces.fault_xneg.eq_srcs.rupture.slip_function]
slip.label = Final slip on -x fault
-slip = spatialdata.spatialdb.UniformDB
-slip.values = [left-lateral-slip,reverse-slip,fault-opening]
-slip.data = [-0.5*m,0.0*m,0.0*m]
+slip = spatialdata.spatialdb.SimpleDB
slip.label = Final slip for fault xneg
+slip.iohandler.filename = finalslip_xneg.spatialdb
slip_time = spatialdata.spatialdb.UniformDB
slip_time.values = [slip-time]
@@ -166,45 +163,38 @@
[pylithapp.problem.formulation.output.domain]
output_freq = time_step
time_step = 0.999*year
-writer.time_format = %04.0f
-writer.time_constant = 0.1*year
+writer = pylith.meshio.DataWriterHDF5Mesh
[pylithapp.problem.formulation.output.subdomain]
label = face_zpos
output_freq = time_step
time_step = 0.999*year
-writer.time_format = %04.0f
-writer.time_constant = 0.1*year
+writer = pylith.meshio.DataWriterHDF5SubMesh
[pylithapp.problem.interfaces.fault_middle.output]
output_freq = time_step
time_step = 0.999*year
-writer.time_format = %04.0f
-writer.time_constant = 0.1*year
+writer = pylith.meshio.DataWriterHDF5SubSubMesh
[pylithapp.problem.interfaces.fault_xpos.output]
output_freq = time_step
time_step = 0.999*year
-writer.time_format = %04.0f
-writer.time_constant = 0.1*year
+writer = pylith.meshio.DataWriterHDF5SubSubMesh
[pylithapp.problem.interfaces.fault_xneg.output]
output_freq = time_step
time_step = 0.999*year
-writer.time_format = %04.0f
-writer.time_constant = 0.1*year
+writer = pylith.meshio.DataWriterHDF5SubSubMesh
[pylithapp.problem.materials.elastic.output]
output_freq = time_step
time_step = 0.999*year
-writer.time_format = %04.0f
-writer.time_constant = 0.1*year
+writer = pylith.meshio.DataWriterHDF5Mesh
[pylithapp.problem.materials.viscoelastic.output]
output_freq = time_step
time_step = 0.999*year
-writer.time_format = %04.0f
-writer.time_constant = 0.1*year
+writer = pylith.meshio.DataWriterHDF5Mesh
# ----------------------------------------------------------------------
# PETSc
@@ -212,7 +202,7 @@
[pylithapp.petsc]
ksp_rtol = 1.0e-8
-ksp_atol = 1.0e-12
+ksp_atol = 1.0e-20
ksp_max_it = 500
ksp_monitor = true
Modified: short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4.cfg 2011-10-05 23:17:27 UTC (rev 19018)
+++ short/3D/PyLith/benchmarks/trunk/static/performance/solvertest/tet4.cfg 2011-10-05 23:19:16 UTC (rev 19019)
@@ -40,25 +40,25 @@
# output
# ----------------------------------------------------------------------
[pylithapp.problem.formulation.output.domain]
-writer.filename = output/tet4.vtk
+writer.filename = output/tet4.h5
[pylithapp.problem.formulation.output.subdomain]
-writer.filename = output/tet4-groundsurf.vtk
+writer.filename = output/tet4-groundsurf.h5
[pylithapp.timedependent.interfaces.fault_middle.output]
-writer.filename = output/tet4-fault-middle.vtk
+writer.filename = output/tet4-fault-middle.h5
[pylithapp.timedependent.interfaces.fault_xpos.output]
-writer.filename = output/tet4-fault-xpos.vtk
+writer.filename = output/tet4-fault-xpos.h5
[pylithapp.timedependent.interfaces.fault_xneg.output]
-writer.filename = output/tet4-fault-xneg.vtk
+writer.filename = output/tet4-fault-xneg.h5
[pylithapp.timedependent.materials.elastic.output]
cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = output/tet4-elastic.vtk
+writer.filename = output/tet4-elastic.h5
[pylithapp.timedependent.materials.viscoelastic.output]
cell_filter = pylith.meshio.CellFilterAvgMesh
-writer.filename = output/tet4-viscoelastic.vtk
+writer.filename = output/tet4-viscoelastic.h5
More information about the CIG-COMMITS
mailing list