犰狳图书馆线性系统的近似解

Question

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

提前致谢

Answer 1

通过检查文件include\armadillo_bits\glue_solve_meat.hpp，似乎对于非方矩阵，Armadillo使用LAPACK方法gelsd（基于SVD的最小二乘解）。如果提供了设置选项solve_opts::fast，它会首先尝试gels（QR / LQ分解），如果失败，它也会转到SVD。