diff --git a/CBLAS/testing/c_cblat2.f b/CBLAS/testing/c_cblat2.f index fce05e472..36b1d37f3 100644 --- a/CBLAS/testing/c_cblat2.f +++ b/CBLAS/testing/c_cblat2.f @@ -269,13 +269,14 @@ N = MIN( 32, NMAX ) DO 120 J = 1, N DO 110 I = 1, N - A( I, J ) = MAX( I - J + 1, 0 ) + A( I, J ) = REAL( MAX( I - J + 1, 0 ) ) 110 CONTINUE - X( J ) = J + X( J ) = REAL( J ) Y( J ) = ZERO 120 CONTINUE DO 130 J = 1, N - YY( J ) = J*( ( J + 1 )*J )/2 - ( ( J + 1 )*J*( J - 1 ) )/3 + YY( J ) = REAL( J*( ( J + 1 )*J )/2 - + $ ( ( J + 1 )*J*( J - 1 ) )/3 ) 130 CONTINUE * YY holds the exact result. On exit from CMVCH YT holds * the result computed by CMVCH. @@ -2749,7 +2750,7 @@ IC = 0 GO TO 10 END IF - CBEG = CMPLX( ( I - 500 )/1001.0, ( J - 500 )/1001.0 ) + CBEG = CMPLX( REAL( I - 500 )/1001.0, REAL( J - 500 )/1001.0 ) RETURN * * End of CBEG. diff --git a/CBLAS/testing/c_cblat3.f b/CBLAS/testing/c_cblat3.f index 3575e6d8f..b0c0d4972 100644 --- a/CBLAS/testing/c_cblat3.f +++ b/CBLAS/testing/c_cblat3.f @@ -221,14 +221,15 @@ N = MIN( 32, NMAX ) DO 100 J = 1, N DO 90 I = 1, N - AB( I, J ) = MAX( I - J + 1, 0 ) + AB( I, J ) = REAL( MAX( I - J + 1, 0 ) ) 90 CONTINUE - AB( J, NMAX + 1 ) = J - AB( 1, NMAX + J ) = J + AB( J, NMAX + 1 ) = REAL( J ) + AB( 1, NMAX + J ) = REAL( J ) C( J, 1 ) = ZERO 100 CONTINUE DO 110 J = 1, N - CC( J ) = J*( ( J + 1 )*J )/2 - ( ( J + 1 )*J*( J - 1 ) )/3 + CC( J ) = REAL( J*( ( J + 1 )*J )/2 - + $ ( ( J + 1 )*J*( J - 1 ) )/3 ) 110 CONTINUE * CC holds the exact result. On exit from CMMCH CT holds * the result computed by CMMCH. @@ -252,12 +253,12 @@ STOP END IF DO 120 J = 1, N - AB( J, NMAX + 1 ) = N - J + 1 - AB( 1, NMAX + J ) = N - J + 1 + AB( J, NMAX + 1 ) = REAL( N - J + 1 ) + AB( 1, NMAX + J ) = REAL( N - J + 1 ) 120 CONTINUE DO 130 J = 1, N - CC( N - J + 1 ) = J*( ( J + 1 )*J )/2 - - $ ( ( J + 1 )*J*( J - 1 ) )/3 + CC( N - J + 1 ) = REAL( J*( ( J + 1 )*J )/2 - + $ ( ( J + 1 )*J*( J - 1 ) )/3 ) 130 CONTINUE TRANSA = 'C' TRANSB = 'N' @@ -2823,7 +2824,7 @@ IC = 0 GO TO 10 END IF - CBEG = CMPLX( ( I - 500 )/1001.0, ( J - 500 )/1001.0 ) + CBEG = CMPLX( REAL( I - 500 )/1001.0, REAL( J - 500 )/1001.0 ) RETURN * * End of CBEG. @@ -2928,7 +2929,7 @@ IF( LDC.GT.NMAX ) $ GO TO 100 LCC = LDC*N - NULL = N.LE.0. + NULL = N.LE.0 * DO 90 IK = 1, NIDIM K = IDIM( IK ) diff --git a/CBLAS/testing/c_sblat1.f b/CBLAS/testing/c_sblat1.f index 44e24e5d0..805dfe83f 100644 --- a/CBLAS/testing/c_sblat1.f +++ b/CBLAS/testing/c_sblat1.f @@ -587,26 +587,26 @@ MWPN(5) = 3 MWPN(10) = 3 DO 160 I = 1, 5 - MWPX(I) = I - MWPY(I) = I - MWPTX(1,I) = I - MWPTY(1,I) = I - MWPTX(2,I) = I - MWPTY(2,I) = -I - MWPTX(3,I) = 6 - I - MWPTY(3,I) = I - 6 - MWPTX(4,I) = I - MWPTY(4,I) = -I - MWPTX(6,I) = 6 - I - MWPTY(6,I) = I - 6 - MWPTX(7,I) = -I - MWPTY(7,I) = I - MWPTX(8,I) = I - 6 - MWPTY(8,I) = 6 - I - MWPTX(9,I) = -I - MWPTY(9,I) = I - MWPTX(11,I) = I - 6 - MWPTY(11,I) = 6 - I + MWPX(I) = REAL( I ) + MWPY(I) = REAL( I ) + MWPTX(1,I) = REAL( I ) + MWPTY(1,I) = REAL( I ) + MWPTX(2,I) = REAL( I ) + MWPTY(2,I) = REAL( -I ) + MWPTX(3,I) = REAL( 6 - I ) + MWPTY(3,I) = REAL( I - 6 ) + MWPTX(4,I) = REAL( I ) + MWPTY(4,I) = REAL( -I ) + MWPTX(6,I) = REAL( 6 - I ) + MWPTY(6,I) = REAL( I - 6 ) + MWPTX(7,I) = REAL( -I ) + MWPTY(7,I) = REAL( I ) + MWPTX(8,I) = REAL( I - 6 ) + MWPTY(8,I) = REAL( 6 - I ) + MWPTX(9,I) = REAL( -I ) + MWPTY(9,I) = REAL( I ) + MWPTX(11,I) = REAL( I - 6 ) + MWPTY(11,I) = REAL( 6 - I ) 160 CONTINUE MWPTX(5,1) = 1 MWPTX(5,2) = 3 diff --git a/CBLAS/testing/c_sblat2.f b/CBLAS/testing/c_sblat2.f index 5e268895f..1ea2203d9 100644 --- a/CBLAS/testing/c_sblat2.f +++ b/CBLAS/testing/c_sblat2.f @@ -279,13 +279,14 @@ N = MIN( 32, NMAX ) DO 120 J = 1, N DO 110 I = 1, N - A( I, J ) = MAX( I - J + 1, 0 ) + A( I, J ) = REAL( MAX( I - J + 1, 0 ) ) 110 CONTINUE - X( J ) = J + X( J ) = REAL( J ) Y( J ) = ZERO 120 CONTINUE DO 130 J = 1, N - YY( J ) = J*( ( J + 1 )*J )/2 - ( ( J + 1 )*J*( J - 1 ) )/3 + YY( J ) = REAL( J*( ( J + 1 )*J )/2 - + $ ( ( J + 1 )*J*( J - 1 ) )/3 ) 130 CONTINUE * YY holds the exact result. On exit from SMVCH YT holds * the result computed by SMVCH. diff --git a/CBLAS/testing/c_sblat3.f b/CBLAS/testing/c_sblat3.f index b02491167..412a1fed5 100644 --- a/CBLAS/testing/c_sblat3.f +++ b/CBLAS/testing/c_sblat3.f @@ -223,14 +223,15 @@ N = MIN( 32, NMAX ) DO 100 J = 1, N DO 90 I = 1, N - AB( I, J ) = MAX( I - J + 1, 0 ) + AB( I, J ) = REAL( MAX( I - J + 1, 0 ) ) 90 CONTINUE - AB( J, NMAX + 1 ) = J - AB( 1, NMAX + J ) = J + AB( J, NMAX + 1 ) = REAL( J ) + AB( 1, NMAX + J ) = REAL( J ) C( J, 1 ) = ZERO 100 CONTINUE DO 110 J = 1, N - CC( J ) = J*( ( J + 1 )*J )/2 - ( ( J + 1 )*J*( J - 1 ) )/3 + CC( J ) = REAL( J*( ( J + 1 )*J )/2 - + $ ( ( J + 1 )*J*( J - 1 ) )/3 ) 110 CONTINUE * CC holds the exact result. On exit from SMMCH CT holds * the result computed by SMMCH. @@ -254,12 +255,12 @@ STOP END IF DO 120 J = 1, N - AB( J, NMAX + 1 ) = N - J + 1 - AB( 1, NMAX + J ) = N - J + 1 + AB( J, NMAX + 1 ) = REAL( N - J + 1 ) + AB( 1, NMAX + J ) = REAL( N - J + 1 ) 120 CONTINUE DO 130 J = 1, N - CC( N - J + 1 ) = J*( ( J + 1 )*J )/2 - - $ ( ( J + 1 )*J*( J - 1 ) )/3 + CC( N - J + 1 ) = REAL( J*( ( J + 1 )*J )/2 - + $ ( ( J + 1 )*J*( J - 1 ) )/3 ) 130 CONTINUE TRANSA = 'T' TRANSB = 'N' @@ -2581,7 +2582,7 @@ IC = 0 GO TO 10 END IF - SBEG = ( I - 500 )/1001.0 + SBEG = REAL( I - 500 )/1001.0 RETURN * * End of SBEG. diff --git a/CBLAS/testing/c_zblat2.f b/CBLAS/testing/c_zblat2.f index 5347864b7..165d4533c 100644 --- a/CBLAS/testing/c_zblat2.f +++ b/CBLAS/testing/c_zblat2.f @@ -2755,7 +2755,7 @@ IC = 0 GO TO 10 END IF - ZBEG = DCMPLX( ( I - 500 )/1001.0, ( J - 500 )/1001.0 ) + ZBEG = DCMPLX( ( I - 500 )/1001.0D0, ( J - 500 )/1001.0D0 ) RETURN * * End of ZBEG. diff --git a/CBLAS/testing/c_zblat3.f b/CBLAS/testing/c_zblat3.f index 8e5be3072..d4b0a67c4 100644 --- a/CBLAS/testing/c_zblat3.f +++ b/CBLAS/testing/c_zblat3.f @@ -2929,7 +2929,7 @@ IF( LDC.GT.NMAX ) $ GO TO 100 LCC = LDC*N - NULL = N.LE.0. + NULL = N.LE.0 * DO 90 IK = 1, NIDIM K = IDIM( IK ) diff --git a/TESTING/EIG/cchkbb.f b/TESTING/EIG/cchkbb.f index eca5e4d4d..6315b0a8d 100644 --- a/TESTING/EIG/cchkbb.f +++ b/TESTING/EIG/cchkbb.f @@ -568,7 +568,7 @@ GO TO 70 * 60 CONTINUE - ANORM = RTUNFL*MAX( M, N )*ULPINV + ANORM = RTUNFL*REAL( MAX( M, N ) )*ULPINV GO TO 70 * 70 CONTINUE diff --git a/TESTING/EIG/cchkbd.f b/TESTING/EIG/cchkbd.f index 84f5d7d9f..7b3ba3bdd 100644 --- a/TESTING/EIG/cchkbd.f +++ b/TESTING/EIG/cchkbd.f @@ -562,7 +562,7 @@ M = MVAL( JSIZE ) N = NVAL( JSIZE ) MNMIN = MIN( M, N ) - AMNINV = ONE / MAX( M, N, 1 ) + AMNINV = ONE / REAL( MAX( M, N, 1 ) ) * IF( NSIZES.NE.1 ) THEN MTYPES = MIN( MAXTYP, NTYPES ) @@ -619,7 +619,7 @@ GO TO 70 * 60 CONTINUE - ANORM = RTUNFL*MAX( M, N )*ULPINV + ANORM = RTUNFL*REAL( MAX( M, N ) )*ULPINV GO TO 70 * 70 CONTINUE diff --git a/TESTING/EIG/cchkgg.f b/TESTING/EIG/cchkgg.f index 5bd766eae..f21c37949 100644 --- a/TESTING/EIG/cchkgg.f +++ b/TESTING/EIG/cchkgg.f @@ -657,7 +657,7 @@ N = NN( JSIZE ) N1 = MAX( 1, N ) RMAGN( 2 ) = SAFMAX*ULP / REAL( N1 ) - RMAGN( 3 ) = SAFMIN*ULPINV*N1 + RMAGN( 3 ) = SAFMIN*ULPINV*REAL( N1 ) * IF( NSIZES.NE.1 ) THEN MTYPES = MIN( MAXTYP, NTYPES ) diff --git a/TESTING/EIG/cchkhb.f b/TESTING/EIG/cchkhb.f index 5d6a4780b..13da2c0bd 100644 --- a/TESTING/EIG/cchkhb.f +++ b/TESTING/EIG/cchkhb.f @@ -489,7 +489,7 @@ GO TO 70 * 60 CONTINUE - ANORM = RTUNFL*N*ULPINV + ANORM = RTUNFL*REAL( N )*ULPINV GO TO 70 * 70 CONTINUE diff --git a/TESTING/EIG/cchkhb2stg.f b/TESTING/EIG/cchkhb2stg.f index a8f7e5895..5bc2da44a 100644 --- a/TESTING/EIG/cchkhb2stg.f +++ b/TESTING/EIG/cchkhb2stg.f @@ -533,7 +533,7 @@ GO TO 70 * 60 CONTINUE - ANORM = RTUNFL*N*ULPINV + ANORM = RTUNFL*REAL( N )*ULPINV GO TO 70 * 70 CONTINUE diff --git a/TESTING/EIG/cchkhs.f b/TESTING/EIG/cchkhs.f index 07eed07e1..71af6f970 100644 --- a/TESTING/EIG/cchkhs.f +++ b/TESTING/EIG/cchkhs.f @@ -612,7 +612,7 @@ GO TO 70 * 60 CONTINUE - ANORM = RTUNFL*N*ULPINV + ANORM = RTUNFL*REAL( N )*ULPINV GO TO 70 * 70 CONTINUE