Changes between Version 9 and Version 10 of mpich/test/cpi

Timestamp:

Apr 15, 2008, 5:13:32 PM (18 years ago)

Author:

wade

Comment:

--

mpich/test/cpi

@@ -4 +4 @@
 #include <stdio.h>
 #include <math.h>
+#include <stdlib.h>
 
 double f( double );
@@ -13 +14 @@
 int main( int argc, char *argv[])
 {
-    int done = 0, n, myid, numprocs, i, count=0;
+    int done = 0, n, myid, numprocs, i=0, count=0;
     double PI25DT = 3.141592653589793238462643;
     double mypi, pi, h, sum, x;
@@ -28 +29 @@
             myid, processor_name);
 
-    n = 0;
+    n = 0; 
     while (!done)
     {
-
-/* 註 1: 我在這邊開始計算 loop 執行的次數 */
         count++;
-
+        printf("Node %d : loop runs  %d  times\n", myid, count);
         if (myid == 0)
         {
@@ -42 +41 @@
 */
             if (n==0) n=100; else n=0;
-
             startwtime = MPI_Wtime();
         }
+        printf("Node %d :  n = %d before MPI_Bcast\n", myid, n);
         MPI_Bcast(&n, 1, MPI_INT, 0, MPI_COMM_WORLD);
+        printf("Node %d :  n = %d before MPI_Bcast\n", myid, n);
         if (n == 0)
             done = 1;
@@ -61 +61 @@
             MPI_Reduce(&mypi, &pi, 1, MPI_DOUBLE, MPI_SUM, 0, MPI_COMM_WORLD);
 
-/* 註 2: 在這邊我把每個 node 所執行的 loop 次數印出 */
             if (myid == 0)
             {
@@ -69 +68 @@
                 printf("wall clock time = %f\n",
                        endwtime-startwtime);           
-                printf("Node 0 loop runs  %d  times\n", count);
-            } else
-            {
-                printf("Node not 0 loop runs  %d  times\n", count);
             }
-
         }
     }
     MPI_Finalize();
+
     return 0;
 }
 }}}
 == 結果 ==
-就算去掉 while loop 也能執行。
-[[Image(cpi-02.png)]]
+去掉 while loop 也能執行，是因為 while loop 的目的為展現出，每個 Node 在第二次進入 while loop 時，Node 0 會將 n = 0 藉由 MPI_Bcast 廣播出去。由下圖可發現，每個 Node 只會執行二次，因為在 MPI_Bcast 中，每個 Node 皆會被 lock 住，直到 Node 0 將 n = 0 送到其它的 Node ，每個 Node 才會接著往下做，不過 lock 的機制，似乎與「群」有關，有些 function 如 sprintf，本身就是一個群，所以可以看見前面所有的 Node 會全部將結果顯示在螢幕上才繼續做下一「群」的程式，而我們也可以發現，在下一「群」的程式中，也是將這「群」的程式執行完後才將結果顯示出。 
+[[Image(cpi_02.png)]]