[cig-commits] r15071 - in short/3D/PyLith/trunk: . examples/3d/tet4 libsrc/meshio pylith/meshio

[brad at geodynamics.org]

Author: brad
Date: 2009-05-27 15:01:24 -0700 (Wed, 27 May 2009)
New Revision: 15071

Modified:
   short/3D/PyLith/trunk/TODO
   short/3D/PyLith/trunk/examples/3d/tet4/mat_elastic.spatialdb
   short/3D/PyLith/trunk/examples/3d/tet4/mat_viscoelastic.spatialdb
   short/3D/PyLith/trunk/examples/3d/tet4/pylithapp.cfg
   short/3D/PyLith/trunk/libsrc/meshio/DataWriterVTK.cc
   short/3D/PyLith/trunk/pylith/meshio/DataWriterVTK.py
   short/3D/PyLith/trunk/pylith/meshio/DataWriterVTKMesh.py
Log:
Fixed normalization of time constant.

Modified: short/3D/PyLith/trunk/TODO
===================================================================
--- short/3D/PyLith/trunk/TODO	2009-05-27 04:18:46 UTC (rev 15070)
+++ short/3D/PyLith/trunk/TODO	2009-05-27 22:01:24 UTC (rev 15071)
@@ -7,9 +7,8 @@
   field split
 
 Brad
-  Test precompute geometry [DONE]
-  Nondimensionalization
-    Material state vars - use dimstatevars and dimproperties to get scales
+  Questions for Matt:
+    How to set symmetric matrix flag?
   cleanup
   full-scale testing
   test uniform refinement

Modified: short/3D/PyLith/trunk/examples/3d/tet4/mat_elastic.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/tet4/mat_elastic.spatialdb	2009-05-27 04:18:46 UTC (rev 15070)
+++ short/3D/PyLith/trunk/examples/3d/tet4/mat_elastic.spatialdb	2009-05-27 22:01:24 UTC (rev 15071)
@@ -7,7 +7,7 @@
 SimpleDB {
   num-values = 3 // number of material property values
   value-names =  density vs vp // names of the material property values
-  value-units =  kg/m^3  m/s  m/s // units (conversions not supported)
+  value-units =  kg/m**3  m/s  m/s // units (conversions not supported)
   num-locs = 1 // number of locations
   data-dim = 0
   space-dim = 3

Modified: short/3D/PyLith/trunk/examples/3d/tet4/mat_viscoelastic.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/tet4/mat_viscoelastic.spatialdb	2009-05-27 04:18:46 UTC (rev 15070)
+++ short/3D/PyLith/trunk/examples/3d/tet4/mat_viscoelastic.spatialdb	2009-05-27 22:01:24 UTC (rev 15071)
@@ -10,7 +10,7 @@
 SimpleDB {
   num-values = 4 // number of material property values
   value-names =  density vs vp viscosity // names of the material property values
-  value-units =  kg/m^3  m/s  m/s Pa-s // units (conversions not supported)
+  value-units =  kg/m**3  m/s  m/s Pa*s // units (conversions not supported)
   num-locs = 1 // number of locations
   data-dim = 0
   space-dim = 3

Modified: short/3D/PyLith/trunk/examples/3d/tet4/pylithapp.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/tet4/pylithapp.cfg	2009-05-27 04:18:46 UTC (rev 15070)
+++ short/3D/PyLith/trunk/examples/3d/tet4/pylithapp.cfg	2009-05-27 22:01:24 UTC (rev 15071)
@@ -65,14 +65,14 @@
 [pylithapp.timedependent.materials.elastic]
 label = Elastic material
 id = 1
-db.iohandler.filename = mat_elastic.spatialdb
+properties_db.iohandler.filename = mat_elastic.spatialdb
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.shape = tetrahedron
 
 [pylithapp.timedependent.materials.viscoelastic]
 label = Viscoelastic material
 id = 2
-db.iohandler.filename = mat_viscoelastic.spatialdb
+properties_db.iohandler.filename = mat_viscoelastic.spatialdb
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.shape = tetrahedron
 

Modified: short/3D/PyLith/trunk/libsrc/meshio/DataWriterVTK.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/meshio/DataWriterVTK.cc	2009-05-27 04:18:46 UTC (rev 15070)
+++ short/3D/PyLith/trunk/libsrc/meshio/DataWriterVTK.cc	2009-05-27 22:01:24 UTC (rev 15071)
@@ -289,7 +289,7 @@
     sprintf(sbuffer, _timeFormat.c_str(), t/_timeConstant);
     std::string timestamp(sbuffer);
     const int pos = timestamp.find(".");
-    if (pos >0 && pos != timestamp.length())
+    if (pos > 0 && pos != timestamp.length())
       timestamp.erase(pos, 1);
     filename
       << std::string(_filename, 0, indexExt) << "_t" << timestamp << ".vtk";

Modified: short/3D/PyLith/trunk/pylith/meshio/DataWriterVTK.py
===================================================================
--- short/3D/PyLith/trunk/pylith/meshio/DataWriterVTK.py	2009-05-27 04:18:46 UTC (rev 15070)
+++ short/3D/PyLith/trunk/pylith/meshio/DataWriterVTK.py	2009-05-27 22:01:24 UTC (rev 15071)
@@ -66,8 +66,8 @@
 
     # Nondimensionalize
     timeScale = normalizer.timeScale()
-    self.timeConstant = normalizer.nondimensionalize(self.timeConstant,
-                                                     timeScale)
+    self.timeConstantN = normalizer.nondimensionalize(self.timeConstant,
+                                                timeScale)
     return
 
 

Modified: short/3D/PyLith/trunk/pylith/meshio/DataWriterVTKMesh.py
===================================================================
--- short/3D/PyLith/trunk/pylith/meshio/DataWriterVTKMesh.py	2009-05-27 04:18:46 UTC (rev 15070)
+++ short/3D/PyLith/trunk/pylith/meshio/DataWriterVTKMesh.py	2009-05-27 22:01:24 UTC (rev 15071)
@@ -46,7 +46,7 @@
     
     ModuleDataWriterVTK.filename(self, self.filename)
     ModuleDataWriterVTK.timeFormat(self, self.timeFormat)
-    ModuleDataWriterVTK.timeConstant(self, self.timeConstant.value)
+    ModuleDataWriterVTK.timeConstant(self, self.timeConstantN)
     return