Fortran Coder

查看: 52|回复: 3

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

[复制链接]

3

帖子

1

主题

0

精华

新人

F 币
17 元
贡献
7 点
发表于 2017-12-6 16:32:54 | 显示全部楼层 |阅读模式
利用一个简单的方程验证自己用的对不对的,但是发现出错了
[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
回复

使用道具 举报

1148

帖子

12

主题

5

精华

论坛跑堂

Fcode跑堂

F 币
1128 元
贡献
892 点

新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

发表于 2017-12-6 17:17:40 | 显示全部楼层
奇怪,我这里是正常的。得到的结果也是正确的

QQ截图20171206171653.png

3

帖子

1

主题

0

精华

新人

F 币
17 元
贡献
7 点
 楼主| 发表于 2017-12-6 18:23:37 | 显示全部楼层
fcode 发表于 2017-12-6 17:17
奇怪,我这里是正常的。得到的结果也是正确的

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

3

帖子

1

主题

0

精华

新人

F 币
17 元
贡献
7 点
 楼主| 发表于 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
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

QQ|捐赠本站|Archiver|关于我们 About Us|QQ群|Fcode

GMT+8, 2017-12-18 15:06

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表