Skip to content

Commit d2a547a

Browse files
angschangsch
committed
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
1 parent b9b3c78 commit d2a547a

File tree

9 files changed

+29
-27
lines changed

9 files changed

+29
-27
lines changed

‎TESTING/EIG/serrbd.f

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -111,7 +111,7 @@ SUBROUTINE SERRBD( PATH, NUNIT )
111111
*
112112
DO 20 J = 1, NMAX
113113
DO 10 I = 1, NMAX
114-
A( I, J ) = 1.D0 / REAL( I+J )
114+
A( I, J ) = 1. / REAL( I+J )
115115
10 CONTINUE
116116
20 CONTINUE
117117
OK = .TRUE.

‎TESTING/LIN/cerrrfp.f

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -89,10 +89,10 @@ SUBROUTINE CERRRFP( NUNIT )
8989
*
9090
NOUT = NUNIT
9191
OK = .TRUE.
92-
A( 1, 1 ) = CMPLX( 1.D0 , 1.D0 )
93-
B( 1, 1 ) = CMPLX( 1.D0 , 1.D0 )
94-
ALPHA = CMPLX( 1.D0 , 1.D0 )
95-
BETA = CMPLX( 1.D0 , 1.D0 )
92+
A( 1, 1 ) = CMPLX( 1.0 , 1.0 )
93+
B( 1, 1 ) = CMPLX( 1.0 , 1.0 )
94+
ALPHA = CMPLX( 1.0 , 1.0 )
95+
BETA = CMPLX( 1.0 , 1.0 )
9696
*
9797
SRNAMT = 'CPFTRF'
9898
INFOT = 1

‎TESTING/LIN/chet01_aa.f

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -159,7 +159,7 @@ SUBROUTINE CHET01_AA( UPLO, N, A, LDA, AFAC, LDAFAC, IPIV, C,
159159
EXTERNAL CLASET, CLAVHE
160160
* ..
161161
* .. Intrinsic Functions ..
162-
INTRINSIC DBLE
162+
INTRINSIC REAL
163163
* ..
164164
* .. Executable Statements ..
165165
*
@@ -255,7 +255,7 @@ SUBROUTINE CHET01_AA( UPLO, N, A, LDA, AFAC, LDAFAC, IPIV, C,
255255
IF( RESID.NE.ZERO )
256256
$ RESID = ONE / EPS
257257
ELSE
258-
RESID = ( ( RESID / DBLE( N ) ) / ANORM ) / EPS
258+
RESID = ( ( RESID / REAL( N ) ) / ANORM ) / EPS
259259
END IF
260260
*
261261
RETURN

‎TESTING/LIN/csyt01_aa.f

Lines changed: 5 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -141,9 +141,10 @@ SUBROUTINE CSYT01_AA( UPLO, N, A, LDA, AFAC, LDAFAC, IPIV, C,
141141
*
142142
* .. Parameters ..
143143
REAL ZERO, ONE
144-
PARAMETER ( ZERO = 0.0D+0, ONE = 1.0D+0 )
144+
PARAMETER ( ZERO = 0.0E+0, ONE = 1.0E+0 )
145145
COMPLEX CZERO, CONE
146-
PARAMETER ( CZERO = 0.0E+0, CONE = 1.0E+0 )
146+
PARAMETER ( CZERO = ( 0.0E+0, 0.0E+0 ),
147+
$ CONE = ( 1.0E+0, 0.0E+0 ) )
147148
* ..
148149
* .. Local Scalars ..
149150
INTEGER I, J
@@ -158,7 +159,7 @@ SUBROUTINE CSYT01_AA( UPLO, N, A, LDA, AFAC, LDAFAC, IPIV, C,
158159
EXTERNAL CLASET, CLAVSY
159160
* ..
160161
* .. Intrinsic Functions ..
161-
INTRINSIC DBLE
162+
INTRINSIC REAL
162163
* ..
163164
* .. Executable Statements ..
164165
*
@@ -250,7 +251,7 @@ SUBROUTINE CSYT01_AA( UPLO, N, A, LDA, AFAC, LDAFAC, IPIV, C,
250251
IF( RESID.NE.ZERO )
251252
$ RESID = ONE / EPS
252253
ELSE
253-
RESID = ( ( RESID / DBLE( N ) ) / ANORM ) / EPS
254+
RESID = ( ( RESID / REAL( N ) ) / ANORM ) / EPS
254255
END IF
255256
*
256257
RETURN

‎TESTING/LIN/serrlqt.f

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -101,11 +101,11 @@ SUBROUTINE SERRLQT( PATH, NUNIT )
101101
*
102102
DO J = 1, NMAX
103103
DO I = 1, NMAX
104-
A( I, J ) = 1.D0 / REAL( I+J )
105-
C( I, J ) = 1.D0 / REAL( I+J )
106-
T( I, J ) = 1.D0 / REAL( I+J )
104+
A( I, J ) = 1. / REAL( I+J )
105+
C( I, J ) = 1. / REAL( I+J )
106+
T( I, J ) = 1. / REAL( I+J )
107107
END DO
108-
W( J ) = 0.D0
108+
W( J ) = 0.
109109
END DO
110110
OK = .TRUE.
111111
*

‎TESTING/LIN/serrlqtp.f

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -101,11 +101,11 @@ SUBROUTINE SERRLQTP( PATH, NUNIT )
101101
*
102102
DO J = 1, NMAX
103103
DO I = 1, NMAX
104-
A( I, J ) = 1.D0 / REAL( I+J )
105-
C( I, J ) = 1.D0 / REAL( I+J )
106-
T( I, J ) = 1.D0 / REAL( I+J )
104+
A( I, J ) = 1. / REAL( I+J )
105+
C( I, J ) = 1. / REAL( I+J )
106+
T( I, J ) = 1. / REAL( I+J )
107107
END DO
108-
W( J ) = 0.0
108+
W( J ) = 0.
109109
END DO
110110
OK = .TRUE.
111111
*

‎TESTING/LIN/serrtsqr.f

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -101,11 +101,11 @@ SUBROUTINE SERRTSQR( PATH, NUNIT )
101101
*
102102
DO J = 1, NMAX
103103
DO I = 1, NMAX
104-
A( I, J ) = 1.D0 / REAL( I+J )
105-
C( I, J ) = 1.D0 / REAL( I+J )
106-
T( I, J ) = 1.D0 / REAL( I+J )
104+
A( I, J ) = 1. / REAL( I+J )
105+
C( I, J ) = 1. / REAL( I+J )
106+
T( I, J ) = 1. / REAL( I+J )
107107
END DO
108-
W( J ) = 0.D0
108+
W( J ) = 0.
109109
END DO
110110
OK = .TRUE.
111111
*

‎TESTING/LIN/ssyt01_aa.f

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -156,7 +156,7 @@ SUBROUTINE SSYT01_AA( UPLO, N, A, LDA, AFAC, LDAFAC, IPIV, C,
156156
EXTERNAL SLASET, SLAVSY, SSWAP, STRMM, SLACPY
157157
* ..
158158
* .. Intrinsic Functions ..
159-
INTRINSIC DBLE
159+
INTRINSIC REAL
160160
* ..
161161
* .. Executable Statements ..
162162
*
@@ -248,7 +248,7 @@ SUBROUTINE SSYT01_AA( UPLO, N, A, LDA, AFAC, LDAFAC, IPIV, C,
248248
IF( RESID.NE.ZERO )
249249
$ RESID = ONE / EPS
250250
ELSE
251-
RESID = ( ( RESID / DBLE( N ) ) / ANORM ) / EPS
251+
RESID = ( ( RESID / REAL( N ) ) / ANORM ) / EPS
252252
END IF
253253
*
254254
RETURN

‎TESTING/LIN/zsyt01_aa.f

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -142,8 +142,9 @@ SUBROUTINE ZSYT01_AA( UPLO, N, A, LDA, AFAC, LDAFAC, IPIV, C,
142142
* .. Parameters ..
143143
DOUBLE PRECISION ZERO, ONE
144144
PARAMETER ( ZERO = 0.0D+0, ONE = 1.0D+0 )
145-
COMPLEX*16 CZERO, CONE
146-
PARAMETER ( CZERO = 0.0E+0, CONE = 1.0E+0 )
145+
COMPLEX*16 CZERO, CONE
146+
PARAMETER ( CZERO = ( 0.0D+0, 0.0D+0 ),
147+
$ CONE = ( 1.0D+0, 0.0D+0 ) )
147148
* ..
148149
* .. Local Scalars ..
149150
INTEGER I, J

0 commit comments

Comments
 (0)