编写一个MPI并行程序,实现矩阵A和矩阵B相乘,结果放在矩阵C中。
要求:
(1)程序不能固定或假定矩阵的行数和列数,也不能事先假定有多少个进程参与运算;
(2)矩阵的行数与列数自行定义,矩阵A和B数值自行定义,计算结果C存入硬盘文件;
(3)请节约各进程的内存,不要随意开辟内存空间,进程0首先读入进程1所需的数据并发送,然后读入进程2所需的数据并发送,依次类推,最后读入本进程所需的数据,个进程接收、计算、并将计算结果发送给进程0。
#include<iostream>
#include<stdio.h>
#include<mpi.h>
#include<math.h>
#include<stdlib.h>
void multi(double *a,int x,int y) //矩阵数组
{
int i;
for(i=0;i<x*y;i++)
a[i]=3+i;
}
void plus(double *a,double *b,double *c,int x,int y,int z) //计算
{
int i,j,k;
double s;
for(i=0;i<x;i++)
{