diff --git a/TESTING/EIG/ddrvst.f b/TESTING/EIG/ddrvst.f index 77e975575..61d6973f2 100644 --- a/TESTING/EIG/ddrvst.f +++ b/TESTING/EIG/ddrvst.f @@ -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, diff --git a/TESTING/EIG/ddrvst2stg.f b/TESTING/EIG/ddrvst2stg.f index 9a24e346d..f1c848d07 100644 --- a/TESTING/EIG/ddrvst2stg.f +++ b/TESTING/EIG/ddrvst2stg.f @@ -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, diff --git a/TESTING/EIG/dlatm4.f b/TESTING/EIG/dlatm4.f index 4a028e5a8..b271dce05 100644 --- a/TESTING/EIG/dlatm4.f +++ b/TESTING/EIG/dlatm4.f @@ -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 diff --git a/TESTING/EIG/sdrvst.f b/TESTING/EIG/sdrvst.f index 22fdf264e..34bd053c3 100644 --- a/TESTING/EIG/sdrvst.f +++ b/TESTING/EIG/sdrvst.f @@ -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, diff --git a/TESTING/EIG/sdrvst2stg.f b/TESTING/EIG/sdrvst2stg.f index 6fe62b071..f20d9fdf3 100644 --- a/TESTING/EIG/sdrvst2stg.f +++ b/TESTING/EIG/sdrvst2stg.f @@ -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, diff --git a/TESTING/EIG/slatm4.f b/TESTING/EIG/slatm4.f index d286f0f8d..c452f5085 100644 --- a/TESTING/EIG/slatm4.f +++ b/TESTING/EIG/slatm4.f @@ -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