[cig-commits] r21852 - in short/3D/PyLith/trunk: libsrc/pylith/meshio unittests/libtests/meshio

[knepley at geodynamics.org]

Author: knepley
Date: 2013-04-12 18:11:13 -0700 (Fri, 12 Apr 2013)
New Revision: 21852

Modified:
   short/3D/PyLith/trunk/libsrc/pylith/meshio/DataWriterVTK.cc
   short/3D/PyLith/trunk/unittests/libtests/meshio/TestDataWriterFaultMesh.cc
Log:
Fixed definition of fault cell field, fixed setting number of timesteps in writer

Modified: short/3D/PyLith/trunk/libsrc/pylith/meshio/DataWriterVTK.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/pylith/meshio/DataWriterVTK.cc	2013-04-13 00:09:52 UTC (rev 21851)
+++ short/3D/PyLith/trunk/libsrc/pylith/meshio/DataWriterVTK.cc	2013-04-13 01:11:13 UTC (rev 21852)
@@ -140,6 +140,8 @@
   _dm = dmMesh;assert(_dm);
   err = PetscObjectReference((PetscObject) _dm);CHECK_PETSC_ERROR(err);
 
+  DataWriter<mesh_type, field_type>::_numTimeSteps = numTimeSteps;
+
   // Create VTK label in DM: Cleared in close().
   if (label) {
     topology::StratumIS cellsIS(dmMesh, label, labelId);

Modified: short/3D/PyLith/trunk/unittests/libtests/meshio/TestDataWriterFaultMesh.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/meshio/TestDataWriterFaultMesh.cc	2013-04-13 00:09:52 UTC (rev 21851)
+++ short/3D/PyLith/trunk/unittests/libtests/meshio/TestDataWriterFaultMesh.cc	2013-04-13 01:11:13 UTC (rev 21852)
@@ -164,7 +164,7 @@
   const int nfields = _data->numCellFields;
 
   PetscDM dmMesh = _faultMesh->dmMesh();CPPUNIT_ASSERT(dmMesh);
-  topology::Stratum cellsStratum(dmMesh, topology::Stratum::HEIGHT, 0); // MATT This seems wrong.
+  topology::Stratum cellsStratum(dmMesh, topology::Stratum::HEIGHT, 1);
   const PetscInt cStart = cellsStratum.begin();
   const PetscInt cEnd = cellsStratum.end();
   PetscInt numCells = cellsStratum.size();
@@ -175,7 +175,7 @@
     const int fiberDim = _data->cellFieldsInfo[i].fiber_dim;
     fields->add(name, name);
     MeshField& field = fields->get(name);
-    field.newSection(topology::FieldBase::CELLS_FIELD, fiberDim);
+    field.newSection(topology::FieldBase::FACES_FIELD, fiberDim);
     field.allocate();
     field.vectorFieldType(_data->cellFieldsInfo[i].field_type);