利用MKL求解矩阵方程出错

wenqiang0606

利用一个简单的方程验证自己用的对不对的，但是发现出错了

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

program main
  use mkl95_precision
  use mkl95_lapack
  implicit none
  integer i,info
  integer ipiv(3)
  real*8 ::a(3,3),b(3,1),X(3,1)
  a=0
  b=1
  
  forall( i=1:3)
  a(i,i)=1
  end forall
  print *,'Matrix A is:'
  do i=1,3
      print *,a(i,:)
  end do
  print *,'Matrix B is:'
  print *,b
  
  call getrf(a,ipiv,info)
  call getrs(a,ipiv,b,'N',info)
  print *,'Solution X is:'
  print *,b
end program

运行的结果，不对，不知道哪里出错了，请教一下，谢谢~

fcode

奇怪，我这里是正常的。得到的结果也是正确的

wenqiang0606

fcode 发表于 2017-12-6 17:17
奇怪，我这里是正常的。得到的结果也是正确的

难道是我的库的链接设置错了吗还是？

wenqiang0606

这个是include的位置：C:\Program Files (x86)\Intel\Composer XE 2013 SP1\mkl\include
linker的library位置：C:\Program Files (x86)\Intel\Composer XE 2013 SP1\mkl\lib\intel64
library那里改成了parallel