DOUBLE PRECISION FUNCTION DCABS1(Z)
* .. Scalar Arguments ..
- DOUBLE COMPLEX Z
+ COMPLEX*16 Z
* ..
* ..
* Purpose
INTEGER INCX,N
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX ZX(*)
+ COMPLEX*16 ZX(*)
* ..
*
* Purpose
INTEGER INCX,N
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX X(*)
+ COMPLEX*16 X(*)
* ..
*
* Purpose
INTEGER INCX,N
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX ZX(*)
+ COMPLEX*16 ZX(*)
* ..
*
* Purpose
SPARAM(5) = SH22
END IF
- 260 CONTINUE
SPARAM(1) = SFLAG
RETURN
END
SUBROUTINE ZAXPY(N,ZA,ZX,INCX,ZY,INCY)
* .. Scalar Arguments ..
- DOUBLE COMPLEX ZA
+ COMPLEX*16 ZA
INTEGER INCX,INCY,N
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX ZX(*),ZY(*)
+ COMPLEX*16 ZX(*),ZY(*)
* ..
*
* Purpose
INTEGER INCX,INCY,N
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX ZX(*),ZY(*)
+ COMPLEX*16 ZX(*),ZY(*)
* ..
*
* Purpose
- DOUBLE COMPLEX FUNCTION ZDOTC(N,ZX,INCX,ZY,INCY)
+ COMPLEX*16 FUNCTION ZDOTC(N,ZX,INCX,ZY,INCY)
* .. Scalar Arguments ..
INTEGER INCX,INCY,N
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX ZX(*),ZY(*)
+ COMPLEX*16 ZX(*),ZY(*)
* ..
*
* Purpose
* =====================================================================
*
* .. Local Scalars ..
- DOUBLE COMPLEX ZTEMP
+ COMPLEX*16 ZTEMP
INTEGER I,IX,IY
* ..
* .. Intrinsic Functions ..
- DOUBLE COMPLEX FUNCTION ZDOTU(N,ZX,INCX,ZY,INCY)
+ COMPLEX*16 FUNCTION ZDOTU(N,ZX,INCX,ZY,INCY)
* .. Scalar Arguments ..
INTEGER INCX,INCY,N
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX ZX(*),ZY(*)
+ COMPLEX*16 ZX(*),ZY(*)
* ..
*
* Purpose
* =====================================================================
*
* .. Local Scalars ..
- DOUBLE COMPLEX ZTEMP
+ COMPLEX*16 ZTEMP
INTEGER I,IX,IY
* ..
ZTEMP = (0.0d0,0.0d0)
INTEGER INCX,N
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX ZX(*)
+ COMPLEX*16 ZX(*)
* ..
*
* Purpose
SUBROUTINE ZGBMV(TRANS,M,N,KL,KU,ALPHA,A,LDA,X,INCX,BETA,Y,INCY)
* .. Scalar Arguments ..
- DOUBLE COMPLEX ALPHA,BETA
+ COMPLEX*16 ALPHA,BETA
INTEGER INCX,INCY,KL,KU,LDA,M,N
CHARACTER TRANS
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),X(*),Y(*)
+ COMPLEX*16 A(LDA,*),X(*),Y(*)
* ..
*
* Purpose
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ONE
+ COMPLEX*16 ONE
PARAMETER (ONE= (1.0D+0,0.0D+0))
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP
+ COMPLEX*16 TEMP
INTEGER I,INFO,IX,IY,J,JX,JY,K,KUP1,KX,KY,LENX,LENY
LOGICAL NOCONJ
* ..
SUBROUTINE ZGEMM(TRANSA,TRANSB,M,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
* .. Scalar Arguments ..
- DOUBLE COMPLEX ALPHA,BETA
+ COMPLEX*16 ALPHA,BETA
INTEGER K,LDA,LDB,LDC,M,N
CHARACTER TRANSA,TRANSB
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),B(LDB,*),C(LDC,*)
+ COMPLEX*16 A(LDA,*),B(LDB,*),C(LDC,*)
* ..
*
* Purpose
INTRINSIC DCONJG,MAX
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP
+ COMPLEX*16 TEMP
INTEGER I,INFO,J,L,NCOLA,NROWA,NROWB
LOGICAL CONJA,CONJB,NOTA,NOTB
* ..
* .. Parameters ..
- DOUBLE COMPLEX ONE
+ COMPLEX*16 ONE
PARAMETER (ONE= (1.0D+0,0.0D+0))
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
*
SUBROUTINE ZGEMV(TRANS,M,N,ALPHA,A,LDA,X,INCX,BETA,Y,INCY)
* .. Scalar Arguments ..
- DOUBLE COMPLEX ALPHA,BETA
+ COMPLEX*16 ALPHA,BETA
INTEGER INCX,INCY,LDA,M,N
CHARACTER TRANS
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),X(*),Y(*)
+ COMPLEX*16 A(LDA,*),X(*),Y(*)
* ..
*
* Purpose
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ONE
+ COMPLEX*16 ONE
PARAMETER (ONE= (1.0D+0,0.0D+0))
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP
+ COMPLEX*16 TEMP
INTEGER I,INFO,IX,IY,J,JX,JY,KX,KY,LENX,LENY
LOGICAL NOCONJ
* ..
SUBROUTINE ZGERC(M,N,ALPHA,X,INCX,Y,INCY,A,LDA)
* .. Scalar Arguments ..
- DOUBLE COMPLEX ALPHA
+ COMPLEX*16 ALPHA
INTEGER INCX,INCY,LDA,M,N
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),X(*),Y(*)
+ COMPLEX*16 A(LDA,*),X(*),Y(*)
* ..
*
* Purpose
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP
+ COMPLEX*16 TEMP
INTEGER I,INFO,IX,J,JY,KX
* ..
* .. External Subroutines ..
SUBROUTINE ZGERU(M,N,ALPHA,X,INCX,Y,INCY,A,LDA)
* .. Scalar Arguments ..
- DOUBLE COMPLEX ALPHA
+ COMPLEX*16 ALPHA
INTEGER INCX,INCY,LDA,M,N
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),X(*),Y(*)
+ COMPLEX*16 A(LDA,*),X(*),Y(*)
* ..
*
* Purpose
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP
+ COMPLEX*16 TEMP
INTEGER I,INFO,IX,J,JY,KX
* ..
* .. External Subroutines ..
SUBROUTINE ZHBMV(UPLO,N,K,ALPHA,A,LDA,X,INCX,BETA,Y,INCY)
* .. Scalar Arguments ..
- DOUBLE COMPLEX ALPHA,BETA
+ COMPLEX*16 ALPHA,BETA
INTEGER INCX,INCY,K,LDA,N
CHARACTER UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),X(*),Y(*)
+ COMPLEX*16 A(LDA,*),X(*),Y(*)
* ..
*
* Purpose
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ONE
+ COMPLEX*16 ONE
PARAMETER (ONE= (1.0D+0,0.0D+0))
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP1,TEMP2
+ COMPLEX*16 TEMP1,TEMP2
INTEGER I,INFO,IX,IY,J,JX,JY,KPLUS1,KX,KY,L
* ..
* .. External Functions ..
SUBROUTINE ZHEMM(SIDE,UPLO,M,N,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
* .. Scalar Arguments ..
- DOUBLE COMPLEX ALPHA,BETA
+ COMPLEX*16 ALPHA,BETA
INTEGER LDA,LDB,LDC,M,N
CHARACTER SIDE,UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),B(LDB,*),C(LDC,*)
+ COMPLEX*16 A(LDA,*),B(LDB,*),C(LDC,*)
* ..
*
* Purpose
INTRINSIC DBLE,DCONJG,MAX
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP1,TEMP2
+ COMPLEX*16 TEMP1,TEMP2
INTEGER I,INFO,J,K,NROWA
LOGICAL UPPER
* ..
* .. Parameters ..
- DOUBLE COMPLEX ONE
+ COMPLEX*16 ONE
PARAMETER (ONE= (1.0D+0,0.0D+0))
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
*
SUBROUTINE ZHEMV(UPLO,N,ALPHA,A,LDA,X,INCX,BETA,Y,INCY)
* .. Scalar Arguments ..
- DOUBLE COMPLEX ALPHA,BETA
+ COMPLEX*16 ALPHA,BETA
INTEGER INCX,INCY,LDA,N
CHARACTER UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),X(*),Y(*)
+ COMPLEX*16 A(LDA,*),X(*),Y(*)
* ..
*
* Purpose
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ONE
+ COMPLEX*16 ONE
PARAMETER (ONE= (1.0D+0,0.0D+0))
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP1,TEMP2
+ COMPLEX*16 TEMP1,TEMP2
INTEGER I,INFO,IX,IY,J,JX,JY,KX,KY
* ..
* .. External Functions ..
CHARACTER UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),X(*)
+ COMPLEX*16 A(LDA,*),X(*)
* ..
*
* Purpose
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP
+ COMPLEX*16 TEMP
INTEGER I,INFO,IX,J,JX,KX
* ..
* .. External Functions ..
SUBROUTINE ZHER2(UPLO,N,ALPHA,X,INCX,Y,INCY,A,LDA)
* .. Scalar Arguments ..
- DOUBLE COMPLEX ALPHA
+ COMPLEX*16 ALPHA
INTEGER INCX,INCY,LDA,N
CHARACTER UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),X(*),Y(*)
+ COMPLEX*16 A(LDA,*),X(*),Y(*)
* ..
*
* Purpose
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP1,TEMP2
+ COMPLEX*16 TEMP1,TEMP2
INTEGER I,INFO,IX,IY,J,JX,JY,KX,KY
* ..
* .. External Functions ..
SUBROUTINE ZHER2K(UPLO,TRANS,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
* .. Scalar Arguments ..
- DOUBLE COMPLEX ALPHA
+ COMPLEX*16 ALPHA
DOUBLE PRECISION BETA
INTEGER K,LDA,LDB,LDC,N
CHARACTER TRANS,UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),B(LDB,*),C(LDC,*)
+ COMPLEX*16 A(LDA,*),B(LDB,*),C(LDC,*)
* ..
*
* Purpose
INTRINSIC DBLE,DCONJG,MAX
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP1,TEMP2
+ COMPLEX*16 TEMP1,TEMP2
INTEGER I,INFO,J,L,NROWA
LOGICAL UPPER
* ..
* .. Parameters ..
DOUBLE PRECISION ONE
PARAMETER (ONE=1.0D+0)
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
*
CHARACTER TRANS,UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),C(LDC,*)
+ COMPLEX*16 A(LDA,*),C(LDC,*)
* ..
*
* Purpose
INTRINSIC DBLE,DCMPLX,DCONJG,MAX
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP
+ COMPLEX*16 TEMP
DOUBLE PRECISION RTEMP
INTEGER I,INFO,J,L,NROWA
LOGICAL UPPER
SUBROUTINE ZHPMV(UPLO,N,ALPHA,AP,X,INCX,BETA,Y,INCY)
* .. Scalar Arguments ..
- DOUBLE COMPLEX ALPHA,BETA
+ COMPLEX*16 ALPHA,BETA
INTEGER INCX,INCY,N
CHARACTER UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX AP(*),X(*),Y(*)
+ COMPLEX*16 AP(*),X(*),Y(*)
* ..
*
* Purpose
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ONE
+ COMPLEX*16 ONE
PARAMETER (ONE= (1.0D+0,0.0D+0))
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP1,TEMP2
+ COMPLEX*16 TEMP1,TEMP2
INTEGER I,INFO,IX,IY,J,JX,JY,K,KK,KX,KY
* ..
* .. External Functions ..
CHARACTER UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX AP(*),X(*)
+ COMPLEX*16 AP(*),X(*)
* ..
*
* Purpose
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP
+ COMPLEX*16 TEMP
INTEGER I,INFO,IX,J,JX,K,KK,KX
* ..
* .. External Functions ..
SUBROUTINE ZHPR2(UPLO,N,ALPHA,X,INCX,Y,INCY,AP)
* .. Scalar Arguments ..
- DOUBLE COMPLEX ALPHA
+ COMPLEX*16 ALPHA
INTEGER INCX,INCY,N
CHARACTER UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX AP(*),X(*),Y(*)
+ COMPLEX*16 AP(*),X(*),Y(*)
* ..
*
* Purpose
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP1,TEMP2
+ COMPLEX*16 TEMP1,TEMP2
INTEGER I,INFO,IX,IY,J,JX,JY,K,KK,KX,KY
* ..
* .. External Functions ..
SUBROUTINE ZROTG(CA,CB,C,S)
* .. Scalar Arguments ..
- DOUBLE COMPLEX CA,CB,S
+ COMPLEX*16 CA,CB,S
DOUBLE PRECISION C
* ..
*
* =====================================================================
*
* .. Local Scalars ..
- DOUBLE COMPLEX ALPHA
+ COMPLEX*16 ALPHA
DOUBLE PRECISION NORM,SCALE
* ..
* .. Intrinsic Functions ..
SUBROUTINE ZSCAL(N,ZA,ZX,INCX)
* .. Scalar Arguments ..
- DOUBLE COMPLEX ZA
+ COMPLEX*16 ZA
INTEGER INCX,N
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX ZX(*)
+ COMPLEX*16 ZX(*)
* ..
*
* Purpose
INTEGER INCX,INCY,N
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX ZX(*),ZY(*)
+ COMPLEX*16 ZX(*),ZY(*)
* ..
*
* Purpose
* =====================================================================
*
* .. Local Scalars ..
- DOUBLE COMPLEX ZTEMP
+ COMPLEX*16 ZTEMP
INTEGER I,IX,IY
* ..
IF (N.LE.0) RETURN
SUBROUTINE ZSYMM(SIDE,UPLO,M,N,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
* .. Scalar Arguments ..
- DOUBLE COMPLEX ALPHA,BETA
+ COMPLEX*16 ALPHA,BETA
INTEGER LDA,LDB,LDC,M,N
CHARACTER SIDE,UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),B(LDB,*),C(LDC,*)
+ COMPLEX*16 A(LDA,*),B(LDB,*),C(LDC,*)
* ..
*
* Purpose
INTRINSIC MAX
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP1,TEMP2
+ COMPLEX*16 TEMP1,TEMP2
INTEGER I,INFO,J,K,NROWA
LOGICAL UPPER
* ..
* .. Parameters ..
- DOUBLE COMPLEX ONE
+ COMPLEX*16 ONE
PARAMETER (ONE= (1.0D+0,0.0D+0))
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
*
SUBROUTINE ZSYR2K(UPLO,TRANS,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
* .. Scalar Arguments ..
- DOUBLE COMPLEX ALPHA,BETA
+ COMPLEX*16 ALPHA,BETA
INTEGER K,LDA,LDB,LDC,N
CHARACTER TRANS,UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),B(LDB,*),C(LDC,*)
+ COMPLEX*16 A(LDA,*),B(LDB,*),C(LDC,*)
* ..
*
* Purpose
INTRINSIC MAX
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP1,TEMP2
+ COMPLEX*16 TEMP1,TEMP2
INTEGER I,INFO,J,L,NROWA
LOGICAL UPPER
* ..
* .. Parameters ..
- DOUBLE COMPLEX ONE
+ COMPLEX*16 ONE
PARAMETER (ONE= (1.0D+0,0.0D+0))
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
*
SUBROUTINE ZSYRK(UPLO,TRANS,N,K,ALPHA,A,LDA,BETA,C,LDC)
* .. Scalar Arguments ..
- DOUBLE COMPLEX ALPHA,BETA
+ COMPLEX*16 ALPHA,BETA
INTEGER K,LDA,LDC,N
CHARACTER TRANS,UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),C(LDC,*)
+ COMPLEX*16 A(LDA,*),C(LDC,*)
* ..
*
* Purpose
INTRINSIC MAX
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP
+ COMPLEX*16 TEMP
INTEGER I,INFO,J,L,NROWA
LOGICAL UPPER
* ..
* .. Parameters ..
- DOUBLE COMPLEX ONE
+ COMPLEX*16 ONE
PARAMETER (ONE= (1.0D+0,0.0D+0))
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
*
CHARACTER DIAG,TRANS,UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),X(*)
+ COMPLEX*16 A(LDA,*),X(*)
* ..
*
* Purpose
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP
+ COMPLEX*16 TEMP
INTEGER I,INFO,IX,J,JX,KPLUS1,KX,L
LOGICAL NOCONJ,NOUNIT
* ..
CHARACTER DIAG,TRANS,UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),X(*)
+ COMPLEX*16 A(LDA,*),X(*)
* ..
*
* Purpose
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP
+ COMPLEX*16 TEMP
INTEGER I,INFO,IX,J,JX,KPLUS1,KX,L
LOGICAL NOCONJ,NOUNIT
* ..
CHARACTER DIAG,TRANS,UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX AP(*),X(*)
+ COMPLEX*16 AP(*),X(*)
* ..
*
* Purpose
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP
+ COMPLEX*16 TEMP
INTEGER I,INFO,IX,J,JX,K,KK,KX
LOGICAL NOCONJ,NOUNIT
* ..
CHARACTER DIAG,TRANS,UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX AP(*),X(*)
+ COMPLEX*16 AP(*),X(*)
* ..
*
* Purpose
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP
+ COMPLEX*16 TEMP
INTEGER I,INFO,IX,J,JX,K,KK,KX
LOGICAL NOCONJ,NOUNIT
* ..
SUBROUTINE ZTRMM(SIDE,UPLO,TRANSA,DIAG,M,N,ALPHA,A,LDA,B,LDB)
* .. Scalar Arguments ..
- DOUBLE COMPLEX ALPHA
+ COMPLEX*16 ALPHA
INTEGER LDA,LDB,M,N
CHARACTER DIAG,SIDE,TRANSA,UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),B(LDB,*)
+ COMPLEX*16 A(LDA,*),B(LDB,*)
* ..
*
* Purpose
INTRINSIC DCONJG,MAX
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP
+ COMPLEX*16 TEMP
INTEGER I,INFO,J,K,NROWA
LOGICAL LSIDE,NOCONJ,NOUNIT,UPPER
* ..
* .. Parameters ..
- DOUBLE COMPLEX ONE
+ COMPLEX*16 ONE
PARAMETER (ONE= (1.0D+0,0.0D+0))
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
*
CHARACTER DIAG,TRANS,UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),X(*)
+ COMPLEX*16 A(LDA,*),X(*)
* ..
*
* Purpose
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP
+ COMPLEX*16 TEMP
INTEGER I,INFO,IX,J,JX,KX
LOGICAL NOCONJ,NOUNIT
* ..
SUBROUTINE ZTRSM(SIDE,UPLO,TRANSA,DIAG,M,N,ALPHA,A,LDA,B,LDB)
* .. Scalar Arguments ..
- DOUBLE COMPLEX ALPHA
+ COMPLEX*16 ALPHA
INTEGER LDA,LDB,M,N
CHARACTER DIAG,SIDE,TRANSA,UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),B(LDB,*)
+ COMPLEX*16 A(LDA,*),B(LDB,*)
* ..
*
* Purpose
INTRINSIC DCONJG,MAX
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP
+ COMPLEX*16 TEMP
INTEGER I,INFO,J,K,NROWA
LOGICAL LSIDE,NOCONJ,NOUNIT,UPPER
* ..
* .. Parameters ..
- DOUBLE COMPLEX ONE
+ COMPLEX*16 ONE
PARAMETER (ONE= (1.0D+0,0.0D+0))
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
*
CHARACTER DIAG,TRANS,UPLO
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A(LDA,*),X(*)
+ COMPLEX*16 A(LDA,*),X(*)
* ..
*
* Purpose
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER (ZERO= (0.0D+0,0.0D+0))
* ..
* .. Local Scalars ..
- DOUBLE COMPLEX TEMP
+ COMPLEX*16 TEMP
INTEGER I,INFO,IX,J,JX,KX
LOGICAL NOCONJ,NOUNIT
* ..
END DO
ELSE IF (ICASE.EQ.13) THEN
* .. DSDOT ..
- CALL TESTDSDOT(SNGL(DSDOT(N,SNGL(SX),INCX,SNGL(SY),INCY)),
- $ SNGL(DT7(KN,KI)),SNGL(SSIZE1(KN)), .3125E-1)
+ CALL TESTDSDOT(REAL(DSDOT(N,REAL(SX),INCX,REAL(SY),INCY)),
+ $ REAL(DT7(KN,KI)),REAL(SSIZE1(KN)), .3125E-1)
ELSE
WRITE (NOUT,*) ' Shouldn''t be here in CHECK2'
STOP
* End of LZERES.
*
END
- DOUBLE COMPLEX FUNCTION ZBEG( RESET )
+ COMPLEX*16 FUNCTION ZBEG( RESET )
*
* Generates complex numbers as pairs of random numbers uniformly
* distributed between -0.5 and 0.5.
* End of LZERES.
*
END
- DOUBLE COMPLEX FUNCTION ZBEG( RESET )
+ COMPLEX*16 FUNCTION ZBEG( RESET )
*
* Generates complex numbers as pairs of random numbers uniformly
* distributed between -0.5 and 0.5.
* ..
* .. Array Arguments ..
INTEGER IPIV( * )
- DOUBLE COMPLEX A( LDA, * ), B( LDB, * ), WORK( * )
+ COMPLEX*16 A( LDA, * ), B( LDB, * ), WORK( * )
* ..
*
* Purpose
* The number of right hand sides, i.e., the number of columns
* of the matrix B. NRHS >= 0.
*
-* A (input) DOUBLE COMPLEX array, dimension (LDA,N)
+* A (input) COMPLEX*16 array, dimension (LDA,N)
* The block diagonal matrix D and the multipliers used to
* obtain the factor U or L as computed by ZHETRF.
*
* Details of the interchanges and the block structure of D
* as determined by ZHETRF.
*
-* B (input/output) DOUBLE COMPLEX array, dimension (LDB,NRHS)
+* B (input/output) COMPLEX*16 array, dimension (LDB,NRHS)
* On entry, the right hand side matrix B.
* On exit, the solution matrix X.
*
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ONE
+ COMPLEX*16 ONE
PARAMETER ( ONE = (1.0D+0,0.0D+0) )
* ..
* .. Local Scalars ..
LOGICAL UPPER
INTEGER I, IINFO, J, K, KP
DOUBLE PRECISION S
- DOUBLE COMPLEX AK, AKM1, AKM1K, BK, BKM1, DENOM
+ COMPLEX*16 AK, AKM1, AKM1K, BK, BKM1, DENOM
* ..
* .. External Functions ..
LOGICAL LSAME
* ..
* .. Array Arguments ..
INTEGER IPIV( * )
- DOUBLE COMPLEX A( LDA, * ), WORK( * )
+ COMPLEX*16 A( LDA, * ), WORK( * )
* ..
*
* Purpose
* N (input) INTEGER
* The order of the matrix A. N >= 0.
*
-* A (input) DOUBLE COMPLEX array, dimension (LDA,N)
+* A (input) COMPLEX*16 array, dimension (LDA,N)
* The block diagonal matrix D and the multipliers used to
* obtain the factor U or L as computed by ZSYTRF.
*
* Details of the interchanges and the block structure of D
* as determined by ZSYTRF.
*
-* WORK (workspace) DOUBLE COMPLEX array, dimension (N)
+* WORK (workspace) COMPLEX*16 array, dimension (N)
*
* LWORK (input) INTEGER
* The length of WORK. LWORK >=1.
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ZERO
+ COMPLEX*16 ZERO
PARAMETER ( ZERO = (0.0D+0,0.0D+0) )
* ..
* .. External Functions ..
* .. Local Scalars ..
LOGICAL UPPER, CONVERT
INTEGER I, IP, J
- DOUBLE COMPLEX TEMP
+ COMPLEX*16 TEMP
* ..
* .. Executable Statements ..
*
INTEGER I1, I2, LDA, N
* ..
* .. Array Arguments ..
- DOUBLE COMPLEX A( LDA, N )
+ COMPLEX*16 A( LDA, N )
*
* Purpose
* =======
* N (input) INTEGER
* The order of the matrix A. N >= 0.
*
-* A (input/output) DOUBLE COMPLEX array, dimension (LDA,N)
+* A (input/output) COMPLEX*16 array, dimension (LDA,N)
* On entry, the NB diagonal matrix D and the multipliers
* used to obtain the factor U or L as computed by ZSYTRF.
*
* .. Local Scalars ..
LOGICAL UPPER
INTEGER I
- DOUBLE COMPLEX TMP
+ COMPLEX*16 TMP
*
* .. External Functions ..
LOGICAL LSAME
* ..
* .. Array Arguments ..
INTEGER IPIV( * )
- DOUBLE COMPLEX A( LDA, * ), WORK( N+NB+1,* )
+ COMPLEX*16 A( LDA, * ), WORK( N+NB+1,* )
* ..
*
* Purpose
* N (input) INTEGER
* The order of the matrix A. N >= 0.
*
-* A (input/output) DOUBLE COMPLEX array, dimension (LDA,N)
+* A (input/output) COMPLEX*16 array, dimension (LDA,N)
* On entry, the NNB diagonal matrix D and the multipliers
* used to obtain the factor U or L as computed by ZSYTRF.
*
* Details of the interchanges and the NNB structure of D
* as determined by ZSYTRF.
*
-* WORK (workspace) DOUBLE COMPLEX array, dimension (N+NNB+1,NNB+3)
+* WORK (workspace) COMPLEX*16 array, dimension (N+NNB+1,NNB+3)
*
* NB (input) INTEGER
* Block size
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ONE, ZERO
+ COMPLEX*16 ONE, ZERO
PARAMETER ( ONE = ( 1.0D+0, 0.0D+0 ),
$ ZERO = ( 0.0D+0, 0.0D+0 ) )
* ..
INTEGER COUNT
INTEGER J, U11, INVD
- DOUBLE COMPLEX AK, AKKP1, AKP1, D, T
- DOUBLE COMPLEX U01_I_J, U01_IP1_J
- DOUBLE COMPLEX U11_I_J, U11_IP1_J
+ COMPLEX*16 AK, AKKP1, AKP1, D, T
+ COMPLEX*16 U01_I_J, U01_IP1_J
+ COMPLEX*16 U11_I_J, U11_IP1_J
* ..
* .. External Functions ..
LOGICAL LSAME
* ..
* .. Array Arguments ..
INTEGER IPIV( * )
- DOUBLE COMPLEX A( LDA, * ), B( LDB, * ), WORK( * )
+ COMPLEX*16 A( LDA, * ), B( LDB, * ), WORK( * )
* ..
*
* Purpose
* The number of right hand sides, i.e., the number of columns
* of the matrix B. NRHS >= 0.
*
-* A (input) DOUBLE COMPLEX array, dimension (LDA,N)
+* A (input) COMPLEX*16 array, dimension (LDA,N)
* The block diagonal matrix D and the multipliers used to
* obtain the factor U or L as computed by ZSYTRF.
*
* Details of the interchanges and the block structure of D
* as determined by ZSYTRF.
*
-* B (input/output) DOUBLE COMPLEX array, dimension (LDB,NRHS)
+* B (input/output) COMPLEX*16 array, dimension (LDB,NRHS)
* On entry, the right hand side matrix B.
* On exit, the solution matrix X.
*
* =====================================================================
*
* .. Parameters ..
- DOUBLE COMPLEX ONE
+ COMPLEX*16 ONE
PARAMETER ( ONE = (1.0D+0,0.0D+0) )
* ..
* .. Local Scalars ..
LOGICAL UPPER
INTEGER I, IINFO, J, K, KP
- DOUBLE COMPLEX AK, AKM1, AKM1K, BK, BKM1, DENOM
+ COMPLEX*16 AK, AKM1, AKM1K, BK, BKM1, DENOM
* ..
* .. External Functions ..
LOGICAL LSAME
- DOUBLE COMPLEX FUNCTION ZLARND( IDIST, ISEED )
+ COMPLEX*16 FUNCTION ZLARND( IDIST, ISEED )
*
* -- LAPACK auxiliary routine (version 3.1) --
* Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd..
- DOUBLE COMPLEX FUNCTION ZLATM2( M, N, I, J, KL, KU, IDIST,
+ COMPLEX*16 FUNCTION ZLATM2( M, N, I, J, KL, KU, IDIST,
$ ISEED, D, IGRADE, DL, DR, IPVTNG, IWORK, SPARSE )
*
* -- LAPACK auxiliary test routine (version 3.1) --
- DOUBLE COMPLEX FUNCTION ZLATM3( M, N, I, J, ISUB, JSUB, KL, KU,
+ COMPLEX*16 FUNCTION ZLATM3( M, N, I, J, ISUB, JSUB, KL, KU,
$ IDIST, ISEED, D, IGRADE, DL, DR, IPVTNG, IWORK,
$ SPARSE )
*