Эйген разреженного матрица резерв NNZ

Question

мне было интересно, если метод reserve в Эйгене делает также выделяет память для внешних индексов и внутренних индексов (в случае 3-CSC.)

Это означает, что если я сделать

n = 1000 
SparseMatrix<float> A; 
A.reserve(n) 

затем

• Алло Cate n память данных
• выделить n памяти для внутренних индексов
• выделить n+1 памяти для внешних индексов (в худшем случае и априори не известно)

документация говорит только

void Eigen::SparseMatrix<...>::reserve (Index reserveSize) 

Предоставляет резервные копии без нулей.

Условие: матрица должна быть в сжатом режиме.

---

Может кто-нибудь разъяснить это? Или сказать что-то о политике резервирования? Спасибо

Answer 1

Нет, размеры матрицы должны быть довольно статичными, поэтому reserve(n) выделяет пространство для не нулей. Если размер неизвестен, вы можете начать с верхней границей:

SparseMatrix<float> A(n,n); 

резерва, например:

A.reserve(10*n); 

и как только вы сделали, и что вы знаете, окончательный размер, уменьшить его:

A.conservativeResize(rows,cols);