4

私には2つのプログラムがあります。いくつかの計算を実行する「ワーカー」を生成する「マスター」。マスターにワーカーから結果を取得して合計を保存してもらいたい。MPI_Reduceを使用してワーカーから結果を収集しようとしていますが、ワーカーはMPI_Reduceを使用してマスターMPI_Commに送信しています。それが正しいかどうかはわかりません。これが私のプログラムです：

マスター：

#include <mpi.h>
#include <iostream>
using namespace std;

int main(int argc, char *argv[]) { 
    int world_size, universe_size, *universe_sizep, flag; 

    int rc, send, recv;

    // intercommunicator
    MPI_Comm everyone;

    MPI_Init(&argc, &argv); 
    MPI_Comm_size(MPI_COMM_WORLD, &world_size); 

    if (world_size != 1) {
        cout << "Top heavy with management" << endl;
    } 

    MPI_Attr_get(MPI_COMM_WORLD, MPI_UNIVERSE_SIZE, &universe_sizep, &flag);  
    if (!flag) { 
        cout << "This MPI does not support UNIVERSE_SIZE. How many processes total?";
        cout << "Enter the universe size: ";
        cin >> universe_size; 
    } else {
        universe_size = *universe_sizep;
    }
    if (universe_size == 1) {
        cout << "No room to start workers" << endl;
    }

    MPI_Comm_spawn("so_worker", MPI_ARGV_NULL, universe_size-1,  
             MPI_INFO_NULL, 0, MPI_COMM_SELF, &everyone,  
             MPI_ERRCODES_IGNORE);

    send = 0;

    rc = MPI_Reduce(&send, &recv, 1, MPI_INT, MPI_SUM, 0, everyone);

    // store result of recv ...
    // other calculations here
    cout << "From spawned workers recv: " << recv << endl;

    MPI_Finalize(); 
    return 0; 
}

ワーカー：

#include <mpi.h>
#include <iostream>
using namespace std;

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

    int rc, send,recv;


    int parent_size, parent_id, my_id, numprocs; 
    // parent intercomm
    MPI_Comm parent; 
    MPI_Init(&argc, &argv); 

    MPI_Comm_get_parent(&parent); 
    if (parent == MPI_COMM_NULL) {
        cout << "No parent!" << endl;
    }
    MPI_Comm_remote_size(parent, &parent_size); 
    MPI_Comm_rank(parent, &parent_id) ; 
    //cout << "Parent is of size: " << size << endl;
    if (parent_size != 1) {
        cout << "Something's wrong with the parent" << endl;
    }

    MPI_Comm_rank(MPI_COMM_WORLD, &my_id) ;     
    MPI_Comm_size(MPI_COMM_WORLD, &numprocs) ;  

    cout << "I'm child process rank "<< my_id << " and we are " << numprocs << endl;
    cout << "The parent process rank "<< parent_id << " and we are " << parent_size << endl;

    // get value of send
    send = 7; // just an example
    recv = 0;

    rc = MPI_Reduce(&send, &recv, 1, MPI_INT, MPI_SUM, parent_id, parent);
    if (rc != MPI_SUCCESS)
        cout << my_id << " failure on mpi_reduce in WORKER" << endl;

    MPI_Finalize(); 
    return 0; 
}

私は両方をコンパイルし、次のように実行します（osxの場合はmpic ++）：

mpic++ so_worker.cpp -o so_worker
mpic++ so_master.cpp -o so_master
mpirun -n 1 so_master

これは、ワーカーを生成するマスターを実行する正しい方法ですか？

マスターでは、MPI_Reduceから常に0が返されます。インターコミュニケーターからMPI_reduceを使用できますか、それともワーカーからMPI_Sendを使用し、マスターからMPI_Recvを使用する必要がありますか？なぜ機能しないのかよくわかりません。

どんな助けでもいただければ幸いです。ありがとう！

4

1 に答える 1

7

MPI_Comm_get_parent元のプロセスと生成されたすべてのプロセスを含む親インターコミュニケーターを返します。この場合、呼び出しMPI_Comm_rank(parent, &parent_id)は親のランクを返すのではなく、インターコミュニケーターのローカルグループ内の現在のプロセスのランクを返します。

I'm child process rank 0 and we are 3
The parent process **rank 0** and we are 1
I'm child process rank 1 and we are 3
The parent process **rank 1** and we are 1
I'm child process rank 2 and we are 3
The parent process **rank 2** and we are 1

（強調表示された値がどのように異なるかを観察してください-親プロセスのランクは同じである必要がありますね？）

これがMPI_Reduce()、すべてのワーカープロセスがルートランクに異なる値を指定するため、呼び出しが成功しない理由です。元々マスタープロセスが1つあったため、のリモートグループでのランクはでparentあり0、したがって、すべてのワーカーは0ルートとして次のように指定する必要がありますMPI_Reduce。

//
// Worker code
//
rc = MPI_Reduce(&send, &recv, 1, MPI_INT, MPI_SUM, 0, parent);

これは問題の半分にすぎません。残りの半分は、根付いた集合的な操作（たとえばMPI_REDUCE）がインターコミュニケーターとは少し異なる動作をすることです。最初に、2つのグループのどちらがルートをホストするかを決定する必要があります。ルートグループが識別されると、ルートプロセスはinMPI_ROOTの値として渡される必要があり、ルートグループ内の他のすべてのプロセスは渡される必要があります。つまり、受信グループのプロセスは、ルート化された集合操作にはまったく参加しません。マスターコードは、マスターのグループに1つのプロセスしか存在できないように記述されているため、マスターコードの呼び出しを次のように変更するだけで十分です。rootMPI_REDUCEMPI_PROC_NULLMPI_Reduce

//
// Master code
//
rc = MPI_Reduce(&send, &recv, 1, MPI_INT, MPI_SUM, MPI_ROOT, everyone);

マスターはリダクション操作自体にも参加しないことに注意してください。たとえば、ルートはリダクション対象のデータを送信しないため、sendbuf（この場合）の値は関係ありません。&sendリモートグループ内のプロセス。

于 2013-02-06T16:30:11.267 に答える