-
-/* Copyright (c) 2007-2015. The SimGrid Team.
- * All rights reserved. */
+/* Copyright (c) 2007-2016. 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 <simgrid/s4u/host.hpp>
#include <xbt/ex.hpp>
#include "private.h"
return retval;
}
+int PMPI_Type_size_x(MPI_Datatype datatype, MPI_Count *size)
+{
+ int retval = 0;
+
+ if (datatype == MPI_DATATYPE_NULL) {
+ retval = MPI_ERR_TYPE;
+ } else if (size == nullptr) {
+ retval = MPI_ERR_ARG;
+ } else {
+ *size = static_cast<MPI_Count>(smpi_datatype_size(datatype));
+ retval = MPI_SUCCESS;
+ }
+ return retval;
+}
+
int PMPI_Type_get_extent(MPI_Datatype datatype, MPI_Aint * lb, MPI_Aint * extent)
{
int retval = 0;
} else {
if (n == 0) {
*newgroup = group;
- if(group!= smpi_comm_group(MPI_COMM_WORLD) && group != MPI_GROUP_NULL
+ if (group != smpi_comm_group(MPI_COMM_WORLD)
&& group != smpi_comm_group(MPI_COMM_SELF) && group != MPI_GROUP_EMPTY)
smpi_group_use(group);
} else if (n == smpi_group_size(group)) {
} else {
if (n == 0) {
*newgroup = group;
- if(group!= smpi_comm_group(MPI_COMM_WORLD) && group != MPI_GROUP_NULL
+ if (group!= smpi_comm_group(MPI_COMM_WORLD)
&& group != smpi_comm_group(MPI_COMM_SELF) && group != MPI_GROUP_EMPTY)
smpi_group_use(group);
} else {
return retval;
}
+int PMPI_Comm_create_group(MPI_Comm comm, MPI_Group group, int, MPI_Comm* comm_out)
+{
+ int retval = 0;
+ smpi_bench_end();
+
+ if (comm_out == nullptr) {
+ retval = MPI_ERR_ARG;
+ } else if (comm == MPI_COMM_NULL) {
+ retval = MPI_ERR_COMM;
+ } else {
+ retval = MPI_Comm_create(comm, group, comm_out);
+ }
+ smpi_bench_begin();
+
+ return retval;
+}
+
int PMPI_Send_init(void *buf, int count, MPI_Datatype datatype, int dst, int tag, MPI_Comm comm, MPI_Request * request)
{
int retval = 0;
smpi_empty_status(status);
retval = MPI_SUCCESS;
} else {
- int rank = (request!=nullptr && (*request)->comm != MPI_COMM_NULL) ? smpi_process_index() : -1;
+ int rank = ((*request)->comm != MPI_COMM_NULL) ? smpi_process_index() : -1;
instr_extra_data extra = xbt_new0(s_instr_extra_data_t,1);
extra->type = TRACING_TEST;
smpi_bench_end();
//save requests information for tracing
int i;
- int *srcs = nullptr, *dsts = nullptr, *recvs = nullptr, *tags = nullptr;
- MPI_Comm* comms = nullptr;
+ typedef struct {
+ int src;
+ int dst;
+ int recv;
+ int tag;
+ MPI_Comm comm;
+ } savedvalstype;
+ savedvalstype* savedvals=nullptr;
if(count>0){
- srcs = xbt_new0(int, count);
- dsts = xbt_new0(int, count);
- recvs = xbt_new0(int, count);
- tags = xbt_new0(int, count);
- comms = xbt_new0(MPI_Comm, count);
+ savedvals = xbt_new0(savedvalstype, count);
}
for (i = 0; i < count; i++) {
MPI_Request req = requests[i]; //already received requests are no longer valid
if (req) {
- srcs[i] = req->src;
- dsts[i] = req->dst;
- recvs[i] = req->recv;
- tags[i] = req->tag;
- comms[i] = req->comm;
+ savedvals[i]=(savedvalstype){req->src, req->dst, req->recv, req->tag, req->comm};
}
}
int rank_traced = smpi_process_index();
*index = smpi_mpi_waitany(count, requests, status);
if(*index!=MPI_UNDEFINED){
- int src_traced = srcs[*index];
+ int src_traced = savedvals[*index].src;
//the src may not have been known at the beginning of the recv (MPI_ANY_SOURCE)
- int dst_traced = dsts[*index];
- int is_wait_for_receive = recvs[*index];
+ int dst_traced = savedvals[*index].dst;
+ int is_wait_for_receive = savedvals[*index].recv;
if (is_wait_for_receive) {
- if(srcs[*index]==MPI_ANY_SOURCE)
+ if(savedvals[*index].src==MPI_ANY_SOURCE)
src_traced = (status!=MPI_STATUSES_IGNORE) ?
- smpi_group_rank(smpi_comm_group(comms[*index]), status->MPI_SOURCE) : srcs[*index];
- TRACE_smpi_recv(rank_traced, src_traced, dst_traced, tags[*index]);
+ smpi_group_rank(smpi_comm_group(savedvals[*index].comm), status->MPI_SOURCE) : savedvals[*index].src;
+ TRACE_smpi_recv(rank_traced, src_traced, dst_traced, savedvals[*index].tag);
}
TRACE_smpi_ptp_out(rank_traced, src_traced, dst_traced, __FUNCTION__);
}
- xbt_free(srcs);
- xbt_free(dsts);
- xbt_free(recvs);
- xbt_free(tags);
- xbt_free(comms);
-
+ xbt_free(savedvals);
smpi_bench_begin();
return MPI_SUCCESS;
smpi_bench_end();
//save information from requests
int i;
- int *srcs = xbt_new0(int, count);
- int *dsts = xbt_new0(int, count);
- int *recvs = xbt_new0(int, count);
- int *tags = xbt_new0(int, count);
- int *valid = xbt_new0(int, count);
- MPI_Comm *comms = xbt_new0(MPI_Comm, count);
+ typedef struct {
+ int src;
+ int dst;
+ int recv;
+ int tag;
+ int valid;
+ MPI_Comm comm;
+ } savedvalstype;
+ savedvalstype* savedvals=xbt_new0(savedvalstype, count);
for (i = 0; i < count; i++) {
MPI_Request req = requests[i];
if(req!=MPI_REQUEST_NULL){
- srcs[i] = req->src;
- dsts[i] = req->dst;
- recvs[i] = req->recv;
- tags[i] = req->tag;
- comms[i] = req->comm;
- valid[i]=1;;
+ savedvals[i]=(savedvalstype){req->src, req->dst, req->recv, req->tag, 1, req->comm};
}else{
- valid[i]=0;
+ savedvals[i].valid=0;
}
}
int rank_traced = smpi_process_index();
int retval = smpi_mpi_waitall(count, requests, status);
for (i = 0; i < count; i++) {
- if(valid[i]){
+ if(savedvals[i].valid){
//the src may not have been known at the beginning of the recv (MPI_ANY_SOURCE)
- int src_traced = srcs[i];
- int dst_traced = dsts[i];
- int is_wait_for_receive = recvs[i];
+ int src_traced = savedvals[i].src;
+ int dst_traced = savedvals[i].dst;
+ int is_wait_for_receive = savedvals[i].recv;
if (is_wait_for_receive) {
if(src_traced==MPI_ANY_SOURCE)
src_traced = (status!=MPI_STATUSES_IGNORE) ?
- smpi_group_rank(smpi_comm_group(comms[i]), status[i].MPI_SOURCE) : srcs[i];
- TRACE_smpi_recv(rank_traced, src_traced, dst_traced,tags[i]);
+ smpi_group_rank(smpi_comm_group(savedvals[i].comm), status[i].MPI_SOURCE) : savedvals[i].src;
+ TRACE_smpi_recv(rank_traced, src_traced, dst_traced,savedvals[i].tag);
}
}
}
TRACE_smpi_ptp_out(rank_traced, -1, -1, __FUNCTION__);
- xbt_free(srcs);
- xbt_free(dsts);
- xbt_free(recvs);
- xbt_free(valid);
- xbt_free(tags);
- xbt_free(comms);
+ xbt_free(savedvals);
smpi_bench_begin();
return retval;
if (comm == MPI_COMM_NULL) {
retval = MPI_ERR_COMM;
} else {
- int rank = comm != MPI_COMM_NULL ? smpi_process_index() : -1;
- instr_extra_data extra = xbt_new0(s_instr_extra_data_t,1);
- extra->type = TRACING_BARRIER;
- TRACE_smpi_collective_in(rank, -1, __FUNCTION__, extra);
+ int rank = comm != MPI_COMM_NULL ? smpi_process_index() : -1;
+ instr_extra_data extra = xbt_new0(s_instr_extra_data_t, 1);
+ extra->type = TRACING_BARRIER;
+ TRACE_smpi_collective_in(rank, -1, __FUNCTION__, extra);
- mpi_coll_barrier_fun(comm);
- retval = MPI_SUCCESS;
+ mpi_coll_barrier_fun(comm);
+ retval = MPI_SUCCESS;
- TRACE_smpi_collective_out(rank, -1, __FUNCTION__);
+ TRACE_smpi_collective_out(rank, -1, __FUNCTION__);
}
smpi_bench_begin();
if(known==0)
dt_size_send = smpi_datatype_size(datatype);
extra->send_size = count*dt_size_send;
+ void * sendtmpbuf = sendbuf;
+ if(sendbuf==MPI_IN_PLACE){
+ sendtmpbuf= static_cast<void*>(xbt_malloc(count*smpi_datatype_size(datatype)));
+ memcpy(sendtmpbuf,recvbuf, count*smpi_datatype_size(datatype));
+ }
TRACE_smpi_collective_in(rank, -1, __FUNCTION__,extra);
- smpi_mpi_exscan(sendbuf, recvbuf, count, datatype, op, comm);
+ smpi_mpi_exscan(sendtmpbuf, recvbuf, count, datatype, op, comm);
retval = MPI_SUCCESS;
TRACE_smpi_collective_out(rank, -1, __FUNCTION__);
+ if(sendbuf==MPI_IN_PLACE)
+ xbt_free(sendtmpbuf);
}
smpi_bench_begin();
dt_size_send = smpi_datatype_size(datatype);
extra->send_size = 0;
extra->recvcounts= xbt_new(int, size);
- for(i=0; i< size; i++)//copy data to avoid bad free
+ int totalcount = 0;
+ for(i=0; i< size; i++){//copy data to avoid bad free
extra->recvcounts[i] = recvcounts[i]*dt_size_send;
- TRACE_smpi_collective_in(rank, -1, __FUNCTION__,extra);
-
+ totalcount+= recvcounts[i];
+ }
void* sendtmpbuf=sendbuf;
- if(sendbuf==MPI_IN_PLACE)
- sendtmpbuf=recvbuf;
+ if(sendbuf==MPI_IN_PLACE){
+ sendtmpbuf= static_cast<void*>(xbt_malloc(totalcount*smpi_datatype_size(datatype)));
+ memcpy(sendtmpbuf,recvbuf, totalcount*smpi_datatype_size(datatype));
+ }
- mpi_coll_reduce_scatter_fun(sendtmpbuf, recvbuf, recvcounts, datatype, op, comm);
- retval = MPI_SUCCESS;
+ TRACE_smpi_collective_in(rank, -1, __FUNCTION__,extra);
+
+ mpi_coll_reduce_scatter_fun(sendtmpbuf, recvbuf, recvcounts, datatype, op, comm);
+ retval = MPI_SUCCESS;
TRACE_smpi_collective_out(rank, -1, __FUNCTION__);
+
+ if(sendbuf==MPI_IN_PLACE)
+ xbt_free(sendtmpbuf);
}
smpi_bench_begin();
} else if (recvcount < 0) {
retval = MPI_ERR_ARG;
} else {
- int count=smpi_comm_size(comm);
+ int count=smpi_comm_size(comm);
int rank = comm != MPI_COMM_NULL ? smpi_process_index() : -1;
instr_extra_data extra = xbt_new0(s_instr_extra_data_t,1);
extra->recvcounts= xbt_new(int, count);
for(i=0; i< count; i++)//copy data to avoid bad free
extra->recvcounts[i] = recvcount*dt_size_send;
+ void* sendtmpbuf=sendbuf;
+ if(sendbuf==MPI_IN_PLACE){
+ sendtmpbuf= static_cast<void*>(xbt_malloc(recvcount*count*smpi_datatype_size(datatype)));
+ memcpy(sendtmpbuf,recvbuf, recvcount*count*smpi_datatype_size(datatype));
+ }
TRACE_smpi_collective_in(rank, -1, __FUNCTION__,extra);
- int* recvcounts=static_cast<int*>(xbt_malloc(count));
- for (i=0; i<count;i++)
- recvcounts[i]=recvcount;
- mpi_coll_reduce_scatter_fun(sendbuf, recvbuf, recvcounts, datatype, op, comm);
- xbt_free(recvcounts);
- retval = MPI_SUCCESS;
+ int* recvcounts=static_cast<int*>(xbt_malloc(count*sizeof(int)));
+ for (i=0; i<count;i++)
+ recvcounts[i]=recvcount;
+ mpi_coll_reduce_scatter_fun(sendtmpbuf, recvbuf, recvcounts, datatype, op, comm);
+ xbt_free(recvcounts);
+ retval = MPI_SUCCESS;
TRACE_smpi_collective_out(rank, -1, __FUNCTION__);
+
+ if(sendbuf==MPI_IN_PLACE)
+ xbt_free(sendtmpbuf);
}
smpi_bench_begin();
if (comm == MPI_COMM_NULL) {
retval = MPI_ERR_COMM;
- } else if (sendtype == MPI_DATATYPE_NULL
+ } else if ((sendbuf != MPI_IN_PLACE && sendtype == MPI_DATATYPE_NULL)
|| recvtype == MPI_DATATYPE_NULL) {
retval = MPI_ERR_TYPE;
} else {
int rank = comm != MPI_COMM_NULL ? smpi_process_index() : -1;
instr_extra_data extra = xbt_new0(s_instr_extra_data_t,1);
extra->type = TRACING_ALLTOALL;
+
+ void* sendtmpbuf = static_cast<char*>(sendbuf);
+ int sendtmpcount = sendcount;
+ MPI_Datatype sendtmptype = sendtype;
+ if( sendbuf == MPI_IN_PLACE ) {
+ sendtmpbuf = static_cast<void*>(xbt_malloc(recvcount*smpi_comm_size(comm)*smpi_datatype_size(recvtype)));
+ memcpy(sendtmpbuf,recvbuf, recvcount*smpi_comm_size(comm)*smpi_datatype_size(recvtype));
+ sendtmpcount = recvcount;
+ sendtmptype = recvtype;
+ }
+
int known=0;
- extra->datatype1 = encode_datatype(sendtype, &known);
+ extra->datatype1 = encode_datatype(sendtmptype, &known);
if(known==0)
- extra->send_size = sendcount*smpi_datatype_size(sendtype);
+ extra->send_size = sendtmpcount*smpi_datatype_size(sendtmptype);
else
- extra->send_size = sendcount;
+ extra->send_size = sendtmpcount;
extra->datatype2 = encode_datatype(recvtype, &known);
if(known==0)
extra->recv_size = recvcount*smpi_datatype_size(recvtype);
else
extra->recv_size = recvcount;
+
TRACE_smpi_collective_in(rank, -1, __FUNCTION__,extra);
- retval = mpi_coll_alltoall_fun(sendbuf, sendcount, sendtype, recvbuf, recvcount, recvtype, comm);
+ retval = mpi_coll_alltoall_fun(sendtmpbuf, sendtmpcount, sendtmptype, recvbuf, recvcount, recvtype, comm);
TRACE_smpi_collective_out(rank, -1, __FUNCTION__);
+
+ if( sendbuf == MPI_IN_PLACE )
+ xbt_free(sendtmpbuf);
}
smpi_bench_begin();
retval = MPI_ERR_COMM;
} else if (sendtype == MPI_DATATYPE_NULL || recvtype == MPI_DATATYPE_NULL) {
retval = MPI_ERR_TYPE;
- } else if (sendcounts == nullptr || senddisps == nullptr || recvcounts == nullptr || recvdisps == nullptr) {
+ } else if ((sendbuf!= MPI_IN_PLACE && (sendcounts == nullptr || senddisps == nullptr)) || recvcounts == nullptr || recvdisps == nullptr) {
retval = MPI_ERR_ARG;
} else {
int rank = comm != MPI_COMM_NULL ? smpi_process_index() : -1;
extra->recvcounts= xbt_new(int, size);
extra->sendcounts= xbt_new(int, size);
int known=0;
- extra->datatype1 = encode_datatype(sendtype, &known);
- int dt_size_send = 1;
- if(known==0)
- dt_size_send = smpi_datatype_size(sendtype);
int dt_size_recv = 1;
extra->datatype2 = encode_datatype(recvtype, &known);
- if(known==0)
- dt_size_recv = smpi_datatype_size(recvtype);
+ dt_size_recv = smpi_datatype_size(recvtype);
+
+ void* sendtmpbuf = static_cast<char*>(sendbuf);
+ int * sendtmpcounts = sendcounts;
+ int *sendtmpdisps = senddisps;
+ MPI_Datatype sendtmptype = sendtype;
+ int maxsize=0;
for(i=0; i< size; i++){//copy data to avoid bad free
- extra->send_size += sendcounts[i]*dt_size_send;
extra->recv_size += recvcounts[i]*dt_size_recv;
-
- extra->sendcounts[i] = sendcounts[i]*dt_size_send;
extra->recvcounts[i] = recvcounts[i]*dt_size_recv;
+ if (((recvdisps[i]+recvcounts[i])*dt_size_recv) > maxsize)
+ maxsize=(recvdisps[i]+recvcounts[i])*dt_size_recv;
+ }
+
+ if( sendbuf == MPI_IN_PLACE ) {
+ sendtmpbuf = static_cast<void*>(xbt_malloc(maxsize));
+ memcpy(sendtmpbuf,recvbuf, maxsize);
+ sendtmpcounts= static_cast<int*>(xbt_malloc(size*sizeof(int)));
+ memcpy(sendtmpcounts,recvcounts, size*sizeof(int));
+ sendtmpdisps= static_cast<int*>(xbt_malloc(size*sizeof(int)));
+ memcpy(sendtmpdisps,recvdisps, size*sizeof(int));
+ sendtmptype=recvtype;
+ }
+
+ extra->datatype1 = encode_datatype(sendtmptype, &known);
+ int dt_size_send = 1;
+ dt_size_send = smpi_datatype_size(sendtmptype);
+
+ for(i=0; i< size; i++){//copy data to avoid bad free
+ extra->send_size += sendtmpcounts[i]*dt_size_send;
+ extra->sendcounts[i] = sendtmpcounts[i]*dt_size_send;
}
extra->num_processes = size;
TRACE_smpi_collective_in(rank, -1, __FUNCTION__,extra);
-
- retval = mpi_coll_alltoallv_fun(sendbuf, sendcounts, senddisps, sendtype, recvbuf, recvcounts, recvdisps, recvtype,
+ retval = mpi_coll_alltoallv_fun(sendtmpbuf, sendtmpcounts, sendtmpdisps, sendtmptype, recvbuf, recvcounts, recvdisps, recvtype,
comm);
TRACE_smpi_collective_out(rank, -1, __FUNCTION__);
+
+ if( sendbuf == MPI_IN_PLACE ) {
+ xbt_free(sendtmpbuf);
+ xbt_free(sendtmpcounts);
+ xbt_free(sendtmpdisps);
+ }
}
smpi_bench_begin();
{
int retval = MPI_SUCCESS;
- strncpy(name, sg_host_get_name(SIMIX_host_self()),
- strlen(sg_host_get_name(SIMIX_host_self())) < MPI_MAX_PROCESSOR_NAME - 1 ?
- strlen(sg_host_get_name(SIMIX_host_self())) +1 : MPI_MAX_PROCESSOR_NAME - 1 );
+ strncpy(name, SIMIX_host_self()->cname(), strlen(SIMIX_host_self()->cname()) < MPI_MAX_PROCESSOR_NAME - 1
+ ? strlen(SIMIX_host_self()->cname()) + 1
+ : MPI_MAX_PROCESSOR_NAME - 1);
*resultlen = strlen(name) > MPI_MAX_PROCESSOR_NAME ? MPI_MAX_PROCESSOR_NAME : strlen(name);
return retval;
} else if (count<0){
retval = MPI_ERR_COUNT;
} else {
- int* blocklens=static_cast<int*>(xbt_malloc(blocklength*count));
+ int* blocklens=static_cast<int*>(xbt_malloc(blocklength*count*sizeof(int)));
for (i=0; i<count;i++)
blocklens[i]=blocklength;
retval = smpi_datatype_indexed(count, blocklens, indices, old_type, new_type);
} else if (count<0){
retval = MPI_ERR_COUNT;
} else {
- int* blocklens=(int*)xbt_malloc(blocklength*count);
+ int* blocklens=(int*)xbt_malloc(blocklength*count*sizeof(int));
for (i=0; i<count;i++)blocklens[i]=blocklength;
retval = smpi_datatype_hindexed(count, blocklens, indices, old_type, new_type);
xbt_free(blocklens);