Lines Matching +full:x +full:- +full:alpha
6 * http://www.netlib.org/lapack/explore-html/
22 *> of the file are read using list-directed input, the last 16 records
29 *> -1 UNIT NUMBER OF SNAPSHOT FILE (NOT USED IF .LT. 0)
39 *> 1 2 -1 -2 VALUES OF INCX AND INCY
40 *> 3 NUMBER OF VALUES OF ALPHA
41 *> 0.0 1.0 0.7 VALUES OF ALPHA
78 *> Street, Suite 100, Downers Grove, Illinois 60515-1263, USA.
81 *> -- Written on 10-August-1987.
85 *> 10-9-00: Change STATUS='NEW' to 'UNKNOWN' so that the testers
105 * -- Reference BLAS test routine (version 3.4.1) --
106 * -- Reference BLAS is a software package provided by Univ. of Tennessee, --
107 * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
136 $ G( NMAX ), X( NMAX ), XS( NMAX*INCMAX ), variable
232 * Values of ALPHA
235 WRITE( NOUT, FMT = 9997 )'ALPHA', NALMAX
292 A( I, J ) = MAX( I - J + 1, 0 )
294 X( J ) = J
298 YY( J ) = J*( ( J + 1 )*J )/2 - ( ( J + 1 )*J*( J - 1 ) )/3
303 CALL SMVCH( TRANS, N, N, ONE, A, NMAX, X, 1, ZERO, Y, 1, YT, G,
311 CALL SMVCH( TRANS, N, N, ONE, A, NMAX, X, -1, ZERO, Y, -1, YT, G,
344 $ X, XX, XS, Y, YY, YS, YT, G )
350 $ X, XX, XS, Y, YY, YS, YT, G )
361 $ NMAX, INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS,
367 $ NMAX, INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS,
373 $ NMAX, INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS,
410 9989 FORMAT( ' FOR ALPHA ', 7F6.1 )
416 9985 FORMAT( ' ERROR IN SMVCH - IN-LINE DOT PRODUCTS ARE BEING EVALU',
422 9983 FORMAT( 1X, A6, ' WAS NOT TESTED' )
424 9981 FORMAT( /' ******* FATAL ERROR - TESTS ABANDONED *******' )
425 9980 FORMAT( ' ERROR-EXITS WILL NOT BE TESTED' )
432 $ BET, NINC, INC, NMAX, INCMAX, A, AA, AS, X, XX, argument
439 * -- Written on 10-August-1987.
455 $ X( NMAX ), XS( NMAX*INCMAX ), local
461 REAL ALPHA, ALS, BETA, BLS, ERR, ERRMAX, TRANSL local
505 $ M = MAX( N - ND, 0 )
517 KL = MAX( KU - 1, 0 )
519 KU = N - 1
520 KL = M - 1
558 * Generate the vector X.
561 CALL SMAKE( 'GE', ' ', ' ', 1, NL, X, 1, XX,
562 $ ABS( INCX ), 0, NL - 1, RESET, TRANSL )
564 X( NL/2 ) = ZERO
565 XX( 1 + ABS( INCX )*( NL/2 - 1 ) ) = ZERO
573 ALPHA = ALF( IA )
582 $ YY, ABS( INCY ), 0, ML - 1,
595 ALS = ALPHA
615 $ TRANS, M, N, ALPHA, LDA, INCX, BETA,
619 CALL SGEMV( TRANS, M, N, ALPHA, AA,
625 $ TRANS, M, N, KL, KU, ALPHA, LDA,
629 CALL SGBMV( TRANS, M, N, KL, KU, ALPHA,
634 * Check if error-exit was taken incorrectly.
648 ISAME( 4 ) = ALS.EQ.ALPHA
665 ISAME( 6 ) = ALS.EQ.ALPHA
699 CALL SMVCH( TRANS, M, N, ALPHA, A,
700 $ NMAX, X, INCX, BETA, Y,
742 WRITE( NOUT, FMT = 9994 )NC, SNAME, TRANS, M, N, ALPHA, LDA,
746 $ ALPHA, LDA, INCX, BETA, INCY
754 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
758 $ ' - SUSPECT *******' )
760 9995 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', 4( I3, ',' ), F4.1,
761 $ ', A,', I3, ', X,', I2, ',', F4.1, ', Y,', I2, ') .' )
762 9994 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', 2( I3, ',' ), F4.1,
763 $ ', A,', I3, ', X,', I2, ',', F4.1, ', Y,', I2,
765 9993 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
773 $ BET, NINC, INC, NMAX, INCMAX, A, AA, AS, X, XX, argument
780 * -- Written on 10-August-1987.
796 $ X( NMAX ), XS( NMAX*INCMAX ), local
802 REAL ALPHA, ALS, BETA, BLS, ERR, ERRMAX, TRANSL local
854 K = N - 1
887 * Generate the vector X.
890 CALL SMAKE( 'GE', ' ', ' ', 1, N, X, 1, XX,
891 $ ABS( INCX ), 0, N - 1, RESET, TRANSL )
893 X( N/2 ) = ZERO
894 XX( 1 + ABS( INCX )*( N/2 - 1 ) ) = ZERO
902 ALPHA = ALF( IA )
911 $ ABS( INCY ), 0, N - 1, RESET,
922 ALS = ALPHA
942 $ UPLO, N, ALPHA, LDA, INCX, BETA, INCY
945 CALL SSYMV( UPLO, N, ALPHA, AA, LDA, XX,
950 $ UPLO, N, K, ALPHA, LDA, INCX, BETA,
954 CALL SSBMV( UPLO, N, K, ALPHA, AA, LDA,
959 $ UPLO, N, ALPHA, INCX, BETA, INCY
962 CALL SSPMV( UPLO, N, ALPHA, AA, XX, INCX,
966 * Check if error-exit was taken incorrectly.
979 ISAME( 3 ) = ALS.EQ.ALPHA
994 ISAME( 4 ) = ALS.EQ.ALPHA
1008 ISAME( 3 ) = ALS.EQ.ALPHA
1040 CALL SMVCH( 'N', N, N, ALPHA, A, NMAX, X,
1080 WRITE( NOUT, FMT = 9993 )NC, SNAME, UPLO, N, ALPHA, LDA, INCX,
1083 WRITE( NOUT, FMT = 9994 )NC, SNAME, UPLO, N, K, ALPHA, LDA,
1086 WRITE( NOUT, FMT = 9995 )NC, SNAME, UPLO, N, ALPHA, INCX,
1095 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
1099 $ ' - SUSPECT *******' )
1101 9995 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',', F4.1, ', AP',
1102 $ ', X,', I2, ',', F4.1, ', Y,', I2, ') .' )
1103 9994 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', 2( I3, ',' ), F4.1,
1104 $ ', A,', I3, ', X,', I2, ',', F4.1, ', Y,', I2,
1106 9993 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',', F4.1, ', A,',
1107 $ I3, ', X,', I2, ',', F4.1, ', Y,', I2, ') .' )
1108 9992 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
1116 $ INCMAX, A, AA, AS, X, XX, XS, XT, G, Z ) argument
1122 * -- Written on 10-August-1987.
1136 $ AS( NMAX*NMAX ), G( NMAX ), X( NMAX ), local
1198 K = N - 1
1237 * Generate the vector X.
1240 CALL SMAKE( 'GE', ' ', ' ', 1, N, X, 1, XX,
1241 $ ABS( INCX ), 0, N - 1, RESET,
1244 X( N/2 ) = ZERO
1245 XX( 1 + ABS( INCX )*( N/2 - 1 ) ) = ZERO
1322 * Check if error-exit was taken incorrectly.
1387 CALL SMVCH( TRANS, N, N, ONE, A, NMAX, X,
1396 Z( I ) = XX( 1 + ( I - 1 )*
1398 XX( 1 + ( I - 1 )*ABS( INCX ) )
1399 $ = X( I )
1402 $ INCX, ZERO, X, INCX, XT, G,
1453 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
1457 $ ' - SUSPECT *******' )
1459 9995 FORMAT( 1X, I6, ': ', A6, '(', 3( '''', A1, ''',' ), I3, ', AP, ',
1460 $ 'X,', I2, ') .' )
1461 9994 FORMAT( 1X, I6, ': ', A6, '(', 3( '''', A1, ''',' ), 2( I3, ',' ),
1462 $ ' A,', I3, ', X,', I2, ') .' )
1463 9993 FORMAT( 1X, I6, ': ', A6, '(', 3( '''', A1, ''',' ), I3, ', A,',
1464 $ I3, ', X,', I2, ') .' )
1465 9992 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
1473 $ INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS, YT, G, argument
1480 * -- Written on 10-August-1987.
1494 $ AS( NMAX*NMAX ), G( NMAX ), X( NMAX ), local
1500 REAL ALPHA, ALS, ERR, ERRMAX, TRANSL local
1534 $ M = MAX( N - ND, 0 )
1552 * Generate the vector X.
1555 CALL SMAKE( 'GE', ' ', ' ', 1, M, X, 1, XX, ABS( INCX ),
1556 $ 0, M - 1, RESET, TRANSL )
1558 X( M/2 ) = ZERO
1559 XX( 1 + ABS( INCX )*( M/2 - 1 ) ) = ZERO
1570 $ ABS( INCY ), 0, N - 1, RESET, TRANSL )
1573 YY( 1 + ABS( INCY )*( N/2 - 1 ) ) = ZERO
1577 ALPHA = ALF( IA )
1583 $ AA, LDA, M - 1, N - 1, RESET, TRANSL )
1591 ALS = ALPHA
1609 $ ALPHA, INCX, INCY, LDA
1612 CALL SGER( M, N, ALPHA, XX, INCX, YY, INCY, AA,
1615 * Check if error-exit was taken incorrectly.
1627 ISAME( 3 ) = ALS.EQ.ALPHA
1659 Z( I ) = X( I )
1663 Z( I ) = X( M - I + 1 )
1670 W( 1 ) = Y( N - J + 1 )
1672 CALL SMVCH( 'N', M, 1, ALPHA, Z, NMAX, W, 1,
1674 $ AA( 1 + ( J - 1 )*LDA ), EPS,
1710 WRITE( NOUT, FMT = 9994 )NC, SNAME, M, N, ALPHA, INCX, INCY, LDA
1717 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
1721 $ ' - SUSPECT *******' )
1724 9994 FORMAT( 1X, I6, ': ', A6, '(', 2( I3, ',' ), F4.1, ', X,', I2,
1726 9993 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
1734 $ INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS, YT, G, argument
1741 * -- Written on 10-August-1987.
1755 $ AS( NMAX*NMAX ), G( NMAX ), X( NMAX ), local
1761 REAL ALPHA, ALS, ERR, ERRMAX, TRANSL local
1821 * Generate the vector X.
1824 CALL SMAKE( 'GE', ' ', ' ', 1, N, X, 1, XX, ABS( INCX ),
1825 $ 0, N - 1, RESET, TRANSL )
1827 X( N/2 ) = ZERO
1828 XX( 1 + ABS( INCX )*( N/2 - 1 ) ) = ZERO
1832 ALPHA = ALF( IA )
1833 NULL = N.LE.0.OR.ALPHA.EQ.ZERO
1839 $ AA, LDA, N - 1, N - 1, RESET, TRANSL )
1847 ALS = ALPHA
1862 $ ALPHA, INCX, LDA
1865 CALL SSYR( UPLO, N, ALPHA, XX, INCX, AA, LDA )
1869 $ ALPHA, INCX
1872 CALL SSPR( UPLO, N, ALPHA, XX, INCX, AA )
1875 * Check if error-exit was taken incorrectly.
1887 ISAME( 3 ) = ALS.EQ.ALPHA
1919 Z( I ) = X( I )
1923 Z( I ) = X( N - I + 1 )
1934 LJ = N - J + 1
1936 CALL SMVCH( 'N', LJ, 1, ALPHA, Z( JJ ), LJ, W,
1983 WRITE( NOUT, FMT = 9993 )NC, SNAME, UPLO, N, ALPHA, INCX, LDA
1985 WRITE( NOUT, FMT = 9994 )NC, SNAME, UPLO, N, ALPHA, INCX
1993 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
1997 $ ' - SUSPECT *******' )
2000 9994 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',', F4.1, ', X,',
2002 9993 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',', F4.1, ', X,',
2004 9992 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
2012 $ INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS, YT, G, argument
2019 * -- Written on 10-August-1987.
2033 $ AS( NMAX*NMAX ), G( NMAX ), X( NMAX ), local
2039 REAL ALPHA, ALS, ERR, ERRMAX, TRANSL local
2100 * Generate the vector X.
2103 CALL SMAKE( 'GE', ' ', ' ', 1, N, X, 1, XX, ABS( INCX ),
2104 $ 0, N - 1, RESET, TRANSL )
2106 X( N/2 ) = ZERO
2107 XX( 1 + ABS( INCX )*( N/2 - 1 ) ) = ZERO
2118 $ ABS( INCY ), 0, N - 1, RESET, TRANSL )
2121 YY( 1 + ABS( INCY )*( N/2 - 1 ) ) = ZERO
2125 ALPHA = ALF( IA )
2126 NULL = N.LE.0.OR.ALPHA.EQ.ZERO
2132 $ NMAX, AA, LDA, N - 1, N - 1, RESET,
2141 ALS = ALPHA
2160 $ ALPHA, INCX, INCY, LDA
2163 CALL SSYR2( UPLO, N, ALPHA, XX, INCX, YY, INCY,
2168 $ ALPHA, INCX, INCY
2171 CALL SSPR2( UPLO, N, ALPHA, XX, INCX, YY, INCY,
2175 * Check if error-exit was taken incorrectly.
2187 ISAME( 3 ) = ALS.EQ.ALPHA
2221 Z( I, 1 ) = X( I )
2225 Z( I, 1 ) = X( N - I + 1 )
2234 Z( I, 2 ) = Y( N - I + 1 )
2246 LJ = N - J + 1
2248 CALL SMVCH( 'N', LJ, 2, ALPHA, Z( JJ, 1 ),
2297 WRITE( NOUT, FMT = 9993 )NC, SNAME, UPLO, N, ALPHA, INCX,
2300 WRITE( NOUT, FMT = 9994 )NC, SNAME, UPLO, N, ALPHA, INCX, INCY
2308 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
2312 $ ' - SUSPECT *******' )
2315 9994 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',', F4.1, ', X,',
2317 9993 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',', F4.1, ', X,',
2319 9992 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
2328 * Requires a special version of the error-handling routine XERBLA.
2329 * ALPHA, BETA, A, X and Y should not need to be defined.
2333 * -- Written on 10-August-1987.
2344 REAL ALPHA, BETA local
2346 REAL A( 1, 1 ), X( 1 ), Y( 1 ) local
2358 * it is called, and is then tested and re-set by CHKXER.
2364 CALL SGEMV( '/', 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2367 CALL SGEMV( 'N', -1, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2370 CALL SGEMV( 'N', 0, -1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2373 CALL SGEMV( 'N', 2, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2376 CALL SGEMV( 'N', 0, 0, ALPHA, A, 1, X, 0, BETA, Y, 1 )
2379 CALL SGEMV( 'N', 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 0 )
2383 CALL SGBMV( '/', 0, 0, 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2386 CALL SGBMV( 'N', -1, 0, 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2389 CALL SGBMV( 'N', 0, -1, 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2392 CALL SGBMV( 'N', 0, 0, -1, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2395 CALL SGBMV( 'N', 2, 0, 0, -1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2398 CALL SGBMV( 'N', 0, 0, 1, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2401 CALL SGBMV( 'N', 0, 0, 0, 0, ALPHA, A, 1, X, 0, BETA, Y, 1 )
2404 CALL SGBMV( 'N', 0, 0, 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 0 )
2408 CALL SSYMV( '/', 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2411 CALL SSYMV( 'U', -1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2414 CALL SSYMV( 'U', 2, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2417 CALL SSYMV( 'U', 0, ALPHA, A, 1, X, 0, BETA, Y, 1 )
2420 CALL SSYMV( 'U', 0, ALPHA, A, 1, X, 1, BETA, Y, 0 )
2424 CALL SSBMV( '/', 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2427 CALL SSBMV( 'U', -1, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2430 CALL SSBMV( 'U', 0, -1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2433 CALL SSBMV( 'U', 0, 1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2436 CALL SSBMV( 'U', 0, 0, ALPHA, A, 1, X, 0, BETA, Y, 1 )
2439 CALL SSBMV( 'U', 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 0 )
2443 CALL SSPMV( '/', 0, ALPHA, A, X, 1, BETA, Y, 1 )
2446 CALL SSPMV( 'U', -1, ALPHA, A, X, 1, BETA, Y, 1 )
2449 CALL SSPMV( 'U', 0, ALPHA, A, X, 0, BETA, Y, 1 )
2452 CALL SSPMV( 'U', 0, ALPHA, A, X, 1, BETA, Y, 0 )
2456 CALL STRMV( '/', 'N', 'N', 0, A, 1, X, 1 )
2459 CALL STRMV( 'U', '/', 'N', 0, A, 1, X, 1 )
2462 CALL STRMV( 'U', 'N', '/', 0, A, 1, X, 1 )
2465 CALL STRMV( 'U', 'N', 'N', -1, A, 1, X, 1 )
2468 CALL STRMV( 'U', 'N', 'N', 2, A, 1, X, 1 )
2471 CALL STRMV( 'U', 'N', 'N', 0, A, 1, X, 0 )
2475 CALL STBMV( '/', 'N', 'N', 0, 0, A, 1, X, 1 )
2478 CALL STBMV( 'U', '/', 'N', 0, 0, A, 1, X, 1 )
2481 CALL STBMV( 'U', 'N', '/', 0, 0, A, 1, X, 1 )
2484 CALL STBMV( 'U', 'N', 'N', -1, 0, A, 1, X, 1 )
2487 CALL STBMV( 'U', 'N', 'N', 0, -1, A, 1, X, 1 )
2490 CALL STBMV( 'U', 'N', 'N', 0, 1, A, 1, X, 1 )
2493 CALL STBMV( 'U', 'N', 'N', 0, 0, A, 1, X, 0 )
2497 CALL STPMV( '/', 'N', 'N', 0, A, X, 1 )
2500 CALL STPMV( 'U', '/', 'N', 0, A, X, 1 )
2503 CALL STPMV( 'U', 'N', '/', 0, A, X, 1 )
2506 CALL STPMV( 'U', 'N', 'N', -1, A, X, 1 )
2509 CALL STPMV( 'U', 'N', 'N', 0, A, X, 0 )
2513 CALL STRSV( '/', 'N', 'N', 0, A, 1, X, 1 )
2516 CALL STRSV( 'U', '/', 'N', 0, A, 1, X, 1 )
2519 CALL STRSV( 'U', 'N', '/', 0, A, 1, X, 1 )
2522 CALL STRSV( 'U', 'N', 'N', -1, A, 1, X, 1 )
2525 CALL STRSV( 'U', 'N', 'N', 2, A, 1, X, 1 )
2528 CALL STRSV( 'U', 'N', 'N', 0, A, 1, X, 0 )
2532 CALL STBSV( '/', 'N', 'N', 0, 0, A, 1, X, 1 )
2535 CALL STBSV( 'U', '/', 'N', 0, 0, A, 1, X, 1 )
2538 CALL STBSV( 'U', 'N', '/', 0, 0, A, 1, X, 1 )
2541 CALL STBSV( 'U', 'N', 'N', -1, 0, A, 1, X, 1 )
2544 CALL STBSV( 'U', 'N', 'N', 0, -1, A, 1, X, 1 )
2547 CALL STBSV( 'U', 'N', 'N', 0, 1, A, 1, X, 1 )
2550 CALL STBSV( 'U', 'N', 'N', 0, 0, A, 1, X, 0 )
2554 CALL STPSV( '/', 'N', 'N', 0, A, X, 1 )
2557 CALL STPSV( 'U', '/', 'N', 0, A, X, 1 )
2560 CALL STPSV( 'U', 'N', '/', 0, A, X, 1 )
2563 CALL STPSV( 'U', 'N', 'N', -1, A, X, 1 )
2566 CALL STPSV( 'U', 'N', 'N', 0, A, X, 0 )
2570 CALL SGER( -1, 0, ALPHA, X, 1, Y, 1, A, 1 )
2573 CALL SGER( 0, -1, ALPHA, X, 1, Y, 1, A, 1 )
2576 CALL SGER( 0, 0, ALPHA, X, 0, Y, 1, A, 1 )
2579 CALL SGER( 0, 0, ALPHA, X, 1, Y, 0, A, 1 )
2582 CALL SGER( 2, 0, ALPHA, X, 1, Y, 1, A, 1 )
2586 CALL SSYR( '/', 0, ALPHA, X, 1, A, 1 )
2589 CALL SSYR( 'U', -1, ALPHA, X, 1, A, 1 )
2592 CALL SSYR( 'U', 0, ALPHA, X, 0, A, 1 )
2595 CALL SSYR( 'U', 2, ALPHA, X, 1, A, 1 )
2599 CALL SSPR( '/', 0, ALPHA, X, 1, A )
2602 CALL SSPR( 'U', -1, ALPHA, X, 1, A )
2605 CALL SSPR( 'U', 0, ALPHA, X, 0, A )
2609 CALL SSYR2( '/', 0, ALPHA, X, 1, Y, 1, A, 1 )
2612 CALL SSYR2( 'U', -1, ALPHA, X, 1, Y, 1, A, 1 )
2615 CALL SSYR2( 'U', 0, ALPHA, X, 0, Y, 1, A, 1 )
2618 CALL SSYR2( 'U', 0, ALPHA, X, 1, Y, 0, A, 1 )
2621 CALL SSYR2( 'U', 2, ALPHA, X, 1, Y, 1, A, 1 )
2625 CALL SSPR2( '/', 0, ALPHA, X, 1, Y, 1, A )
2628 CALL SSPR2( 'U', -1, ALPHA, X, 1, Y, 1, A )
2631 CALL SSPR2( 'U', 0, ALPHA, X, 0, Y, 1, A )
2634 CALL SSPR2( 'U', 0, ALPHA, X, 1, Y, 0, A )
2644 9999 FORMAT( ' ', A6, ' PASSED THE TESTS OF ERROR-EXITS' )
2645 9998 FORMAT( ' ******* ', A6, ' FAILED THE TESTS OF ERROR-EXITS *****',
2663 * -- Written on 10-August-1987.
2671 PARAMETER ( ROGUE = -1.0E10 )
2702 IF( ( I.LE.J.AND.J - I.LE.KU ).OR.
2703 $ ( I.GE.J.AND.I - J.LE.KL ) )THEN
2728 AA( I + ( J - 1 )*LDA ) = A( I, J )
2731 AA( I + ( J - 1 )*LDA ) = ROGUE
2736 DO 60 I1 = 1, KU + 1 - J
2737 AA( I1 + ( J - 1 )*LDA ) = ROGUE
2739 DO 70 I2 = I1, MIN( KL + KU + 1, KU + 1 + M - J )
2740 AA( I2 + ( J - 1 )*LDA ) = A( I2 + J - KU - 1, J )
2743 AA( I3 + ( J - 1 )*LDA ) = ROGUE
2751 IEND = J - 1
2763 DO 100 I = 1, IBEG - 1
2764 AA( I + ( J - 1 )*LDA ) = ROGUE
2767 AA( I + ( J - 1 )*LDA ) = A( I, J )
2770 AA( I + ( J - 1 )*LDA ) = ROGUE
2777 IBEG = MAX( 1, KL + 2 - J )
2790 IEND = MIN( KL + 1, 1 + M - J )
2792 DO 140 I = 1, IBEG - 1
2793 AA( I + ( J - 1 )*LDA ) = ROGUE
2796 AA( I + ( J - 1 )*LDA ) = A( I + J - KK, J )
2799 AA( I + ( J - 1 )*LDA ) = ROGUE
2827 SUBROUTINE SMVCH( TRANS, M, N, ALPHA, A, NMAX, X, INCX, BETA, Y, argument
2834 * -- Written on 10-August-1987.
2842 REAL ALPHA, BETA, EPS, ERR local
2847 REAL A( NMAX, * ), G( * ), X( * ), Y( * ), YT( * ), local
2866 INCXL = -1
2873 INCYL = -1
2879 * Compute expected result in YT using data in A, X and Y.
2889 YT( IY ) = YT( IY ) + A( J, I )*X( JX )
2890 G( IY ) = G( IY ) + ABS( A( J, I )*X( JX ) )
2895 YT( IY ) = YT( IY ) + A( I, J )*X( JX )
2896 G( IY ) = G( IY ) + ABS( A( I, J )*X( JX ) )
2900 YT( IY ) = ALPHA*YT( IY ) + BETA*Y( IY )
2901 G( IY ) = ABS( ALPHA )*G( IY ) + ABS( BETA*Y( IY ) )
2909 ERRI = ABS( YT( I ) - YY( 1 + ( I - 1 )*ABS( INCY ) ) )/EPS
2926 $ YY( 1 + ( I - 1 )*ABS( INCY ) )
2929 $ YY( 1 + ( I - 1 )*ABS( INCY ) ), YT(I)
2936 9999 FORMAT( ' ******* FATAL ERROR - COMPUTED RESULT IS LESS THAN HAL',
2939 9998 FORMAT( 1X, I7, 2G18.6 )
2950 * -- Written on 10-August-1987.
2982 * -- Written on 10-August-1987.
3013 DO 30 I = 1, IBEG - 1
3035 * Generates random numbers uniformly distributed between -0.5 and 0.5.
3039 * -- Written on 10-August-1987.
3068 I = I - 1000*( I/1000 )
3073 SBEG = REAL( I - 500 )/1001.0
3079 REAL FUNCTION SDIFF( X, Y )
3083 * -- Written on 10-August-1987.
3087 REAL X, Y local
3089 SDIFF = X - Y
3101 * -- Written on 10-August-1987.
3136 * -- Written on 10-August-1987.