矩阵求逆的并行算法.doc

《矩阵求逆的并行算法.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、.矩阵求逆的并行算法：#include"stdio.h"#include"stdlib.h"#include"mpi.h"#defineintsizesizeof(int)#definefloatsizesizeof(float)#defineA(x,y)A[x*N+y]#defineQ(x,y)Q[x*N+y]#definea(x,y)a[x*N+y]#definef(x)f[x]float*A,*Q;float*a,*f;intN,v,m;intp;intmyid;MPI_Statusstatus;FILE*

2、dataFile;doublestarttime,endtime,time1;voidreadData(){inti,j;starttime=MPI_Wtime();dataFile=fopen("dataIn.txt","r");fscanf(dataFile,"%d",&N);A=(float*)malloc(floatsize*N*N);for(i=0;i

3、ntf("Inputoffile"dataIn.txt"");printf("%d",N);for(i=0;i

4、Q(i,j));}printf("");}endtime=MPI_Wtime();printf("");printf("Wholerunningtime=%fseconds",endtime-starttime);printf("Distributedatatime=%fseconds",time1-starttime);printf("Parallelcomputetime=%fseconds",endtime-time1);dataFile=fopen("dataOut.txt","w")

5、;fprintf(dataFile,"Inputoffile"dataIn.txt"");fprintf(dataFile,"%d",N);for(i=0;i

6、(dataFile,"%ft",Q(i,j));}fprintf(dataFile,"");}/*0号进程将时间统计写入目标文件*/fprintf(dataFile,"");fprintf(dataFile,"Wholerunningtime=%fseconds",endtime-starttime);fprintf(dataFile,"Distributedatatime=%fseconds",time1-starttime);fprintf(dataFile,"Parallelcompute

7、time=%fseconds",endtime-time1);fprintf(dataFile,"ParallelProcessnumber=%d",p);fclose(dataFile);}voidbroadcast(inti,intj,intv){intk;if(myid==j){a(i,v)=(float)(1/a(i,v));for(k=0;k

8、t(&f(0),N,MPI_FLOAT,myid,MPI_COMM_WORLD);}else{MPI_Bcast(f,N,MPI_FLOAT,j,MPI_COMM_WORLD);}}voidtransform(inti,intj,intv){intk,w;if(myid!=j){for(k=0;k