我只是想知道，在计算的情况下B=inv(A)，假设I是单位矩阵（对角线），
A_factorization通过计算得到因式分解后getrf(A)，进行以下求解步骤在数值上是否相同？
getrs(A_factorization, piv, I)
getri(A_factorization, piv)
lapack getri 在数值上与单位矩阵为 RHS 的 getrs 相同吗？
计算科学
数值分析
线性求解器
拉帕克
2021-12-19 09:59:02
1个回答
基本上我的LAPACK文档说明了getri
此方法反转 U，然后通过求解系统 inv(A)*L = inv(U) 来计算 inv(A) 以获得 inv(A)。
似乎数值上与求解L*inv(L)=I然后是不同的过程U*inv(A)=inv(L)。我的理解是getri应该比getrs. 否则就没有理由创建它，除了可能节省一些内存工作空间。
当然，两种情况下的结果是相同的。
其它你可能感兴趣的问题