c++ - 在 MPI 中使用 Armadillo 稀疏矩阵

Question

我正在尝试初始化 Armadillo 稀疏矩阵 sp_mat在 MPI 中如下:

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 也无济于事(我遇到 undefined reference 错误)。
我如何初始化 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);