[cig-commits] r14709 - in short/3D/PyLith/branches/pylith-swig: libsrc/meshio unittests/libtests/meshio

[brad at geodynamics.org]

Author: brad
Date: 2009-04-14 17:53:20 -0700 (Tue, 14 Apr 2009)
New Revision: 14709

Modified:
   short/3D/PyLith/branches/pylith-swig/libsrc/meshio/DataWriterVTK.cc
   short/3D/PyLith/branches/pylith-swig/unittests/libtests/meshio/TestDataWriterVTKBCMesh.cc
Log:
Fixed bug in output unit test.

Modified: short/3D/PyLith/branches/pylith-swig/libsrc/meshio/DataWriterVTK.cc
===================================================================
--- short/3D/PyLith/branches/pylith-swig/libsrc/meshio/DataWriterVTK.cc	2009-04-15 00:19:58 UTC (rev 14708)
+++ short/3D/PyLith/branches/pylith-swig/libsrc/meshio/DataWriterVTK.cc	2009-04-15 00:53:20 UTC (rev 14709)
@@ -221,17 +221,13 @@
   typedef typename field_type::Mesh::RealSection RealSection;
 
   try {
-    const ALE::Obj<SieveMesh>& sieveMesh = field.mesh().sieveMesh();
-    assert(!sieveMesh.isNull());
-    const ALE::Obj<typename SieveMesh::label_sequence>& cells = (0 == label) ?
-      sieveMesh->heightStratum(0) :
-      sieveMesh->getLabelStratum(label, labelId);
-    assert(!cells.isNull());
     int rank = 0;
     MPI_Comm_rank(field.mesh().comm(), &rank);
 
     // Correctly handle boundary and fault meshes
     //   Cannot just use mesh->depth() because boundaries report the wrong thing
+    const ALE::Obj<SieveMesh>& sieveMesh = field.mesh().sieveMesh();
+    assert(!sieveMesh.isNull());
     const int cellDepth = (sieveMesh->depth() == -1) ? -1 : 1;
     const int depth = (0 == label) ? cellDepth : labelId;
     const std::string labelName = (0 == label) ?

Modified: short/3D/PyLith/branches/pylith-swig/unittests/libtests/meshio/TestDataWriterVTKBCMesh.cc
===================================================================
--- short/3D/PyLith/branches/pylith-swig/unittests/libtests/meshio/TestDataWriterVTKBCMesh.cc	2009-04-15 00:19:58 UTC (rev 14708)
+++ short/3D/PyLith/branches/pylith-swig/unittests/libtests/meshio/TestDataWriterVTKBCMesh.cc	2009-04-15 00:53:20 UTC (rev 14709)
@@ -302,7 +302,7 @@
       const int fiberDim = _data->cellFieldsInfo[i].fiber_dim;
       fields->add(name, name);
       SubMeshField& field = fields->get(name);
-      field.newSection(topology::FieldBase::CELLS_FIELD, fiberDim);
+      field.newSection(topology::FieldBase::CELLS_FIELD, fiberDim, 1);
       field.allocate();
       field.vectorFieldType(_data->cellFieldsInfo[i].field_type);