利用MKL求解矩阵方程出错

wenqiang0606 · 发表于 2017-12-6 16:32:54

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

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

01program main
02  use mkl95_precision
03  use mkl95_lapack
04  implicit none
05  integer i,info
06  integer ipiv(3)
07  real*8 ::a(3,3),b(3,1),X(3,1)
08  a=0
09  b=1
10   
11  forall( i=1:3)
12  a(i,i)=1
13  end forall
14  print *,'Matrix A is:'
15  do i=1,3
16      print *,a(i,:)
17  end do
18  print *,'Matrix B is:'
19  print *,b
20   
21  call getrf(a,ipiv,info)
22  call getrs(a,ipiv,b,'N',info)
23  print *,'Solution X is:'
24  print *,b
25end program

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

fcode · 发表于 2017-12-6 17:17:40

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

wenqiang0606 · 发表于 2017-12-6 18:23:37

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

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

wenqiang0606 · 发表于 2017-12-6 18:29:49

这个是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

[数学库] 利用MKL求解矩阵方程出错

美女勋章

热心勋章

星光勋章

新人勋章

贡献勋章

管理勋章

帅哥勋章

爱心勋章

规矩勋章

元老勋章

水王勋章