src/smpi/bindings/smpi_pmpi_coll.cpp

   1 /* Copyright (c) 2007-2019. The SimGrid Team. All rights reserved.          */
   2
   3 /* This program is free software; you can redistribute it and/or modify it
   4  * under the terms of the license (GNU LGPL) which comes with this package. */
   5
   6 #include "private.hpp"
   7 #include "smpi_coll.hpp"
   8 #include "smpi_comm.hpp"
   9 #include "smpi_request.hpp"
  10 #include "smpi_datatype_derived.hpp"
  11 #include "smpi_op.hpp"
  12 #include "src/smpi/include/smpi_actor.hpp"
  13
  14 XBT_LOG_EXTERNAL_DEFAULT_CATEGORY(smpi_pmpi);
  15
  16 /* PMPI User level calls */
  17
  18 int PMPI_Barrier(MPI_Comm comm)
  19 {
  20   return PMPI_Ibarrier(comm, MPI_REQUEST_IGNORED);
  21 }
  22
  23 int PMPI_Ibarrier(MPI_Comm comm, MPI_Request *request)
  24 {
  25   if (comm == MPI_COMM_NULL)
  26     return MPI_ERR_COMM;
  27   if (request == nullptr)
  28     return MPI_ERR_ARG;
  29
  30   smpi_bench_end();
  31   int rank = simgrid::s4u::this_actor::get_pid();
  32   TRACE_smpi_comm_in(rank, request == MPI_REQUEST_IGNORED ? "PMPI_Barrier" : "PMPI_Ibarrier",
  33                      new simgrid::instr::NoOpTIData(request == MPI_REQUEST_IGNORED ? "barrier" : "ibarrier"));
  34   if (request == MPI_REQUEST_IGNORED) {
  35     simgrid::smpi::Colls::barrier(comm);
  36     // Barrier can be used to synchronize RMA calls. Finish all requests from comm before.
  37     comm->finish_rma_calls();
  38   } else
  39     simgrid::smpi::Colls::ibarrier(comm, request);
  40
  41   TRACE_smpi_comm_out(rank);
  42   smpi_bench_begin();
  43   return MPI_SUCCESS;
  44 }
  45
  46 int PMPI_Bcast(void *buf, int count, MPI_Datatype datatype, int root, MPI_Comm comm)
  47 {
  48   return PMPI_Ibcast(buf, count, datatype, root, comm, MPI_REQUEST_IGNORED);
  49 }
  50
  51 int PMPI_Ibcast(void *buf, int count, MPI_Datatype datatype,
  52                    int root, MPI_Comm comm, MPI_Request* request)
  53 {
  54   if (comm == MPI_COMM_NULL)
  55     return MPI_ERR_COMM;
  56   if (buf == nullptr && count > 0)
  57     return MPI_ERR_BUFFER;
  58   if (datatype == MPI_DATATYPE_NULL || not datatype->is_valid())
  59     return MPI_ERR_TYPE;
  60   if (count < 0)
  61     return MPI_ERR_COUNT;
  62   if (root < 0 || root >= comm->size())
  63     return MPI_ERR_ROOT;
  64   if (request == nullptr)
  65     return MPI_ERR_ARG;
  66
  67   smpi_bench_end();
  68   int rank = simgrid::s4u::this_actor::get_pid();
  69   TRACE_smpi_comm_in(rank, request == MPI_REQUEST_IGNORED ? "PMPI_Bcast" : "PMPI_Ibcast",
  70                      new simgrid::instr::CollTIData(request == MPI_REQUEST_IGNORED ? "bcast" : "ibcast", root, -1.0,
  71                                                     datatype->is_replayable() ? count : count * datatype->size(), -1,
  72                                                     simgrid::smpi::Datatype::encode(datatype), ""));
  73   if (comm->size() > 1) {
  74     if (request == MPI_REQUEST_IGNORED)
  75       simgrid::smpi::Colls::bcast(buf, count, datatype, root, comm);
  76     else
  77       simgrid::smpi::Colls::ibcast(buf, count, datatype, root, comm, request);
  78   } else {
  79     if (request != MPI_REQUEST_IGNORED)
  80       *request = MPI_REQUEST_NULL;
  81   }
  82
  83   TRACE_smpi_comm_out(rank);
  84   smpi_bench_begin();
  85   return MPI_SUCCESS;
  86 }
  87
  88 int PMPI_Gather(const void *sendbuf, int sendcount, MPI_Datatype sendtype,void *recvbuf, int recvcount, MPI_Datatype recvtype,
  89                 int root, MPI_Comm comm){
  90   return PMPI_Igather(sendbuf, sendcount, sendtype, recvbuf, recvcount, recvtype, root, comm, MPI_REQUEST_IGNORED);
  91 }
  92
  93 int PMPI_Igather(const void* sendbuf, int sendcount, MPI_Datatype sendtype, void* recvbuf, int recvcount,
  94                  MPI_Datatype recvtype, int root, MPI_Comm comm, MPI_Request* request)
  95 {
  96   if (comm == MPI_COMM_NULL)
  97     return MPI_ERR_COMM;
  98   if ((sendbuf == nullptr) || ((comm->rank() == root) && recvbuf == nullptr))
  99     return MPI_ERR_BUFFER;
 100   if (((sendbuf != MPI_IN_PLACE && sendcount > 0) && (sendtype == MPI_DATATYPE_NULL)) ||
 101       ((comm->rank() == root) && (recvtype == MPI_DATATYPE_NULL)))
 102     return MPI_ERR_TYPE;
 103   if (((sendbuf != MPI_IN_PLACE) && (sendcount < 0)) || ((comm->rank() == root) && (recvcount < 0)))
 104     return MPI_ERR_COUNT;
 105   if (root < 0 || root >= comm->size())
 106     return MPI_ERR_ROOT;
 107   if (request == nullptr)
 108     return MPI_ERR_ARG;
 109
 110   smpi_bench_end();
 111   const void* real_sendbuf   = sendbuf;
 112   int real_sendcount         = sendcount;
 113   MPI_Datatype real_sendtype = sendtype;
 114   if ((comm->rank() == root) && (sendbuf == MPI_IN_PLACE)) {
 115     real_sendcount = 0;
 116     real_sendtype  = recvtype;
 117   }
 118   int rank = simgrid::s4u::this_actor::get_pid();
 119
 120   TRACE_smpi_comm_in(rank, request == MPI_REQUEST_IGNORED ? "PMPI_Gather" : "PMPI_Igather",
 121                      new simgrid::instr::CollTIData(
 122                          request == MPI_REQUEST_IGNORED ? "gather" : "igather", root, -1.0,
 123                          real_sendtype->is_replayable() ? real_sendcount : real_sendcount * real_sendtype->size(),
 124                          (comm->rank() != root || recvtype->is_replayable()) ? recvcount : recvcount * recvtype->size(),
 125                          simgrid::smpi::Datatype::encode(real_sendtype), simgrid::smpi::Datatype::encode(recvtype)));
 126   if (request == MPI_REQUEST_IGNORED)
 127     simgrid::smpi::Colls::gather(real_sendbuf, real_sendcount, real_sendtype, recvbuf, recvcount, recvtype, root, comm);
 128   else
 129     simgrid::smpi::Colls::igather(real_sendbuf, real_sendcount, real_sendtype, recvbuf, recvcount, recvtype, root, comm,
 130                                   request);
 131
 132   TRACE_smpi_comm_out(rank);
 133   smpi_bench_begin();
 134   return MPI_SUCCESS;
 135 }
 136
 137 int PMPI_Gatherv(const void *sendbuf, int sendcount, MPI_Datatype sendtype, void *recvbuf, const int *recvcounts, const int *displs,
 138                 MPI_Datatype recvtype, int root, MPI_Comm comm){
 139   return PMPI_Igatherv(sendbuf, sendcount, sendtype, recvbuf, recvcounts, displs, recvtype, root, comm, MPI_REQUEST_IGNORED);
 140 }
 141
 142 int PMPI_Igatherv(const void* sendbuf, int sendcount, MPI_Datatype sendtype, void* recvbuf, const int* recvcounts, const int* displs,
 143                   MPI_Datatype recvtype, int root, MPI_Comm comm, MPI_Request* request)
 144 {
 145   if (comm == MPI_COMM_NULL)
 146     return MPI_ERR_COMM;
 147   if ((sendbuf == nullptr && sendcount > 0) || ((comm->rank() == root) && recvbuf == nullptr))
 148     return MPI_ERR_BUFFER;
 149   if (((sendbuf != MPI_IN_PLACE) && (sendtype == MPI_DATATYPE_NULL)) ||
 150       ((comm->rank() == root) && (recvtype == MPI_DATATYPE_NULL)))
 151     return MPI_ERR_TYPE;
 152   if ((sendbuf != MPI_IN_PLACE) && (sendcount < 0))
 153     return MPI_ERR_COUNT;
 154   if ((comm->rank() == root) && (recvcounts == nullptr || displs == nullptr))
 155     return MPI_ERR_ARG;
 156   if (root < 0 || root >= comm->size())
 157     return MPI_ERR_ROOT;
 158   if (request == nullptr)
 159     return MPI_ERR_ARG;
 160
 161   for (int i = 0; i < comm->size(); i++) {
 162     if ((comm->rank() == root) && (recvcounts[i] < 0))
 163       return MPI_ERR_COUNT;
 164   }
 165
 166   smpi_bench_end();
 167   const void* real_sendbuf   = sendbuf;
 168   int real_sendcount         = sendcount;
 169   MPI_Datatype real_sendtype = sendtype;
 170   if ((comm->rank() == root) && (sendbuf == MPI_IN_PLACE)) {
 171     real_sendcount = 0;
 172     real_sendtype  = recvtype;
 173   }
 174
 175   int rank         = simgrid::s4u::this_actor::get_pid();
 176   int dt_size_recv = recvtype->is_replayable() ? 1 : recvtype->size();
 177
 178   std::vector<int>* trace_recvcounts = new std::vector<int>;
 179   if (comm->rank() == root) {
 180     for (int i = 0; i < comm->size(); i++) // copy data to avoid bad free
 181       trace_recvcounts->push_back(recvcounts[i] * dt_size_recv);
 182   }
 183
 184   TRACE_smpi_comm_in(rank, request == MPI_REQUEST_IGNORED ? "PMPI_Gatherv" : "PMPI_Igatherv",
 185                      new simgrid::instr::VarCollTIData(
 186                          request == MPI_REQUEST_IGNORED ? "gatherv" : "igatherv", root,
 187                          real_sendtype->is_replayable() ? real_sendcount : real_sendcount * real_sendtype->size(),
 188                          nullptr, dt_size_recv, trace_recvcounts, simgrid::smpi::Datatype::encode(real_sendtype),
 189                          simgrid::smpi::Datatype::encode(recvtype)));
 190   if (request == MPI_REQUEST_IGNORED)
 191     simgrid::smpi::Colls::gatherv(real_sendbuf, real_sendcount, real_sendtype, recvbuf, recvcounts, displs, recvtype,
 192                                   root, comm);
 193   else
 194     simgrid::smpi::Colls::igatherv(real_sendbuf, real_sendcount, real_sendtype, recvbuf, recvcounts, displs, recvtype,
 195                                    root, comm, request);
 196
 197   TRACE_smpi_comm_out(rank);
 198   smpi_bench_begin();
 199   return MPI_SUCCESS;
 200 }
 201
 202 int PMPI_Allgather(const void *sendbuf, int sendcount, MPI_Datatype sendtype,
 203                    void *recvbuf, int recvcount, MPI_Datatype recvtype, MPI_Comm comm){
 204   return PMPI_Iallgather(sendbuf, sendcount, sendtype, recvbuf, recvcount, recvtype, comm, MPI_REQUEST_IGNORED);
 205 }
 206
 207 int PMPI_Iallgather(const void* sendbuf, int sendcount, MPI_Datatype sendtype, void* recvbuf, int recvcount,
 208                     MPI_Datatype recvtype, MPI_Comm comm, MPI_Request* request)
 209 {
 210   if (comm == MPI_COMM_NULL)
 211     return MPI_ERR_COMM;
 212   if ((sendbuf == nullptr && sendcount > 0) || (recvbuf == nullptr))
 213     return MPI_ERR_BUFFER;
 214   if (((sendbuf != MPI_IN_PLACE) && (sendtype == MPI_DATATYPE_NULL)) || (recvtype == MPI_DATATYPE_NULL))
 215     return MPI_ERR_TYPE;
 216   if (((sendbuf != MPI_IN_PLACE) && (sendcount < 0)) || (recvcount < 0))
 217     return MPI_ERR_COUNT;
 218   if (request == nullptr)
 219     return MPI_ERR_ARG;
 220
 221   smpi_bench_end();
 222   if (sendbuf == MPI_IN_PLACE) {
 223     sendbuf   = static_cast<char*>(recvbuf) + recvtype->get_extent() * recvcount * comm->rank();
 224     sendcount = recvcount;
 225     sendtype  = recvtype;
 226   }
 227   int rank = simgrid::s4u::this_actor::get_pid();
 228
 229   TRACE_smpi_comm_in(rank, request == MPI_REQUEST_IGNORED ? "PMPI_Allgather" : "PMPI_Iallggather",
 230                      new simgrid::instr::CollTIData(
 231                          request == MPI_REQUEST_IGNORED ? "allgather" : "iallgather", -1, -1.0,
 232                          sendtype->is_replayable() ? sendcount : sendcount * sendtype->size(),
 233                          recvtype->is_replayable() ? recvcount : recvcount * recvtype->size(),
 234                          simgrid::smpi::Datatype::encode(sendtype), simgrid::smpi::Datatype::encode(recvtype)));
 235   if (request == MPI_REQUEST_IGNORED)
 236     simgrid::smpi::Colls::allgather(sendbuf, sendcount, sendtype, recvbuf, recvcount, recvtype, comm);
 237   else
 238     simgrid::smpi::Colls::iallgather(sendbuf, sendcount, sendtype, recvbuf, recvcount, recvtype, comm, request);
 239
 240   TRACE_smpi_comm_out(rank);
 241   smpi_bench_begin();
 242   return MPI_SUCCESS;
 243 }
 244
 245 int PMPI_Allgatherv(const void *sendbuf, int sendcount, MPI_Datatype sendtype,
 246                    void *recvbuf, const int *recvcounts, const int *displs, MPI_Datatype recvtype, MPI_Comm comm){
 247   return PMPI_Iallgatherv(sendbuf, sendcount, sendtype, recvbuf, recvcounts, displs, recvtype, comm, MPI_REQUEST_IGNORED);
 248 }
 249
 250 int PMPI_Iallgatherv(const void* sendbuf, int sendcount, MPI_Datatype sendtype, void* recvbuf, const int* recvcounts, const int* displs,
 251                      MPI_Datatype recvtype, MPI_Comm comm, MPI_Request* request)
 252 {
 253   if (comm == MPI_COMM_NULL)
 254     return MPI_ERR_COMM;
 255   if ((sendbuf == nullptr && sendcount > 0) || (recvbuf == nullptr))
 256     return MPI_ERR_BUFFER;
 257   if (((sendbuf != MPI_IN_PLACE) && (sendtype == MPI_DATATYPE_NULL)) || (recvtype == MPI_DATATYPE_NULL))
 258     return MPI_ERR_TYPE;
 259   if ((sendbuf != MPI_IN_PLACE) && (sendcount < 0))
 260     return MPI_ERR_COUNT;
 261   if (recvcounts == nullptr || displs == nullptr)
 262     return MPI_ERR_ARG;
 263   if (request == nullptr)
 264     return MPI_ERR_ARG;
 265
 266   for (int i = 0; i < comm->size(); i++) { // copy data to avoid bad free
 267     if (recvcounts[i] < 0)
 268       return MPI_ERR_COUNT;
 269   }
 270
 271   smpi_bench_end();
 272   if (sendbuf == MPI_IN_PLACE) {
 273     sendbuf   = static_cast<char*>(recvbuf) + recvtype->get_extent() * displs[comm->rank()];
 274     sendcount = recvcounts[comm->rank()];
 275     sendtype  = recvtype;
 276   }
 277   int rank         = simgrid::s4u::this_actor::get_pid();
 278   int dt_size_recv = recvtype->is_replayable() ? 1 : recvtype->size();
 279
 280   std::vector<int>* trace_recvcounts = new std::vector<int>;
 281   for (int i = 0; i < comm->size(); i++) { // copy data to avoid bad free
 282     trace_recvcounts->push_back(recvcounts[i] * dt_size_recv);
 283   }
 284
 285   TRACE_smpi_comm_in(
 286       rank, request == MPI_REQUEST_IGNORED ? "PMPI_Allgatherv" : "PMPI_Iallgatherv",
 287       new simgrid::instr::VarCollTIData(request == MPI_REQUEST_IGNORED ? "allgatherv" : "iallgatherv", -1,
 288                                         sendtype->is_replayable() ? sendcount : sendcount * sendtype->size(), nullptr,
 289                                         dt_size_recv, trace_recvcounts, simgrid::smpi::Datatype::encode(sendtype),
 290                                         simgrid::smpi::Datatype::encode(recvtype)));
 291   if (request == MPI_REQUEST_IGNORED)
 292     simgrid::smpi::Colls::allgatherv(sendbuf, sendcount, sendtype, recvbuf, recvcounts, displs, recvtype, comm);
 293   else
 294     simgrid::smpi::Colls::iallgatherv(sendbuf, sendcount, sendtype, recvbuf, recvcounts, displs, recvtype, comm,
 295                                       request);
 296
 297   TRACE_smpi_comm_out(rank);
 298   smpi_bench_begin();
 299   return MPI_SUCCESS;
 300 }
 301
 302 int PMPI_Scatter(const void *sendbuf, int sendcount, MPI_Datatype sendtype,
 303                 void *recvbuf, int recvcount, MPI_Datatype recvtype, int root, MPI_Comm comm){
 304   return PMPI_Iscatter(sendbuf, sendcount, sendtype, recvbuf, recvcount, recvtype, root, comm, MPI_REQUEST_IGNORED);
 305 }
 306
 307 int PMPI_Iscatter(const void* sendbuf, int sendcount, MPI_Datatype sendtype, void* recvbuf, int recvcount,
 308                   MPI_Datatype recvtype, int root, MPI_Comm comm, MPI_Request* request)
 309 {
 310   if (comm == MPI_COMM_NULL)
 311     return MPI_ERR_COMM;
 312   if (((comm->rank() == root) && (sendtype == MPI_DATATYPE_NULL || not sendtype->is_valid())) ||
 313       ((recvbuf != MPI_IN_PLACE) && (recvtype == MPI_DATATYPE_NULL || not recvtype->is_valid())))
 314     return MPI_ERR_TYPE;
 315   if (((comm->rank() == root) && (sendcount < 0)) || ((recvbuf != MPI_IN_PLACE) && (recvcount < 0)))
 316     return MPI_ERR_COUNT;
 317   if ((sendbuf == recvbuf) || ((comm->rank() == root) && sendcount > 0 && (sendbuf == nullptr)) ||
 318       (recvcount > 0 && recvbuf == nullptr))
 319     return MPI_ERR_BUFFER;
 320   if (root < 0 || root >= comm->size())
 321     return MPI_ERR_ROOT;
 322   if (request == nullptr)
 323     return MPI_ERR_ARG;
 324
 325   smpi_bench_end();
 326   if (recvbuf == MPI_IN_PLACE) {
 327     recvtype  = sendtype;
 328     recvcount = sendcount;
 329   }
 330   int rank = simgrid::s4u::this_actor::get_pid();
 331
 332   TRACE_smpi_comm_in(rank, request == MPI_REQUEST_IGNORED ? "PMPI_Scatter" : "PMPI_Iscatter",
 333                      new simgrid::instr::CollTIData(
 334                          request == MPI_REQUEST_IGNORED ? "scatter" : "iscatter", root, -1.0,
 335                          (comm->rank() != root || sendtype->is_replayable()) ? sendcount : sendcount * sendtype->size(),
 336                          recvtype->is_replayable() ? recvcount : recvcount * recvtype->size(),
 337                          simgrid::smpi::Datatype::encode(sendtype), simgrid::smpi::Datatype::encode(recvtype)));
 338   if (request == MPI_REQUEST_IGNORED)
 339     simgrid::smpi::Colls::scatter(sendbuf, sendcount, sendtype, recvbuf, recvcount, recvtype, root, comm);
 340   else
 341     simgrid::smpi::Colls::iscatter(sendbuf, sendcount, sendtype, recvbuf, recvcount, recvtype, root, comm, request);
 342
 343   TRACE_smpi_comm_out(rank);
 344   smpi_bench_begin();
 345   return MPI_SUCCESS;
 346 }
 347
 348 int PMPI_Scatterv(const void *sendbuf, const int *sendcounts, const int *displs,
 349                  MPI_Datatype sendtype, void *recvbuf, int recvcount, MPI_Datatype recvtype, int root, MPI_Comm comm){
 350   return PMPI_Iscatterv(sendbuf, sendcounts, displs, sendtype, recvbuf, recvcount, recvtype, root, comm, MPI_REQUEST_IGNORED);
 351 }
 352
 353 int PMPI_Iscatterv(const void* sendbuf, const int* sendcounts, const int* displs, MPI_Datatype sendtype, void* recvbuf, int recvcount,
 354                    MPI_Datatype recvtype, int root, MPI_Comm comm, MPI_Request* request)
 355 {
 356   if (comm == MPI_COMM_NULL)
 357     return MPI_ERR_COMM;
 358   if (sendcounts == nullptr || displs == nullptr)
 359     return MPI_ERR_ARG;
 360   if (((comm->rank() == root) && (sendtype == MPI_DATATYPE_NULL)) ||
 361       ((recvbuf != MPI_IN_PLACE) && (recvtype == MPI_DATATYPE_NULL)))
 362     return MPI_ERR_TYPE;
 363   if (request == nullptr)
 364     return MPI_ERR_ARG;
 365   if (recvbuf != MPI_IN_PLACE && recvcount < 0)
 366     return MPI_ERR_COUNT;
 367   if (root < 0 || root >= comm->size())
 368     return MPI_ERR_ROOT;
 369
 370   if (comm->rank() == root) {
 371     if (recvbuf == MPI_IN_PLACE) {
 372       recvtype  = sendtype;
 373       recvcount = sendcounts[comm->rank()];
 374     }
 375     for (int i = 0; i < comm->size(); i++) {
 376       if (sendcounts[i] < 0)
 377         return MPI_ERR_COUNT;
 378     }
 379   }
 380
 381   smpi_bench_end();
 382
 383   int rank         = simgrid::s4u::this_actor::get_pid();
 384   int dt_size_send = sendtype->is_replayable() ? 1 : sendtype->size();
 385
 386   std::vector<int>* trace_sendcounts = new std::vector<int>;
 387   if (comm->rank() == root) {
 388     for (int i = 0; i < comm->size(); i++) { // copy data to avoid bad free
 389       trace_sendcounts->push_back(sendcounts[i] * dt_size_send);
 390     }
 391   }
 392
 393   TRACE_smpi_comm_in(rank, request == MPI_REQUEST_IGNORED ? "PMPI_Scatterv" : "PMPI_Iscatterv",
 394                      new simgrid::instr::VarCollTIData(
 395                          request == MPI_REQUEST_IGNORED ? "scatterv" : "iscatterv", root, dt_size_send,
 396                          trace_sendcounts, recvtype->is_replayable() ? recvcount : recvcount * recvtype->size(),
 397                          nullptr, simgrid::smpi::Datatype::encode(sendtype),
 398                          simgrid::smpi::Datatype::encode(recvtype)));
 399   if (request == MPI_REQUEST_IGNORED)
 400     simgrid::smpi::Colls::scatterv(sendbuf, sendcounts, displs, sendtype, recvbuf, recvcount, recvtype, root, comm);
 401   else
 402     simgrid::smpi::Colls::iscatterv(sendbuf, sendcounts, displs, sendtype, recvbuf, recvcount, recvtype, root, comm,
 403                                     request);
 404
 405   TRACE_smpi_comm_out(rank);
 406   smpi_bench_begin();
 407   return MPI_SUCCESS;
 408 }
 409
 410 int PMPI_Reduce(const void *sendbuf, void *recvbuf, int count, MPI_Datatype datatype, MPI_Op op, int root, MPI_Comm comm)
 411 {
 412   return PMPI_Ireduce(sendbuf, recvbuf, count, datatype, op, root, comm, MPI_REQUEST_IGNORED);
 413 }
 414
 415 int PMPI_Ireduce(const void *sendbuf, void *recvbuf, int count, MPI_Datatype datatype, MPI_Op op, int root, MPI_Comm comm, MPI_Request* request)
 416 {
 417   if (comm == MPI_COMM_NULL)
 418     return MPI_ERR_COMM;
 419   if ((sendbuf == nullptr && count > 0) || ((comm->rank() == root) && recvbuf == nullptr))
 420     return MPI_ERR_BUFFER;
 421   if (datatype == MPI_DATATYPE_NULL || not datatype->is_valid())
 422     return MPI_ERR_TYPE;
 423   if (op == MPI_OP_NULL)
 424     return MPI_ERR_OP;
 425   if (request == nullptr)
 426     return MPI_ERR_ARG;
 427   if (root < 0 || root >= comm->size())
 428     return MPI_ERR_ROOT;
 429   if (count < 0)
 430     return MPI_ERR_COUNT;
 431
 432   smpi_bench_end();
 433   int rank = simgrid::s4u::this_actor::get_pid();
 434
 435   TRACE_smpi_comm_in(rank, request == MPI_REQUEST_IGNORED ? "PMPI_Reduce" : "PMPI_Ireduce",
 436                      new simgrid::instr::CollTIData(request == MPI_REQUEST_IGNORED ? "reduce" : "ireduce", root, 0,
 437                                                     datatype->is_replayable() ? count : count * datatype->size(), -1,
 438                                                     simgrid::smpi::Datatype::encode(datatype), ""));
 439   if (request == MPI_REQUEST_IGNORED)
 440     simgrid::smpi::Colls::reduce(sendbuf, recvbuf, count, datatype, op, root, comm);
 441   else
 442     simgrid::smpi::Colls::ireduce(sendbuf, recvbuf, count, datatype, op, root, comm, request);
 443
 444   TRACE_smpi_comm_out(rank);
 445   smpi_bench_begin();
 446   return MPI_SUCCESS;
 447 }
 448
 449 int PMPI_Reduce_local(const void* inbuf, void* inoutbuf, int count, MPI_Datatype datatype, MPI_Op op)
 450 {
 451   if (datatype == MPI_DATATYPE_NULL || not datatype->is_valid())
 452     return MPI_ERR_TYPE;
 453   if (op == MPI_OP_NULL)
 454     return MPI_ERR_OP;
 455   if (count < 0)
 456     return MPI_ERR_COUNT;
 457
 458   smpi_bench_end();
 459   op->apply(inbuf, inoutbuf, &count, datatype);
 460   smpi_bench_begin();
 461   return MPI_SUCCESS;
 462 }
 463
 464 int PMPI_Allreduce(const void *sendbuf, void *recvbuf, int count, MPI_Datatype datatype, MPI_Op op, MPI_Comm comm)
 465 {
 466   return PMPI_Iallreduce(sendbuf, recvbuf, count, datatype, op, comm, MPI_REQUEST_IGNORED);
 467 }
 468
 469 int PMPI_Iallreduce(const void *sendbuf, void *recvbuf, int count, MPI_Datatype datatype, MPI_Op op, MPI_Comm comm, MPI_Request *request)
 470 {
 471   if (comm == MPI_COMM_NULL)
 472     return MPI_ERR_COMM;
 473   if ((sendbuf == nullptr && count > 0) || (recvbuf == nullptr))
 474     return MPI_ERR_BUFFER;
 475   if (datatype == MPI_DATATYPE_NULL || not datatype->is_valid())
 476     return MPI_ERR_TYPE;
 477   if (count < 0)
 478     return MPI_ERR_COUNT;
 479   if (op == MPI_OP_NULL)
 480     return MPI_ERR_OP;
 481   if (request == nullptr)
 482     return MPI_ERR_ARG;
 483
 484   smpi_bench_end();
 485   const void* real_sendbuf = sendbuf;
 486   std::unique_ptr<unsigned char[]> tmp_sendbuf;
 487   if (sendbuf == MPI_IN_PLACE) {
 488     tmp_sendbuf.reset(new unsigned char[count * datatype->get_extent()]);
 489     simgrid::smpi::Datatype::copy(recvbuf, count, datatype, tmp_sendbuf.get(), count, datatype);
 490     real_sendbuf = tmp_sendbuf.get();
 491   }
 492   int rank = simgrid::s4u::this_actor::get_pid();
 493
 494   TRACE_smpi_comm_in(rank, request == MPI_REQUEST_IGNORED ? "PMPI_Allreduce" : "PMPI_Iallreduce",
 495                      new simgrid::instr::CollTIData(request == MPI_REQUEST_IGNORED ? "allreduce" : "iallreduce", -1, 0,
 496                                                     datatype->is_replayable() ? count : count * datatype->size(), -1,
 497                                                     simgrid::smpi::Datatype::encode(datatype), ""));
 498
 499   if (request == MPI_REQUEST_IGNORED)
 500     simgrid::smpi::Colls::allreduce(real_sendbuf, recvbuf, count, datatype, op, comm);
 501   else
 502     simgrid::smpi::Colls::iallreduce(real_sendbuf, recvbuf, count, datatype, op, comm, request);
 503
 504   TRACE_smpi_comm_out(rank);
 505   smpi_bench_begin();
 506   return MPI_SUCCESS;
 507 }
 508
 509 int PMPI_Scan(const void *sendbuf, void *recvbuf, int count, MPI_Datatype datatype, MPI_Op op, MPI_Comm comm)
 510 {
 511   return PMPI_Iscan(sendbuf, recvbuf, count, datatype, op, comm, MPI_REQUEST_IGNORED);
 512 }
 513
 514 int PMPI_Iscan(const void *sendbuf, void *recvbuf, int count, MPI_Datatype datatype, MPI_Op op, MPI_Comm comm, MPI_Request* request)
 515 {
 516   if (comm == MPI_COMM_NULL)
 517     return MPI_ERR_COMM;
 518   if (datatype == MPI_DATATYPE_NULL || not datatype->is_valid())
 519     return MPI_ERR_TYPE;
 520   if (op == MPI_OP_NULL)
 521     return MPI_ERR_OP;
 522   if (request == nullptr)
 523     return MPI_ERR_ARG;
 524   if (count < 0)
 525     return MPI_ERR_COUNT;
 526   if (sendbuf == nullptr || recvbuf == nullptr)
 527     return MPI_ERR_BUFFER;
 528
 529   smpi_bench_end();
 530   int rank         = simgrid::s4u::this_actor::get_pid();
 531   const void* real_sendbuf = sendbuf;
 532   std::unique_ptr<unsigned char[]> tmp_sendbuf;
 533   if (sendbuf == MPI_IN_PLACE) {
 534     tmp_sendbuf.reset(new unsigned char[count * datatype->size()]);
 535     real_sendbuf = memcpy(tmp_sendbuf.get(), recvbuf, count * datatype->size());
 536   }
 537   TRACE_smpi_comm_in(rank, request == MPI_REQUEST_IGNORED ? "PMPI_Scan" : "PMPI_Iscan",
 538                      new simgrid::instr::Pt2PtTIData(request == MPI_REQUEST_IGNORED ? "scan" : "iscan", -1,
 539                                                      datatype->is_replayable() ? count : count * datatype->size(),
 540                                                      simgrid::smpi::Datatype::encode(datatype)));
 541
 542   int retval;
 543   if (request == MPI_REQUEST_IGNORED)
 544     retval = simgrid::smpi::Colls::scan(real_sendbuf, recvbuf, count, datatype, op, comm);
 545   else
 546     retval = simgrid::smpi::Colls::iscan(real_sendbuf, recvbuf, count, datatype, op, comm, request);
 547
 548   TRACE_smpi_comm_out(rank);
 549   smpi_bench_begin();
 550   return retval;
 551 }
 552
 553 int PMPI_Exscan(const void *sendbuf, void *recvbuf, int count, MPI_Datatype datatype, MPI_Op op, MPI_Comm comm)
 554 {
 555   return PMPI_Iexscan(sendbuf, recvbuf, count, datatype, op, comm, MPI_REQUEST_IGNORED);
 556 }
 557
 558 int PMPI_Iexscan(const void *sendbuf, void *recvbuf, int count, MPI_Datatype datatype, MPI_Op op, MPI_Comm comm, MPI_Request* request){
 559   if (comm == MPI_COMM_NULL)
 560     return MPI_ERR_COMM;
 561   if (not datatype->is_valid())
 562     return MPI_ERR_TYPE;
 563   if (op == MPI_OP_NULL)
 564     return MPI_ERR_OP;
 565   if (request == nullptr)
 566     return MPI_ERR_ARG;
 567   if (count < 0)
 568     return MPI_ERR_COUNT;
 569   if (sendbuf == nullptr || recvbuf == nullptr)
 570     return MPI_ERR_BUFFER;
 571
 572   smpi_bench_end();
 573   int rank         = simgrid::s4u::this_actor::get_pid();
 574   const void* real_sendbuf = sendbuf;
 575   std::unique_ptr<unsigned char[]> tmp_sendbuf;
 576   if (sendbuf == MPI_IN_PLACE) {
 577     tmp_sendbuf.reset(new unsigned char[count * datatype->size()]);
 578     real_sendbuf = memcpy(tmp_sendbuf.get(), recvbuf, count * datatype->size());
 579   }
 580
 581   TRACE_smpi_comm_in(rank, request == MPI_REQUEST_IGNORED ? "PMPI_Exscan" : "PMPI_Iexscan",
 582                      new simgrid::instr::Pt2PtTIData(request == MPI_REQUEST_IGNORED ? "exscan" : "iexscan", -1,
 583                                                      datatype->is_replayable() ? count : count * datatype->size(),
 584                                                      simgrid::smpi::Datatype::encode(datatype)));
 585
 586   int retval;
 587   if (request == MPI_REQUEST_IGNORED)
 588     retval = simgrid::smpi::Colls::exscan(real_sendbuf, recvbuf, count, datatype, op, comm);
 589   else
 590     retval = simgrid::smpi::Colls::iexscan(real_sendbuf, recvbuf, count, datatype, op, comm, request);
 591
 592   TRACE_smpi_comm_out(rank);
 593   smpi_bench_begin();
 594   return retval;
 595 }
 596
 597 int PMPI_Reduce_scatter(const void *sendbuf, void *recvbuf, const int *recvcounts, MPI_Datatype datatype, MPI_Op op, MPI_Comm comm)
 598 {
 599   return PMPI_Ireduce_scatter(sendbuf, recvbuf, recvcounts, datatype, op, comm, MPI_REQUEST_IGNORED);
 600 }
 601
 602 int PMPI_Ireduce_scatter(const void *sendbuf, void *recvbuf, const int *recvcounts, MPI_Datatype datatype, MPI_Op op, MPI_Comm comm, MPI_Request *request)
 603 {
 604   if (comm == MPI_COMM_NULL)
 605     return MPI_ERR_COMM;
 606   if ((sendbuf == nullptr) || (recvbuf == nullptr))
 607     return MPI_ERR_BUFFER;
 608   if (datatype == MPI_DATATYPE_NULL || not datatype->is_valid())
 609     return MPI_ERR_TYPE;
 610   if (op == MPI_OP_NULL)
 611     return MPI_ERR_OP;
 612   if (recvcounts == nullptr)
 613     return MPI_ERR_ARG;
 614   if (request == nullptr)
 615     return MPI_ERR_ARG;
 616
 617   for (int i = 0; i < comm->size(); i++) { // copy data to avoid bad free
 618     if (recvcounts[i] < 0)
 619       return MPI_ERR_COUNT;
 620   }
 621
 622   smpi_bench_end();
 623   int rank                           = simgrid::s4u::this_actor::get_pid();
 624   std::vector<int>* trace_recvcounts = new std::vector<int>;
 625   int dt_send_size                   = datatype->is_replayable() ? 1 : datatype->size();
 626   int totalcount                     = 0;
 627
 628   for (int i = 0; i < comm->size(); i++) { // copy data to avoid bad free
 629     trace_recvcounts->push_back(recvcounts[i] * dt_send_size);
 630     totalcount += recvcounts[i];
 631   }
 632
 633   const void* real_sendbuf = sendbuf;
 634   std::unique_ptr<unsigned char[]> tmp_sendbuf;
 635   if (sendbuf == MPI_IN_PLACE) {
 636     tmp_sendbuf.reset(new unsigned char[totalcount * datatype->size()]);
 637     real_sendbuf = memcpy(tmp_sendbuf.get(), recvbuf, totalcount * datatype->size());
 638   }
 639
 640   TRACE_smpi_comm_in(rank, request == MPI_REQUEST_IGNORED ? "PMPI_Reduce_scatter" : "PMPI_Ireduce_scatter",
 641                      new simgrid::instr::VarCollTIData(
 642                          request == MPI_REQUEST_IGNORED ? "reducescatter" : "ireducescatter", -1, dt_send_size, nullptr,
 643                          -1, trace_recvcounts, simgrid::smpi::Datatype::encode(datatype), ""));
 644
 645   if (request == MPI_REQUEST_IGNORED)
 646     simgrid::smpi::Colls::reduce_scatter(real_sendbuf, recvbuf, recvcounts, datatype, op, comm);
 647   else
 648     simgrid::smpi::Colls::ireduce_scatter(real_sendbuf, recvbuf, recvcounts, datatype, op, comm, request);
 649
 650   TRACE_smpi_comm_out(rank);
 651   smpi_bench_begin();
 652   return MPI_SUCCESS;
 653 }
 654
 655 int PMPI_Reduce_scatter_block(const void *sendbuf, void *recvbuf, int recvcount,
 656                               MPI_Datatype datatype, MPI_Op op, MPI_Comm comm)
 657 {
 658   return PMPI_Ireduce_scatter_block(sendbuf, recvbuf, recvcount, datatype, op, comm, MPI_REQUEST_IGNORED);
 659 }
 660
 661 int PMPI_Ireduce_scatter_block(const void* sendbuf, void* recvbuf, int recvcount, MPI_Datatype datatype, MPI_Op op,
 662                                MPI_Comm comm, MPI_Request* request)
 663 {
 664   if (comm == MPI_COMM_NULL)
 665     return MPI_ERR_COMM;
 666   if (not datatype->is_valid())
 667     return MPI_ERR_TYPE;
 668   if (op == MPI_OP_NULL)
 669     return MPI_ERR_OP;
 670   if (recvcount < 0)
 671     return MPI_ERR_ARG;
 672   if (request == nullptr)
 673     return MPI_ERR_ARG;
 674
 675   smpi_bench_end();
 676   int count = comm->size();
 677
 678   int rank                           = simgrid::s4u::this_actor::get_pid();
 679   int dt_send_size                   = datatype->is_replayable() ? 1 : datatype->size();
 680   std::vector<int>* trace_recvcounts = new std::vector<int>(recvcount * dt_send_size); // copy data to avoid bad free
 681
 682   const void* real_sendbuf = sendbuf;
 683   std::unique_ptr<unsigned char[]> tmp_sendbuf;
 684   if (sendbuf == MPI_IN_PLACE) {
 685     tmp_sendbuf.reset(new unsigned char[recvcount * count * datatype->size()]);
 686     real_sendbuf = memcpy(tmp_sendbuf.get(), recvbuf, recvcount * count * datatype->size());
 687   }
 688
 689   TRACE_smpi_comm_in(
 690       rank, request == MPI_REQUEST_IGNORED ? "PMPI_Reduce_scatter_block" : "PMPI_Ireduce_scatter_block",
 691       new simgrid::instr::VarCollTIData(request == MPI_REQUEST_IGNORED ? "reducescatter" : "ireducescatter", -1, 0,
 692                                         nullptr, -1, trace_recvcounts, simgrid::smpi::Datatype::encode(datatype), ""));
 693
 694   int* recvcounts = new int[count];
 695   for (int i      = 0; i < count; i++)
 696     recvcounts[i] = recvcount;
 697   if (request == MPI_REQUEST_IGNORED)
 698     simgrid::smpi::Colls::reduce_scatter(real_sendbuf, recvbuf, recvcounts, datatype, op, comm);
 699   else
 700     simgrid::smpi::Colls::ireduce_scatter(real_sendbuf, recvbuf, recvcounts, datatype, op, comm, request);
 701   delete[] recvcounts;
 702
 703   TRACE_smpi_comm_out(rank);
 704   smpi_bench_begin();
 705   return MPI_SUCCESS;
 706 }
 707
 708 int PMPI_Alltoall(const void* sendbuf, int sendcount, MPI_Datatype sendtype, void* recvbuf, int recvcount,
 709                   MPI_Datatype recvtype, MPI_Comm comm){
 710   return PMPI_Ialltoall(sendbuf, sendcount, sendtype, recvbuf, recvcount, recvtype, comm, MPI_REQUEST_IGNORED);
 711 }
 712
 713 int PMPI_Ialltoall(const void* sendbuf, int sendcount, MPI_Datatype sendtype, void* recvbuf, int recvcount,
 714                    MPI_Datatype recvtype, MPI_Comm comm, MPI_Request* request)
 715 {
 716   if (comm == MPI_COMM_NULL)
 717     return MPI_ERR_COMM;
 718   if ((sendbuf == nullptr && sendcount > 0) || (recvbuf == nullptr && recvcount > 0))
 719     return MPI_ERR_BUFFER;
 720   if ((sendbuf != MPI_IN_PLACE && sendtype == MPI_DATATYPE_NULL) || recvtype == MPI_DATATYPE_NULL)
 721     return MPI_ERR_TYPE;
 722   if ((sendbuf != MPI_IN_PLACE && sendcount < 0) || recvcount < 0)
 723     return MPI_ERR_COUNT;
 724   if (request == nullptr)
 725     return MPI_ERR_ARG;
 726
 727   smpi_bench_end();
 728   int rank                 = simgrid::s4u::this_actor::get_pid();
 729   const void* real_sendbuf = sendbuf;
 730   int real_sendcount         = sendcount;
 731   MPI_Datatype real_sendtype = sendtype;
 732   std::unique_ptr<unsigned char[]> tmp_sendbuf;
 733   if (sendbuf == MPI_IN_PLACE) {
 734     tmp_sendbuf.reset(new unsigned char[recvcount * comm->size() * recvtype->size()]);
 735     // memcpy(??,nullptr,0) is actually undefined behavor, even if harmless.
 736     if (recvbuf != nullptr)
 737       memcpy(tmp_sendbuf.get(), recvbuf, recvcount * comm->size() * recvtype->size());
 738     real_sendbuf = tmp_sendbuf.get();
 739     real_sendcount = recvcount;
 740     real_sendtype  = recvtype;
 741   }
 742
 743   TRACE_smpi_comm_in(rank, request == MPI_REQUEST_IGNORED ? "PMPI_Alltoall" : "PMPI_Ialltoall",
 744                      new simgrid::instr::CollTIData(
 745                          request == MPI_REQUEST_IGNORED ? "alltoall" : "ialltoall", -1, -1.0,
 746                          real_sendtype->is_replayable() ? real_sendcount : real_sendcount * real_sendtype->size(),
 747                          recvtype->is_replayable() ? recvcount : recvcount * recvtype->size(),
 748                          simgrid::smpi::Datatype::encode(real_sendtype), simgrid::smpi::Datatype::encode(recvtype)));
 749   int retval;
 750   if (request == MPI_REQUEST_IGNORED)
 751     retval =
 752         simgrid::smpi::Colls::alltoall(real_sendbuf, real_sendcount, real_sendtype, recvbuf, recvcount, recvtype, comm);
 753   else
 754     retval = simgrid::smpi::Colls::ialltoall(real_sendbuf, real_sendcount, real_sendtype, recvbuf, recvcount, recvtype,
 755                                              comm, request);
 756
 757   TRACE_smpi_comm_out(rank);
 758   smpi_bench_begin();
 759   return retval;
 760 }
 761
 762 int PMPI_Alltoallv(const void* sendbuf, const int* sendcounts, const int* senddisps, MPI_Datatype sendtype, void* recvbuf,
 763                    const int* recvcounts, const int* recvdisps, MPI_Datatype recvtype, MPI_Comm comm)
 764 {
 765   return PMPI_Ialltoallv(sendbuf, sendcounts, senddisps, sendtype, recvbuf, recvcounts, recvdisps, recvtype, comm, MPI_REQUEST_IGNORED);
 766 }
 767
 768 int PMPI_Ialltoallv(const void* sendbuf, const int* sendcounts, const int* senddisps, MPI_Datatype sendtype, void* recvbuf,
 769                     const int* recvcounts, const int* recvdisps, MPI_Datatype recvtype, MPI_Comm comm, MPI_Request* request)
 770 {
 771   if (comm == MPI_COMM_NULL)
 772     return MPI_ERR_COMM;
 773   if (sendbuf == nullptr || recvbuf == nullptr)
 774     return MPI_ERR_BUFFER;
 775   if ((sendbuf != MPI_IN_PLACE && sendtype == MPI_DATATYPE_NULL) || recvtype == MPI_DATATYPE_NULL)
 776     return MPI_ERR_TYPE;
 777   if ((sendbuf != MPI_IN_PLACE && (sendcounts == nullptr || senddisps == nullptr)) || recvcounts == nullptr ||
 778       recvdisps == nullptr)
 779     return MPI_ERR_ARG;
 780   if (request == nullptr)
 781     return MPI_ERR_ARG;
 782
 783   int rank = simgrid::s4u::this_actor::get_pid();
 784   int size = comm->size();
 785   for (int i = 0; i < size; i++) {
 786     if (recvcounts[i] < 0 || (sendbuf != MPI_IN_PLACE && sendcounts[i] < 0))
 787       return MPI_ERR_COUNT;
 788   }
 789
 790   smpi_bench_end();
 791   int send_size                      = 0;
 792   int recv_size                      = 0;
 793   std::vector<int>* trace_sendcounts = new std::vector<int>;
 794   std::vector<int>* trace_recvcounts = new std::vector<int>;
 795   int dt_size_recv                   = recvtype->size();
 796
 797   const void* real_sendbuf   = sendbuf;
 798   const int* real_sendcounts = sendcounts;
 799   const int* real_senddisps  = senddisps;
 800   MPI_Datatype real_sendtype = sendtype;
 801   int maxsize              = 0;
 802   for (int i = 0; i < size; i++) { // copy data to avoid bad free
 803     recv_size += recvcounts[i] * dt_size_recv;
 804     trace_recvcounts->push_back(recvcounts[i] * dt_size_recv);
 805     if (((recvdisps[i] + recvcounts[i]) * dt_size_recv) > maxsize)
 806       maxsize = (recvdisps[i] + recvcounts[i]) * dt_size_recv;
 807   }
 808
 809   std::unique_ptr<unsigned char[]> tmp_sendbuf;
 810   std::unique_ptr<int[]> tmp_sendcounts;
 811   std::unique_ptr<int[]> tmp_senddisps;
 812   if (sendbuf == MPI_IN_PLACE) {
 813     tmp_sendbuf.reset(new unsigned char[maxsize]);
 814     real_sendbuf = memcpy(tmp_sendbuf.get(), recvbuf, maxsize);
 815     tmp_sendcounts.reset(new int[size]);
 816     std::copy(recvcounts, recvcounts + size, tmp_sendcounts.get());
 817     real_sendcounts = tmp_sendcounts.get();
 818     tmp_senddisps.reset(new int[size]);
 819     std::copy(recvdisps, recvdisps + size, tmp_senddisps.get());
 820     real_senddisps = tmp_senddisps.get();
 821     real_sendtype  = recvtype;
 822   }
 823
 824   int dt_size_send = real_sendtype->size();
 825
 826   for (int i = 0; i < size; i++) { // copy data to avoid bad free
 827     send_size += real_sendcounts[i] * dt_size_send;
 828     trace_sendcounts->push_back(real_sendcounts[i] * dt_size_send);
 829   }
 830
 831   TRACE_smpi_comm_in(rank, request == MPI_REQUEST_IGNORED ? "PMPI_Alltoallv" : "PMPI_Ialltoallv",
 832                      new simgrid::instr::VarCollTIData(request == MPI_REQUEST_IGNORED ? "alltoallv" : "ialltoallv", -1,
 833                                                        send_size, trace_sendcounts, recv_size, trace_recvcounts,
 834                                                        simgrid::smpi::Datatype::encode(real_sendtype),
 835                                                        simgrid::smpi::Datatype::encode(recvtype)));
 836
 837   int retval;
 838   if (request == MPI_REQUEST_IGNORED)
 839     retval = simgrid::smpi::Colls::alltoallv(real_sendbuf, real_sendcounts, real_senddisps, real_sendtype, recvbuf,
 840                                              recvcounts, recvdisps, recvtype, comm);
 841   else
 842     retval = simgrid::smpi::Colls::ialltoallv(real_sendbuf, real_sendcounts, real_senddisps, real_sendtype, recvbuf,
 843                                               recvcounts, recvdisps, recvtype, comm, request);
 844
 845   TRACE_smpi_comm_out(rank);
 846   smpi_bench_begin();
 847   return retval;
 848 }
 849
 850 int PMPI_Alltoallw(const void* sendbuf, const int* sendcounts, const int* senddisps, const MPI_Datatype* sendtypes, void* recvbuf,
 851                    const int* recvcounts, const int* recvdisps, const MPI_Datatype* recvtypes, MPI_Comm comm)
 852 {
 853   return PMPI_Ialltoallw(sendbuf, sendcounts, senddisps, sendtypes, recvbuf, recvcounts, recvdisps, recvtypes, comm, MPI_REQUEST_IGNORED);
 854 }
 855
 856 int PMPI_Ialltoallw(const void* sendbuf, const int* sendcounts, const int* senddisps, const MPI_Datatype* sendtypes, void* recvbuf,
 857                     const int* recvcounts, const int* recvdisps, const MPI_Datatype* recvtypes, MPI_Comm comm, MPI_Request* request)
 858 {
 859   if (comm == MPI_COMM_NULL)
 860     return MPI_ERR_COMM;
 861   if (sendbuf == nullptr || recvbuf == nullptr)
 862     return MPI_ERR_BUFFER;
 863   if ((sendbuf != MPI_IN_PLACE && sendtypes == nullptr) || recvtypes == nullptr)
 864     return MPI_ERR_TYPE;
 865   if ((sendbuf != MPI_IN_PLACE && (sendcounts == nullptr || senddisps == nullptr)) || recvcounts == nullptr ||
 866       recvdisps == nullptr)
 867     return MPI_ERR_ARG;
 868   if (request == nullptr)
 869     return MPI_ERR_ARG;
 870
 871   smpi_bench_end();
 872   int rank = simgrid::s4u::this_actor::get_pid();
 873   int size = comm->size();
 874   for (int i = 0; i < size; i++) { // copy data to avoid bad free
 875     if (recvcounts[i] < 0 || (sendbuf != MPI_IN_PLACE && sendcounts[i] < 0))
 876       return MPI_ERR_COUNT;
 877   }
 878   int send_size                      = 0;
 879   int recv_size                      = 0;
 880   std::vector<int>* trace_sendcounts = new std::vector<int>;
 881   std::vector<int>* trace_recvcounts = new std::vector<int>;
 882
 883   const void* real_sendbuf           = sendbuf;
 884   const int* real_sendcounts         = sendcounts;
 885   const int* real_senddisps          = senddisps;
 886   const MPI_Datatype* real_sendtypes = sendtypes;
 887   unsigned long maxsize      = 0;
 888   for (int i = 0; i < size; i++) { // copy data to avoid bad free
 889     if (recvtypes[i] == MPI_DATATYPE_NULL) {
 890       delete trace_recvcounts;
 891       delete trace_sendcounts;
 892       return MPI_ERR_TYPE;
 893     }
 894     recv_size += recvcounts[i] * recvtypes[i]->size();
 895     trace_recvcounts->push_back(recvcounts[i] * recvtypes[i]->size());
 896     if ((recvdisps[i] + (recvcounts[i] * recvtypes[i]->size())) > maxsize)
 897       maxsize = recvdisps[i] + (recvcounts[i] * recvtypes[i]->size());
 898   }
 899
 900   std::unique_ptr<unsigned char[]> tmp_sendbuf;
 901   std::unique_ptr<int[]> tmp_sendcounts;
 902   std::unique_ptr<int[]> tmp_senddisps;
 903   std::unique_ptr<MPI_Datatype[]> tmp_sendtypes;
 904   if (sendbuf == MPI_IN_PLACE) {
 905     tmp_sendbuf.reset(new unsigned char[maxsize]);
 906     real_sendbuf = memcpy(tmp_sendbuf.get(), recvbuf, maxsize);
 907     tmp_sendcounts.reset(new int[size]);
 908     std::copy(recvcounts, recvcounts + size, tmp_sendcounts.get());
 909     real_sendcounts = tmp_sendcounts.get();
 910     tmp_senddisps.reset(new int[size]);
 911     std::copy(recvdisps, recvdisps + size, tmp_senddisps.get());
 912     real_senddisps = tmp_senddisps.get();
 913     tmp_sendtypes.reset(new MPI_Datatype[size]);
 914     std::copy(recvtypes, recvtypes + size, tmp_sendtypes.get());
 915     real_sendtypes = tmp_sendtypes.get();
 916   }
 917
 918   for (int i = 0; i < size; i++) { // copy data to avoid bad free
 919     send_size += real_sendcounts[i] * real_sendtypes[i]->size();
 920     trace_sendcounts->push_back(real_sendcounts[i] * real_sendtypes[i]->size());
 921   }
 922
 923   TRACE_smpi_comm_in(rank, request == MPI_REQUEST_IGNORED ? "PMPI_Alltoallw" : "PMPI_Ialltoallw",
 924                      new simgrid::instr::VarCollTIData(request == MPI_REQUEST_IGNORED ? "alltoallv" : "ialltoallv", -1,
 925                                                        send_size, trace_sendcounts, recv_size, trace_recvcounts,
 926                                                        simgrid::smpi::Datatype::encode(real_sendtypes[0]),
 927                                                        simgrid::smpi::Datatype::encode(recvtypes[0])));
 928
 929   int retval;
 930   if (request == MPI_REQUEST_IGNORED)
 931     retval = simgrid::smpi::Colls::alltoallw(real_sendbuf, real_sendcounts, real_senddisps, real_sendtypes, recvbuf,
 932                                              recvcounts, recvdisps, recvtypes, comm);
 933   else
 934     retval = simgrid::smpi::Colls::ialltoallw(real_sendbuf, real_sendcounts, real_senddisps, real_sendtypes, recvbuf,
 935                                               recvcounts, recvdisps, recvtypes, comm, request);
 936
 937   TRACE_smpi_comm_out(rank);
 938   smpi_bench_begin();
 939   return retval;
 940 }