MPI_ludepOrder

for(i=0;i<m-1;i++)
MPI_Recv(&a(i,0),M,MPI_FLOAT,0,i+1,MPI_COMM_WORLD,&status);
else
for(i=0;i<m;i++)
MPI_Send(&A(i,0),M,MPI_FLOAT,i1,i2,MPI_COMM_WORLD);
}
}
}//my_rank = 0
else//my_rank != 0
{
if (left != 0)
if (my_rank < left)
if(left == 0)
v = j*m + i;
else if (j < left)
v = j * m + i;
else
v = left * m + (j - left)*mm + i;
if (my_rank>j)
for(k=0;k<mm;k++)
{
a(k,v)=a(k,v)/f[v];
for(w=v+1;w<M;w++)
a(k,w)=a(k,w)-f[w]*a(k,v);
}
}//m==1
}//left !=0
else//left == 0
{
i2 = 0;
for(i=m;i< M;i++)
{
i1=i/m;
i2=i2%m+1;
#define floatsize sizeof(float)
#define intsize sizeof(int)
int M,N;
int m;
int mm;
int mmm;
float *A;
int my_rank;
int p;
MPI_Status status;
int left=0;
MPI_Bcast(f,M,MPI_FLOAT,j,MPI_COMM_WORLD);
}
/*编号等于my_rank的进程利用主行对其第i+1,…,m-1行数据做行变换*/
if (my_rank==j)
for(k=i+1;k<mm;k++)
{
a(k,v)=a(k,v)/f[v];
for(w=v+1;w<M;w++)
a(k,w)=a(k,w)-f[w]*a(k,v);
}
/*编号大于my_rank的进程利用主行对其各行数据做行变换*/
{
for(i = m;i < M;i++)
{
i1=i;
i2=1;
MPI_Send(&A(i,0),M,MPI_FLOAT,i1,i2,MPI_COMM_WORLD);
if(left != 0)
{
i2 = 0;
for(i=m;i< left * m;i++)
{
i1=i/m;
i2=i2%m+1;
MPI_Send(&A(i,0),M,MPI_FLOAT,i1,i2,MPI_COMM_WORLD);
i2=i2%(m-1)+1;
MPI_Send(&A(i,0),M,MPI_FLOAT,i1,i2,MPI_COMM_WORLD);
}
}//m!=1
else//m==1
for(i=0;i<m;i++)
MPI_Recv(&a(i,0),M,MPI_FLOAT,0,i+1,MPI_COMM_WORLD,&status);
else
{
int i,j,k,my_rank,group_size;
int i1,i2;
int v,w;
float *a,*f,*l,*u;
FILE *fdA;
double time1,time2;
MPI_Init(&argc,&argv);
for (k=v;k<M;k++)
f[k]=a(i,k);
MPI_Bcast(f,M,MPI_FLOAT,my_rank,MPI_COMM_WORLD);
}
else
{
if(left == 0)
v = j*m + i;
else if (j < left)
v = j * m + i;
else
v = left * m + (j - left)*mm + i;
void fatal(char *message)
{
printf("%s\n",message);
exit(1);
}
void Environment_Finalize(float *a,float *f)
{
free(a);
free(f);
}
int main(int argc, char **argv)
if (my_rank==0)
{
l=(float*)malloc(floatsize*M*M);
u=(float*)malloc(floatsize*M*M);
}
/*0号进程采用顺序行划分将矩阵A划分为大小m*M的p块子矩阵，依次发送给1至p-1号进程*/
}
/*分配至各进程的子矩阵大小为m*M*/
a=(float*)malloc(floatsize*m*M);
/*各进程为主行元素建立发送和接收缓冲区*/
f=(float*)malloc(floatsize*M);
/*0号进程为l和u矩阵分配内存，以分离出经过变换后的A矩阵中的l和u矩阵*/
}
}
}
/*0号进程从其余各进程中接收子矩阵a，得到经过变换的矩阵A*/
if (my_rank==0)
{
for(i=0;i<m;i++)
for(j=0;j<M;j++)
A(i,j)=a(i,j);
}
if (my_rank!=0)
{
if (left == 0)
mm = m;
else if (my_rank < left)
mm = m;
else if(j < left)
mm = m;
else
mm = m - 1;
for(i=0; i<mm; i++)
{
/*j号进程负责广播主行元素*/
if (my_rank==j)
{
//v=i*p+j;
exit(0);
}
A=(float *)malloc(floatsize*M*M);
for(i = 0; i < M; i ++)
for(j = 0; j < M; j ++)
fscanf(fdA, "%f", A+i*M+j);
fclose(fdA);
}
/*0号进程将M广播给所有进程*/
MPI_Bcast(&M,1,MPI_INT,0,MPI_COMM_WORLD);
m=M/p;
if (M%p!=0) Leabharlann m++;
left = M%p;
}
if (m != 1)
{
i2 = 0;
for(i = left * m; i < M; i++)
{
i1=left + (I – left*m) / (m - 1);
if (a==NULL) fatal("allocate error\n");
if (my_rank==0)
{
for(i=0;i<m;i++)
for(j=0;j<M;j++)
a(i,j)=A(i,j);
#include "stdio.h"
#include "stdlib.h"
#include "mpi.h"
#define a(x,y) a[x*M+y]
/*A为M*M矩阵*/
#define A(x,y) A[x*M+y]
#define l(x,y) l[x*M+y]
#define u(x,y) u[x*M+y]
MPI_Comm_size(MPI_COMM_WORLD,&group_size);
MPI_Comm_rank(MPI_COMM_WORLD,&my_rank);
p=group_size;
time1=MPI_Wtime();
if (my_rank==0)
MPI_Recv(&a(i,0),M,MPI_FLOAT,0,i+1,MPI_COMM_WORLD,&status);
}
for(j=0;j<p;j++){
if (left == 0)