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的接口呢

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

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

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

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

看如下代码之注释

[Fortran] 查看源码 复制源码

      PROGRAM MAIN
    
      USE MKL95_PRECISION
      USE MKL95_LAPACK
      
      IMPLICIT NONE

      INTEGER :: I, J, N, INFO
      
!      INTEGER,ALLOCATABLE:: IPIV(:)
      COMPLEX(8), ALLOCATABLE :: A(:,:), B(:,:), IPIV(:)


      N = 3
      ALLOCATE( IPIV(N) )
      ALLOCATE( A(N,N),  B(N,1) )

      A = 0.0
      Open (5, File='A.TXT', Status='OLD') ! INPUT FILE

      DO J=1,N
         DO I=1,N
            READ(5,*) A(I,J)
         ENDDO
      ENDDO

      B = 0.0
      Open (6, File='B.TXT', Status='OLD') ! INPUT FILE

      DO I=1,N
            READ(6,*) B(I,1)
      ENDDO
      
      WRITE(*,*)
      WRITE(*,*)'Matrix A : '
      DO I=1,N;
         WRITE(*,*) (A(I,:));
      ENDDO      

      WRITE(*,*)
      WRITE(*,*)'Matrix B : '
      DO I=1,N;
         WRITE(*,*) (B(I,1));
      ENDDO


      CALL ZGETRS( 'N', N, 1, A, N, IPIV, B, N, INFO )

!      CALL GETRS(A, IPIV, B, 'N', INFO)
   

      WRITE(*,*)
      WRITE(*,*)'Matrix B : '
      DO I=1,N;
         WRITE(*,*) (B(I,1));
      ENDDO
      
      WRITE(*,*) INFO



      END PROGRAM MAIN    

试试把 A 转置一下

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

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

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