Boost 的 y=Ax 线性代数解

boost有吗？ 其中 A、y 和 x 分别是矩阵（稀疏且可以很大）和向量。 y 或 x 可能是未知的。

我在这里似乎找不到它：http://www.boost.org/doc/libs/1_39_0/libs/numeric/ublas/doc/index.htm

是的，您可以使用 boost 的 ublas 库求解线性方程。这是使用 LU 因式分解和回代来获得逆矩阵的一种简短方法：

using namespace boost::ublas;

Ainv = identity_matrix<float>(A.size1());
permutation_matrix<size_t> pm(A.size1());
lu_factorize(A,pm)
lu_substitute(A, pm, Ainv);

因此，要求解线性系统 Ax=y，您可以通过取 (trans(A)A)^-1 的倒数来求解方程 trans(A)Ax=trans(A)y 以获得 x： x = (trans( A)A)^-1Ay。