1

我想解决以下步骤:

  • 填充 boost::numeric::ublas::compressed_matrix;

  • 现在,我需要应用 Cholesky 分解。

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

提前致谢。

葡聚糖

4

0 回答 0