安裝OpenMPI來(lái)配合C語(yǔ)言程序進(jìn)行并行計(jì)算

來(lái)源：本站原創(chuàng)|時(shí)間：2020-01-10|欄目：C語(yǔ)言|點(diǎn)擊：次

安裝OPENMPI
由于是實(shí)驗(yàn)，也不進(jìn)行多機(jī)的配置了，只在虛擬機(jī)里安裝吧。多個(gè)機(jī)器的配置可以參考此文
最簡(jiǎn)單的方法，apt安裝

sudo apt-get install libcr-dev mpich2 mpich2-doc

測(cè)試
hello.c

/* C Example */
#include <mpi.h>
#include <stdio.h>

int main (int argc, char* argv[])
{
 int rank, size;

 MPI_Init (&argc, &argv);   /* starts MPI */
 MPI_Comm_rank (MPI_COMM_WORLD, &rank);    /* get current process id */
 MPI_Comm_size (MPI_COMM_WORLD, &size);    /* get number of processes */
 printf( "Hello world from process %d of %d\n", rank, size );
 MPI_Finalize();
 return 0;
}

編譯運(yùn)行及顯示結(jié)果

mpicc mpi_hello.c -o hello
mpirun -np 2 ./hello
Hello world from process 0 of 2
Hello world from process 1 of 2

正常出現(xiàn)結(jié)果表明沒(méi)有問(wèn)題，
看下openmpi的版本

mpirun --version

mpirun (Open MPI) 1.6.5
Report bugs to http://www.open-mpi.org/community/help/

MPI計(jì)算矩陣乘法

通過(guò)opemMPI加速矩陣乘法運(yùn)算。采用主從模式，0號(hào)是master，其他是child（或者叫worker，as you wish）。

基本思路
兩個(gè)矩陣A，B進(jìn)行乘法運(yùn)算，則A的行 i 乘以B的列 j 得出的數(shù)是新矩陣(i,j)坐標(biāo)的數(shù)值。A(MN) B(NK)最后矩陣是M*K的，實(shí)驗(yàn)中M=N=K=1000，我也就沒(méi)有明確區(qū)分MNK，全部用MATRIX_SIZE定義的。
最簡(jiǎn)單的思路就是每個(gè)worker分配(MATRIX_SIZE/(numprocess-1))個(gè)，然后如果有余下的，就分給余數(shù)對(duì)應(yīng)的worker。比如MATRIX_SIZE=10，numprocess=4 則實(shí)際的worker有3個(gè)，每個(gè)人分3行，最后的一行給id是1的。可以很簡(jiǎn)單的利用循環(huán)類分配。最后Master收集所有的結(jié)果，并按照順序組裝起來(lái)就行。
每個(gè)worker的工作就是接收來(lái)自master的一行，和B矩陣運(yùn)算，得出新一行的結(jié)果，然后發(fā)送回master

代碼
多加了很多注釋來(lái)解釋，函數(shù)的說(shuō)明下一節(jié)解釋下。

#include <mpi.h>
#include <stdio.h>
#define MATRIX_SIZE 10
#define FROM_MASTER 1 //這里的類型可以區(qū)分消息的種類，以便區(qū)分worker發(fā)送來(lái)的結(jié)果
#define FROM_CHILD 2 
#define MASTER 0
MPI_Status status;
int myid,numprocess;
//最終保存的結(jié)果
int ans [MATRIX_SIZE*MATRIX_SIZE];
int A[MATRIX_SIZE*MATRIX_SIZE],B[MATRIX_SIZE*MATRIX_SIZE];
//讀取文件，注意讀取文件要放在master里，不然會(huì)讀兩遍，出現(xiàn)錯(cuò)誤
void readFile(){
  FILE* fina,*finb;
  fina=fopen("a.txt","r");
  int i;
  for (i = 0; i < MATRIX_SIZE*MATRIX_SIZE ; ++i)
  {
    fscanf(fina,"%d ",&A[i]);
  }
  fclose(fina);
  finb=fopen("b.txt","r");
  for(i=0;i<MATRIX_SIZE*MATRIX_SIZE;i++)
    fscanf(finb,"%d ",&B[i]);
  fclose(finb);
  printf("read file ok\n");
}
int master(){
  int workid,dest,i,j;
  printf("numprocess %d\n",numprocess );
  //給每個(gè)worker發(fā)送B矩陣過(guò)去
  for(i=0;i<numprocess-1;i++){
    //send B matrix
    MPI_Send(&B,MATRIX_SIZE*MATRIX_SIZE,MPI_INT,i+1,FROM_MASTER,MPI_COMM_WORLD);
  }
  //開(kāi)始給每個(gè)worker分配任務(wù)，取模即可
  for (i = 0; i < MATRIX_SIZE; i++)
  {
    //attention: num of workers is numprocess-1
    workid=i%(numprocess-1)+1;
    //send single line in A
    MPI_Send(&A[i*MATRIX_SIZE],MATRIX_SIZE,MPI_INT,workid,FROM_MASTER,MPI_COMM_WORLD);
  }

  //等待從worker發(fā)送來(lái)的數(shù)據(jù)
  int tempLine[MATRIX_SIZE];
  for (i = 0; i < MATRIX_SIZE*MATRIX_SIZE; i++)
  {
    ans[i]=0;
  }
  for (i = 0; i < MATRIX_SIZE; ++i)
  {
    int myprocess=i%(numprocess-1)+1;
    printf("Master is waiting %d\n",myprocess);
    //receive every line from every process
    MPI_Recv(&tempLine,MATRIX_SIZE,MPI_INT,myprocess,FROM_CHILD,MPI_COMM_WORLD,&status);
    //發(fā)送過(guò)來(lái)的都是計(jì)算好了的一行的數(shù)據(jù)，直接組裝到ans里就行
    for(j=0;j<MATRIX_SIZE;j++){
      ans[MATRIX_SIZE*i+j]=tempLine[j];
    }
    printf("Master gets %d\n",i);
  }

  for(i=0;i<MATRIX_SIZE*MATRIX_SIZE;i++){
    printf("%d ",ans[i] );
    if(i%MATRIX_SIZE==(MATRIX_SIZE-1))printf("\n");
  }
  printf("The Master is out\n");

}
int worker(){
  int mA[MATRIX_SIZE],mB[MATRIX_SIZE*MATRIX_SIZE],mC[MATRIX_SIZE];
  int i,j,bi;
  MPI_Recv(&mB,MATRIX_SIZE*MATRIX_SIZE,MPI_INT,MASTER,FROM_MASTER,MPI_COMM_WORLD,&status);
  //接收來(lái)自master的A的行
  for(i=0;i<MATRIX_SIZE/(numprocess-1);i++){
    MPI_Recv(&mA,MATRIX_SIZE,MPI_INT,MASTER,FROM_MASTER,MPI_COMM_WORLD,&status);
    //矩陣乘法，A 的一行和B矩陣相乘
    for(bi=0;bi<MATRIX_SIZE;bi++){
      mC[bi]=0;
      for(j=0;j<MATRIX_SIZE;j++){
        mC[bi]+=mA[j]*mB[bi*MATRIX_SIZE+j];
      }
    }
    MPI_Send(&mC,MATRIX_SIZE,MPI_INT,MASTER,FROM_CHILD,MPI_COMM_WORLD);
  }
  //如果處于余數(shù)范圍內(nèi)，則需要多計(jì)算一行
  if(MATRIX_SIZE%(numprocess-1)!=0){
    if (myid<=(MATRIX_SIZE%(numprocess-1)))
    {
      MPI_Recv(&mA,MATRIX_SIZE,MPI_INT,MASTER,FROM_MASTER,MPI_COMM_WORLD,&status);
      for(bi=0;bi<MATRIX_SIZE;bi++){
        mC[bi]=0;
        for(j=0;j<MATRIX_SIZE;j++){
          mC[bi]+=mA[j]*mB[bi*MATRIX_SIZE+j];
        }
      }
      MPI_Send(&mC,MATRIX_SIZE,MPI_INT,MASTER,FROM_CHILD,MPI_COMM_WORLD);
    }
  }
  printf("The worker %d is out\n",myid);
}
int main(int argc, char **argv)
{
  MPI_Init (&argc, &argv); 
  MPI_Comm_rank(MPI_COMM_WORLD,&myid);
  MPI_Comm_size(MPI_COMM_WORLD,&numprocess);

  if(myid==MASTER){
    readFile();
    master();
  }
  if(myid>MASTER){
    worker();
  }

  MPI_Finalize();
  return 0;
}

OPENMPI簡(jiǎn)單函數(shù)介紹
針對(duì)實(shí)驗(yàn)用到的幾個(gè)函數(shù)進(jìn)行說(shuō)明。
MPI為程序員提供一個(gè)并行環(huán)境庫(kù)，程序員通過(guò)調(diào)用MPI的庫(kù)程序來(lái)達(dá)到程序員所要達(dá)到的并行目的，可以只使用其中的6個(gè)最基本的函數(shù)就能編寫一個(gè)完整的MPI程序去求解很多問(wèn)題。這6個(gè)基本函數(shù)，包括啟動(dòng)和結(jié)束MPI環(huán)境，識(shí)別進(jìn)程以及發(fā)送和接收消息：
理論上說(shuō)，MPI所有的通信功能可以用它的六個(gè)基本的調(diào)用來(lái)實(shí)現(xiàn)：

MPI_INIT 啟動(dòng)MPI環(huán)境
MPI_COMM_SIZE 確定進(jìn)程數(shù)
MPI_COMM_RANK 確定自己的進(jìn)程標(biāo)識(shí)符
MPI_SEND 發(fā)送一條消息
MPI_RECV 接收一條消息
MPI_FINALIZE 結(jié)束MPI環(huán)境

初始化和結(jié)束

MPI初始化：通過(guò)MPI_Init函數(shù)進(jìn)入MPI環(huán)境并完成所有的初始化工作。

int MPI_Init( int *argc, char * * * argv )

MPI結(jié)束：通過(guò)MPI_Finalize函數(shù)從MPI環(huán)境中退出。

int MPI_Finalize(void)

獲取進(jìn)程的編號(hào)
調(diào)用MPI_Comm_rank函數(shù)獲得當(dāng)前進(jìn)程在指定通信域中的編號(hào)，將自身與其他程序區(qū)分。

int MPI_Comm_rank(MPI_Comm comm, int *rank)

獲取指定通信域的進(jìn)程數(shù)
調(diào)用MPI_Comm_size函數(shù)獲取指定通信域的進(jìn)程個(gè)數(shù)，確定自身完成任務(wù)比例。

int MPI_Comm_size(MPI_Comm comm, int *size)

MPI消息
一個(gè)消息好比一封信
消息的內(nèi)容的內(nèi)容即信的內(nèi)容，在MPI中成為消息緩沖(Message Buffer)
消息的接收發(fā)送者即信的地址，在MPI中成為消息封裝(Message Envelop)
MPI中，消息緩沖由三元組<起始地址，數(shù)據(jù)個(gè)數(shù)，數(shù)據(jù)類型>標(biāo)識(shí)
消息信封由三元組<源/目標(biāo)進(jìn)程，消息標(biāo)簽，通信域>標(biāo)識(shí)

消息發(fā)送

MPI_Send函數(shù)用于發(fā)送一個(gè)消息到目標(biāo)進(jìn)程。

int MPI_Send(void *buf, int count, MPI_Datatype dataytpe, int dest, int tag, MPI_Comm comm)

buf是要發(fā)送數(shù)據(jù)的指針，比如一個(gè)A數(shù)組，可以直接&A，count則是數(shù)據(jù)長(zhǎng)度，datatype都要改成MPI的type。dest就是worker的id了。tag則可以通過(guò)不同的type來(lái)區(qū)分消息類型，比如是master發(fā)送的還是worker發(fā)送的。

消息接收

MPI_Recv函數(shù)用于從指定進(jìn)程接收一個(gè)消息

 int MPI_Recv(void *buf, int count, MPI_Datatype datatyepe,int source, int tag, MPI_Comm comm, MPI_Status *status)

編譯和執(zhí)行
生成執(zhí)行文件data

mpicc -o programname programname.c

一個(gè)MPI并行程序由若干個(gè)并發(fā)進(jìn)程組成，這些進(jìn)程可以相同也可以不同。MPI只支持靜態(tài)進(jìn)程創(chuàng)建，即：每個(gè)進(jìn)程在執(zhí)行前必須在MPI環(huán)境中登記，且它們必須一起啟動(dòng)。通常啟動(dòng)可執(zhí)行的MPI程序是通過(guò)命令行來(lái)實(shí)現(xiàn)的。啟動(dòng)方法由具體實(shí)現(xiàn)確定。例如在MPICH實(shí)現(xiàn)中通過(guò)下列命令行可同時(shí)在獨(dú)立的機(jī)器上啟動(dòng)相同的可執(zhí)行程序：

 mpirun –np N programname

其中N是同時(shí)運(yùn)行的進(jìn)程的個(gè)數(shù)，programname是可執(zhí)行的MPI程序的程序名。

上一篇：C語(yǔ)言在屏幕上輸出楊輝三角

欄目：C語(yǔ)言

下一篇：C++中strcpy函數(shù)的實(shí)現(xiàn)

本文標(biāo)題：安裝OpenMPI來(lái)配合C語(yǔ)言程序進(jìn)行并行計(jì)算

本文地址：http://mengdiqiu.com.cn/a1/Cyuyan/2659.html

更多C語(yǔ)言

欧美大屁股bbbbxxxx,狼人大香伊蕉国产www亚洲,男ji大巴进入女人的视频小说,男人把ji大巴放进女人免费视频,免费情侣作爱视频

C語(yǔ)言

安裝OpenMPI來(lái)配合C語(yǔ)言程序進(jìn)行并行計(jì)算

您可能感興趣的文章

閱讀排行

本欄相關(guān)

隨機(jī)閱讀