Fortran Coder

查看: 85|回复: 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 | 显示全部楼层
已经解决了。
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

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

GMT+8, 2018-7-22 12:55

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

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