-/* Copyright (c) 2013-2017. The SimGrid Team.
+/* Copyright (c) 2013-2019. The SimGrid Team.
* All rights reserved. */
/* This program is free software; you can redistribute it and/or modify it
* under the terms of the license (GNU LGPL) which comes with this package. */
-#include "../colls_private.h"
+#include "../colls_private.hpp"
/*****************************************************************************
static int is_2dmesh(int num, int *i, int *j)
{
int x, max = num / 2;
- x = sqrt(num);
+ x = sqrt(double(num));
while (x <= max) {
if ((num % x) == 0) {
* comm: communication
* Descrp: Function realizes the allgather operation using the 2dmesh
* algorithm. Allgather ommunication occurs first in the x dimension then in
- * the y dimension. The communication in each dimension follows
+ * the y dimension. The communication in each dimension follows
* "simple"
* Auther: Ahmad Faraj
****************************************************************************/
int
-Coll_allgather_2dmesh::allgather(void *send_buff, int send_count, MPI_Datatype
+Coll_allgather_2dmesh::allgather(const void *send_buff, int send_count, MPI_Datatype
send_type, void *recv_buff, int recv_count,
MPI_Datatype recv_type, MPI_Comm comm)
{
- MPI_Request *req, *req_ptr;
MPI_Aint extent;
int i, src, dst, rank, num_procs;
block_size = extent * send_count;
if (not is_2dmesh(num_procs, &X, &Y))
- THROWF(arg_error,0, "allgather_2dmesh algorithm can't be used with this number of processes! ");
+ throw std::invalid_argument("allgather_2dmesh algorithm can't be used with this number of processes!");
my_row_base = (rank / Y) * Y;
my_col_base = rank % Y;
if (Y > X)
num_reqs = Y;
- req = (MPI_Request *) xbt_malloc(num_reqs * sizeof(MPI_Request));
+ MPI_Request* req = new MPI_Request[num_reqs];
+ MPI_Request* req_ptr = req;
- req_ptr = req;
-
- // do local allgather/local copy
+ // do local allgather/local copy
recv_offset = rank * block_size;
Datatype::copy(send_buff, send_count, send_type, (char *)recv_buff + recv_offset,
recv_count, recv_type);
Request::waitall(X - 1, req, MPI_STATUSES_IGNORE);
- free(req);
+ delete[] req;
return MPI_SUCCESS;
}