[cig-commits] r11633 - short/3D/PyLith/trunk/libsrc/meshio

[knepley at geodynamics.org]

Author: knepley
Date: 2008-03-29 10:45:44 -0700 (Sat, 29 Mar 2008)
New Revision: 11633

Modified:
   short/3D/PyLith/trunk/libsrc/meshio/DataWriterVTK.cc
Log:
I think I fixed the problem with VTK output for fields


Modified: short/3D/PyLith/trunk/libsrc/meshio/DataWriterVTK.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/meshio/DataWriterVTK.cc	2008-03-29 17:05:05 UTC (rev 11632)
+++ short/3D/PyLith/trunk/libsrc/meshio/DataWriterVTK.cc	2008-03-29 17:45:44 UTC (rev 11633)
@@ -159,8 +159,10 @@
       mesh->getFactory()->getNumbering(mesh, labelName, 0);
     assert(!numbering.isNull());
 
-    const int fiberDim = 
+    const int localFiberDim = 
       field->getFiberDimension(*mesh->getLabelStratum(labelName, 0)->begin());
+    int fiberDim;
+    MPI_Allreduce((void *) &localFiberDim, (void *) &fiberDim, 1, MPI_INT, MPI_MAX, mesh->comm());
     const int enforceDim = (fieldType != VECTOR_FIELD) ? fiberDim : 3;
 
     PetscErrorCode err = 0;
@@ -224,8 +226,10 @@
     const ALE::Obj<Mesh::numbering_type>& numbering = 
       mesh->getFactory()->getNumbering(mesh, labelName, depth);
     assert(!numbering.isNull());
-    const int fiberDim = 
+    const int localFiberDim = 
       field->getFiberDimension(*mesh->getLabelStratum(labelName, depth)->begin());
+    int fiberDim;
+    MPI_Allreduce((void *) &localFiberDim, (void *) &fiberDim, 1, MPI_INT, MPI_MAX, mesh->comm());
     const int enforceDim = (fieldType != VECTOR_FIELD) ? fiberDim : 3;
 
     PetscErrorCode err = 0;