Skip to content
9 changes: 5 additions & 4 deletions BLAS/TESTING/cblat2.f
Original file line number Diff line number Diff line change
Expand Up @@ -292,13 +292,14 @@ PROGRAM CBLAT2
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.
Expand Down Expand Up @@ -3212,7 +3213,7 @@ COMPLEX FUNCTION CBEG( RESET )
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
Expand Down
19 changes: 10 additions & 9 deletions BLAS/TESTING/cblat3.f
Original file line number Diff line number Diff line change
Expand Up @@ -243,14 +243,15 @@ PROGRAM CBLAT3
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.
Expand All @@ -274,12 +275,12 @@ PROGRAM CBLAT3
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'
Expand Down Expand Up @@ -3581,7 +3582,7 @@ COMPLEX FUNCTION CBEG( RESET )
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
Expand Down
40 changes: 20 additions & 20 deletions BLAS/TESTING/sblat1.f
Original file line number Diff line number Diff line change
Expand Up @@ -916,26 +916,26 @@ SUBROUTINE CHECK3(SFAC)
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
Expand Down
7 changes: 4 additions & 3 deletions BLAS/TESTING/sblat2.f
Original file line number Diff line number Diff line change
Expand Up @@ -294,13 +294,14 @@ PROGRAM SBLAT2
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.
Expand Down
19 changes: 10 additions & 9 deletions BLAS/TESTING/sblat3.f
Original file line number Diff line number Diff line change
Expand Up @@ -241,14 +241,15 @@ PROGRAM SBLAT3
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.
Expand All @@ -272,12 +273,12 @@ PROGRAM SBLAT3
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'
Expand Down Expand Up @@ -3058,7 +3059,7 @@ REAL FUNCTION SBEG( RESET )
IC = 0
GO TO 10
END IF
SBEG = ( I - 500 )/1001.0
SBEG = REAL( I - 500 )/1001.0
RETURN
*
* End of SBEG
Expand Down
9 changes: 5 additions & 4 deletions CBLAS/testing/c_cblat2.f
Original file line number Diff line number Diff line change
Expand Up @@ -269,13 +269,14 @@ PROGRAM CBLAT2
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.
Expand Down Expand Up @@ -2749,7 +2750,7 @@ COMPLEX FUNCTION CBEG( RESET )
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.
Expand Down
21 changes: 11 additions & 10 deletions CBLAS/testing/c_cblat3.f
Original file line number Diff line number Diff line change
Expand Up @@ -221,14 +221,15 @@ PROGRAM CBLAT3
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.
Expand All @@ -252,12 +253,12 @@ PROGRAM CBLAT3
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'
Expand Down Expand Up @@ -2823,7 +2824,7 @@ COMPLEX FUNCTION CBEG( RESET )
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.
Expand Down Expand Up @@ -2928,7 +2929,7 @@ SUBROUTINE CCHK6( SNAME, EPS, THRESH, NOUT, NTRA, TRACE, REWI,
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 )
Expand Down
40 changes: 20 additions & 20 deletions CBLAS/testing/c_sblat1.f
Original file line number Diff line number Diff line change
Expand Up @@ -587,26 +587,26 @@ SUBROUTINE CHECK3(SFAC)
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
Expand Down
7 changes: 4 additions & 3 deletions CBLAS/testing/c_sblat2.f
Original file line number Diff line number Diff line change
Expand Up @@ -279,13 +279,14 @@ PROGRAM SBLAT2
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.
Expand Down
19 changes: 10 additions & 9 deletions CBLAS/testing/c_sblat3.f
Original file line number Diff line number Diff line change
Expand Up @@ -223,14 +223,15 @@ PROGRAM SBLAT3
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.
Expand All @@ -254,12 +255,12 @@ PROGRAM SBLAT3
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'
Expand Down Expand Up @@ -2581,7 +2582,7 @@ REAL FUNCTION SBEG( RESET )
IC = 0
GO TO 10
END IF
SBEG = ( I - 500 )/1001.0
SBEG = REAL( I - 500 )/1001.0
RETURN
*
* End of SBEG.
Expand Down
2 changes: 1 addition & 1 deletion CBLAS/testing/c_zblat2.f
Original file line number Diff line number Diff line change
Expand Up @@ -2755,7 +2755,7 @@ COMPLEX*16 FUNCTION ZBEG( RESET )
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.
Expand Down
2 changes: 1 addition & 1 deletion CBLAS/testing/c_zblat3.f
Original file line number Diff line number Diff line change
Expand Up @@ -2929,7 +2929,7 @@ SUBROUTINE ZCHK6( SNAME, EPS, THRESH, NOUT, NTRA, TRACE, REWI,
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 )
Expand Down
Loading
Loading