]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Actually fix compute_local_to_global_vertex_index_map. 6507/head
authorDavid Wells <wellsd2@rpi.edu>
Sun, 6 May 2018 19:57:18 +0000 (15:57 -0400)
committerMatthias Maier <tamiko@43-1.org>
Sun, 6 May 2018 23:32:43 +0000 (18:32 -0500)
The MPI types do not match the deal.II types, which is causes chaos when we use
32 bit indices.

source/grid/grid_tools.cc

index a38f743fda84c8b8b35b0da67732d71183bce063..fab5e1a1b5b412bc6aaead2a56fab816f72cfedb 100644 (file)
@@ -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<types::global_vertex_index> 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());

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.