我浏览了 LAPACK 用户指南,但找不到 LAPACK 是否为基于 Krylov 子空间的方法(例如 CG 或 BiCGSTAB 等)和基于牛顿方法的非线性求解器提供例程。此外,LAPACK 中是否提供迭代方法?
LAPACK 是否为基于 Krylov 子空间的求解器和非线性求解器提供例程?
计算科学
拉帕克
克雷洛夫法
2021-12-07 17:48:45
2个回答
据我所知,LAPACK 中没有这样的方法。由于 LAPACK 是线性代数包,因此不包括非线性求解器。
但是,您可以使用底层 BLAS 来实现迭代方法。对于使用雅可比矩阵的非线性求解器(例如,牛顿法),LAPACK 的矩阵分解可能会派上用场。
不过,您可能想看看PETSc。该库实现了许多迭代方法(例如,CG、GMRES、BiCG)和非线性求解器(例如,Quasi-Newton、信任域 Newton、Richardson / Picard 迭代、Anderson 混合)。
LAPACK 不包括任何迭代求解器。LAPACK 中的例程用于特征值、矩阵分解和涉及密集矩阵的方程组的解,而迭代方法通常用于大型和稀疏矩阵。
其它你可能感兴趣的问题