在MPI中使用Armadillo稀疏矩阵

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) );
}