Две функции MPI_Allreduce() не работают, выдает ошибку NULL Communicator

Question

Я использую пример кода из книги MPI [скоро сообщу название].

Что он делает, это следующее:

а) создает два коммуникатора world = MPI_COMM_WORLD содержащий все процессы и worker что исключает random number generator server (последний процесс ранга).

б) Итак, server генерирует случайные числа и передает их workers на requests от workers,

в) Что workers сделать, они считают отдельно количество образцов падения inside а также outside unit circle внутри unit square,

г) после достаточного уровня точности подсчеты inside а также outside являются Allreduce d, чтобы вычислить значение PI как их соотношение.

** Код хорошо компилируется. Однако при запуске с помощью следующей команды (на самом деле с любым значением n) **

>mpiexec -n 2 apple.exe 0.0001

Я получаю следующие ошибки:

Fatal error in MPI_Allreduce: Invalid communicator, error stack:
MPI_Allreduce(855): MPI_Allreduce(sbuf=000000000022EDCC, rbuf=000000000022EDDC,
count=1, MPI_INT, MPI_SUM, MPI_COMM_NULL) failed
MPI_Allreduce(780): Null communicator
pi =  0.00000000000000000000
job aborted:
rank: node: exit code[: error message]
0: PC: 1: process 0 exited without calling finalize
1: PC: 123

Редактировать: (((Удалено: Но когда я удаляю любой из двух MPI_Allreduce() функции, он работает без каких-либо ошибок во время выполнения, хотя и с неправильным ответом.))

Код:

#include <mpi.h>
#include <mpe.h>
#include <stdlib.h>

#define CHUNKSIZE 1000
/* message tags */
#define REQUEST 1
#define REPLY 2

int main(int argc, char *argv[])
{
    int iter;
    int in, out, i, iters, max, ix, iy, ranks [1], done, temp;
    double x, y, Pi, error, epsilon;
    int numprocs, myid, server, totalin, totalout, workerid;
    int rands[CHUNKSIZE], request;
    MPI_Comm world, workers;
    MPI_Group world_group, worker_group;
    MPI_Status status;
    MPI_Init(&argc,&argv);
    world = MPI_COMM_WORLD;
    MPI_Comm_size(world,&numprocs);
    MPI_Comm_rank(world,&myid);
    server = numprocs-1; /* last proc is server */
    if(myid==0) sscanf(argv[1], "%lf", &epsilon);
    MPI_Bcast(&epsilon, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD);
    MPI_Comm_group(world, &world_group);
    ranks[0] = server;
    MPI_Group_excl(world_group, 1, ranks, &worker_group);
    MPI_Comm_create(world, worker_group, &workers);
    MPI_Group_free(&worker_group);
    if(myid==server)   /* I am the rand server */
    {
        srand(time(NULL));
        do
        {
            MPI_Recv(&request, 1, MPI_INT, MPI_ANY_SOURCE, REQUEST, world, &status);
            if(request)
            {
                for(i=0; i<CHUNKSIZE;)
                {
                    rands[i] = rand();
                    if(rands[i]<=INT_MAX) ++i;
                }
                MPI_Send(rands, CHUNKSIZE, MPI_INT,status.MPI_SOURCE, REPLY, world);
            }
        }
        while(request>0);
    }
    else   /* I am a worker process */
    {
        request = 1;
        done = in = out = 0;
        max = INT_MAX; /* max int, for normalization */
        MPI_Send(&request, 1, MPI_INT, server, REQUEST, world);
        MPI_Comm_rank(workers, &workerid);
        iter = 0;
        while(!done)
        {
            ++iter;
            request = 1;
            MPI_Recv(rands, CHUNKSIZE, MPI_INT, server, REPLY, world, &status);
            for(i=0; i<CHUNKSIZE;)
            {
                x = (((double) rands[i++])/max)*2-1;
                y = (((double) rands[i++])/max)*2-1;
                if(x*x+y*y<1.0) ++in;
                else ++out;
            }

            /* ** see error here ** */
            MPI_Allreduce(&in, &totalin, 1, MPI_INT, MPI_SUM, workers);
            MPI_Allreduce(&out, &totalout, 1, MPI_INT, MPI_SUM, workers);
            /* only one of the above two MPI_Allreduce() functions working */

            Pi = (4.0*totalin)/(totalin+totalout);
            error = fabs( Pi-3.141592653589793238462643);
            done = (error<epsilon||(totalin+totalout)>1000000);
            request = (done)?0:1;
            if(myid==0)
            {
                printf("\rpi = %23.20f", Pi);
                MPI_Send(&request, 1, MPI_INT, server, REQUEST, world);
            }
            else
            {
                if(request)
                    MPI_Send(&request, 1, MPI_INT, server, REQUEST, world);
            }
            MPI_Comm_free(&workers);
        }
    }
    if(myid==0)
    {
        printf("\npoints: %d\nin: %d, out: %d, <ret> to exit\n", totalin+totalout, totalin, totalout);
        getchar();
    }
    MPI_Finalize();
}

В чем здесь ошибка? Я что-то пропустил? Любая помощь или указатель будет высоко оценен.

Answer 1

Вы освобождаете workers коммуникатор, прежде чем вы закончите с его использованием. Переместить MPI_Comm_free(&workers) позвонить после while(!done) { ... } петля.