Replace DOUBLE COMPLEX by COMPLEX*16

Replace SNGL by REAL

50 files changed, 158 insertions, 159 deletions

diff --git a/BLAS/SRC/dcabs1.f b/BLAS/SRC/dcabs1.f
index 4d5fb697..73a51fbe 100644
--- a/BLAS/SRC/dcabs1.f
+++ b/BLAS/SRC/dcabs1.f
@@ -1,6 +1,6 @@
       DOUBLE PRECISION FUNCTION DCABS1(Z)
 *     .. Scalar Arguments ..
-      DOUBLE COMPLEX Z
+      COMPLEX*16 Z
 *     ..
 *     ..
 *  Purpose
diff --git a/BLAS/SRC/dzasum.f b/BLAS/SRC/dzasum.f
index a5c440c5..992b3a5c 100644
--- a/BLAS/SRC/dzasum.f
+++ b/BLAS/SRC/dzasum.f
@@ -3,7 +3,7 @@
       INTEGER INCX,N
 *     ..
 *     .. Array Arguments ..
-      DOUBLE COMPLEX ZX(*)
+      COMPLEX*16 ZX(*)
 *     ..
 *
 *  Purpose
diff --git a/BLAS/SRC/dznrm2.f b/BLAS/SRC/dznrm2.f
index 71b5674b..a5e3da08 100644
--- a/BLAS/SRC/dznrm2.f
+++ b/BLAS/SRC/dznrm2.f
@@ -3,7 +3,7 @@
       INTEGER INCX,N
 *     ..
 *     .. Array Arguments ..
-      DOUBLE COMPLEX X(*)
+      COMPLEX*16 X(*)
 *     ..
 *
 *  Purpose
diff --git a/BLAS/SRC/izamax.f b/BLAS/SRC/izamax.f
index af14fb16..9ab4b45d 100644
--- a/BLAS/SRC/izamax.f
+++ b/BLAS/SRC/izamax.f
@@ -3,7 +3,7 @@
       INTEGER INCX,N
 *     ..
 *     .. Array Arguments ..
-      DOUBLE COMPLEX ZX(*)
+      COMPLEX*16 ZX(*)
 *     ..
 *
 *  Purpose
diff --git a/BLAS/SRC/srotmg.f b/BLAS/SRC/srotmg.f
index 4a64e687..aae4f284 100644
--- a/BLAS/SRC/srotmg.f
+++ b/BLAS/SRC/srotmg.f
@@ -188,7 +188,6 @@
          SPARAM(5) = SH22
       END IF
 
-  260 CONTINUE
       SPARAM(1) = SFLAG
       RETURN
       END
diff --git a/BLAS/SRC/zaxpy.f b/BLAS/SRC/zaxpy.f
index fc32a86e..07afb5f3 100644
--- a/BLAS/SRC/zaxpy.f
+++ b/BLAS/SRC/zaxpy.f
@@ -1,10 +1,10 @@
       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
diff --git a/BLAS/SRC/zcopy.f b/BLAS/SRC/zcopy.f
index 49bf4aca..1038cb96 100644
--- a/BLAS/SRC/zcopy.f
+++ b/BLAS/SRC/zcopy.f
@@ -3,7 +3,7 @@
       INTEGER INCX,INCY,N
 *     ..
 *     .. Array Arguments ..
-      DOUBLE COMPLEX ZX(*),ZY(*)
+      COMPLEX*16 ZX(*),ZY(*)
 *     ..
 *
 *  Purpose
diff --git a/BLAS/SRC/zdotc.f b/BLAS/SRC/zdotc.f
index 506ee9ef..2f6b6532 100644
--- a/BLAS/SRC/zdotc.f
+++ b/BLAS/SRC/zdotc.f
@@ -1,9 +1,9 @@
-      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
@@ -20,7 +20,7 @@
 *  =====================================================================
 *
 *     .. Local Scalars ..
-      DOUBLE COMPLEX ZTEMP
+      COMPLEX*16 ZTEMP
       INTEGER I,IX,IY
 *     ..
 *     .. Intrinsic Functions ..
diff --git a/BLAS/SRC/zdotu.f b/BLAS/SRC/zdotu.f
index 816b2014..1722c8af 100644
--- a/BLAS/SRC/zdotu.f
+++ b/BLAS/SRC/zdotu.f
@@ -1,9 +1,9 @@
-      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
@@ -20,7 +20,7 @@
 *  =====================================================================
 *
 *     .. Local Scalars ..
-      DOUBLE COMPLEX ZTEMP
+      COMPLEX*16 ZTEMP
       INTEGER I,IX,IY
 *     ..
       ZTEMP = (0.0d0,0.0d0)
diff --git a/BLAS/SRC/zdscal.f b/BLAS/SRC/zdscal.f
index 23030e10..d4d96c75 100644
--- a/BLAS/SRC/zdscal.f
+++ b/BLAS/SRC/zdscal.f
@@ -4,7 +4,7 @@
       INTEGER INCX,N
 *     ..
 *     .. Array Arguments ..
-      DOUBLE COMPLEX ZX(*)
+      COMPLEX*16 ZX(*)
 *     ..
 *
 *  Purpose
diff --git a/BLAS/SRC/zgbmv.f b/BLAS/SRC/zgbmv.f
index 04a29d46..3000be42 100644
--- a/BLAS/SRC/zgbmv.f
+++ b/BLAS/SRC/zgbmv.f
@@ -1,11 +1,11 @@
       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
@@ -133,13 +133,13 @@
 *  =====================================================================
 *
 *     .. 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
 *     ..
diff --git a/BLAS/SRC/zgemm.f b/BLAS/SRC/zgemm.f
index f4f98981..8954cbfb 100644
--- a/BLAS/SRC/zgemm.f
+++ b/BLAS/SRC/zgemm.f
@@ -1,11 +1,11 @@
       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
@@ -142,14 +142,14 @@
       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))
 *     ..
 *
diff --git a/BLAS/SRC/zgemv.f b/BLAS/SRC/zgemv.f
index bb2ae4fc..9415940e 100644
--- a/BLAS/SRC/zgemv.f
+++ b/BLAS/SRC/zgemv.f
@@ -1,11 +1,11 @@
       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
@@ -106,13 +106,13 @@
 *  =====================================================================
 *
 *     .. 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
 *     ..
diff --git a/BLAS/SRC/zgerc.f b/BLAS/SRC/zgerc.f
index 06fdf597..a8872b49 100644
--- a/BLAS/SRC/zgerc.f
+++ b/BLAS/SRC/zgerc.f
@@ -1,10 +1,10 @@
       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
@@ -81,11 +81,11 @@
 *  =====================================================================
 *
 *     .. 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 ..
diff --git a/BLAS/SRC/zgeru.f b/BLAS/SRC/zgeru.f
index dd43b7c5..3ece8018 100644
--- a/BLAS/SRC/zgeru.f
+++ b/BLAS/SRC/zgeru.f
@@ -1,10 +1,10 @@
       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
@@ -81,11 +81,11 @@
 *  =====================================================================
 *
 *     .. 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 ..
diff --git a/BLAS/SRC/zhbmv.f b/BLAS/SRC/zhbmv.f
index 6f6ef0b1..90988684 100644
--- a/BLAS/SRC/zhbmv.f
+++ b/BLAS/SRC/zhbmv.f
@@ -1,11 +1,11 @@
       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
@@ -135,13 +135,13 @@
 *  =====================================================================
 *
 *     .. 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 ..
diff --git a/BLAS/SRC/zhemm.f b/BLAS/SRC/zhemm.f
index 5c8d9500..0537e919 100644
--- a/BLAS/SRC/zhemm.f
+++ b/BLAS/SRC/zhemm.f
@@ -1,11 +1,11 @@
       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
@@ -148,14 +148,14 @@
       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))
 *     ..
 *
diff --git a/BLAS/SRC/zhemv.f b/BLAS/SRC/zhemv.f
index d2e7d132..8e61e93b 100644
--- a/BLAS/SRC/zhemv.f
+++ b/BLAS/SRC/zhemv.f
@@ -1,11 +1,11 @@
       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
@@ -104,13 +104,13 @@
 *  =====================================================================
 *
 *     .. 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 ..
diff --git a/BLAS/SRC/zher.f b/BLAS/SRC/zher.f
index 2f82e019..898a1e03 100644
--- a/BLAS/SRC/zher.f
+++ b/BLAS/SRC/zher.f
@@ -5,7 +5,7 @@
       CHARACTER UPLO
 *     ..
 *     .. Array Arguments ..
-      DOUBLE COMPLEX A(LDA,*),X(*)
+      COMPLEX*16 A(LDA,*),X(*)
 *     ..
 *
 *  Purpose
@@ -91,11 +91,11 @@
 *  =====================================================================
 *
 *     .. 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 ..
diff --git a/BLAS/SRC/zher2.f b/BLAS/SRC/zher2.f
index a2cc1164..bb414e9e 100644
--- a/BLAS/SRC/zher2.f
+++ b/BLAS/SRC/zher2.f
@@ -1,11 +1,11 @@
       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
@@ -102,11 +102,11 @@
 *  =====================================================================
 *
 *     .. 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 ..
diff --git a/BLAS/SRC/zher2k.f b/BLAS/SRC/zher2k.f
index e20f2040..e5f7be83 100644
--- a/BLAS/SRC/zher2k.f
+++ b/BLAS/SRC/zher2k.f
@@ -1,12 +1,12 @@
       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
@@ -154,14 +154,14 @@
       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))
 *     ..
 *
diff --git a/BLAS/SRC/zherk.f b/BLAS/SRC/zherk.f
index d0684707..af04413a 100644
--- a/BLAS/SRC/zherk.f
+++ b/BLAS/SRC/zherk.f
@@ -5,7 +5,7 @@
       CHARACTER TRANS,UPLO
 *     ..
 *     .. Array Arguments ..
-      DOUBLE COMPLEX A(LDA,*),C(LDC,*)
+      COMPLEX*16 A(LDA,*),C(LDC,*)
 *     ..
 *
 *  Purpose
@@ -134,7 +134,7 @@
       INTRINSIC DBLE,DCMPLX,DCONJG,MAX
 *     ..
 *     .. Local Scalars ..
-      DOUBLE COMPLEX TEMP
+      COMPLEX*16 TEMP
       DOUBLE PRECISION RTEMP
       INTEGER I,INFO,J,L,NROWA
       LOGICAL UPPER
diff --git a/BLAS/SRC/zhpmv.f b/BLAS/SRC/zhpmv.f
index c76c46a3..514cd20a 100644
--- a/BLAS/SRC/zhpmv.f
+++ b/BLAS/SRC/zhpmv.f
@@ -1,11 +1,11 @@
       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
@@ -101,13 +101,13 @@
 *  =====================================================================
 *
 *     .. 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 ..
diff --git a/BLAS/SRC/zhpr.f b/BLAS/SRC/zhpr.f
index 577d52cc..859805d2 100644
--- a/BLAS/SRC/zhpr.f
+++ b/BLAS/SRC/zhpr.f
@@ -5,7 +5,7 @@
       CHARACTER UPLO
 *     ..
 *     .. Array Arguments ..
-      DOUBLE COMPLEX AP(*),X(*)
+      COMPLEX*16 AP(*),X(*)
 *     ..
 *
 *  Purpose
@@ -88,11 +88,11 @@
 *  =====================================================================
 *
 *     .. 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 ..
diff --git a/BLAS/SRC/zhpr2.f b/BLAS/SRC/zhpr2.f
index 024fc88f..48b3591f 100644
--- a/BLAS/SRC/zhpr2.f
+++ b/BLAS/SRC/zhpr2.f
@@ -1,11 +1,11 @@
       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
@@ -99,11 +99,11 @@
 *  =====================================================================
 *
 *     .. 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 ..
diff --git a/BLAS/SRC/zrotg.f b/BLAS/SRC/zrotg.f
index 6656302f..c25fa654 100644
--- a/BLAS/SRC/zrotg.f
+++ b/BLAS/SRC/zrotg.f
@@ -1,6 +1,6 @@
       SUBROUTINE ZROTG(CA,CB,C,S)
 *     .. Scalar Arguments ..
-      DOUBLE COMPLEX CA,CB,S
+      COMPLEX*16 CA,CB,S
       DOUBLE PRECISION C
 *     ..
 *
@@ -12,7 +12,7 @@
 *  =====================================================================
 *
 *     .. Local Scalars ..
-      DOUBLE COMPLEX ALPHA
+      COMPLEX*16 ALPHA
       DOUBLE PRECISION NORM,SCALE
 *     ..
 *     .. Intrinsic Functions ..
diff --git a/BLAS/SRC/zscal.f b/BLAS/SRC/zscal.f
index a84c618a..d5acf683 100644
--- a/BLAS/SRC/zscal.f
+++ b/BLAS/SRC/zscal.f
@@ -1,10 +1,10 @@
       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
diff --git a/BLAS/SRC/zswap.f b/BLAS/SRC/zswap.f
index 47894581..d3bb047e 100644
--- a/BLAS/SRC/zswap.f
+++ b/BLAS/SRC/zswap.f
@@ -3,7 +3,7 @@
       INTEGER INCX,INCY,N
 *     ..
 *     .. Array Arguments ..
-      DOUBLE COMPLEX ZX(*),ZY(*)
+      COMPLEX*16 ZX(*),ZY(*)
 *     ..
 *
 *  Purpose
@@ -20,7 +20,7 @@
 *  =====================================================================
 *
 *     .. Local Scalars ..
-      DOUBLE COMPLEX ZTEMP
+      COMPLEX*16 ZTEMP
       INTEGER I,IX,IY
 *     ..
       IF (N.LE.0) RETURN
diff --git a/BLAS/SRC/zsymm.f b/BLAS/SRC/zsymm.f
index f5d07e1c..e0192e9f 100644
--- a/BLAS/SRC/zsymm.f
+++ b/BLAS/SRC/zsymm.f
@@ -1,11 +1,11 @@
       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
@@ -146,14 +146,14 @@
       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))
 *     ..
 *
diff --git a/BLAS/SRC/zsyr2k.f b/BLAS/SRC/zsyr2k.f
index 808b7ede..a9bd0f4d 100644
--- a/BLAS/SRC/zsyr2k.f
+++ b/BLAS/SRC/zsyr2k.f
@@ -1,11 +1,11 @@
       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
@@ -145,14 +145,14 @@
       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))
 *     ..
 *
diff --git a/BLAS/SRC/zsyrk.f b/BLAS/SRC/zsyrk.f
index ecc1285d..b2acf31e 100644
--- a/BLAS/SRC/zsyrk.f
+++ b/BLAS/SRC/zsyrk.f
@@ -1,11 +1,11 @@
       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
@@ -128,14 +128,14 @@
       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))
 *     ..
 *
diff --git a/BLAS/SRC/ztbmv.f b/BLAS/SRC/ztbmv.f
index c7fdb8f3..2d2fbe17 100644
--- a/BLAS/SRC/ztbmv.f
+++ b/BLAS/SRC/ztbmv.f
@@ -4,7 +4,7 @@
       CHARACTER DIAG,TRANS,UPLO
 *     ..
 *     .. Array Arguments ..
-      DOUBLE COMPLEX A(LDA,*),X(*)
+      COMPLEX*16 A(LDA,*),X(*)
 *     ..
 *
 *  Purpose
@@ -140,11 +140,11 @@
 *  =====================================================================
 *
 *     .. 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
 *     ..
diff --git a/BLAS/SRC/ztbsv.f b/BLAS/SRC/ztbsv.f
index 9f3bf510..80ebb740 100644
--- a/BLAS/SRC/ztbsv.f
+++ b/BLAS/SRC/ztbsv.f
@@ -4,7 +4,7 @@
       CHARACTER DIAG,TRANS,UPLO
 *     ..
 *     .. Array Arguments ..
-      DOUBLE COMPLEX A(LDA,*),X(*)
+      COMPLEX*16 A(LDA,*),X(*)
 *     ..
 *
 *  Purpose
@@ -143,11 +143,11 @@
 *  =====================================================================
 *
 *     .. 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
 *     ..
diff --git a/BLAS/SRC/ztpmv.f b/BLAS/SRC/ztpmv.f
index bf74ec04..33f4b52f 100644
--- a/BLAS/SRC/ztpmv.f
+++ b/BLAS/SRC/ztpmv.f
@@ -4,7 +4,7 @@
       CHARACTER DIAG,TRANS,UPLO
 *     ..
 *     .. Array Arguments ..
-      DOUBLE COMPLEX AP(*),X(*)
+      COMPLEX*16 AP(*),X(*)
 *     ..
 *
 *  Purpose
@@ -100,11 +100,11 @@
 *  =====================================================================
 *
 *     .. 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
 *     ..
diff --git a/BLAS/SRC/ztpsv.f b/BLAS/SRC/ztpsv.f
index d8aae2ed..55eaf9e1 100644
--- a/BLAS/SRC/ztpsv.f
+++ b/BLAS/SRC/ztpsv.f
@@ -4,7 +4,7 @@
       CHARACTER DIAG,TRANS,UPLO
 *     ..
 *     .. Array Arguments ..
-      DOUBLE COMPLEX AP(*),X(*)
+      COMPLEX*16 AP(*),X(*)
 *     ..
 *
 *  Purpose
@@ -102,11 +102,11 @@
 *  =====================================================================
 *
 *     .. 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
 *     ..
diff --git a/BLAS/SRC/ztrmm.f b/BLAS/SRC/ztrmm.f
index 9ea6b9b9..e02efd43 100644
--- a/BLAS/SRC/ztrmm.f
+++ b/BLAS/SRC/ztrmm.f
@@ -1,11 +1,11 @@
       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
@@ -138,14 +138,14 @@
       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))
 *     ..
 *
diff --git a/BLAS/SRC/ztrmv.f b/BLAS/SRC/ztrmv.f
index 50708e1a..345c8cd9 100644
--- a/BLAS/SRC/ztrmv.f
+++ b/BLAS/SRC/ztrmv.f
@@ -4,7 +4,7 @@
       CHARACTER DIAG,TRANS,UPLO
 *     ..
 *     .. Array Arguments ..
-      DOUBLE COMPLEX A(LDA,*),X(*)
+      COMPLEX*16 A(LDA,*),X(*)
 *     ..
 *
 *  Purpose
@@ -103,11 +103,11 @@
 *  =====================================================================
 *
 *     .. 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
 *     ..
diff --git a/BLAS/SRC/ztrsm.f b/BLAS/SRC/ztrsm.f
index e8bee0eb..4cfe33b8 100644
--- a/BLAS/SRC/ztrsm.f
+++ b/BLAS/SRC/ztrsm.f
@@ -1,11 +1,11 @@
       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
@@ -140,14 +140,14 @@
       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))
 *     ..
 *
diff --git a/BLAS/SRC/ztrsv.f b/BLAS/SRC/ztrsv.f
index d5d54837..2527445f 100644
--- a/BLAS/SRC/ztrsv.f
+++ b/BLAS/SRC/ztrsv.f
@@ -4,7 +4,7 @@
       CHARACTER DIAG,TRANS,UPLO
 *     ..
 *     .. Array Arguments ..
-      DOUBLE COMPLEX A(LDA,*),X(*)
+      COMPLEX*16 A(LDA,*),X(*)
 *     ..
 *
 *  Purpose
@@ -105,11 +105,11 @@
 *  =====================================================================
 *
 *     .. 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
 *     ..
diff --git a/BLAS/TESTING/dblat1.f b/BLAS/TESTING/dblat1.f
index d2435042..249b6abf 100644
--- a/BLAS/TESTING/dblat1.f
+++ b/BLAS/TESTING/dblat1.f
@@ -625,8 +625,8 @@
                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
diff --git a/BLAS/TESTING/zblat2.f b/BLAS/TESTING/zblat2.f
index 8e49cabf..e7829fdf 100644
--- a/BLAS/TESTING/zblat2.f
+++ b/BLAS/TESTING/zblat2.f
@@ -3100,7 +3100,7 @@
 *     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.
diff --git a/BLAS/TESTING/zblat3.f b/BLAS/TESTING/zblat3.f
index fb471255..379d65b8 100644
--- a/BLAS/TESTING/zblat3.f
+++ b/BLAS/TESTING/zblat3.f
@@ -3305,7 +3305,7 @@
 *     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.
diff --git a/SRC/zhetrs2.f b/SRC/zhetrs2.f
index d0536136..6557e3e3 100644
--- a/SRC/zhetrs2.f
+++ b/SRC/zhetrs2.f
@@ -15,7 +15,7 @@
 *     ..
 *     .. Array Arguments ..
       INTEGER            IPIV( * )
-      DOUBLE COMPLEX   A( LDA, * ), B( LDB, * ), WORK( * )
+      COMPLEX*16       A( LDA, * ), B( LDB, * ), WORK( * )
 *     ..
 *
 *  Purpose
@@ -41,7 +41,7 @@
 *          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.
 *
@@ -52,7 +52,7 @@
 *          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.
 *
@@ -68,14 +68,14 @@
 *  =====================================================================
 *
 *     .. 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
diff --git a/SRC/zsyconv.f b/SRC/zsyconv.f
index ae01e61c..43cae75c 100644
--- a/SRC/zsyconv.f
+++ b/SRC/zsyconv.f
@@ -14,7 +14,7 @@
 *     ..
 *     .. Array Arguments ..
       INTEGER            IPIV( * )
-      DOUBLE COMPLEX     A( LDA, * ), WORK( * )
+      COMPLEX*16         A( LDA, * ), WORK( * )
 *     ..
 *
 *  Purpose
@@ -40,7 +40,7 @@
 *  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.
 *
@@ -51,7 +51,7 @@
 *          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. 
@@ -69,7 +69,7 @@
 *  =====================================================================
 *
 *     .. Parameters ..
-      DOUBLE COMPLEX     ZERO
+      COMPLEX*16         ZERO
       PARAMETER          ( ZERO = (0.0D+0,0.0D+0) )
 *     ..
 *     .. External Functions ..
@@ -81,7 +81,7 @@
 *     .. Local Scalars ..
       LOGICAL            UPPER, CONVERT
       INTEGER            I, IP, J
-      DOUBLE COMPLEX     TEMP
+      COMPLEX*16         TEMP
 *     ..
 *     .. Executable Statements ..
 *
diff --git a/SRC/zsyswapr.f b/SRC/zsyswapr.f
index 62cace36..40239000 100644
--- a/SRC/zsyswapr.f
+++ b/SRC/zsyswapr.f
@@ -10,7 +10,7 @@
       INTEGER          I1, I2, LDA, N
 *     ..
 *     .. Array Arguments ..
-      DOUBLE COMPLEX   A( LDA, N )
+      COMPLEX*16       A( LDA, N )
 *
 *  Purpose
 *  =======
@@ -30,7 +30,7 @@
 *  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.
 *
@@ -56,7 +56,7 @@
 *     .. Local Scalars ..
       LOGICAL            UPPER
       INTEGER            I
-      DOUBLE COMPLEX     TMP
+      COMPLEX*16         TMP
 *
 *     .. External Functions ..
       LOGICAL            LSAME
diff --git a/SRC/zsytri2x.f b/SRC/zsytri2x.f
index 584a9fea..51645389 100644
--- a/SRC/zsytri2x.f
+++ b/SRC/zsytri2x.f
@@ -13,7 +13,7 @@
 *     ..
 *     .. Array Arguments ..
       INTEGER            IPIV( * )
-      DOUBLE COMPLEX     A( LDA, * ), WORK( N+NB+1,* )
+      COMPLEX*16         A( LDA, * ), WORK( N+NB+1,* )
 *     ..
 *
 *  Purpose
@@ -35,7 +35,7 @@
 *  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.
 *
@@ -53,7 +53,7 @@
 *          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
@@ -67,7 +67,7 @@
 *  =====================================================================
 *
 *     .. Parameters ..
-      DOUBLE COMPLEX     ONE, ZERO
+      COMPLEX*16         ONE, ZERO
       PARAMETER          ( ONE = ( 1.0D+0, 0.0D+0 ),
      $                   ZERO = ( 0.0D+0, 0.0D+0 ) )
 *     ..
@@ -77,9 +77,9 @@
       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
diff --git a/SRC/zsytrs2.f b/SRC/zsytrs2.f
index b2df5e6d..3d01aa96 100644
--- a/SRC/zsytrs2.f
+++ b/SRC/zsytrs2.f
@@ -15,7 +15,7 @@
 *     ..
 *     .. Array Arguments ..
       INTEGER            IPIV( * )
-      DOUBLE COMPLEX   A( LDA, * ), B( LDB, * ), WORK( * )
+      COMPLEX*16       A( LDA, * ), B( LDB, * ), WORK( * )
 *     ..
 *
 *  Purpose
@@ -41,7 +41,7 @@
 *          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.
 *
@@ -52,7 +52,7 @@
 *          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.
 *
@@ -68,13 +68,13 @@
 *  =====================================================================
 *
 *     .. 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
diff --git a/TESTING/MATGEN/zlarnd.f b/TESTING/MATGEN/zlarnd.f
index b5b413fe..09a20c4f 100644
--- a/TESTING/MATGEN/zlarnd.f
+++ b/TESTING/MATGEN/zlarnd.f
@@ -1,4 +1,4 @@
-      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..
diff --git a/TESTING/MATGEN/zlatm2.f b/TESTING/MATGEN/zlatm2.f
index 26afabb1..c5c2d950 100644
--- a/TESTING/MATGEN/zlatm2.f
+++ b/TESTING/MATGEN/zlatm2.f
@@ -1,4 +1,4 @@
-      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) --
diff --git a/TESTING/MATGEN/zlatm3.f b/TESTING/MATGEN/zlatm3.f
index cf8f1fdd..e07ae49e 100644
--- a/TESTING/MATGEN/zlatm3.f
+++ b/TESTING/MATGEN/zlatm3.f
@@ -1,4 +1,4 @@
-      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 )
 *