[cig-commits] r16996 - in short/3D/PyLith/trunk: . applications/utilities
brad at geodynamics.org
brad at geodynamics.org
Fri Jun 11 14:36:53 PDT 2010
Author: brad
Date: 2010-06-11 14:36:53 -0700 (Fri, 11 Jun 2010)
New Revision: 16996
Removed:
short/3D/PyLith/trunk/applications/utilities/pylith_vtktohdf5
Modified:
short/3D/PyLith/trunk/TODO
Log:
Updated TODO.
Modified: short/3D/PyLith/trunk/TODO
===================================================================
--- short/3D/PyLith/trunk/TODO 2010-06-11 20:41:23 UTC (rev 16995)
+++ short/3D/PyLith/trunk/TODO 2010-06-11 21:36:53 UTC (rev 16996)
@@ -2,57 +2,6 @@
CURRENT ISSUES/PRIORITIES
======================================================================
-BINBOTS
-
- need cygwin binbot
-
- need Darwin 10.6 binary
-
-MANUAL
-
- Tutorials
- 3d/hex8 [Charles]
-
- Benchmarks
- Update performance?
-
-TUTORIALS/EXAMPLES
-
- 3d/hex8
- directory structure
- 1. Dirichlet BC (static)
- simple shear w/axial compression (max horiz at 45 degrees)
- 2. Neumann BC (static)
- remove Dirichlet BC on two sides and replace w/tractions
- 3. Earthquake rupture (static)
- add fault slip w/greater slip at depth
- 4. Dirichlet BC (quasi-static)
- constant velocity BC + MaxwellIsotropic3D (lower crust)
- 5. Dirichlet + Neumann (quasi-static)
- constant velocity + time-dependent traction BC
- 6. Multiple earthquake rupture + creep (quasi-static)
- fixed boundaries, w/multiple eq ruptures
- 7. Earthquake rupture + creep + Dirichlet BC (quasi-static)
- mini savage-prescott
- 8. Add power-law rheology (quasi-static)
- nonlinear solver
- 9. Replace power-law with Drucker-Prager elastoplastic
- nonlinear solver
- 10. Static friction + Dirichlet BC (static)
- stick, simple shear w/axial compression
- 11. Static friction + Dirichlet BC (static)
- slip, simple shear w/axial compression
- 12. Static friction + Dirichlet BC (quasi-static)
- slip, simple shear w/axial compression
- 13. Slip-weakening friction + Dirichlet BC (quasi-static)
- slip, simple shear w/axial compression
- 14. Rate- and state-friction + Dirichlet BC (quasi-static)
- slip, simple shear w/axial compression
- 15. Gravity
- 16. Gravity w/prestress
- 17. Gravity w/small strain
-
-
----------------------------------------------------------------------
SECONDARY PRIORITIES
----------------------------------------------------------------------
@@ -61,7 +10,7 @@
Need field split working for both SolverLinear and SolverNonlinear.
-* Output to VTK XML files w/data values as binary
+* Output to HDF5 files.
* Uniform global refinement for tets with faults
Deleted: short/3D/PyLith/trunk/applications/utilities/pylith_vtktohdf5
===================================================================
--- short/3D/PyLith/trunk/applications/utilities/pylith_vtktohdf5 2010-06-11 20:41:23 UTC (rev 16995)
+++ short/3D/PyLith/trunk/applications/utilities/pylith_vtktohdf5 2010-06-11 21:36:53 UTC (rev 16996)
@@ -1,142 +0,0 @@
-#!/usr/bin/env nemesis
-#
-# ======================================================================
-#
-# Brad T. Aagaard
-# U.S. Geological Survey
-#
-# {LicenseText}
-#
-# ======================================================================
-#
-
-# This script combines VTK files written by PyLith into a single HDF5
-# file. The filenames must be of the pattern [name]_t[number].vtk and
-# follow the usual VTK conventions for files with time slices. Usage:
-# pylith_vtktohdf5 [--verbose] [--filename=FILE]
-
-__requires__ = "PyLith"
-
-
-# ======================================================================
-class ConvertApp(object):
- """
- Application for combining VTK files into a single HDF5 file.
- """
-
- # PUBLIC METHODS /////////////////////////////////////////////////////
-
- def __init__(self, name="pylith_vtktohdf5"):
- """
- Constructor.
- """
- self.verbose = False
- self.filename = ""
- return
-
-
- def main(self, *args, **kwds):
- """
- Main entry point for application.
- """
-
- reader = VTKDataReader()
- filenames = self._getFilenames()
- self.h5 = ()
-
- nsteps = len(filenames)
- assert(nsteps > 0)
-
- # Initial time step
- if self.verbose:
- print "Reading data from VTK file '%s'." % filenames[0]
- data = reader.read(filenames[0])
- self._writeVertices(data['vertices'])
- self._writeCells(data['cells'])
- self._setupGroups(nsteps, data['vertex_fields'], data['cell_fields'])
- self._writeVertexFields(data['vertex_fields'])
- self._writeVellFields(data['cell_fields'])
-
- # Time steps
- for istep in xrange(nsteps):
- if self.verbose:
- print "Reading data from VTK file '%s'." % filenames[0]
- data = reader.read(filenames[istep])
- self._writeVertexFields(data['vertex_fields'])
- self._writeCellFields(data['cell_fields'])
-
- h5.close()
-
- return
-
-
- # PRIVATE METHODS ////////////////////////////////////////////////////
-
- def _getFilenames(self):
- """
- Get filenames with same pattern as provided filename but different time stamps.
- """
- filenames = []
- return filenames
-
-
- def _setupGroups(self, data):
- """
- Setup groups in HDF5 file.
- """
- return
-
-
- def _writeVertices(self, data):
- """
- Write vertices to HDF5 file.
- """
- return
-
-
- def _writeCells(self, data):
- """
- Write cells to HDF5 file.
- """
- return
-
-
- def _writeVertexFields(self, data):
- """
- Add vertex fields at current time slice to HDF5 file.
- """
- return
-
-
- def _writeCellFields(self, data):
- """
- Add cell fields at current time slice to HDF5 file.
- """
- return
-
-
-# ----------------------------------------------------------------------
-if __name__ == "__main__":
-
- usage = "%prog [--verbose] [--file=FILE]"
- from optparse import OptionParser
- parser = OptionParser(usage=usage)
- parser.add_option("-v", "--verbose", dest="verbose",
- action="store_true", default=False,
- help="Print verbose output.")
- parser.add_option("-o", "--file", dest="filename",
- type="string", metavar="FILE",
- help="Read data file files matching FILE.",
- default="")
- (options, args) = parser.parse_args()
- if options.filename == "":
- raise IOError("No input provided. Specify a file using the --file command "
- "line argument.")
-
- app = ConvertApp()
- app.verbose = options.verbose
- app.filename = options.filename
- app.main()
-
-
-# End of file
More information about the CIG-COMMITS
mailing list