Fortran Coder

标题: 利用MKL求解矩阵方程出错 [打印本页]

作者: wenqiang0606 时间: 2017-12-6 16:32
标题: 利用MKL求解矩阵方程出错
利用一个简单的方程验证自己用的对不对的，但是发现出错了

[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

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

捕获.PNG (7.34 KB, 下载次数: 222)

作者: fcode 时间: 2017-12-6 17:17
奇怪，我这里是正常的。得到的结果也是正确的

QQ截图20171206171653.png (12.8 KB, 下载次数: 214)

作者: wenqiang0606 时间: 2017-12-6 18:23

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

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

作者: wenqiang0606 时间: 2017-12-6 18:29
这个是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

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)