/*
mpi规约函数
MPI_Reduce
时间:15.7.26
Jason Zhou
热爱你所写下的程序,他是你的伙伴,而不是工具.
*/
#include"mpi.h"
#include<stdio.h>
int main(int argc,char * argv[])
{
int np,rank;
MPI_Init(&argc,&argv);
MPI_Comm_rank(MPI_COMM_WORLD,&rank);
MPI_Comm_size(MPI_COMM_WORLD,&np);
printf("Hello world! rank:%d of np:%d\n",rank,np);
int a[6][6];
int b[6][6];
for(int i=0;i<6;i++)
for(int j=0;j<6;j++)
{
a[i][j]=rank;
b[i][j]=0;
}
MPI_Reduce(a,b,6,MPI_INT,MPI_SUM,0,MPI_COMM_WORLD);
if(0==rank)
{
for(int i=0;i<6;i++)
{
for(int j=0;j<6;j++)
{
fprintf(stderr,"%d ",b[i][j]);
}
fprintf(stderr,"\n");
}
}
MPI_Finalize();
return 1;
}
MPI_Reduce mpi规约函数
最新推荐文章于 2024-02-05 17:27:59 发布