разреженная матрица из триплета

Question

С собственному (C++) я хочу, чтобы иметь такой же результат, как Matlab в одном:

% Matlab code: Create a 1500-by-1500 sparse matrix from the triplets i, j, 
 
      and v 
 
i = [900 1000]; 
 
    j = [900 1000]; 
 
    v = [10 100]; 
 
    S = sparse(i,j,v,1500,1500) 
 

 
    %result 
 
    S = 
 

 
(900,900)  10 
 
(1000,1000) 100

Я хочу сделать то же самое в C++, я видел в разных сообщения, которые мы можем обрабатывать разреженные матрицы с Eigen Library, но я действительно не понимал, как это сделать. Вы можете мне помочь ?

Благодаря

Answer 1

Создать вектор троек

typedef Eigen::Triplet<int> Trip; 
std::vector<Trip> trp; 

Создать триплеты

trp.push_back(Trip(900,900,10)); // (index, index, value) 
trp.push_back(Trip(1000,1000,100)); 

Назначьте их в разреженную матрицу Эйгеном

int rows, cols; 
rows = cols = 1500; 
Eigen::SparseMatrix<int> A(rows,cols); 

A.setFromTriplets(trp.begin(), trp.end()); 

Распечатать все (плотная) матрица

std::cout << A << std::endl; 

или печатать только индексы и значения

std::cout << "Row\tCol\tVal" <<std::endl; 
for (int k=0; k < A.outerSize(); ++k) 
{ 
    for (SparseMatrix<int>::InnerIterator it(A,k); it; ++it) 
    { 
     std::cout << it.row() << "\t"; // row index 
     std::cout << it.col() << "\t"; // col index (here it is equal to k) 
     std::cout << it.value() << std::endl; 
    } 
}