MPI, нахождение простых чисел

@Ide · Регистрация: 12.06.2013

Здравсвуйте!
Пишу программу для нахождения простых чисел. Пишу в Visual Studio 2008 на C++.
Программу написал и она работает на нескольких ПК. Но проблема в том, что результаты в главный процесс возвращаются неправильно.
Допустим запускаем:
компьютер с главным процессом: 5 процессов
компьютер раб: 5 процессов.

В результате в главный файл запишутся результаты от первых 5 и только от ОДНОГО от компьютера раба. Хотя все 5 там отработали и результаты получили.

Главный процесс остается висеть в памяти, и 4 процесса на другом ПК - тоже. Т.е. возникает тупик. Видимо первый ждет от них результата, а из тех отсылает только первый.

Исходник:

C++

    MPI_Init(&argc, &argv);
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);
    MPI_Comm_size(MPI_COMM_WORLD, &size);
 
if (rank == 0) 
    {  
 
    FILE *f;
    double *a;
    char name[6]="0.txt";
    name[0]=rank+48;
    f=fopen(name,"w");
  int position = 0;
    int how_many=0;
  fprintf(f,"N %i\n",N);
        for(int i=1; i<size; i++) 
       { 
        position=0;
        fprintf(f," %i\n",i);
       
        MPI_Recv(&how_many,1,MPI_INT,i,77,MPI_COMM_WORLD,&status);
        int*  my_array = (int*)malloc(how_many * sizeof(int));
        fprintf(f,"Array size %i\n",how_many);
            MPI_Recv(my_array,how_many,MPI_INT,i,78,MPI_COMM_WORLD,&status);
 
 
                for(int j=0; j<(int)how_many; j++) 
                 {
                fprintf(f,"My rank %i I got : %i from %i\n",rank,my_array[j],i);
                 }
 
 
}
 
 
printf("My rank %d\n", rank);
fclose(f);
    }
else
    {   
 
 
    FILE *f;
    double *a;
    char name[6]="0.txt";
    name[0]=rank+48;
    f=fopen(name,"w");
 
        fprintf(f,"A: %i  \n",min);
        fprintf(f,"B: %i  \n",max);
        
    
         myAL= IsSimple(min, max);
        for (int i=0; i<myAL->Count; i++)
        {
            fprintf(f,"Counted: %i  \n",(int)myAL[i]);
            }
 
 
    int kol = myAL->Count;
    int*  my_array = new int[myAL->Count];
 
            for (int i=0; i<myAL->Count; i++)
        {
            my_array[i] = (int)myAL[i];
        }
 
   fprintf(f,"Kol: %i  \n",(int)kol);
   fclose(f);
 
  int position = 0;
 
   MPI_Send(&kol, 1, MPI_INT, 0, 77, MPI_COMM_WORLD);
   MPI_Send(my_array, myAL->Count, MPI_INT, 0, 78, MPI_COMM_WORLD);
        
}

Пожалуйста помогите, бьюсь уже 2 недели.

Опции темы