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

[brad at geodynamics.org]

Author: brad
Date: 2008-03-29 19:54:05 -0700 (Sat, 29 Mar 2008)
New Revision: 11635

Modified:
   short/3D/PyLith/trunk/libsrc/meshio/DataWriterVTK.cc
Log:
Readded assert to verify fiber dimension for output is positive.

Modified: short/3D/PyLith/trunk/libsrc/meshio/DataWriterVTK.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/meshio/DataWriterVTK.cc	2008-03-30 02:52:51 UTC (rev 11634)
+++ short/3D/PyLith/trunk/libsrc/meshio/DataWriterVTK.cc	2008-03-30 02:54:05 UTC (rev 11635)
@@ -162,7 +162,9 @@
     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());
+    MPI_Allreduce((void *) &localFiberDim, (void *) &fiberDim, 1, MPI_INT, MPI_MAX,
+		  mesh->comm());
+    assert(fiberDim > 0);
     const int enforceDim = (fieldType != VECTOR_FIELD) ? fiberDim : 3;
 
     PetscErrorCode err = 0;
@@ -202,8 +204,8 @@
 				       const ALE::Obj<real_section_type>& field,
 				       const VectorFieldEnum fieldType,
 				       const ALE::Obj<ALE::Mesh>& mesh,
-                       const char* label,
-                       const int labelId)
+				       const char* label,
+				       const int labelId)
 { // writeCellField
   assert(0 != name);
   assert(!mesh.isNull());
@@ -229,7 +231,9 @@
     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());
+    MPI_Allreduce((void *) &localFiberDim, (void *) &fiberDim, 1, MPI_INT, MPI_MAX, 
+		  mesh->comm());
+    assert(fiberDim > 0);
     const int enforceDim = (fieldType != VECTOR_FIELD) ? fiberDim : 3;
 
     PetscErrorCode err = 0;