我将来自每个进程的一组数组元素加在一起:
double rho[1024];
//Some operation to calculate rho for each process;
MPI_Allreduce(rho,rho,1024,MPI_DOUBLE,MPI_SUM,MPI_COMM_WORLD);
将 rho 作为 sendbuf 和 recvbuf 工作吗?
你检查过MPI_IN_PLACE
吗?根据MPI_AllReduce 手册页和MPI 文档,只要您在同一个组中工作,它就可以用于为 sendbuf 和 recvbuf 指定相同的缓冲区。
调用看起来像:
MPI_Allreduce(MPI_IN_PLACE,rho,1024,MPI_DOUBLE,MPI_SUM,MPI_COMM_WORLD);