[cig-commits] r19228 - in short/3D/PyLith/branches/v1.6-stable: examples/3d/hex8 libsrc/pylith/meshio

[brad at geodynamics.org]

Author: brad
Date: 2011-11-21 15:19:41 -0800 (Mon, 21 Nov 2011)
New Revision: 19228

Modified:
   short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8/pylithapp.cfg
   short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/meshio/DataWriterVTK.cc
Log:
Small change to cleanup data writer.

Modified: short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8/pylithapp.cfg
===================================================================
--- short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8/pylithapp.cfg	2011-11-21 23:18:00 UTC (rev 19227)
+++ short/3D/PyLith/branches/v1.6-stable/examples/3d/hex8/pylithapp.cfg	2011-11-21 23:19:41 UTC (rev 19228)
@@ -25,7 +25,7 @@
 # mesh_generator
 # ----------------------------------------------------------------------
 [pylithapp.mesh_generator]
-#debug = 1   ; uncomment to get very verbose mesh information
+debug = 1   ; uncomment to get very verbose mesh information
 
 # Change the default mesh reader to the CUBIT reader.
 reader = pylith.meshio.MeshIOCubit

Modified: short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/meshio/DataWriterVTK.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/meshio/DataWriterVTK.cc	2011-11-21 23:18:00 UTC (rev 19227)
+++ short/3D/PyLith/branches/v1.6-stable/libsrc/pylith/meshio/DataWriterVTK.cc	2011-11-21 23:19:41 UTC (rev 19228)
@@ -266,13 +266,15 @@
     const ALE::Obj<RealSection>& section = field.section();
     assert(!section.isNull());
 
-    const int localFiberDim = 
+    assert(!sieveMesh->getLabelStratum(labelName, depth).isNull());
+    int fiberDimLocal = 
       (sieveMesh->getLabelStratum(labelName, depth)->size() > 0) ? 
       section->getFiberDimension(*sieveMesh->getLabelStratum(labelName, depth)->begin()) : 0;
     int fiberDim = 0;
-    MPI_Allreduce((void *) &localFiberDim, (void *) &fiberDim, 1, 
-		  MPI_INT, MPI_MAX, field.mesh().comm());
+    MPI_Allreduce(&fiberDimLocal, &fiberDim, 1, MPI_INT, MPI_MAX,
+		  field.mesh().comm());
     assert(fiberDim > 0);
+
     const int enforceDim =
       (field.vectorFieldType() != topology::FieldBase::VECTOR) ? fiberDim : 3;