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 DMVCH( TRANS, N, N, ONE, A, NMAX, X, 1, ZERO, Y, 1, YT, G,
311 CALL DMVCH( 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 DMVCH - 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 DOUBLE PRECISION 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 DMAKE( '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 DGEMV( TRANS, M, N, ALPHA, AA,
625 $ TRANS, M, N, KL, KU, ALPHA, LDA,
629 CALL DGBMV( 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 DMVCH( 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 DOUBLE PRECISION ALPHA, ALS, BETA, BLS, ERR, ERRMAX, TRANSL local
854 K = N - 1
887 * Generate the vector X.
890 CALL DMAKE( '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 DSYMV( UPLO, N, ALPHA, AA, LDA, XX,
950 $ UPLO, N, K, ALPHA, LDA, INCX, BETA,
954 CALL DSBMV( UPLO, N, K, ALPHA, AA, LDA,
959 $ UPLO, N, ALPHA, INCX, BETA, INCY
962 CALL DSPMV( 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 DMVCH( '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 DMAKE( '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 DMVCH( 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 DOUBLE PRECISION ALPHA, ALS, ERR, ERRMAX, TRANSL local
1534 $ M = MAX( N - ND, 0 )
1552 * Generate the vector X.
1555 CALL DMAKE( '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 DGER( 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 DMVCH( '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 DOUBLE PRECISION ALPHA, ALS, ERR, ERRMAX, TRANSL local
1821 * Generate the vector X.
1824 CALL DMAKE( '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 DSYR( UPLO, N, ALPHA, XX, INCX, AA, LDA )
1869 $ ALPHA, INCX
1872 CALL DSPR( 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 DMVCH( '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 DOUBLE PRECISION ALPHA, ALS, ERR, ERRMAX, TRANSL local
2100 * Generate the vector X.
2103 CALL DMAKE( '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 DSYR2( UPLO, N, ALPHA, XX, INCX, YY, INCY,
2168 $ ALPHA, INCX, INCY
2171 CALL DSPR2( 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 DMVCH( '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 DOUBLE PRECISION ALPHA, BETA local
2346 DOUBLE PRECISION A( 1, 1 ), X( 1 ), Y( 1 ) local
2358 * it is called, and is then tested and re-set by CHKXER.
2364 CALL DGEMV( '/', 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2367 CALL DGEMV( 'N', -1, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2370 CALL DGEMV( 'N', 0, -1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2373 CALL DGEMV( 'N', 2, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2376 CALL DGEMV( 'N', 0, 0, ALPHA, A, 1, X, 0, BETA, Y, 1 )
2379 CALL DGEMV( 'N', 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 0 )
2383 CALL DGBMV( '/', 0, 0, 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2386 CALL DGBMV( 'N', -1, 0, 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2389 CALL DGBMV( 'N', 0, -1, 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2392 CALL DGBMV( 'N', 0, 0, -1, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2395 CALL DGBMV( 'N', 2, 0, 0, -1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2398 CALL DGBMV( 'N', 0, 0, 1, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2401 CALL DGBMV( 'N', 0, 0, 0, 0, ALPHA, A, 1, X, 0, BETA, Y, 1 )
2404 CALL DGBMV( 'N', 0, 0, 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 0 )
2408 CALL DSYMV( '/', 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2411 CALL DSYMV( 'U', -1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2414 CALL DSYMV( 'U', 2, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2417 CALL DSYMV( 'U', 0, ALPHA, A, 1, X, 0, BETA, Y, 1 )
2420 CALL DSYMV( 'U', 0, ALPHA, A, 1, X, 1, BETA, Y, 0 )
2424 CALL DSBMV( '/', 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2427 CALL DSBMV( 'U', -1, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2430 CALL DSBMV( 'U', 0, -1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2433 CALL DSBMV( 'U', 0, 1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2436 CALL DSBMV( 'U', 0, 0, ALPHA, A, 1, X, 0, BETA, Y, 1 )
2439 CALL DSBMV( 'U', 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 0 )
2443 CALL DSPMV( '/', 0, ALPHA, A, X, 1, BETA, Y, 1 )
2446 CALL DSPMV( 'U', -1, ALPHA, A, X, 1, BETA, Y, 1 )
2449 CALL DSPMV( 'U', 0, ALPHA, A, X, 0, BETA, Y, 1 )
2452 CALL DSPMV( 'U', 0, ALPHA, A, X, 1, BETA, Y, 0 )
2456 CALL DTRMV( '/', 'N', 'N', 0, A, 1, X, 1 )
2459 CALL DTRMV( 'U', '/', 'N', 0, A, 1, X, 1 )
2462 CALL DTRMV( 'U', 'N', '/', 0, A, 1, X, 1 )
2465 CALL DTRMV( 'U', 'N', 'N', -1, A, 1, X, 1 )
2468 CALL DTRMV( 'U', 'N', 'N', 2, A, 1, X, 1 )
2471 CALL DTRMV( 'U', 'N', 'N', 0, A, 1, X, 0 )
2475 CALL DTBMV( '/', 'N', 'N', 0, 0, A, 1, X, 1 )
2478 CALL DTBMV( 'U', '/', 'N', 0, 0, A, 1, X, 1 )
2481 CALL DTBMV( 'U', 'N', '/', 0, 0, A, 1, X, 1 )
2484 CALL DTBMV( 'U', 'N', 'N', -1, 0, A, 1, X, 1 )
2487 CALL DTBMV( 'U', 'N', 'N', 0, -1, A, 1, X, 1 )
2490 CALL DTBMV( 'U', 'N', 'N', 0, 1, A, 1, X, 1 )
2493 CALL DTBMV( 'U', 'N', 'N', 0, 0, A, 1, X, 0 )
2497 CALL DTPMV( '/', 'N', 'N', 0, A, X, 1 )
2500 CALL DTPMV( 'U', '/', 'N', 0, A, X, 1 )
2503 CALL DTPMV( 'U', 'N', '/', 0, A, X, 1 )
2506 CALL DTPMV( 'U', 'N', 'N', -1, A, X, 1 )
2509 CALL DTPMV( 'U', 'N', 'N', 0, A, X, 0 )
2513 CALL DTRSV( '/', 'N', 'N', 0, A, 1, X, 1 )
2516 CALL DTRSV( 'U', '/', 'N', 0, A, 1, X, 1 )
2519 CALL DTRSV( 'U', 'N', '/', 0, A, 1, X, 1 )
2522 CALL DTRSV( 'U', 'N', 'N', -1, A, 1, X, 1 )
2525 CALL DTRSV( 'U', 'N', 'N', 2, A, 1, X, 1 )
2528 CALL DTRSV( 'U', 'N', 'N', 0, A, 1, X, 0 )
2532 CALL DTBSV( '/', 'N', 'N', 0, 0, A, 1, X, 1 )
2535 CALL DTBSV( 'U', '/', 'N', 0, 0, A, 1, X, 1 )
2538 CALL DTBSV( 'U', 'N', '/', 0, 0, A, 1, X, 1 )
2541 CALL DTBSV( 'U', 'N', 'N', -1, 0, A, 1, X, 1 )
2544 CALL DTBSV( 'U', 'N', 'N', 0, -1, A, 1, X, 1 )
2547 CALL DTBSV( 'U', 'N', 'N', 0, 1, A, 1, X, 1 )
2550 CALL DTBSV( 'U', 'N', 'N', 0, 0, A, 1, X, 0 )
2554 CALL DTPSV( '/', 'N', 'N', 0, A, X, 1 )
2557 CALL DTPSV( 'U', '/', 'N', 0, A, X, 1 )
2560 CALL DTPSV( 'U', 'N', '/', 0, A, X, 1 )
2563 CALL DTPSV( 'U', 'N', 'N', -1, A, X, 1 )
2566 CALL DTPSV( 'U', 'N', 'N', 0, A, X, 0 )
2570 CALL DGER( -1, 0, ALPHA, X, 1, Y, 1, A, 1 )
2573 CALL DGER( 0, -1, ALPHA, X, 1, Y, 1, A, 1 )
2576 CALL DGER( 0, 0, ALPHA, X, 0, Y, 1, A, 1 )
2579 CALL DGER( 0, 0, ALPHA, X, 1, Y, 0, A, 1 )
2582 CALL DGER( 2, 0, ALPHA, X, 1, Y, 1, A, 1 )
2586 CALL DSYR( '/', 0, ALPHA, X, 1, A, 1 )
2589 CALL DSYR( 'U', -1, ALPHA, X, 1, A, 1 )
2592 CALL DSYR( 'U', 0, ALPHA, X, 0, A, 1 )
2595 CALL DSYR( 'U', 2, ALPHA, X, 1, A, 1 )
2599 CALL DSPR( '/', 0, ALPHA, X, 1, A )
2602 CALL DSPR( 'U', -1, ALPHA, X, 1, A )
2605 CALL DSPR( 'U', 0, ALPHA, X, 0, A )
2609 CALL DSYR2( '/', 0, ALPHA, X, 1, Y, 1, A, 1 )
2612 CALL DSYR2( 'U', -1, ALPHA, X, 1, Y, 1, A, 1 )
2615 CALL DSYR2( 'U', 0, ALPHA, X, 0, Y, 1, A, 1 )
2618 CALL DSYR2( 'U', 0, ALPHA, X, 1, Y, 0, A, 1 )
2621 CALL DSYR2( 'U', 2, ALPHA, X, 1, Y, 1, A, 1 )
2625 CALL DSPR2( '/', 0, ALPHA, X, 1, Y, 1, A )
2628 CALL DSPR2( 'U', -1, ALPHA, X, 1, Y, 1, A )
2631 CALL DSPR2( 'U', 0, ALPHA, X, 0, Y, 1, A )
2634 CALL DSPR2( '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.0D10 )
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 DMVCH( TRANS, M, N, ALPHA, A, NMAX, X, INCX, BETA, Y, argument
2834 * -- Written on 10-August-1987.
2842 DOUBLE PRECISION ALPHA, BETA, EPS, ERR local
2847 DOUBLE PRECISION 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 DBEG = DBLE( I - 500 )/1001.0D0
3079 DOUBLE PRECISION FUNCTION DDIFF( X, Y )
3083 * -- Written on 10-August-1987.
3087 DOUBLE PRECISION X, Y local
3089 DDIFF = X - Y
3101 * -- Written on 10-August-1987.
3136 * -- Written on 10-August-1987.