使用MKL库函数ZGETRS计算线性复方程组结果不对，小白求指导

warriorVIP · 发表于 2018-9-13 13:31:08

问题：采用Fortran77接口计算输出的复数解不正确，在采用Fortran95接口时，一直显示说没有这个函数
环境：计算环境完全按照这个网站的设置 https://blog.csdn.net/liuyijie1990/article/details/49923853，采用vs2013+IVF2013

下面把源代码、错误结果还有正确结果发一下，希望各位有空给看下，这是老师给安排的任务，把CVF编译器上的程序转到IVF截面，主要是把 IMSL 库函数 LSACG 用 MKL 中的库函数代替下
我选的是GETRS，如果有更好的意见可告诉我，本人工科博一小白一个

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

01      PROGRAM MAIN
02     
03      USE MKL95_PRECISION
04      USE MKL95_LAPACK
05       
06      IMPLICIT NONE
07 
08      INTEGER :: I, J, N, INFO
09       
10!      INTEGER,ALLOCATABLE:: IPIV(:)
11      COMPLEX(8), ALLOCATABLE :: A(:,:), B(:,:), IPIV(:)
12 
13 
14      N = 3
15      ALLOCATE( IPIV(N) )
16      ALLOCATE( A(N,N),  B(N,1) )
17 
18      A = 0.0
19      Open (5, File='A.TXT', Status='OLD') ! INPUT FILE
20 
21      DO J=1,N
22         DO I=1,N
23            READ(5,*) A(I,J)
24         ENDDO
25      ENDDO
26 
27      B = 0.0
28      Open (6, File='B.TXT', Status='OLD') ! INPUT FILE
29 
30      DO I=1,N
31            READ(6,*) B(I,1)
32      ENDDO
33       
34      WRITE(*,*)
35      WRITE(*,*)'Matrix A : '
36      DO I=1,N;
37         WRITE(*,*) (A(I,:));
38      ENDDO     
39 
40      WRITE(*,*)
41      WRITE(*,*)'Matrix B : '
42      DO I=1,N;
43         WRITE(*,*) (B(I,1));
44      ENDDO
45 
46 
47      CALL ZGETRS( 'N', N, 1, A, N, IPIV, B, N, INFO )
48 
49!      CALL GETRS(A, IPIV, B, 'N', INFO)
50    
51 
52      WRITE(*,*)
53      WRITE(*,*)'Matrix B : '
54      DO I=1,N;
55         WRITE(*,*) (B(I,1));
56      ENDDO
57       
58      WRITE(*,*) INFO
59 
60 
61 
62      END PROGRAM MAIN

错误结果：

正确结果

pasuka · 发表于 2018-9-13 13:54:52

lz不熟悉fortran和MKL很正常，可是怎么连正确答案也算不出来呢？

warriorVIP · 发表于 2018-9-13 14:05:01

pasuka 发表于 2018-9-13 13:54
lz不熟悉fortran和MKL很正常，可是怎么连正确答案也算不出来呢？

哈哈，我的天，您这是已经看穿了吗，指导一下嘛

百事可乐 · 发表于 2018-9-13 14:09:25

试试把 A 转置一下

vvt · 发表于 2018-9-13 15:12:07

看如下代码之注释

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

01program main
02  use f95_precision
03  use lapack95
04  implicit none
05  integer :: i, j, n, info
06  integer,allocatable:: ipiv(:)!// ipiv 为整数
07  complex(8), allocatable :: a(:,:), b(:,:)
08  n = 3
09  allocate( ipiv(n) )
10  allocate( a(n,n),  b(n,1) )
11  a = 0.0
12  open (5, file='a.txt', status='old') ! input file
13  do j=1,n
14    do i=1,n
15      read(5,*) a(j,i)!//a转置
16    enddo
17  enddo
18  !a = transpose(a)!//或者这样转置
19  b = 0.0
20  open (6, file='b.txt', status='old') ! input file
21  do i=1,n
22    read(6,*) b(i,1)
23  enddo
24  write(*,*)
25  write(*,*)'matrix a : '
26  do i=1,n;
27    write(*,*) (a(i,:));
28  enddo
29  write(*,*)
30  write(*,*)'matrix b : '
31  do i=1,n;
32    write(*,*) (b(i,1));
33  enddo
34  !call zgetrs( 'n', n, 1, a, n, ipiv, b, n, info )
35  call getrf(a,ipiv,info)!//先 getrf
36  call getrs(a,ipiv,b,'n',info) !//后 getrs
37  write(*,*)
38  write(*,*)'matrix b : '
39  do i=1,n;
40    write(*,*) (b(i,1));
41  enddo
42  write(*,*) info
43end program main

warriorVIP · 发表于 2018-9-13 15:43:24

百事可乐发表于 2018-9-13 14:09
试试把 A 转置一下

不行啊，我确实是按列读入，按行输出的

warriorVIP · 发表于 2018-9-13 16:47:47

vvt 发表于 2018-9-13 15:12
看如下代码之注释
[mw_shl_code=fortran,true]program main
use f95_precision

我试试，真的好优秀啊，谢谢

warriorVIP · 发表于 2018-9-13 18:13:24

vvt 发表于 2018-9-13 15:12
看如下代码之注释
[mw_shl_code=fortran,true]program main
use f95_precision

非常感谢啊，终于解决了，问下MKL具体的设置您能告诉我一下吗，包括这个USE什么什么的，我之前写的 USE MKL95_PRECISION，USE MKL95_LAPACK为何不能用95的接口呢

warriorVIP · 发表于 2018-9-15 10:30:36

[数学库] 使用MKL库函数ZGETRS计算线性复方程组结果不对，小白求指导

水王勋章

元老勋章

热心勋章

规矩勋章

新人勋章