0

いくつかの単純な練習アプリ以外で初めて MPI を使用しましたが、何かがうまくいきません。

次のメンバーで定義されたクラスがあります (読みやすさと画面スペースの節約のためにメソッドは省略されています)。

class particle
{
    public:
        double _lastUpdate;
    float _x, _y, _xvel, _yvel;
    bool _isStatic;
        bool _isForeign;
        float _size;

    private:
        int _isStaticInt;           // integer copy of _isStatic to be sent over MPI ( since there's no MPI_BOOL :C )
};

各粒子のいくつかの重要なメンバーの値を送信し、その場で他のメンバーを複製することにより、プロセス間で一連の粒子を送信したいと考えています。そのために、MPI データ型をそのように定義します。ご覧のとおり、メンバー _lastUpdate、_isStatic、および _isForeign は含まれていません。

MPI_Datatype types[] = { MPI_FLOAT, MPI_FLOAT, MPI_FLOAT, MPI_FLOAT, MPI_INTEGER, MPI_FLOAT };
std::vector<int> len(6, 1);
std::vector<MPI_Aint> disp(6, 0);
particle temp;
MPI_Aint base;
MPI_Address(&temp, &base);
MPI_Address(&temp._x, &disp[0]);
MPI_Address(&temp._y, &disp[1]);
MPI_Address(&temp._xvel, &disp[2]);
MPI_Address(&temp._yvel, &disp[3]);
MPI_Address(&temp._isStaticInt, &disp[4]);
MPI_Address(&temp._size, &disp[5]);
for (int i=0; i<6; ++i)
{
    disp[i] = disp[i] - base;
}
MPI_Type_struct(6, &len[0], &disp[0], types, &_particleType);
MPI_Type_commit(&_particleType);

これがパーティクルを送る方法です。'parts' は、送信したいパーティクル オブジェクトへのポインタが含まれる、particle* のベクトルで、'size' は parts.size() です。

std::vector<int> len(size, 1);
std::vector<MPI_Aint> disp(size, 0);
MPI_Aint base;
MPI_Address(parts[0], &base);                   // datatype begins at the first selected object
for (int select = 1; select < size; ++select)
{
    MPI_Address(parts[select], &disp[select]);
    disp[select] = disp[select] - base;
}
MPI_Type_hindexed(size, &len[0], &disp[0], _particleType, &_sendType);
MPI_Type_commit(&_sendType);

MPI_Request payload_req;
MPI_Isend(parts[0], 1, _sendType, ngb, 0, _cartesian_comm, &payload_req);

受信も同様に行われますが、この場合のみ、'parts' は、以前に作成された「空白の」パーティクル オブジェクトを指す、particle* のベクトルであり、そのメンバーは受信したデータで満たされます。

std::vector<int> len(size, 1);
std::vector<MPI_Aint> disp(size, 0);
MPI_Aint base;
MPI_Address(parts[0], &base);                       // datatype begins at the first newly inserted object
for (int part = 1; part < size; ++part)
{
    MPI_Address(parts[part], &disp[part]);
    disp[part] = disp[part] - base;
}
MPI_Type_hindexed(size, &len[0], &disp[0], _particleType, &_recvType);
MPI_Type_commit(&_recvType);
MPI_Status status;
MPI_Recv(parts[0], size, _particleType, ngb, 0, _cartesian_comm, &status);

問題は、最初のパーティクルを除くすべての受信パーティクルのメンバーにデフォルトの「空白」値があることです。私は以前に、いくつかの単純な値しか送信しなかったにもかかわらず、同様のことを行う小さなテスト アプリを作成しました。これは、私がここで気付いていないコーディングの間違いがない限り (完全に可能性があります)、この種のデータ型の策略が機能することを保証するものではなく、機能する小さなものは偶然にしか機能しないと私に信じさせます。

このタイプのメモリ操作が安全であり、信頼できるかどうかを確認/拒否できる人はいますか?

4

1 に答える 1

1

気にしないでください。質問を入力するだけで、何が問題なのかがわかったようです。

受信コマンドは次のとおりです。

MPI_Recv(parts[0], 1, _recvType, ngb, 0, _cartesian_comm, &status);
于 2009-07-23T11:46:08.167 に答える