Replace DOUBLE COMPLEX by COMPLEX*16

Replace SNGL by REAL

42 files changed, 130 insertions, 131 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.