From: David Wells Date: Sun, 6 May 2018 19:57:18 +0000 (-0400) Subject: Actually fix compute_local_to_global_vertex_index_map. X-Git-Tag: v9.0.0-rc3~4^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F6507%2Fhead;p=dealii.git Actually fix compute_local_to_global_vertex_index_map. The MPI types do not match the deal.II types, which is causes chaos when we use 32 bit indices. --- diff --git a/source/grid/grid_tools.cc b/source/grid/grid_tools.cc index a38f743fda..fab5e1a1b5 100644 --- a/source/grid/grid_tools.cc +++ b/source/grid/grid_tools.cc @@ -2139,8 +2139,8 @@ next_cell: // processors and shifting the indices accordingly const unsigned int n_cpu = Utilities::MPI::n_mpi_processes(triangulation.get_communicator()); std::vector indices(n_cpu); - int ierr = MPI_Allgather(&next_index, 1, DEAL_II_DOF_INDEX_MPI_TYPE, indices.data(), - 1, DEAL_II_DOF_INDEX_MPI_TYPE, triangulation.get_communicator()); + int ierr = MPI_Allgather(&next_index, 1, DEAL_II_VERTEX_INDEX_MPI_TYPE, indices.data(), + 1, DEAL_II_VERTEX_INDEX_MPI_TYPE, triangulation.get_communicator()); AssertThrowMPI(ierr); Assert(indices.begin() + triangulation.locally_owned_subdomain() < indices.end(), ExcInternalError());