Fix precision mismatches in tests

* For REAL tests, replace 1.D0 with 1.0
* For CMPLX tests, replace CMPLX(1.D0, 1.D0) with CMPLX(1.0, 1.0)
* [sz]het01_aa.f: replace dimension cast DBLE(N) with REAL(N)
* [cz]syt01_aa.f: replace dimension cast DBLE(N) with REAL(N)
  and init CONE and CZERO as complex numbers

Author: angsch

TESTING/EIG/serrbd.f

@@ -111,7 +111,7 @@ SUBROUTINE SERRBD( PATH, NUNIT )
 *
       DO 20 J = 1, NMAX
          DO 10 I = 1, NMAX
-            A( I, J ) = 1.D0 / REAL( I+J )
+            A( I, J ) = 1. / REAL( I+J )
    10    CONTINUE
    20 CONTINUE
       OK = .TRUE.

TESTING/LIN/cerrrfp.f

@@ -89,10 +89,10 @@ SUBROUTINE CERRRFP( NUNIT )
 *
       NOUT = NUNIT
       OK = .TRUE.
-      A( 1, 1 ) = CMPLX( 1.D0 , 1.D0  )
-      B( 1, 1 ) = CMPLX( 1.D0 , 1.D0  )
-      ALPHA     = CMPLX( 1.D0 , 1.D0  )
-      BETA      = CMPLX( 1.D0 , 1.D0  )
+      A( 1, 1 ) = CMPLX( 1.0 , 1.0 )
+      B( 1, 1 ) = CMPLX( 1.0 , 1.0  )
+      ALPHA     = CMPLX( 1.0 , 1.0  )
+      BETA      = CMPLX( 1.0 , 1.0  )
 *
       SRNAMT = 'CPFTRF'
       INFOT = 1

TESTING/LIN/chet01_aa.f

@@ -159,7 +159,7 @@ SUBROUTINE CHET01_AA( UPLO, N, A, LDA, AFAC, LDAFAC, IPIV, C,
       EXTERNAL           CLASET, CLAVHE
 *     ..
 *     .. Intrinsic Functions ..
-      INTRINSIC          DBLE
+      INTRINSIC          REAL
 *     ..
 *     .. Executable Statements ..
 *
@@ -255,7 +255,7 @@ SUBROUTINE CHET01_AA( UPLO, N, A, LDA, AFAC, LDAFAC, IPIV, C,
          IF( RESID.NE.ZERO )
      $      RESID = ONE / EPS
       ELSE
-         RESID = ( ( RESID / DBLE( N ) ) / ANORM ) / EPS
+         RESID = ( ( RESID / REAL( N ) ) / ANORM ) / EPS
       END IF
 *
       RETURN

TESTING/LIN/csyt01_aa.f

@@ -141,9 +141,10 @@ SUBROUTINE CSYT01_AA( UPLO, N, A, LDA, AFAC, LDAFAC, IPIV, C,
 *
 *     .. Parameters ..
       REAL               ZERO, ONE
-      PARAMETER          ( ZERO = 0.0D+0, ONE = 1.0D+0 )
+      PARAMETER          ( ZERO = 0.0E+0, ONE = 1.0E+0 )
       COMPLEX            CZERO, CONE
-      PARAMETER          ( CZERO = 0.0E+0, CONE = 1.0E+0 )
+      PARAMETER          ( CZERO = ( 0.0E+0, 0.0E+0 ),
+     $                     CONE  = ( 1.0E+0, 0.0E+0 ) )
 *     ..
 *     .. Local Scalars ..
       INTEGER            I, J
@@ -158,7 +159,7 @@ SUBROUTINE CSYT01_AA( UPLO, N, A, LDA, AFAC, LDAFAC, IPIV, C,
       EXTERNAL           CLASET, CLAVSY
 *     ..
 *     .. Intrinsic Functions ..
-      INTRINSIC          DBLE
+      INTRINSIC          REAL
 *     ..
 *     .. Executable Statements ..
 *
@@ -250,7 +251,7 @@ SUBROUTINE CSYT01_AA( UPLO, N, A, LDA, AFAC, LDAFAC, IPIV, C,
          IF( RESID.NE.ZERO )
      $      RESID = ONE / EPS
       ELSE
-         RESID = ( ( RESID / DBLE( N ) ) / ANORM ) / EPS
+         RESID = ( ( RESID / REAL( N ) ) / ANORM ) / EPS
       END IF
 *
       RETURN

TESTING/LIN/serrlqt.f

@@ -101,11 +101,11 @@ SUBROUTINE SERRLQT( PATH, NUNIT )
 *
       DO J = 1, NMAX
          DO I = 1, NMAX
-            A( I, J ) = 1.D0 / REAL( I+J )
-            C( I, J ) = 1.D0 / REAL( I+J )
-            T( I, J ) = 1.D0 / REAL( I+J )
+            A( I, J ) = 1. / REAL( I+J )
+            C( I, J ) = 1. / REAL( I+J )
+            T( I, J ) = 1. / REAL( I+J )
          END DO
-         W( J ) = 0.D0
+         W( J ) = 0.
       END DO
       OK = .TRUE.
 *

TESTING/LIN/serrlqtp.f

@@ -101,11 +101,11 @@ SUBROUTINE SERRLQTP( PATH, NUNIT )
 *
       DO J = 1, NMAX
          DO I = 1, NMAX
-            A( I, J ) = 1.D0 / REAL( I+J )
-            C( I, J ) = 1.D0 / REAL( I+J )
-            T( I, J ) = 1.D0 / REAL( I+J )
+            A( I, J ) = 1. / REAL( I+J )
+            C( I, J ) = 1. / REAL( I+J )
+            T( I, J ) = 1. / REAL( I+J )
          END DO
-         W( J ) = 0.0
+         W( J ) = 0.
       END DO
       OK = .TRUE.
 *

TESTING/LIN/serrtsqr.f

@@ -101,11 +101,11 @@ SUBROUTINE SERRTSQR( PATH, NUNIT )
 *
       DO J = 1, NMAX
          DO I = 1, NMAX
-            A( I, J ) = 1.D0 / REAL( I+J )
-            C( I, J ) = 1.D0 / REAL( I+J )
-            T( I, J ) = 1.D0 / REAL( I+J )
+            A( I, J ) = 1. / REAL( I+J )
+            C( I, J ) = 1. / REAL( I+J )
+            T( I, J ) = 1. / REAL( I+J )
          END DO
-         W( J ) = 0.D0
+         W( J ) = 0.
       END DO
       OK = .TRUE.
 *

TESTING/LIN/ssyt01_aa.f

@@ -156,7 +156,7 @@ SUBROUTINE SSYT01_AA( UPLO, N, A, LDA, AFAC, LDAFAC, IPIV, C,
       EXTERNAL           SLASET, SLAVSY, SSWAP, STRMM, SLACPY
 *     ..
 *     .. Intrinsic Functions ..
-      INTRINSIC          DBLE
+      INTRINSIC          REAL
 *     ..
 *     .. Executable Statements ..
 *
@@ -248,7 +248,7 @@ SUBROUTINE SSYT01_AA( UPLO, N, A, LDA, AFAC, LDAFAC, IPIV, C,
          IF( RESID.NE.ZERO )
      $      RESID = ONE / EPS
       ELSE
-         RESID = ( ( RESID / DBLE( N ) ) / ANORM ) / EPS
+         RESID = ( ( RESID / REAL( N ) ) / ANORM ) / EPS
       END IF
 *
       RETURN

TESTING/LIN/zsyt01_aa.f

@@ -142,8 +142,9 @@ SUBROUTINE ZSYT01_AA( UPLO, N, A, LDA, AFAC, LDAFAC, IPIV, C,
 *     .. Parameters ..
       DOUBLE PRECISION   ZERO, ONE
       PARAMETER          ( ZERO = 0.0D+0, ONE = 1.0D+0 )
-      COMPLEX*16         CZERO, CONE
-      PARAMETER          ( CZERO = 0.0E+0, CONE = 1.0E+0 )
+      COMPLEX*16          CZERO, CONE
+      PARAMETER          ( CZERO = ( 0.0D+0, 0.0D+0 ),
+     $                     CONE  = ( 1.0D+0, 0.0D+0 ) )
 *     ..
 *     .. Local Scalars ..
       INTEGER            I, J