Лучший способ обнулить коэффициенты sparsematrix в Eigen?

Question

Типа разреженной матрицы в библиотеке Эйген C++, скажем

SparseMatrix<double, RowMajor> A; 

, кажется, отсутствует coeffFind (I, J) метод с InnerIterator или двойным * в качестве типа возвращаемого значения. Таким образом, чтобы гарантировать, что не имеет ненулевого коэффициента в позиции (I, J), я могу либо сделать

A.coeffRef(i,j) = 0.0; // Adds an entry with zero value! 

или

if (A.coeff(i,j) != 0.0) { 
    A.coeffRef(i,j) = 0.0; // Double look-up cost! 
} 

или написать свой собственный взгляд вверх метод, который будет дублирующим бинарный поиск, используемый внутри Eigen. Что мне делать?

Answer 1

Я бы пошел писать самостоятельно, если производительность критическая (достаточно важная). В противном случае, на основе только сложности:

1. просмотрова представляет собой О (журнал (NNZ))
2. вставка (O (1) или О (NNZ), в зависимости от того, где она вставлена ​​- также оставляет больше памяти если матрица сжимается)

Таким образом, затраты на O (nnz) всегда для второго варианта или O (log (nnz)) для первого, если нет элемента. Если есть элемент, стоимость будет в два раза выше (log (nnz)), так что ваше сохранение O (log (nnz)) не будет значительным относительно вставки во втором методе. Если это действительно важно, тогда напишите свое.