MKL中getrf省略参数ipiv的问题

fcode

[Fortran] 纯文本查看 复制代码

call getrf(a,ipiv,info)!矩阵分解成LU
!call getrf(a)!矩阵分解成LU
call getri(a,ipiv,info)!求逆阵

我这里两行代码计算结果是一样的。
  0.1764706     -5.8823526E-02  0.2941177
  0.2941176      0.2352941     -0.1764706
-0.2352941      0.4117647     -5.8823526E-02
你再检查一下看是不是别的问题吧

fcode

getrf 函数进行LU分解，使用了部分选主元。也就是会进行部分的行变换。变换结果会返回给 ipiv，在 getri 中使用。

如果你不返回，那么 getri 中的就是错误的。自然结果就不对了。（其实只是进行了行变换（交换）而已）

[Fortran] 纯文本查看 复制代码

ipiv=3
ipiv(1)=2 !// 加上这行就对了
call getrf(a)
call getri(a,ipiv=ipiv,info=info)

从上面的代码可以看出，ipiv不一定非要从getrf里返回。如果你事先知道，你自己生成也可以。所以与 getrf 无关。

fcode

是的，就是如此。如果你为了求逆，那么最好保留。如果只是做LU分解，则可以不必保留。

PS：对管理员和版主，你不必评分，会从你的分数中扣除。管理员和版主不缺分。