Использование Армадилла разреженных матриц в MPI

Question

Я пытаюсь инициализирует Armadillo разреженный matriix 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) 

extern sp_mat X Использование Перед if block также не помогает (я получил неопределенную ссылочную ошибку).
Как инициализировать 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

Вместо того, с помощью указателя трюков, намного чище std::move() от C++ 11 могут быть использованы:

sp_mat X; 

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