调用MKL的dgesv函数时出现错误

Yowai · 发表于 2018-5-9 19:44:06

请问大佬们这个代码哪里有问题。
info运行结果等于4
谢谢！

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

01program main
02use lapack95
03    implicit none
04 
05integer,parameter::size=5
06integer::i,j
07integer::n,lda,ldb
08integer::nrhs
09integer::info
10integer::ipiv(size)
11real::a(size,size),b(size,size)
12n=size;lda=size;ldb=size
13nrhs=size
14 
15 
16forall(i=1:size,j=1:size,i>j)     a(i,j)=2
17forall(i=1:size,j=1:size,i==j)    a(i,j)=3
18forall(i=1:size,j=1:size,i<j)     a(i,j)=4
19 
20forall(i=1:size,j=1:size,i>j)     b(i,j)=0
21forall(i=1:size,j=1:size,i==j)    b(i,j)=1
22forall(i=1:size,j=1:size,i<j)     b(i,j)=0
23 
24call dgesv(n,nrhs,a,lda,ipiv,b,ldb,info)
25write(*,*)'solution:'
26write(*,*)b
27write(*,*)'info=',info
28 
29end program

Yowai · 发表于 2018-5-9 20:03:52

已经解决了。

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