MPIを使ってEigen :: MatrixXdでデータを送信する方法

Question

マシン間でmatrcesを送信するためにMPIを使用したいと思います。以下は私のテストコードです

#include <iostream> 
#include <Eigen/Dense> 
#include <mpi.h> 
using std::cin; 
using std::cout; 
using std::endl; 
using namespace Eigen; 

int main(int argc, char ** argv) 
{ 
    MatrixXd a = MatrixXd::Ones(3, 4); 
    int myrank; 
    MPI_Init(&argc, &argv); 
    MPI_Comm_rank(MPI_COMM_WORLD, &myrank); 
    MPI_Status status; 
    if (0 == myrank) 
    { 
     MPI_Send(&a, 96, MPI_BYTE, 1, 99, MPI_COMM_WORLD); 
    } 
    else if (1 == myrank) 
    { 
     MPI_Recv(&a, 96, MPI_BYTE, 0, 99, MPI_COMM_WORLD, &status); 
     cout << "RANK " << myrank << endl; 
     cout << a << endl; 
    } 
    MPI_Finalize(); 
    return 0; 
} 

エラーなしで正常にコンパイルされましたが、起動したときに次のエラーが返されました。

$ MPI mpiexec -n 2 ./sendMatrixTest 
RANK 1 
[HPNotebook:11633] *** Process received signal *** 
[HPNotebook:11633] Signal: Segmentation fault (11) 
[HPNotebook:11633] Signal code: Address not mapped (1) 
[HPNotebook:11633] Failing at address: 0xf4ba40 
-------------------------------------------------------------------------- 
mpiexec noticed that process rank 1 with PID 11633 on node HPNotebook exited on signal 11 (Segmentation fault). 
-------------------------------------------------------------------------- 

どうすればよいですか？ありがとうございました！

Answer 1

@Mattが指摘しているように、MatrixXdコンテナにはデータだけでなく、しかし、ここにいるので、あなたは、これが機能するように、あなただけの、data() methodを使用し、昔ながらのデータへのポインタを取得することができ、マトリックスのサイズと種類を知っている：

#include <iostream> 
#include <Eigen/Dense> 
#include <mpi.h> 
using std::cin; 
using std::cout; 
using std::endl; 
using namespace Eigen; 

int main(int argc, char ** argv) 
{ 
    MatrixXd a = MatrixXd::Ones(3, 4); 
    int myrank; 
    MPI_Init(&argc, &argv); 
    MPI_Comm_rank(MPI_COMM_WORLD, &myrank); 
    MPI_Status status; 
    if (0 == myrank) 
    { 
     MPI_Send(a.data(), 12, MPI_DOUBLE, 1, 99, MPI_COMM_WORLD); 
    } 
    else if (1 == myrank) 
    { 
     MPI_Recv(a.data(), 12, MPI_DOUBLE, 0, 99, MPI_COMM_WORLD, &status); 
     cout << "RANK " << myrank << endl; 
     cout << a << endl; 
    } 
    MPI_Finalize(); 
    return 0; 
}