Fortran Coder

查看: 5797|回复: 1
打印 上一主题 下一主题

[数学库] 调用MKL的dgesv函数时出现错误

[复制链接]

21

帖子

4

主题

0

精华

入门

F 币
87 元
贡献
50 点
跳转到指定楼层
楼主
发表于 2018-5-9 19:44:06 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
请问大佬们这个代码哪里有问题。
info运行结果等于4
谢谢!

[Fortran] 纯文本查看 复制代码
program main
use lapack95
    implicit none

integer,parameter::size=5
integer::i,j
integer::n,lda,ldb
integer::nrhs
integer::info
integer::ipiv(size)
real::a(size,size),b(size,size)
n=size;lda=size;ldb=size
nrhs=size


forall(i=1:size,j=1:size,i>j)     a(i,j)=2
forall(i=1:size,j=1:size,i==j)    a(i,j)=3
forall(i=1:size,j=1:size,i<j)     a(i,j)=4

forall(i=1:size,j=1:size,i>j)     b(i,j)=0
forall(i=1:size,j=1:size,i==j)    b(i,j)=1
forall(i=1:size,j=1:size,i<j)     b(i,j)=0

call dgesv(n,nrhs,a,lda,ipiv,b,ldb,info)
write(*,*)'solution:'
write(*,*)b
write(*,*)'info=',info

end program
分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

21

帖子

4

主题

0

精华

入门

F 币
87 元
贡献
50 点
沙发
 楼主| 发表于 2018-5-9 20:03:52 | 只看该作者
已经解决了。
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-11-15 06:50

Powered by Tencent X3.4

© 2013-2024 Tencent

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