Armadillo库中线性系统的近似解

Question

有谁知道 Armadillo 库如何找到具有不可逆 A 矩阵 (Ax=b) 的线性系统的近似解。正在考虑的函数是 solve(x,A,b) .

提前致谢

Answer 1

通过检查文件 include\armadillo_bits\glue_solve_meat.hpp，Armadillo 似乎对非方阵使用了 LAPACK 方法 gelsd (least-squares solution based on SVD). If the settings option solve_opts::fast is supplied, it tries first gels（QR/LQ 分解），如果失败，它也求助于 SVD。