PVM上模擬分散式系統

摘要:將一個10*10的矩陣，分給5台電腦去做 並算出時間，要算出分給5台 4台 3台 2台 1台的時間

 

心得:這是我第一個程式.....也花了我很多時間去寫，大概花了10幾個小時吧!真的很辛苦，過程中也很@#% 但寫出來後很有出就感^^

 

  1 /*
  2 (c) mysterywho 
  3 compiler : pvm
  4 date : 15/10/2008
  5 */
  6 
  7 //master端
  8 #include "pvm3.h"
  9 #include <stdio.h>
 10 #include <time.h>
 11 #include <sys/times.h>
 12 
 13 #define MAX 10
 14 
 15 
 16  main()
 17   {
 18 
 19       int mytid,number,task_ids[10];
 20       int i,j,NPROCS=4;
 21       int row,index,count,count1;
 22       int A[MAX][MAX],B[MAX][MAX];
 23       int result[MAX][MAX];
 24     
 25       int gettimeofday(struct timeval *tv);
 26 
 27 
 28       struct timeval tpstart,tpend;
 29       float  timeuse;
 30 
 31 
 32 
 33 
 34       // enroll in PVM //
 35          mytid = pvm_mytid();
 36 
 37      gettimeofday(&tpstart,NULL);
 38   
 39 
 40       // define A[][]
 41          printf("A matrix");
 42          printf("\n");
 43 
 44          for( i=0 ; i<10 ; i++)
 45           {
 46              for( j=0 ; j<10 ; j++)
 47                 {
 48                      A[i][j]= rand() %10;
 49 
 50                      printf("%2d  " ,A[i][j] ) ;
 51                 }
 52                    printf("\n");
 53          }
 54 
 55 
 56 
 57      // define B[][]
 58          printf("B matrix");
 59           printf("\n");
 60 
 61            for( i=0 ; i<10 ; i++)
 62              {
 63              for( j=0 ; j<10 ; j++)
 64                {
 65                        B[i][j]= rand() % 10;
 66                        printf("%2d  ",B[i][j]);
 67                }
 68 
 69                printf("\n");
 70 
 71             }
 72 
 73                printf("\n");
 74 
 75 
 76 
 77 
 78 
 79 
 80     // spawn worker tasks //
 81        number = pvm_spawn("worker_mystery",NULL,PvmTaskDefault,"", NPROCS, task_ids);
 82        printf("The number of spawn workers is %d\n",number);
 83 
 84        //send data to worker tasks //
 85 
 86          index = 0;
 87          row = 0;
 88           for (i=0;i<NPROCS;i++)
 89              {
 90                if(10% NPROCS > i)
 91                    {
 92                 row = 10 / NPROCS + 1;
 93                 count = row*10;
 94                    }
 95 
 96                 else
 97                   {
 98                   row = 10 / NPROCS;
 99                   count = row*10;
                   }
 
                   pvm_initsend(PvmDataDefault);
                   pvm_pkint(&row,1,1);
                   pvm_pkint(&index,1,1);
                   pvm_pkint(&count,1,1);
                   pvm_pkint(&A[index][0], count, 1);
                   pvm_pkint(&B[0][0], 100, 1);
                   pvm_send(task_ids[i], 4);
 
                   index = row + index ;
 
                   }
 
 
         // wait and gather results //
 
           int  msgtype=7;
 
           for (i = 0; i < NPROCS; i++)
             {
                pvm_recv(task_ids[i],msgtype);
                pvm_upkint(&index, 1, 1);
           //   pvm_upkint(&row, 1, 1);
                pvm_upkint(&count1, 1, 1);
           //   pvm_upkint(&A[index][0],count1,1);
                pvm_upkint(&result[index][0], count1, 1);
 
             }
 
         printf("C matrix");
         printf("\n");
 
         for(i=0;i< 10;i++)
           {
           for(j=0;j< 10;j++)
              {
                printf(" %6d  ", result[i][j]);
 
              }
            printf("\n");
 
           }
          printf("\n");
 
    gettimeofday(&tpend,NULL);
      timeuse=1000000*(tpend.tv_sec-tpstart.tv_sec)+  tpend.tv_usec-tpstart.tv_usec;
      timeuse/=1000000;
     
      printf("Used Time:%f\n",timeuse);
 
            printf("\n");
 
             pvm_exit();
                               }
 
 
 
 //worker端
 #include "pvm3.h"
 #include <stdio.h>
 #include <stdlib.h>
 
 #define MAX 10
 
 
  main()
  {
    int A[MAX][MAX],B[MAX][MAX],C[MAX][MAX];
    int mytid,master;
    int i,j,k,m=0,count,count1;
    int row,index;
 
 
    // enroll in PVM //
     mytid = pvm_mytid();
  
    // receive portion of array to be summed //
 
      pvm_recv(-1, -1);
      pvm_upkint(&row,1,1);
      pvm_upkint(&index,1,1);
      pvm_upkint(&count,1,1);
      pvm_upkint(&A[index][0], count, 1);
      pvm_upkint(&B[0][0], 100, 1);
 
      //compute //
 
     for(i=index;i<row+index;i++)
       {
        for(j=0;j<10;j++)
         {
           C[i][j]=0;
           for(k=0;k<10;k++)
             {
              C[i][j]=A[i][k]*B[k][j]+C[i][j];
             }
          }
        }
 
  
     count1 = count;
 
     // send computed sum back to master //
        master = pvm_parent();
 
        pvm_initsend(PvmDataRaw);
        pvm_pkint(&index, 1, 1);
        pvm_pkint(&count1, 1,1);
        pvm_pkint(&C[index][0],count1 ,1);
        pvm_send(master,7);
 
          pvm_exit();
          }