SUBROUTINE CHESVX1_F95(A, B, X, UPLO, AF, IPIV, FACT, & FERR, BERR, RCOND, INFO) ! ! -- LAPACK95 interface driver routine (version 3.0) -- ! UNI-C, Denmark; Univ. of Tennessee, USA; NAG Ltd., UK ! September, 2000 ! ! .. USE STATEMENTS .. USE LA_PRECISION, ONLY: WP => SP USE LA_AUXMOD, ONLY: LSAME, ERINFO USE F77_LAPACK, ONLY: HESVX_F77 => LA_HESVX, ILAENV_F77 => ILAENV ! .. IMPLICIT STATEMENT .. IMPLICIT NONE ! .. SCALAR ARGUMENTS .. CHARACTER(LEN=1), INTENT(IN), OPTIONAL :: UPLO, FACT INTEGER, INTENT(OUT), OPTIONAL :: INFO REAL(WP), INTENT(OUT), OPTIONAL :: RCOND, FERR, BERR ! .. ARRAY ARGUMENTS .. COMPLEX(WP), INTENT(IN) :: A(:,:), B(:) COMPLEX(WP), INTENT(OUT) :: X(:) INTEGER, INTENT(INOUT), OPTIONAL, TARGET :: IPIV(:) COMPLEX(WP), INTENT(INOUT), OPTIONAL, TARGET :: AF(:,:) ! .. PARAMETERS .. CHARACTER(LEN=8), PARAMETER :: SRNAME = 'LA_HESVX' CHARACTER(LEN=6), PARAMETER :: BSNAME = 'CHETRF' ! .. LOCAL SCALARS .. CHARACTER(LEN=1) :: LFACT, LUPLO INTEGER :: LINFO, N, NB, LWORK, ISTAT, ISTAT1, SIPIV, S1AF, S2AF REAL(WP) :: LRCOND, LFERR, LBERR ! .. LOCAL POINTERS .. INTEGER, POINTER :: LPIV(:) REAL(WP), POINTER :: RWORK(:) COMPLEX(WP), POINTER :: WORK(:), LAF(:, :) ! .. INTRINSIC FUNCTIONS .. INTRINSIC PRESENT, SIZE, MAX ! .. EXECUTABLE STATEMENTS .. LINFO = 0; ISTAT = 0; N = SIZE(A, 1) IF( PRESENT(RCOND) ) RCOND = 1.0_WP IF( PRESENT(FACT) )THEN; LFACT = FACT; ELSE; LFACT='N'; END IF IF( PRESENT(UPLO) ) THEN; LUPLO = UPLO; ELSE; LUPLO = 'U'; END IF IF( PRESENT(IPIV) )THEN; SIPIV = SIZE(IPIV); ELSE; SIPIV = N; END IF IF( PRESENT(AF) )THEN; S1AF = SIZE(AF,1); S2AF = SIZE(AF,2) ELSE; S1AF = N; S2AF = N; END IF ! .. TEST THE ARGUMENTS IF( SIZE(A, 2) /= N .OR. N < 0 )THEN; LINFO = -1 ELSE IF( SIZE(B) /= N )THEN; LINFO = -2 ELSE IF( SIZE(X) /= N )THEN; LINFO = -3 ELSE IF( .NOT.LSAME(LUPLO,'U') .AND. .NOT.LSAME(LUPLO,'L') )THEN; LINFO = -4 ELSE IF( S1AF /= N .OR. S2AF /= N ) THEN; LINFO = -5 ELSE IF( SIPIV /= N )THEN; LINFO = -6 ELSE IF( ( .NOT. LSAME(LFACT,'F') .AND. .NOT. LSAME(LFACT,'N') ) .OR. & ( LSAME(LFACT,'F') .AND. .NOT.( PRESENT(AF) .AND. PRESENT(IPIV) ) ) )THEN; LINFO = -7 ELSE IF ( N > 0 )THEN IF( .NOT.PRESENT(AF) ) THEN; ALLOCATE( LAF(N,N), STAT=ISTAT ) ELSE; LAF => AF; END IF IF( ISTAT == 0 )THEN IF( .NOT.PRESENT(IPIV) )THEN; ALLOCATE( LPIV(N), STAT=ISTAT ) ELSE; LPIV => IPIV; END IF END IF IF( ISTAT == 0 )THEN NB = ILAENV_F77( 1, BSNAME, LUPLO, N, -1, -1, -1 ) IF( NB <= 1 .OR. NB >= N ) NB = 1; LWORK = MAX(1,2*N,N*NB) ALLOCATE(WORK(LWORK), RWORK(N), STAT=ISTAT) IF( ISTAT /= 0 )THEN DEALLOCATE(WORK, RWORK, STAT=ISTAT1); LWORK = MAX(1,2*N) ALLOCATE(WORK(LWORK), RWORK(N), STAT=ISTAT) IF( ISTAT /= 0 ) THEN; LINFO = - 100 ELSE; CALL ERINFO( -200, SRNAME, LINFO ); ENDIF ENDIF END IF IF( ISTAT == 0 )THEN ! .. CALL LAPACK77 ROUTINE CALL HESVX_F77( LFACT, LUPLO, N, 1, A, N, LAF, N, LPIV, B, N, X, N, & LRCOND, LFERR, LBERR, WORK, LWORK, RWORK, LINFO ) ELSE; LINFO = -100; END IF IF( .NOT.PRESENT(AF) ) DEALLOCATE( LAF, STAT=ISTAT1 ) IF( .NOT.PRESENT(IPIV) ) DEALLOCATE( LPIV, STAT=ISTAT1 ) IF( PRESENT(FERR) ) FERR = LFERR IF( PRESENT(BERR) ) BERR = LBERR IF( PRESENT(RCOND) ) RCOND=LRCOND DEALLOCATE( WORK, RWORK, STAT=ISTAT1 ) END IF CALL ERINFO( LINFO, SRNAME, INFO, ISTAT ) END SUBROUTINE CHESVX1_F95