[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