matrix3

最新推荐文章于 2024-11-03 21:41:14 发布

原创最新推荐文章于 2024-11-03 21:41:14 发布 · 495 阅读

0 ·

CC 4.0 BY-SA版权

C 同时被 2 个专栏收录

16 篇文章

订阅专栏

MPI与矩阵运算

14 篇文章

订阅专栏

本文展示了一个使用MPI实现的简单并行通信程序示例，该程序涉及两个进程之间的数据发送与接收。通过定义特定的数据类型，实现了矩阵数据的有效传输。

//Writed by YaoPu, Oct 10,2016
//if any problem,send email yaopu@iccas.ac.cn
//only for 2 processors !!!!!
#include"mpi.h"
#include<stdio.h>
int main(int argc,char **argv)
{
int rows=6,cols=5;
int tag=0;
MPI_Init(&argc,&argv);
int rank;
MPI_Comm_rank(MPI_COMM_WORLD,&rank);
int A[rows][cols];
int B[cols][rows];
MPI_Datatype stype, ntype;
MPI_Type_vector(rows,1,cols,MPI_INT,&stype);
MPI_Type_hvector(cols,1,sizeof(int),stype,&ntype);
MPI_Type_commit(&ntype);
MPI_Status sta;
if (rank==0){

for(int i=0;i<rows;i++){
for(int j=0;j<cols;j++){

A[i][j]=i*cols+j;
printf("%d\t",A[i][j]);
}
printf("\n");
}
printf("\n");
MPI_Send(A,1,ntype,1,tag,MPI_COMM_WORLD);
}
else{
MPI_Recv(B,rows*cols,MPI_INT,0,tag,MPI_COMM_WORLD,&sta);
for(int i=0;i<cols;i++){
for(int j=0;j<rows;j++)
printf("%d\t",B[i][j]);
printf("\n");
}
}
MPI_Type_free(&ntype);
MPI_Finalize();
}