我正在尝试使用 Eigen Library SparseQR 线性求解器(矩阵不是正方形而不保证是PSD)来解决Ax = b问题。麻烦的是,即使我使用双精度条目并确保A矩阵中的条目大致相同的顺序,我也没有得到像 Matlab 那样准确的答案。同样的问题。我检查了
给出的残差r = norm(Ax-b)/norm(b)
最终成为0.00174
。我也检查了排名,矩阵是满级。
有没有办法为sparseQR设置容差或残留阈值?或者也许是减少错误的更好方法?
谢谢!
亚洲时报Siddharth