我想解决以下步骤:
填充 boost::numeric::ublas::compressed_matrix;
现在,我需要应用 Cholesky 分解。
- 但是,没有这种boost的功能。所以我开始寻找另一个库——我在考虑 BLAS 或一种 LAPACK 库。但是,虽然我需要一些算法,但如何将 boost::numeric::ublas::compressed_matrix 与 BLAS 或 LAPACK 的函数/算法绑定?有没有这样的方法?
最后,我需要解决“Ax=b”,其中“A”是由 Cholesky 算法分解的 boost 的压缩矩阵。那么如何使用 boost 和/或 LAPACK 的算法或函数来解决“x”?
提前致谢。
葡聚糖