我正在编写一个向后消除算法。在每次迭代中,我需要从 SparseMatrix 的列中消除一些系数并更新其他非零系数。
但是,将系数的引用更改为零并不会解除分配它,因此非零系数的数量是相同的。如何删除引用?我尝试使用 makeCompressed() 无济于事,并且编译器不知道 pruned 。
基本代码如下。
我怎么解决这个问题?
#include <Eigen/SparseCore>
void nukeit(){
Eigen::SparseMatrix<double> A(4, 3);
cout << "non zeros of empty: " << A.nonZeros() << "\n" << endl;
A.insert(0, 0) = 1;
A.insert(2, 1) = 5;
cout << "non zeros are two: " << A.nonZeros() << "\n" << endl;
A.coeffRef(0, 0) = 0;
cout << "non zeros should be one but it's 2: " << A.nonZeros() << "\n" << endl;
cout << "However the matrix has only one non zero element\n" << A << endl;
}
输出
non zeros of empty: 0
non zeros are two: 2
non zeros should be one but it's 2: 2
However the matrix has only one non zero element
0 0 0
0 0 0
0 5 0
0 0 0
最佳答案
将当前列的一些系数设置为零后,您可以通过调用 A.prune(0.0)
显式删除它们。请参阅相应的 doc 。
但是,请注意,这将触发剩余列条目的昂贵内存拷贝。对于稀疏矩阵,我们通常不会就地工作。
关于c++ - 特征:如何从稀疏矩阵中删除初始化系数,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/33474981/