MPI_Gather не работает

Мне нужно было использовать функцию MPI_Gather в моей программе умножения матриц, но последние пару дней у меня с ней проблемы.

Поэтому я написал простую программу MPI, используя только функцию сбора, и пытался заставить ее работать... Для этого я сослался на книгу «Введение в параллельное программирование» Питера Пачеко.

Программа сразу завершает работу, не давая мне абсолютно никаких результатов или ошибок... Я до сих пор не смог понять ошибку.

   /******************NOTE**********************

      The program simply uses the MPI_Gather() function. 
      The program is exiting directly.
      I have written it for only TWO processes.
      (-np 2)

   ******************************************/

    #include<stdio.h>
    #include"mpi.h"

    int main()
    {
     int i,j,proc,rank;
     double d[4];
     double local_a[2];


     MPI_Init(NULL,NULL);
     MPI_Comm_size(MPI_COMM_WORLD, &proc);
     MPI_Comm_rank(MPI_COMM_WORLD, &rank);

     if(rank==0)
     { 

       local_a[0]=1.0;
       local_a[1]=2.0;
     }

     else
     {
       local_a[0]=3.0;
       local_a[1]=4.0;
     }

     int local=2;

     if(rank==0)
     {   

       MPI_Gather(local_a,local,MPI_DOUBLE,d,local,MPI_DOUBLE,0,MPI_COMM_WORLD);
     //MPI_Gather(&local_a,local,MPI_DOUBLE,&d,local,MPI_DOUBLE,0,MPI_COMM_WORLD);
     //also tried the above line just to be certain.


       printf("\n");
       for(j=0;j<4;j++)
         printf("\t%f",d[j]);
     }
     else
     {
       MPI_Gather(local_a,local,MPI_DOUBLE,d,local,MPI_DOUBLE,0,MPI_COMM_WORLD);
     //MPI_Gather(&local_a,local,MPI_DOUBLE,&d,local,MPI_DOUBLE,0,MPI_COMM_WORLD); 
     }

     MPI_Finalize();

     return 0;

    }

Может кто-нибудь, пожалуйста, помогите мне.

Спасибо.

Анага Мадхусуданан


c mpi
person user1715985    schedule 03.10.2012    source источник
comment
в чем твоя ошибка? Как мы можем помочь, не зная, какую ошибку вы получаете для кода, который компилируется и работает нормально в моей системе   -  person pyCthon    schedule 04.10.2012

Ответы (1)


arrow_upward
2
arrow_downward

Ваша программа отлично работает для меня, выдавая на выходе:

    1.000000    2.000000    3.000000    4.000000

Можете ли вы поделиться дополнительной информацией о том, как вы запускаете и компилируете исполняемый файл, чтобы я попытался воспроизвести ошибку и, в случае необходимости, отредактировал ответ?

К вашему сведению, ниже вы можете найти немного измененную версию программы, из которой видно, что приемный буфер может быть выделен только в корне:

#include <stdio.h>
#include <stdlib.h>
#include "mpi.h"

int main() {

  int rank     = -1;
  int commsize = -1;
  double sendbuffer[2];

  MPI_Init(NULL,NULL);
  MPI_Comm_size(MPI_COMM_WORLD, &commsize);
  MPI_Comm_rank(MPI_COMM_WORLD, &rank);

  sendbuffer[0]=2.0*rank;
  sendbuffer[1]=2.0*rank + 1;

  int count=2;

  if(rank==0) {   

    // Recvbuffer is significant only at root
    double * recvbuffer = malloc(2*commsize*sizeof(double));
    // Gather values at root
    MPI_Gather(sendbuffer,count,MPI_DOUBLE,recvbuffer,count,MPI_DOUBLE,0,MPI_COMM_WORLD);
    // Print to screen
    printf("\n");
    for(int jj=0; jj<2*commsize ;jj++)
      printf("%f\n",recvbuffer[jj]);
    // Free recvbuffer
    free(recvbuffer);

  } else {
    MPI_Gather(sendbuffer,count,MPI_DOUBLE,NULL,0,MPI_DOUBLE,0,MPI_COMM_WORLD);
  }

  MPI_Finalize();

  return 0;

}
person Massimiliano    schedule 03.10.2012
comment
Можно также использовать сбор на месте, чтобы использовать только один буфер в основном ранге. - person Hristo Iliev; 04.10.2012