在 MPI 中使用犰狳稀疏矩阵

Question

我正尝试在 MPI 中初始化 Armadillo 稀疏矩阵 sp_mat，如下所示：

if(rank==0)
{  // some code for locations, values
    sp_mat X(locations,values)
}
// this is where I want to use X
if(rank==0)
   some_fun(X)

如您所见，Armadillo 构造函数是 if block 的本地构造函数，因此在 if block 之后无法使用它。换一种方式提出同样的问题：

extern sp_mat X
if(rank==0)
{  // some code for locations, values
    sp_mat X(locations,values)
}
// this is where I want to use X
 if(rank==0)
   some_fun(X)

在 if block 之前使用 extern sp_mat X 也无济于事（我遇到未定义的引用错误）。
如何初始化 X 并在之后重新使用它？

Answer 1

使用（智能）指针：

std::unique_ptr<sp_mat> X; // or std::shared_ptr<sp_mat> or sp_mat*
if (rank == 0) {
   // some code for locations and values     
   X = std::unique_ptr<sp_mat>(new sp_mat(locations, values));
}
...
if (rank == 0)
   some_fun(*X);

Answer 2

与其使用指针技巧，不如使用 C++11 中更简洁的 std::move()：

sp_mat X;

if(rank==0)
{  // some code for locations, values
    X = std::move( sp_mat(locations,values) );
}

在 MPI 中使用犰狳稀疏矩阵

Using Armadillo Sparse Matrices in MPI

c++

mpi

armadillo