mpi-program例子c++.docVIP

/**************************************************************** File name: mpipi.c Description: Caculate PI, which uses parallel algorithm. Author: lijianjun(ict) Time: 2007-6-5 *****************************************************************/ #include mpi.h #include stdio.h #include math.h int main(int argc, char **argv) { int myid, np; int n, i, precision, is, ie; double pi, h, eps, xi, s; MPI_Comm mycomm; MPI_Init(argc, argv); MPI_Comm_dup(MPI_COMM_WORLD,mycomm); MPI_Comm_rank(mycomm, myid); MPI_Comm_size(mycomm, np); printf(\nThe process %d of %d is running.\n, myid, np); if(myid == 0) { printf(Number of digits(1-16)(precision):); scanf(%d, precision); } MPI_Bcast(precision, 1, MPI_INT, 0, mycomm); eps = 1; for(i=1; i=precision; i++) eps = eps * 0.1; n = sqrt(4.0/(3.0*eps)); h = 1.0/n; precision = n/np; if(myid == 0) { s = 2.0; xi = 0; is = 0; ie = precision; } else if(myid == np-1) { s = 1.0; is = myid*precision; ie = n-1; xi = is*h; } else { s = 0.0; is = myid*precision; ie = is + precision; xi = is*h; } for(i=is+1; i=ie; i++) { xi = xi + h; s = s + 4.0/(1.0+xi*xi); } MPI_Reduce(s, pi, 1, MPI_DOUBLE, MPI_SUM, 0, mycomm); if(myid == 0) { pi = h*pi; printf(The pi = %16.15f\n, pi); } MPI_Comm_free(mycomm); MPI_Finalize(); return 1; } /**************************************************************** File name: mpicannon.c Description: Caculate C = A × B, which uses cannon algorithm. Author: lijianjun(ict) Time: 2007-6-5 *****************************************************************/ #include stdio.h #include math.h #include mpi.h #define lda 50 #define ldb 50 #define ldc 50 #define ldt 50 #define maxnp 4 void mat_init_a(int npsr, double a[lda][lda], int m, int ma, int k, int ka, int rowid, int colid) { int i, j, mt, mr, ist, jst, kt, kr, offside; mt = m / npsr; mr = m % npsr; if(rowid mr) ist = (mt + 1) * rowid; else ist = mt * rowid + mr; kt = k / npsr; kr = k % npsr; if