[cig-commits] [commit] knepley/fix-parallel-mult-faults: MeshIO: Broadcast the space dimension so that block size gets set correctly (9e687af)

[cig_noreply at geodynamics.org]

Repository : ssh://shell.geodynamics.org/pylith

On branch  : knepley/fix-parallel-mult-faults
Link       : https://github.com/geodynamics/pylith/compare/f04210a9939f366b499dc580318f1e8f0eec0a9c...a68f689441f9cea02053b4daf857d29c93e00fdb

>---------------------------------------------------------------

commit 9e687af4510d5e63fe67578ac8a3731eca95a9ff
Author: Matthew G. Knepley <knepley at gmail.com>
Date:   Mon Mar 31 10:40:02 2014 -0500

    MeshIO: Broadcast the space dimension so that block size gets set correctly


>---------------------------------------------------------------

9e687af4510d5e63fe67578ac8a3731eca95a9ff
 libsrc/pylith/meshio/MeshIOCubit.cc | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/libsrc/pylith/meshio/MeshIOCubit.cc b/libsrc/pylith/meshio/MeshIOCubit.cc
index 72b4846..91fb95c 100644
--- a/libsrc/pylith/meshio/MeshIOCubit.cc
+++ b/libsrc/pylith/meshio/MeshIOCubit.cc
@@ -77,6 +77,7 @@ pylith::meshio::MeshIOCubit::_read(void)
   scalar_array coordinates;
   int_array cells;
   int_array materialIds;
+  PetscErrorCode err = 0;
 
   if (0 == commRank) {
     try {
@@ -85,6 +86,7 @@ pylith::meshio::MeshIOCubit::_read(void)
       const int meshDim = exofile.getDim("num_dim");
 
       _readVertices(exofile, &coordinates, &numVertices, &spaceDim);
+      err = MPI_Bcast(&spaceDim, 1, MPI_INT, 0, _mesh->comm());PYLITH_CHECK_ERROR(err);
       _readCells(exofile, &cells, &materialIds, &numCells, &numCorners);
       _orientCells(&cells, numCells, numCorners, meshDim);
       MeshBuilder::buildMesh(_mesh, &coordinates, numVertices, spaceDim,
@@ -105,6 +107,7 @@ pylith::meshio::MeshIOCubit::_read(void)
       throw std::runtime_error(msg.str());
     } // try/catch
   } else {
+    err = MPI_Bcast(&spaceDim, 1, MPI_INT, 0, _mesh->comm());PYLITH_CHECK_ERROR(err);
     MeshBuilder::buildMesh(_mesh, &coordinates, numVertices, spaceDim,
 			   cells, numCells, numCorners, meshDim,
 			   _interpolate);