phynexis-deps/lapack
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Remove redundant DBLE()/REAL() type conversions in TESTING

Browse Source 
In ddrvst.f, ddrvst2stg.f, and dlatm4.f, DBLE() was applied to
DOUBLE PRECISION array elements assigned to DOUBLE PRECISION
variables ― a no-op conversion. Similarly, in sdrvst.f,
sdrvst2stg.f, and slatm4.f, REAL() was applied to REAL array
elements assigned to REAL variables.

Remove these 168 redundant calls. Integer-to-real conversions
(e.g., DBLE(N), DBLE(KLEN-1), REAL(MAX(1,N))) are retained as
they are semantically necessary.
This commit is contained in:
NAKATA Maho
2026-03-24 11:37:16 +09:00
parent 7f18edae13
commit 8b1494be8c
6 changed files with 168 additions and 168 deletions
Download Patch File Download Diff File Expand all files Collapse all files
TESTING/EIG/ddrvst.f
+40 -40
View File
@@ -773,10 +773,10 @@ c LIWEDC = 12
IF( JTYPE.LE.7 ) THEN
NTEST = 1
DO 120 I = 1, N
D1( I ) = DBLE( A( I, I ) )
D1( I ) = A( I, I )
120 CONTINUE
DO 130 I = 1, N - 1
D2( I ) = DBLE( A( I+1, I ) )
D2( I ) = A( I+1, I )
130 CONTINUE
SRNAMT = 'DSTEV'
CALL DSTEV( 'V', N, D1, D2, Z, LDU, WORK, IINFO )
@@ -797,17 +797,17 @@ c LIWEDC = 12
* Do tests 1 and 2.
*
DO 140 I = 1, N
D3( I ) = DBLE( A( I, I ) )
D3( I ) = A( I, I )
140 CONTINUE
DO 150 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
150 CONTINUE
CALL DSTT21( N, 0, D3, D4, D1, D2, Z, LDU, WORK,
$ RESULT( 1 ) )
*
NTEST = 3
DO 160 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
160 CONTINUE
SRNAMT = 'DSTEV'
CALL DSTEV( 'N', N, D3, D4, Z, LDU, WORK, IINFO )
@@ -839,10 +839,10 @@ c LIWEDC = 12
NTEST = 4
DO 190 I = 1, N
EVEIGS( I ) = D3( I )
D1( I ) = DBLE( A( I, I ) )
D1( I ) = A( I, I )
190 CONTINUE
DO 200 I = 1, N - 1
D2( I ) = DBLE( A( I+1, I ) )
D2( I ) = A( I+1, I )
200 CONTINUE
SRNAMT = 'DSTEVX'
CALL DSTEVX( 'V', 'A', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -870,17 +870,17 @@ c LIWEDC = 12
* Do tests 4 and 5.
*
DO 210 I = 1, N
D3( I ) = DBLE( A( I, I ) )
D3( I ) = A( I, I )
210 CONTINUE
DO 220 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
220 CONTINUE
CALL DSTT21( N, 0, D3, D4, WA1, D2, Z, LDU, WORK,
$ RESULT( 4 ) )
*
NTEST = 6
DO 230 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
230 CONTINUE
SRNAMT = 'DSTEVX'
CALL DSTEVX( 'N', 'A', N, D3, D4, VL, VU, IL, IU, ABSTOL,
@@ -914,10 +914,10 @@ c LIWEDC = 12
*
NTEST = 7
DO 260 I = 1, N
D1( I ) = DBLE( A( I, I ) )
D1( I ) = A( I, I )
260 CONTINUE
DO 270 I = 1, N - 1
D2( I ) = DBLE( A( I+1, I ) )
D2( I ) = A( I+1, I )
270 CONTINUE
SRNAMT = 'DSTEVR'
CALL DSTEVR( 'V', 'A', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -944,17 +944,17 @@ c LIWEDC = 12
* Do tests 7 and 8.
*
DO 280 I = 1, N
D3( I ) = DBLE( A( I, I ) )
D3( I ) = A( I, I )
280 CONTINUE
DO 290 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
290 CONTINUE
CALL DSTT21( N, 0, D3, D4, WA1, D2, Z, LDU, WORK,
$ RESULT( 7 ) )
*
NTEST = 9
DO 300 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
300 CONTINUE
SRNAMT = 'DSTEVR'
CALL DSTEVR( 'N', 'A', N, D3, D4, VL, VU, IL, IU, ABSTOL,
@@ -989,10 +989,10 @@ c LIWEDC = 12
*
NTEST = 10
DO 330 I = 1, N
D1( I ) = DBLE( A( I, I ) )
D1( I ) = A( I, I )
330 CONTINUE
DO 340 I = 1, N - 1
D2( I ) = DBLE( A( I+1, I ) )
D2( I ) = A( I+1, I )
340 CONTINUE
SRNAMT = 'DSTEVX'
CALL DSTEVX( 'V', 'I', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -1015,10 +1015,10 @@ c LIWEDC = 12
* Do tests 10 and 11.
*
DO 350 I = 1, N
D3( I ) = DBLE( A( I, I ) )
D3( I ) = A( I, I )
350 CONTINUE
DO 360 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
360 CONTINUE
CALL DSTT22( N, M2, 0, D3, D4, WA2, D2, Z, LDU, WORK,
$ MAX( 1, M2 ), RESULT( 10 ) )
@@ -1026,7 +1026,7 @@ c LIWEDC = 12
*
NTEST = 12
DO 370 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
370 CONTINUE
SRNAMT = 'DSTEVX'
CALL DSTEVX( 'N', 'I', N, D3, D4, VL, VU, IL, IU, ABSTOL,
@@ -1076,10 +1076,10 @@ c LIWEDC = 12
END IF
*
DO 390 I = 1, N
D1( I ) = DBLE( A( I, I ) )
D1( I ) = A( I, I )
390 CONTINUE
DO 400 I = 1, N - 1
D2( I ) = DBLE( A( I+1, I ) )
D2( I ) = A( I+1, I )
400 CONTINUE
SRNAMT = 'DSTEVX'
CALL DSTEVX( 'V', 'V', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -1109,17 +1109,17 @@ c LIWEDC = 12
* Do tests 13 and 14.
*
DO 410 I = 1, N
D3( I ) = DBLE( A( I, I ) )
D3( I ) = A( I, I )
410 CONTINUE
DO 420 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
420 CONTINUE
CALL DSTT22( N, M2, 0, D3, D4, WA2, D2, Z, LDU, WORK,
$ MAX( 1, M2 ), RESULT( 13 ) )
*
NTEST = 15
DO 430 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
430 CONTINUE
SRNAMT = 'DSTEVX'
CALL DSTEVX( 'N', 'V', N, D3, D4, VL, VU, IL, IU, ABSTOL,
@@ -1147,10 +1147,10 @@ c LIWEDC = 12
*
NTEST = 16
DO 450 I = 1, N
D1( I ) = DBLE( A( I, I ) )
D1( I ) = A( I, I )
450 CONTINUE
DO 460 I = 1, N - 1
D2( I ) = DBLE( A( I+1, I ) )
D2( I ) = A( I+1, I )
460 CONTINUE
SRNAMT = 'DSTEVD'
CALL DSTEVD( 'V', N, D1, D2, Z, LDU, WORK, LWEDC, IWORK,
@@ -1172,17 +1172,17 @@ c LIWEDC = 12
* Do tests 16 and 17.
*
DO 470 I = 1, N
D3( I ) = DBLE( A( I, I ) )
D3( I ) = A( I, I )
470 CONTINUE
DO 480 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
480 CONTINUE
CALL DSTT21( N, 0, D3, D4, D1, D2, Z, LDU, WORK,
$ RESULT( 16 ) )
*
NTEST = 18
DO 490 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
490 CONTINUE
SRNAMT = 'DSTEVD'
CALL DSTEVD( 'N', N, D3, D4, Z, LDU, WORK, LWEDC, IWORK,
@@ -1215,10 +1215,10 @@ c LIWEDC = 12
*
NTEST = 19
DO 520 I = 1, N
D1( I ) = DBLE( A( I, I ) )
D1( I ) = A( I, I )
520 CONTINUE
DO 530 I = 1, N - 1
D2( I ) = DBLE( A( I+1, I ) )
D2( I ) = A( I+1, I )
530 CONTINUE
SRNAMT = 'DSTEVR'
CALL DSTEVR( 'V', 'I', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -1241,10 +1241,10 @@ c LIWEDC = 12
* DO tests 19 and 20.
*
DO 540 I = 1, N
D3( I ) = DBLE( A( I, I ) )
D3( I ) = A( I, I )
540 CONTINUE
DO 550 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
550 CONTINUE
CALL DSTT22( N, M2, 0, D3, D4, WA2, D2, Z, LDU, WORK,
$ MAX( 1, M2 ), RESULT( 19 ) )
@@ -1252,7 +1252,7 @@ c LIWEDC = 12
*
NTEST = 21
DO 560 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
560 CONTINUE
SRNAMT = 'DSTEVR'
CALL DSTEVR( 'N', 'I', N, D3, D4, VL, VU, IL, IU, ABSTOL,
@@ -1302,10 +1302,10 @@ c LIWEDC = 12
END IF
*
DO 580 I = 1, N
D1( I ) = DBLE( A( I, I ) )
D1( I ) = A( I, I )
580 CONTINUE
DO 590 I = 1, N - 1
D2( I ) = DBLE( A( I+1, I ) )
D2( I ) = A( I+1, I )
590 CONTINUE
SRNAMT = 'DSTEVR'
CALL DSTEVR( 'V', 'V', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -1335,17 +1335,17 @@ c LIWEDC = 12
* Do tests 22 and 23.
*
DO 600 I = 1, N
D3( I ) = DBLE( A( I, I ) )
D3( I ) = A( I, I )
600 CONTINUE
DO 610 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
610 CONTINUE
CALL DSTT22( N, M2, 0, D3, D4, WA2, D2, Z, LDU, WORK,
$ MAX( 1, M2 ), RESULT( 22 ) )
*
NTEST = 24
DO 620 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
620 CONTINUE
SRNAMT = 'DSTEVR'
CALL DSTEVR( 'N', 'V', N, D3, D4, VL, VU, IL, IU, ABSTOL,
TESTING/EIG/ddrvst2stg.f
+40 -40
View File
@@ -776,10 +776,10 @@ c LIWEDC = 12
IF( JTYPE.LE.7 ) THEN
NTEST = 1
DO 120 I = 1, N
D1( I ) = DBLE( A( I, I ) )
D1( I ) = A( I, I )
120 CONTINUE
DO 130 I = 1, N - 1
D2( I ) = DBLE( A( I+1, I ) )
D2( I ) = A( I+1, I )
130 CONTINUE
SRNAMT = 'DSTEV'
CALL DSTEV( 'V', N, D1, D2, Z, LDU, WORK, IINFO )
@@ -800,17 +800,17 @@ c LIWEDC = 12
* Do tests 1 and 2.
*
DO 140 I = 1, N
D3( I ) = DBLE( A( I, I ) )
D3( I ) = A( I, I )
140 CONTINUE
DO 150 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
150 CONTINUE
CALL DSTT21( N, 0, D3, D4, D1, D2, Z, LDU, WORK,
$ RESULT( 1 ) )
*
NTEST = 3
DO 160 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
160 CONTINUE
SRNAMT = 'DSTEV'
CALL DSTEV( 'N', N, D3, D4, Z, LDU, WORK, IINFO )
@@ -842,10 +842,10 @@ c LIWEDC = 12
NTEST = 4
DO 190 I = 1, N
EVEIGS( I ) = D3( I )
D1( I ) = DBLE( A( I, I ) )
D1( I ) = A( I, I )
190 CONTINUE
DO 200 I = 1, N - 1
D2( I ) = DBLE( A( I+1, I ) )
D2( I ) = A( I+1, I )
200 CONTINUE
SRNAMT = 'DSTEVX'
CALL DSTEVX( 'V', 'A', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -873,17 +873,17 @@ c LIWEDC = 12
* Do tests 4 and 5.
*
DO 210 I = 1, N
D3( I ) = DBLE( A( I, I ) )
D3( I ) = A( I, I )
210 CONTINUE
DO 220 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
220 CONTINUE
CALL DSTT21( N, 0, D3, D4, WA1, D2, Z, LDU, WORK,
$ RESULT( 4 ) )
*
NTEST = 6
DO 230 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
230 CONTINUE
SRNAMT = 'DSTEVX'
CALL DSTEVX( 'N', 'A', N, D3, D4, VL, VU, IL, IU, ABSTOL,
@@ -917,10 +917,10 @@ c LIWEDC = 12
*
NTEST = 7
DO 260 I = 1, N
D1( I ) = DBLE( A( I, I ) )
D1( I ) = A( I, I )
260 CONTINUE
DO 270 I = 1, N - 1
D2( I ) = DBLE( A( I+1, I ) )
D2( I ) = A( I+1, I )
270 CONTINUE
SRNAMT = 'DSTEVR'
CALL DSTEVR( 'V', 'A', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -947,17 +947,17 @@ c LIWEDC = 12
* Do tests 7 and 8.
*
DO 280 I = 1, N
D3( I ) = DBLE( A( I, I ) )
D3( I ) = A( I, I )
280 CONTINUE
DO 290 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
290 CONTINUE
CALL DSTT21( N, 0, D3, D4, WA1, D2, Z, LDU, WORK,
$ RESULT( 7 ) )
*
NTEST = 9
DO 300 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
300 CONTINUE
SRNAMT = 'DSTEVR'
CALL DSTEVR( 'N', 'A', N, D3, D4, VL, VU, IL, IU, ABSTOL,
@@ -992,10 +992,10 @@ c LIWEDC = 12
*
NTEST = 10
DO 330 I = 1, N
D1( I ) = DBLE( A( I, I ) )
D1( I ) = A( I, I )
330 CONTINUE
DO 340 I = 1, N - 1
D2( I ) = DBLE( A( I+1, I ) )
D2( I ) = A( I+1, I )
340 CONTINUE
SRNAMT = 'DSTEVX'
CALL DSTEVX( 'V', 'I', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -1018,10 +1018,10 @@ c LIWEDC = 12
* Do tests 10 and 11.
*
DO 350 I = 1, N
D3( I ) = DBLE( A( I, I ) )
D3( I ) = A( I, I )
350 CONTINUE
DO 360 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
360 CONTINUE
CALL DSTT22( N, M2, 0, D3, D4, WA2, D2, Z, LDU, WORK,
$ MAX( 1, M2 ), RESULT( 10 ) )
@@ -1029,7 +1029,7 @@ c LIWEDC = 12
*
NTEST = 12
DO 370 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
370 CONTINUE
SRNAMT = 'DSTEVX'
CALL DSTEVX( 'N', 'I', N, D3, D4, VL, VU, IL, IU, ABSTOL,
@@ -1079,10 +1079,10 @@ c LIWEDC = 12
END IF
*
DO 390 I = 1, N
D1( I ) = DBLE( A( I, I ) )
D1( I ) = A( I, I )
390 CONTINUE
DO 400 I = 1, N - 1
D2( I ) = DBLE( A( I+1, I ) )
D2( I ) = A( I+1, I )
400 CONTINUE
SRNAMT = 'DSTEVX'
CALL DSTEVX( 'V', 'V', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -1112,17 +1112,17 @@ c LIWEDC = 12
* Do tests 13 and 14.
*
DO 410 I = 1, N
D3( I ) = DBLE( A( I, I ) )
D3( I ) = A( I, I )
410 CONTINUE
DO 420 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
420 CONTINUE
CALL DSTT22( N, M2, 0, D3, D4, WA2, D2, Z, LDU, WORK,
$ MAX( 1, M2 ), RESULT( 13 ) )
*
NTEST = 15
DO 430 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
430 CONTINUE
SRNAMT = 'DSTEVX'
CALL DSTEVX( 'N', 'V', N, D3, D4, VL, VU, IL, IU, ABSTOL,
@@ -1150,10 +1150,10 @@ c LIWEDC = 12
*
NTEST = 16
DO 450 I = 1, N
D1( I ) = DBLE( A( I, I ) )
D1( I ) = A( I, I )
450 CONTINUE
DO 460 I = 1, N - 1
D2( I ) = DBLE( A( I+1, I ) )
D2( I ) = A( I+1, I )
460 CONTINUE
SRNAMT = 'DSTEVD'
CALL DSTEVD( 'V', N, D1, D2, Z, LDU, WORK, LWEDC, IWORK,
@@ -1175,17 +1175,17 @@ c LIWEDC = 12
* Do tests 16 and 17.
*
DO 470 I = 1, N
D3( I ) = DBLE( A( I, I ) )
D3( I ) = A( I, I )
470 CONTINUE
DO 480 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
480 CONTINUE
CALL DSTT21( N, 0, D3, D4, D1, D2, Z, LDU, WORK,
$ RESULT( 16 ) )
*
NTEST = 18
DO 490 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
490 CONTINUE
SRNAMT = 'DSTEVD'
CALL DSTEVD( 'N', N, D3, D4, Z, LDU, WORK, LWEDC, IWORK,
@@ -1218,10 +1218,10 @@ c LIWEDC = 12
*
NTEST = 19
DO 520 I = 1, N
D1( I ) = DBLE( A( I, I ) )
D1( I ) = A( I, I )
520 CONTINUE
DO 530 I = 1, N - 1
D2( I ) = DBLE( A( I+1, I ) )
D2( I ) = A( I+1, I )
530 CONTINUE
SRNAMT = 'DSTEVR'
CALL DSTEVR( 'V', 'I', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -1244,10 +1244,10 @@ c LIWEDC = 12
* DO tests 19 and 20.
*
DO 540 I = 1, N
D3( I ) = DBLE( A( I, I ) )
D3( I ) = A( I, I )
540 CONTINUE
DO 550 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
550 CONTINUE
CALL DSTT22( N, M2, 0, D3, D4, WA2, D2, Z, LDU, WORK,
$ MAX( 1, M2 ), RESULT( 19 ) )
@@ -1255,7 +1255,7 @@ c LIWEDC = 12
*
NTEST = 21
DO 560 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
560 CONTINUE
SRNAMT = 'DSTEVR'
CALL DSTEVR( 'N', 'I', N, D3, D4, VL, VU, IL, IU, ABSTOL,
@@ -1305,10 +1305,10 @@ c LIWEDC = 12
END IF
*
DO 580 I = 1, N
D1( I ) = DBLE( A( I, I ) )
D1( I ) = A( I, I )
580 CONTINUE
DO 590 I = 1, N - 1
D2( I ) = DBLE( A( I+1, I ) )
D2( I ) = A( I+1, I )
590 CONTINUE
SRNAMT = 'DSTEVR'
CALL DSTEVR( 'V', 'V', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -1338,17 +1338,17 @@ c LIWEDC = 12
* Do tests 22 and 23.
*
DO 600 I = 1, N
D3( I ) = DBLE( A( I, I ) )
D3( I ) = A( I, I )
600 CONTINUE
DO 610 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
610 CONTINUE
CALL DSTT22( N, M2, 0, D3, D4, WA2, D2, Z, LDU, WORK,
$ MAX( 1, M2 ), RESULT( 22 ) )
*
NTEST = 24
DO 620 I = 1, N - 1
D4( I ) = DBLE( A( I+1, I ) )
D4( I ) = A( I+1, I )
620 CONTINUE
SRNAMT = 'DSTEVR'
CALL DSTEVR( 'N', 'V', N, D3, D4, VL, VU, IL, IU, ABSTOL,
TESTING/EIG/dlatm4.f
+4 -4
View File
@@ -343,10 +343,10 @@
* Scale by AMAGN
*
DO 230 JD = KBEG, KEND
A( JD, JD ) = AMAGN*DBLE( A( JD, JD ) )
A( JD, JD ) = AMAGN*A( JD, JD )
230 CONTINUE
DO 240 JD = ISDB, ISDE
A( JD+1, JD ) = AMAGN*DBLE( A( JD+1, JD ) )
A( JD+1, JD ) = AMAGN*A( JD+1, JD )
240 CONTINUE
*
* If ISIGN = 1 or 2, assign random signs to diagonal and
@@ -354,13 +354,13 @@
*
IF( ISIGN.GT.0 ) THEN
DO 250 JD = KBEG, KEND
IF( DBLE( A( JD, JD ) ).NE.ZERO ) THEN
IF( A( JD, JD ).NE.ZERO ) THEN
IF( DLARAN( ISEED ).GT.HALF )
$ A( JD, JD ) = -A( JD, JD )
END IF
250 CONTINUE
DO 260 JD = ISDB, ISDE
IF( DBLE( A( JD+1, JD ) ).NE.ZERO ) THEN
IF( A( JD+1, JD ).NE.ZERO ) THEN
IF( DLARAN( ISEED ).GT.HALF )
$ A( JD+1, JD ) = -A( JD+1, JD )
END IF
TESTING/EIG/sdrvst.f
+40 -40
View File
@@ -773,10 +773,10 @@ c LIWEDC = 12
IF( JTYPE.LE.7 ) THEN
NTEST = 1
DO 120 I = 1, N
D1( I ) = REAL( A( I, I ) )
D1( I ) = A( I, I )
120 CONTINUE
DO 130 I = 1, N - 1
D2( I ) = REAL( A( I+1, I ) )
D2( I ) = A( I+1, I )
130 CONTINUE
SRNAMT = 'SSTEV'
CALL SSTEV( 'V', N, D1, D2, Z, LDU, WORK, IINFO )
@@ -797,17 +797,17 @@ c LIWEDC = 12
* Do tests 1 and 2.
*
DO 140 I = 1, N
D3( I ) = REAL( A( I, I ) )
D3( I ) = A( I, I )
140 CONTINUE
DO 150 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
150 CONTINUE
CALL SSTT21( N, 0, D3, D4, D1, D2, Z, LDU, WORK,
$ RESULT( 1 ) )
*
NTEST = 3
DO 160 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
160 CONTINUE
SRNAMT = 'SSTEV'
CALL SSTEV( 'N', N, D3, D4, Z, LDU, WORK, IINFO )
@@ -839,10 +839,10 @@ c LIWEDC = 12
NTEST = 4
DO 190 I = 1, N
EVEIGS( I ) = D3( I )
D1( I ) = REAL( A( I, I ) )
D1( I ) = A( I, I )
190 CONTINUE
DO 200 I = 1, N - 1
D2( I ) = REAL( A( I+1, I ) )
D2( I ) = A( I+1, I )
200 CONTINUE
SRNAMT = 'SSTEVX'
CALL SSTEVX( 'V', 'A', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -870,17 +870,17 @@ c LIWEDC = 12
* Do tests 4 and 5.
*
DO 210 I = 1, N
D3( I ) = REAL( A( I, I ) )
D3( I ) = A( I, I )
210 CONTINUE
DO 220 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
220 CONTINUE
CALL SSTT21( N, 0, D3, D4, WA1, D2, Z, LDU, WORK,
$ RESULT( 4 ) )
*
NTEST = 6
DO 230 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
230 CONTINUE
SRNAMT = 'SSTEVX'
CALL SSTEVX( 'N', 'A', N, D3, D4, VL, VU, IL, IU, ABSTOL,
@@ -914,10 +914,10 @@ c LIWEDC = 12
*
NTEST = 7
DO 260 I = 1, N
D1( I ) = REAL( A( I, I ) )
D1( I ) = A( I, I )
260 CONTINUE
DO 270 I = 1, N - 1
D2( I ) = REAL( A( I+1, I ) )
D2( I ) = A( I+1, I )
270 CONTINUE
SRNAMT = 'SSTEVR'
CALL SSTEVR( 'V', 'A', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -944,17 +944,17 @@ c LIWEDC = 12
* Do tests 7 and 8.
*
DO 280 I = 1, N
D3( I ) = REAL( A( I, I ) )
D3( I ) = A( I, I )
280 CONTINUE
DO 290 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
290 CONTINUE
CALL SSTT21( N, 0, D3, D4, WA1, D2, Z, LDU, WORK,
$ RESULT( 7 ) )
*
NTEST = 9
DO 300 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
300 CONTINUE
SRNAMT = 'SSTEVR'
CALL SSTEVR( 'N', 'A', N, D3, D4, VL, VU, IL, IU, ABSTOL,
@@ -989,10 +989,10 @@ c LIWEDC = 12
*
NTEST = 10
DO 330 I = 1, N
D1( I ) = REAL( A( I, I ) )
D1( I ) = A( I, I )
330 CONTINUE
DO 340 I = 1, N - 1
D2( I ) = REAL( A( I+1, I ) )
D2( I ) = A( I+1, I )
340 CONTINUE
SRNAMT = 'SSTEVX'
CALL SSTEVX( 'V', 'I', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -1015,10 +1015,10 @@ c LIWEDC = 12
* Do tests 10 and 11.
*
DO 350 I = 1, N
D3( I ) = REAL( A( I, I ) )
D3( I ) = A( I, I )
350 CONTINUE
DO 360 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
360 CONTINUE
CALL SSTT22( N, M2, 0, D3, D4, WA2, D2, Z, LDU, WORK,
$ MAX( 1, M2 ), RESULT( 10 ) )
@@ -1026,7 +1026,7 @@ c LIWEDC = 12
*
NTEST = 12
DO 370 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
370 CONTINUE
SRNAMT = 'SSTEVX'
CALL SSTEVX( 'N', 'I', N, D3, D4, VL, VU, IL, IU, ABSTOL,
@@ -1076,10 +1076,10 @@ c LIWEDC = 12
END IF
*
DO 390 I = 1, N
D1( I ) = REAL( A( I, I ) )
D1( I ) = A( I, I )
390 CONTINUE
DO 400 I = 1, N - 1
D2( I ) = REAL( A( I+1, I ) )
D2( I ) = A( I+1, I )
400 CONTINUE
SRNAMT = 'SSTEVX'
CALL SSTEVX( 'V', 'V', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -1109,17 +1109,17 @@ c LIWEDC = 12
* Do tests 13 and 14.
*
DO 410 I = 1, N
D3( I ) = REAL( A( I, I ) )
D3( I ) = A( I, I )
410 CONTINUE
DO 420 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
420 CONTINUE
CALL SSTT22( N, M2, 0, D3, D4, WA2, D2, Z, LDU, WORK,
$ MAX( 1, M2 ), RESULT( 13 ) )
*
NTEST = 15
DO 430 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
430 CONTINUE
SRNAMT = 'SSTEVX'
CALL SSTEVX( 'N', 'V', N, D3, D4, VL, VU, IL, IU, ABSTOL,
@@ -1147,10 +1147,10 @@ c LIWEDC = 12
*
NTEST = 16
DO 450 I = 1, N
D1( I ) = REAL( A( I, I ) )
D1( I ) = A( I, I )
450 CONTINUE
DO 460 I = 1, N - 1
D2( I ) = REAL( A( I+1, I ) )
D2( I ) = A( I+1, I )
460 CONTINUE
SRNAMT = 'SSTEVD'
CALL SSTEVD( 'V', N, D1, D2, Z, LDU, WORK, LWEDC, IWORK,
@@ -1172,17 +1172,17 @@ c LIWEDC = 12
* Do tests 16 and 17.
*
DO 470 I = 1, N
D3( I ) = REAL( A( I, I ) )
D3( I ) = A( I, I )
470 CONTINUE
DO 480 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
480 CONTINUE
CALL SSTT21( N, 0, D3, D4, D1, D2, Z, LDU, WORK,
$ RESULT( 16 ) )
*
NTEST = 18
DO 490 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
490 CONTINUE
SRNAMT = 'SSTEVD'
CALL SSTEVD( 'N', N, D3, D4, Z, LDU, WORK, LWEDC, IWORK,
@@ -1215,10 +1215,10 @@ c LIWEDC = 12
*
NTEST = 19
DO 520 I = 1, N
D1( I ) = REAL( A( I, I ) )
D1( I ) = A( I, I )
520 CONTINUE
DO 530 I = 1, N - 1
D2( I ) = REAL( A( I+1, I ) )
D2( I ) = A( I+1, I )
530 CONTINUE
SRNAMT = 'SSTEVR'
CALL SSTEVR( 'V', 'I', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -1241,10 +1241,10 @@ c LIWEDC = 12
* DO tests 19 and 20.
*
DO 540 I = 1, N
D3( I ) = REAL( A( I, I ) )
D3( I ) = A( I, I )
540 CONTINUE
DO 550 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
550 CONTINUE
CALL SSTT22( N, M2, 0, D3, D4, WA2, D2, Z, LDU, WORK,
$ MAX( 1, M2 ), RESULT( 19 ) )
@@ -1252,7 +1252,7 @@ c LIWEDC = 12
*
NTEST = 21
DO 560 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
560 CONTINUE
SRNAMT = 'SSTEVR'
CALL SSTEVR( 'N', 'I', N, D3, D4, VL, VU, IL, IU, ABSTOL,
@@ -1302,10 +1302,10 @@ c LIWEDC = 12
END IF
*
DO 580 I = 1, N
D1( I ) = REAL( A( I, I ) )
D1( I ) = A( I, I )
580 CONTINUE
DO 590 I = 1, N - 1
D2( I ) = REAL( A( I+1, I ) )
D2( I ) = A( I+1, I )
590 CONTINUE
SRNAMT = 'SSTEVR'
CALL SSTEVR( 'V', 'V', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -1335,17 +1335,17 @@ c LIWEDC = 12
* Do tests 22 and 23.
*
DO 600 I = 1, N
D3( I ) = REAL( A( I, I ) )
D3( I ) = A( I, I )
600 CONTINUE
DO 610 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
610 CONTINUE
CALL SSTT22( N, M2, 0, D3, D4, WA2, D2, Z, LDU, WORK,
$ MAX( 1, M2 ), RESULT( 22 ) )
*
NTEST = 24
DO 620 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
620 CONTINUE
SRNAMT = 'SSTEVR'
CALL SSTEVR( 'N', 'V', N, D3, D4, VL, VU, IL, IU, ABSTOL,
TESTING/EIG/sdrvst2stg.f
+40 -40
View File
@@ -776,10 +776,10 @@ c LIWEDC = 12
IF( JTYPE.LE.7 ) THEN
NTEST = 1
DO 120 I = 1, N
D1( I ) = REAL( A( I, I ) )
D1( I ) = A( I, I )
120 CONTINUE
DO 130 I = 1, N - 1
D2( I ) = REAL( A( I+1, I ) )
D2( I ) = A( I+1, I )
130 CONTINUE
SRNAMT = 'SSTEV'
CALL SSTEV( 'V', N, D1, D2, Z, LDU, WORK, IINFO )
@@ -800,17 +800,17 @@ c LIWEDC = 12
* Do tests 1 and 2.
*
DO 140 I = 1, N
D3( I ) = REAL( A( I, I ) )
D3( I ) = A( I, I )
140 CONTINUE
DO 150 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
150 CONTINUE
CALL SSTT21( N, 0, D3, D4, D1, D2, Z, LDU, WORK,
$ RESULT( 1 ) )
*
NTEST = 3
DO 160 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
160 CONTINUE
SRNAMT = 'SSTEV'
CALL SSTEV( 'N', N, D3, D4, Z, LDU, WORK, IINFO )
@@ -842,10 +842,10 @@ c LIWEDC = 12
NTEST = 4
DO 190 I = 1, N
EVEIGS( I ) = D3( I )
D1( I ) = REAL( A( I, I ) )
D1( I ) = A( I, I )
190 CONTINUE
DO 200 I = 1, N - 1
D2( I ) = REAL( A( I+1, I ) )
D2( I ) = A( I+1, I )
200 CONTINUE
SRNAMT = 'SSTEVX'
CALL SSTEVX( 'V', 'A', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -873,17 +873,17 @@ c LIWEDC = 12
* Do tests 4 and 5.
*
DO 210 I = 1, N
D3( I ) = REAL( A( I, I ) )
D3( I ) = A( I, I )
210 CONTINUE
DO 220 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
220 CONTINUE
CALL SSTT21( N, 0, D3, D4, WA1, D2, Z, LDU, WORK,
$ RESULT( 4 ) )
*
NTEST = 6
DO 230 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
230 CONTINUE
SRNAMT = 'SSTEVX'
CALL SSTEVX( 'N', 'A', N, D3, D4, VL, VU, IL, IU, ABSTOL,
@@ -917,10 +917,10 @@ c LIWEDC = 12
*
NTEST = 7
DO 260 I = 1, N
D1( I ) = REAL( A( I, I ) )
D1( I ) = A( I, I )
260 CONTINUE
DO 270 I = 1, N - 1
D2( I ) = REAL( A( I+1, I ) )
D2( I ) = A( I+1, I )
270 CONTINUE
SRNAMT = 'SSTEVR'
CALL SSTEVR( 'V', 'A', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -947,17 +947,17 @@ c LIWEDC = 12
* Do tests 7 and 8.
*
DO 280 I = 1, N
D3( I ) = REAL( A( I, I ) )
D3( I ) = A( I, I )
280 CONTINUE
DO 290 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
290 CONTINUE
CALL SSTT21( N, 0, D3, D4, WA1, D2, Z, LDU, WORK,
$ RESULT( 7 ) )
*
NTEST = 9
DO 300 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
300 CONTINUE
SRNAMT = 'SSTEVR'
CALL SSTEVR( 'N', 'A', N, D3, D4, VL, VU, IL, IU, ABSTOL,
@@ -992,10 +992,10 @@ c LIWEDC = 12
*
NTEST = 10
DO 330 I = 1, N
D1( I ) = REAL( A( I, I ) )
D1( I ) = A( I, I )
330 CONTINUE
DO 340 I = 1, N - 1
D2( I ) = REAL( A( I+1, I ) )
D2( I ) = A( I+1, I )
340 CONTINUE
SRNAMT = 'SSTEVX'
CALL SSTEVX( 'V', 'I', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -1018,10 +1018,10 @@ c LIWEDC = 12
* Do tests 10 and 11.
*
DO 350 I = 1, N
D3( I ) = REAL( A( I, I ) )
D3( I ) = A( I, I )
350 CONTINUE
DO 360 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
360 CONTINUE
CALL SSTT22( N, M2, 0, D3, D4, WA2, D2, Z, LDU, WORK,
$ MAX( 1, M2 ), RESULT( 10 ) )
@@ -1029,7 +1029,7 @@ c LIWEDC = 12
*
NTEST = 12
DO 370 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
370 CONTINUE
SRNAMT = 'SSTEVX'
CALL SSTEVX( 'N', 'I', N, D3, D4, VL, VU, IL, IU, ABSTOL,
@@ -1079,10 +1079,10 @@ c LIWEDC = 12
END IF
*
DO 390 I = 1, N
D1( I ) = REAL( A( I, I ) )
D1( I ) = A( I, I )
390 CONTINUE
DO 400 I = 1, N - 1
D2( I ) = REAL( A( I+1, I ) )
D2( I ) = A( I+1, I )
400 CONTINUE
SRNAMT = 'SSTEVX'
CALL SSTEVX( 'V', 'V', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -1112,17 +1112,17 @@ c LIWEDC = 12
* Do tests 13 and 14.
*
DO 410 I = 1, N
D3( I ) = REAL( A( I, I ) )
D3( I ) = A( I, I )
410 CONTINUE
DO 420 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
420 CONTINUE
CALL SSTT22( N, M2, 0, D3, D4, WA2, D2, Z, LDU, WORK,
$ MAX( 1, M2 ), RESULT( 13 ) )
*
NTEST = 15
DO 430 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
430 CONTINUE
SRNAMT = 'SSTEVX'
CALL SSTEVX( 'N', 'V', N, D3, D4, VL, VU, IL, IU, ABSTOL,
@@ -1150,10 +1150,10 @@ c LIWEDC = 12
*
NTEST = 16
DO 450 I = 1, N
D1( I ) = REAL( A( I, I ) )
D1( I ) = A( I, I )
450 CONTINUE
DO 460 I = 1, N - 1
D2( I ) = REAL( A( I+1, I ) )
D2( I ) = A( I+1, I )
460 CONTINUE
SRNAMT = 'SSTEVD'
CALL SSTEVD( 'V', N, D1, D2, Z, LDU, WORK, LWEDC, IWORK,
@@ -1175,17 +1175,17 @@ c LIWEDC = 12
* Do tests 16 and 17.
*
DO 470 I = 1, N
D3( I ) = REAL( A( I, I ) )
D3( I ) = A( I, I )
470 CONTINUE
DO 480 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
480 CONTINUE
CALL SSTT21( N, 0, D3, D4, D1, D2, Z, LDU, WORK,
$ RESULT( 16 ) )
*
NTEST = 18
DO 490 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
490 CONTINUE
SRNAMT = 'SSTEVD'
CALL SSTEVD( 'N', N, D3, D4, Z, LDU, WORK, LWEDC, IWORK,
@@ -1218,10 +1218,10 @@ c LIWEDC = 12
*
NTEST = 19
DO 520 I = 1, N
D1( I ) = REAL( A( I, I ) )
D1( I ) = A( I, I )
520 CONTINUE
DO 530 I = 1, N - 1
D2( I ) = REAL( A( I+1, I ) )
D2( I ) = A( I+1, I )
530 CONTINUE
SRNAMT = 'SSTEVR'
CALL SSTEVR( 'V', 'I', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -1244,10 +1244,10 @@ c LIWEDC = 12
* DO tests 19 and 20.
*
DO 540 I = 1, N
D3( I ) = REAL( A( I, I ) )
D3( I ) = A( I, I )
540 CONTINUE
DO 550 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
550 CONTINUE
CALL SSTT22( N, M2, 0, D3, D4, WA2, D2, Z, LDU, WORK,
$ MAX( 1, M2 ), RESULT( 19 ) )
@@ -1255,7 +1255,7 @@ c LIWEDC = 12
*
NTEST = 21
DO 560 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
560 CONTINUE
SRNAMT = 'SSTEVR'
CALL SSTEVR( 'N', 'I', N, D3, D4, VL, VU, IL, IU, ABSTOL,
@@ -1305,10 +1305,10 @@ c LIWEDC = 12
END IF
*
DO 580 I = 1, N
D1( I ) = REAL( A( I, I ) )
D1( I ) = A( I, I )
580 CONTINUE
DO 590 I = 1, N - 1
D2( I ) = REAL( A( I+1, I ) )
D2( I ) = A( I+1, I )
590 CONTINUE
SRNAMT = 'SSTEVR'
CALL SSTEVR( 'V', 'V', N, D1, D2, VL, VU, IL, IU, ABSTOL,
@@ -1338,17 +1338,17 @@ c LIWEDC = 12
* Do tests 22 and 23.
*
DO 600 I = 1, N
D3( I ) = REAL( A( I, I ) )
D3( I ) = A( I, I )
600 CONTINUE
DO 610 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
610 CONTINUE
CALL SSTT22( N, M2, 0, D3, D4, WA2, D2, Z, LDU, WORK,
$ MAX( 1, M2 ), RESULT( 22 ) )
*
NTEST = 24
DO 620 I = 1, N - 1
D4( I ) = REAL( A( I+1, I ) )
D4( I ) = A( I+1, I )
620 CONTINUE
SRNAMT = 'SSTEVR'
CALL SSTEVR( 'N', 'V', N, D3, D4, VL, VU, IL, IU, ABSTOL,
TESTING/EIG/slatm4.f
+4 -4
View File
@@ -343,10 +343,10 @@
* Scale by AMAGN
*
DO 230 JD = KBEG, KEND
A( JD, JD ) = AMAGN*REAL( A( JD, JD ) )
A( JD, JD ) = AMAGN*A( JD, JD )
230 CONTINUE
DO 240 JD = ISDB, ISDE
A( JD+1, JD ) = AMAGN*REAL( A( JD+1, JD ) )
A( JD+1, JD ) = AMAGN*A( JD+1, JD )
240 CONTINUE
*
* If ISIGN = 1 or 2, assign random signs to diagonal and
@@ -354,13 +354,13 @@
*
IF( ISIGN.GT.0 ) THEN
DO 250 JD = KBEG, KEND
IF( REAL( A( JD, JD ) ).NE.ZERO ) THEN
IF( A( JD, JD ).NE.ZERO ) THEN
IF( SLARAN( ISEED ).GT.HALF )
$ A( JD, JD ) = -A( JD, JD )
END IF
250 CONTINUE
DO 260 JD = ISDB, ISDE
IF( REAL( A( JD+1, JD ) ).NE.ZERO ) THEN
IF( A( JD+1, JD ).NE.ZERO ) THEN
IF( SLARAN( ISEED ).GT.HALF )
$ A( JD+1, JD ) = -A( JD+1, JD )
END IF