解决本征中的稀疏系统

Question

我正在尝试实现this paper第4.3节中描述的方法，并希望验证我对Eigen API的使用是否正确。我尝试使用不完整的Cholesky作为前置条件进行共轭梯度下降

ConjugateGradient<SparseMatrix<float>, Lower | Upper, IncompleteCholesky<f32, Lower | Upper, NaturalOrdering<int>>> solver;
solver.compute(miBlockSum);
optimalCoeffs = solver.solve(subOptimalCoeffs);

无论上下标志如何，即使更改为默认的对角前置条件，我都收到完全相同的结果，因此问题可能出在代码的其他地方。如果是这样，很高兴知道不是我使用Eigen引起了任何问题。