c - MPI_Finalize() がファイナライズしない

Question

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

int main(int argc, char *argv[])
{
    int i, done = 0, n;
    double PI25DT = 3.141592653589793238462643;
    double pi, tmp, h, sum, x;

    int numprocs, rank;
    MPI_Status status;

    MPI_Init(&argc, &argv);
    MPI_Comm_size(MPI_COMM_WORLD, &numprocs);
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);

    if (numprocs < 1)
         printf("numprocs = %d, should be run with numprocs > 1\n", numprocs);
    else {
        while (!done)
        {
            if (rank == 0) {
                printf("Enter the number of intervals: (0 quits) \n");
                scanf("%d",&n);
            }
            if (n == 0) break;

            MPI_Bcast(&n, 1, MPI_INT, 0, MPI_COMM_WORLD);

            h   = 1.0 / (double) n;
            sum = 0.0;
            for (i = 1 + rank; i <= n; i+=numprocs) {
                x = h * ((double)i - 0.5);
                sum += 4.0 / (1.0 + x*x);
            }

            MPI_Reduce(&sum, &tmp, 1, MPI_DOUBLE, MPI_SUM, 0, MPI_COMM_WORLD);
            if (rank == 0) {
                pi = h * tmp;
                printf("pi is approximately %.16f, Error is %.16f\n", pi, fabs(pi - PI25DT));
            }
        }
    }
    MPI_Finalize();
}

プロセス数 = 10 でこのプログラムを実行し、scanf("%d", &n); に 0 を入力した後。プログラムがファイナライズされていないようで、Ctrl+C で閉じる必要があります。

合計 10 個のプロセスが強制終了されました (一部はクリーンアップ中に mpirun によって実行された可能性があります) mpirun: クリーンな終了が完了しました

プロセスを強制終了した後にコンソールに表示されます。

Answer 1

そうですね。ここで何が起こっているか見てください：

        if (rank == 0) {
            printf("Enter the number of intervals: (0 quits) \n");
            scanf("%d",&n);
        }
        if (n == 0) break;

        MPI_Bcast(&n, 1, MPI_INT, 0, MPI_COMM_WORLD);

ランク 0 で '0' を入力すると、ランク 0 が壊れますが、他のすべてのランクは、前回 n で持っていたもの (または、最初の反復の場合、そのメモリ位置にあったランダムなもの) をすべて持ちます)。したがって、ランク 1-(N-1) は に入り、MPI_Bcast何があるかを知るために放送に参加するのを辛抱強く待っていますn。その間、ランク 0 はループを終了し、MPI_Finalize他の全員がどこにいるのか不思議に思っています。

これらの行を反転するだけです：

        MPI_Bcast(&n, 1, MPI_INT, 0, MPI_COMM_WORLD);
        if (n == 0) break;

今では、ループから抜け出すかどうかを全員が同時に決定します。

Answer 2

無限の while ループがあるため、MPI_Finalize() を実行できません。