MKL函数库中ZGETRF和ZGETRI函数的使用

li913 · 发表于 2019-11-10 14:36:24

本帖最后由 li913 于 2019-11-10 14:38 编辑

放宽心，代码没问题，是误差问题。A不是对角占优的，导致逆矩阵B某些值特别大(15次方)，B的误差可能在0.1-1之间，最终导致C非单位阵。你试试让对角元素大一些，比如 if(i==j) a(i,j)=2*a(i,j)，结果就ok。
另外，推荐使用f95接口，更加简单好使。CALL GETRF(b,lu_ipiv,lu_info)
call getri( b, lu_ipiv,in_info)

li913 · 发表于 2019-11-12 13:27:15

1、mkl库自带f95接口，无需额外安装；
2、徐士良程序集，你自己写代码，先不说正确性，效率肯定不如mkl。仅可作为验证手段。

li913 · 发表于 2019-11-14 14:55:46

没有详细研究过。

[通用算法] MKL函数库中ZGETRF和ZGETRI函数的使用