Приближенное решение линейной системы в библиотеке Armadillo

Question

Кто-нибудь знает, как библиотека Armadillo находит приближенное решение для линейной системы с необратимой матрицей A (Ax = b). Рассматриваемая функция: решить (х, А, б) .

Заранее спасибо

Answer 1

При проверке файла include\armadillo_bits\glue_solve_meat.hpp кажется, что для неквадратной матрицы Армадилло использует метод LAPACK gelsd (решение наименьших квадратов на основе SVD). Если указан параметр настройки solve_opts::fast, сначала он пытается gels (разложение QR / LQ), а если это не удается, он также обращается к SVD.