1、目錄一.系統安裝11.下載地址12.安裝步驟13.Visual Stdio設置1二.實驗程序21.簡單的MPI編程示例22.消息傳遞MPI編程示例133.消息傳遞MPI編程示例244.Monte Carlo方法計算圓周率65.計算積分8三.心得體會9一. 系統安裝1. 下載地址FTP匿名登陸，在pub/mpi/nt文件夾中2. 安裝步驟1) 在安裝有MPI的計算機上要建立一個有管理員權限的賬戶，不可以沒有密碼；2) 雙擊exe文件，按默認設置安裝MPI；3) 注冊MPI賬戶，調用MPIRegister.exe，用戶名和密碼即為第一步中的賬戶。3. Visual Stdio設置為避免每新建一個項

2、目都要設置一次，可以對它進行通用設置。打開視圖-其他窗口-屬性管理器，點擊Debug|Win32目錄下的Microsoft，在VC+目錄下的包含目錄中添加MPICH的Include路徑，庫目錄中添加MPICH的Lib路徑；在鏈接器-輸入目錄下的附加依賴項中添加mpich.lib、mpe.lib、mped.lib、mpichd.lib。二. 實驗程序1. 簡單的MPI編程示例1) 源代碼#include <mpi.h>#include <stdio.h>int main(int argc, char* argv)int num, rk;MPI_Init(&argc

3、, &argv);MPI_Comm_size(MPI_COMM_WORLD, &num);MPI_Comm_rank(MPI_COMM_WORLD, &rk);printf("Hello world from Process %d of %dn", rk, num);MPI_Finalize();return 0;2) 運行截圖2. 消息傳遞MPI編程示例11) 源代碼#include <mpi.h>#include <stdio.h>int main(int argc, char* argv)int myid, numpro

4、cs,source;MPI_Status status;char messages100;MPI_Init(&argc, &argv);MPI_Comm_size(MPI_COMM_WORLD, &numprocs);MPI_Comm_rank(MPI_COMM_WORLD, &myid);if (myid != 0)strcpy(messages, "Hello World!");MPI_Send(messages, strlen(messages) + 1, MPI_CHAR, 0, 99, MPI_COMM_WORLD);elsefor

5、 (source = 1; source < numprocs; source+)MPI_Recv(messages, 100, MPI_CHAR, source, 99, MPI_COMM_WORLD,&status);printf("I am process %d I recv string '%s' from process %dn", myid, messages, source);MPI_Finalize();return 0;2) 運行截圖3. 消息傳遞MPI編程示例21) 源代碼#include <stdio.h>#i

6、nclude <mpi.h> #include<math.h>#define N 1002int main(int argc, char* argv)int myid, numprocs, C = 0;double dataN, SqrtSum=0.0;_double d;for (int i = 0; i < N; i+)datai = i;MPI_Status status;MPI_Init(&argc, &argv);_MPI_Comm_rank(MPI_COMM_WORLD, &myid);MPI_Comm_size(MPI_COM

7、M_WORLD, &numprocs); -numprocs;if (myid = 0)for (int i = 0; i < N; +i) MPI_Send(&datai, 1, MPI_DOUBLE, i%numprocs + 1, 1, MPI_COMM_WORLD);for (int source = 1; source <= numprocs; +source)MPI_Recv(&d, 1, MPI_DOUBLE, source, 99, MPI_COMM_WORLD, &status);SqrtSum += d;elsefor (int

8、i = myid; i < N; i = i + numprocs)MPI_Recv(&d, 1, MPI_DOUBLE, 0, 1, MPI_COMM_WORLD, &status);SqrtSum += sqrt(d);+C;MPI_Send(&SqrtSum, 1, MPI_DOUBLE, 0, 99, MPI_COMM_WORLD); printf("I am process %d. I recv total %d from process 0, and SqrtSum=%f.n", myid, C, SqrtSum);_MPI_Fin

9、alize();_return 0;2) 運行截圖4. Monte Carlo方法計算圓周率1) 源代碼#include "mpi.h"#include <stdio.h>#include <stdlib.h>int main(int argc, char *argv)int myid, numprocs;int source;long count = 1000000;MPI_Status status;MPI_Init(&argc, &argv);MPI_Comm_rank(MPI_COMM_WORLD, &myid);MP

10、I_Comm_size(MPI_COMM_WORLD, &numprocs);srand(myid);double y, x, pi = 0.0, n = count;long m = 0, m1 = 0, i = 0, p = 0;for (i = 0; i<count; i+) x = (double)rand() / (double)RAND_MAX;y = (double)rand() / (double)RAND_MAX;if (x - 0.5)*(x - 0.5) + (y - 0.5)*(y - 0.5)<0.25)+m; pi = 4.0*m / n;pri

11、ntf("Process %d of % pi = %fn", myid, numprocs, pi);if (myid != 0) MPI_Send(&m, 1, MPI_DOUBLE, 0, 1, MPI_COMM_WORLD); else p = m;for (source = 1; source<numprocs; source+) MPI_Recv(&m1, 1, MPI_DOUBLE, source, 1, MPI_COMM_WORLD, &status); p = p + m1; printf("pi = %fn&quo

12、t;, 4.0*p / (n*numprocs); MPI_Finalize();return 0;2) 運行截圖5. 計算積分1) 源代碼#define N 100000000#define a 0#define b 10#include <stdio.h>#include <stdlib.h>#include <time.h>#include "mpi.h"int main(int argc, char* argv)int myid, numprocs; int i; _ double local = 0.0, dx = (doubl

13、e)b - (double)a) / (double)N;double inte, x;MPI_Init(&argc, &argv);_ MPI_Comm_rank(MPI_COMM_WORLD, &myid);_ MPI_Comm_size(MPI_COMM_WORLD, &numprocs);for (i = myid; i<N; i = i + numprocs) x = a + i*dx + dx / 2;local += x*x*dx;MPI_Reduce(&local, &inte, 1, MPI_DOUBLE, MPI_SUM, 0, MPI_COMM_WORLD);_ if (myid = 0) printf("The integal of x*x in region %d,%d =%16