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 17 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,0.0) (1.0,0.0) (0.7,-0.9) VALUES OF ALPHA
43 *> (0.0,0.0) (1.0,0.0) (1.3,-1.1) VALUES OF BETA
79 *> Street, Suite 100, Downers Grove, Illinois 60515-1263, USA.
82 *> -- Written on 10-August-1987.
86 *> 10-9-00: Change STATUS='NEW' to 'UNKNOWN' so that the testers
106 * -- Reference BLAS test routine (version 3.4.1) --
107 * -- Reference BLAS is a software package provided by Univ. of Tennessee, --
108 * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
140 $ X( NMAX ), XS( NMAX*INCMAX ), variable
238 * Values of ALPHA
241 WRITE( NOUT, FMT = 9997 )'ALPHA', NALMAX
298 A( I, J ) = MAX( I - J + 1, 0 )
300 X( J ) = J
304 YY( J ) = J*( ( J + 1 )*J )/2 - ( ( J + 1 )*J*( J - 1 ) )/3
309 CALL ZMVCH( TRANS, N, N, ONE, A, NMAX, X, 1, ZERO, Y, 1, YT, G,
317 CALL ZMVCH( TRANS, N, N, ONE, A, NMAX, X, -1, ZERO, Y, -1, YT, G,
350 $ X, XX, XS, Y, YY, YS, YT, G )
356 $ X, XX, XS, Y, YY, YS, YT, G )
367 $ NMAX, INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS,
373 $ NMAX, INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS,
379 $ NMAX, INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS,
416 9989 FORMAT( ' FOR ALPHA ',
424 9985 FORMAT( ' ERROR IN ZMVCH - IN-LINE DOT PRODUCTS ARE BEING EVALU',
430 9983 FORMAT( 1X, A6, ' WAS NOT TESTED' )
432 9981 FORMAT( /' ******* FATAL ERROR - TESTS ABANDONED *******' )
433 9980 FORMAT( ' ERROR-EXITS WILL NOT BE TESTED' )
440 $ BET, NINC, INC, NMAX, INCMAX, A, AA, AS, X, XX, argument
447 * -- Written on 10-August-1987.
465 $ AS( NMAX*NMAX ), BET( NBET ), X( NMAX ), local
472 COMPLEX*16 ALPHA, ALS, BETA, BLS, TRANSL local
517 $ M = MAX( N - ND, 0 )
529 KL = MAX( KU - 1, 0 )
531 KU = N - 1
532 KL = M - 1
570 * Generate the vector X.
573 CALL ZMAKE( 'GE', ' ', ' ', 1, NL, X, 1, XX,
574 $ ABS( INCX ), 0, NL - 1, RESET, TRANSL )
576 X( NL/2 ) = ZERO
577 XX( 1 + ABS( INCX )*( NL/2 - 1 ) ) = ZERO
585 ALPHA = ALF( IA )
594 $ YY, ABS( INCY ), 0, ML - 1,
607 ALS = ALPHA
627 $ TRANS, M, N, ALPHA, LDA, INCX, BETA,
631 CALL ZGEMV( TRANS, M, N, ALPHA, AA,
637 $ TRANS, M, N, KL, KU, ALPHA, LDA,
641 CALL ZGBMV( TRANS, M, N, KL, KU, ALPHA,
646 * Check if error-exit was taken incorrectly.
660 ISAME( 4 ) = ALS.EQ.ALPHA
677 ISAME( 6 ) = ALS.EQ.ALPHA
711 CALL ZMVCH( TRANS, M, N, ALPHA, A,
712 $ NMAX, X, INCX, BETA, Y,
754 WRITE( NOUT, FMT = 9994 )NC, SNAME, TRANS, M, N, ALPHA, LDA,
758 $ ALPHA, LDA, INCX, BETA, INCY
766 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
770 $ ' - SUSPECT *******' )
772 9995 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', 4( I3, ',' ), '(',
773 $ F4.1, ',', F4.1, '), A,', I3, ', X,', I2, ',(', F4.1, ',',
775 9994 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', 2( I3, ',' ), '(',
776 $ F4.1, ',', F4.1, '), A,', I3, ', X,', I2, ',(', F4.1, ',',
778 9993 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
786 $ BET, NINC, INC, NMAX, INCMAX, A, AA, AS, X, XX, argument
793 * -- Written on 10-August-1987.
811 $ AS( NMAX*NMAX ), BET( NBET ), X( NMAX ), local
818 COMPLEX*16 ALPHA, ALS, BETA, BLS, TRANSL local
871 K = N - 1
904 * Generate the vector X.
907 CALL ZMAKE( 'GE', ' ', ' ', 1, N, X, 1, XX,
908 $ ABS( INCX ), 0, N - 1, RESET, TRANSL )
910 X( N/2 ) = ZERO
911 XX( 1 + ABS( INCX )*( N/2 - 1 ) ) = ZERO
919 ALPHA = ALF( IA )
928 $ ABS( INCY ), 0, N - 1, RESET,
939 ALS = ALPHA
959 $ UPLO, N, ALPHA, LDA, INCX, BETA, INCY
962 CALL ZHEMV( UPLO, N, ALPHA, AA, LDA, XX,
967 $ UPLO, N, K, ALPHA, LDA, INCX, BETA,
971 CALL ZHBMV( UPLO, N, K, ALPHA, AA, LDA,
976 $ UPLO, N, ALPHA, INCX, BETA, INCY
979 CALL ZHPMV( UPLO, N, ALPHA, AA, XX, INCX,
983 * Check if error-exit was taken incorrectly.
996 ISAME( 3 ) = ALS.EQ.ALPHA
1011 ISAME( 4 ) = ALS.EQ.ALPHA
1025 ISAME( 3 ) = ALS.EQ.ALPHA
1057 CALL ZMVCH( 'N', N, N, ALPHA, A, NMAX, X,
1097 WRITE( NOUT, FMT = 9993 )NC, SNAME, UPLO, N, ALPHA, LDA, INCX,
1100 WRITE( NOUT, FMT = 9994 )NC, SNAME, UPLO, N, K, ALPHA, LDA,
1103 WRITE( NOUT, FMT = 9995 )NC, SNAME, UPLO, N, ALPHA, INCX,
1112 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
1116 $ ' - SUSPECT *******' )
1118 9995 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',(', F4.1, ',',
1119 $ F4.1, '), AP, X,', I2, ',(', F4.1, ',', F4.1, '), Y,', I2,
1121 9994 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', 2( I3, ',' ), '(',
1122 $ F4.1, ',', F4.1, '), A,', I3, ', X,', I2, ',(', F4.1, ',',
1124 9993 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',(', F4.1, ',',
1125 $ F4.1, '), A,', I3, ', X,', I2, ',(', F4.1, ',', F4.1, '), ',
1127 9992 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
1135 $ INCMAX, A, AA, AS, X, XX, XS, XT, G, Z ) argument
1141 * -- Written on 10-August-1987.
1159 $ AS( NMAX*NMAX ), X( NMAX ), XS( NMAX*INCMAX ), local
1222 K = N - 1
1261 * Generate the vector X.
1264 CALL ZMAKE( 'GE', ' ', ' ', 1, N, X, 1, XX,
1265 $ ABS( INCX ), 0, N - 1, RESET,
1268 X( N/2 ) = ZERO
1269 XX( 1 + ABS( INCX )*( N/2 - 1 ) ) = ZERO
1346 * Check if error-exit was taken incorrectly.
1411 CALL ZMVCH( TRANS, N, N, ONE, A, NMAX, X,
1420 Z( I ) = XX( 1 + ( I - 1 )*
1422 XX( 1 + ( I - 1 )*ABS( INCX ) )
1423 $ = X( I )
1426 $ INCX, ZERO, X, INCX, XT, G,
1477 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
1481 $ ' - SUSPECT *******' )
1483 9995 FORMAT( 1X, I6, ': ', A6, '(', 3( '''', A1, ''',' ), I3, ', AP, ',
1484 $ 'X,', I2, ') .' )
1485 9994 FORMAT( 1X, I6, ': ', A6, '(', 3( '''', A1, ''',' ), 2( I3, ',' ),
1486 $ ' A,', I3, ', X,', I2, ') .' )
1487 9993 FORMAT( 1X, I6, ': ', A6, '(', 3( '''', A1, ''',' ), I3, ', A,',
1488 $ I3, ', X,', I2, ') .' )
1489 9992 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
1497 $ INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS, YT, G, argument
1504 * -- Written on 10-August-1987.
1522 $ AS( NMAX*NMAX ), X( NMAX ), XS( NMAX*INCMAX ), local
1529 COMPLEX*16 ALPHA, ALS, TRANSL local
1565 $ M = MAX( N - ND, 0 )
1583 * Generate the vector X.
1586 CALL ZMAKE( 'GE', ' ', ' ', 1, M, X, 1, XX, ABS( INCX ),
1587 $ 0, M - 1, RESET, TRANSL )
1589 X( M/2 ) = ZERO
1590 XX( 1 + ABS( INCX )*( M/2 - 1 ) ) = ZERO
1601 $ ABS( INCY ), 0, N - 1, RESET, TRANSL )
1604 YY( 1 + ABS( INCY )*( N/2 - 1 ) ) = ZERO
1608 ALPHA = ALF( IA )
1614 $ AA, LDA, M - 1, N - 1, RESET, TRANSL )
1622 ALS = ALPHA
1640 $ ALPHA, INCX, INCY, LDA
1644 CALL ZGERC( M, N, ALPHA, XX, INCX, YY, INCY, AA,
1649 CALL ZGERU( M, N, ALPHA, XX, INCX, YY, INCY, AA,
1653 * Check if error-exit was taken incorrectly.
1665 ISAME( 3 ) = ALS.EQ.ALPHA
1697 Z( I ) = X( I )
1701 Z( I ) = X( M - I + 1 )
1708 W( 1 ) = Y( N - J + 1 )
1712 CALL ZMVCH( 'N', M, 1, ALPHA, Z, NMAX, W, 1,
1714 $ AA( 1 + ( J - 1 )*LDA ), EPS,
1750 WRITE( NOUT, FMT = 9994 )NC, SNAME, M, N, ALPHA, INCX, INCY, LDA
1757 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
1761 $ ' - SUSPECT *******' )
1764 9994 FORMAT( 1X, I6, ': ', A6, '(', 2( I3, ',' ), '(', F4.1, ',', F4.1,
1765 $ '), X,', I2, ', Y,', I2, ', A,', I3, ') ',
1767 9993 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
1775 $ INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS, YT, G, argument
1782 * -- Written on 10-August-1987.
1800 $ AS( NMAX*NMAX ), X( NMAX ), XS( NMAX*INCMAX ), local
1807 COMPLEX*16 ALPHA, TRANSL local
1868 * Generate the vector X.
1871 CALL ZMAKE( 'GE', ' ', ' ', 1, N, X, 1, XX, ABS( INCX ),
1872 $ 0, N - 1, RESET, TRANSL )
1874 X( N/2 ) = ZERO
1875 XX( 1 + ABS( INCX )*( N/2 - 1 ) ) = ZERO
1880 ALPHA = DCMPLX( RALPHA, RZERO )
1887 $ AA, LDA, N - 1, N - 1, RESET, TRANSL )
1923 * Check if error-exit was taken incorrectly.
1967 Z( I ) = X( I )
1971 Z( I ) = X( N - I + 1 )
1982 LJ = N - J + 1
1984 CALL ZMVCH( 'N', LJ, 1, ALPHA, Z( JJ ), LJ, W,
2041 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
2045 $ ' - SUSPECT *******' )
2048 9994 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',', F4.1, ', X,',
2050 9993 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',', F4.1, ', X,',
2052 9992 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
2060 $ INCMAX, A, AA, AS, X, XX, XS, Y, YY, YS, YT, G, argument
2067 * -- Written on 10-August-1987.
2085 $ AS( NMAX*NMAX ), X( NMAX ), XS( NMAX*INCMAX ), local
2092 COMPLEX*16 ALPHA, ALS, TRANSL local
2154 * Generate the vector X.
2157 CALL ZMAKE( 'GE', ' ', ' ', 1, N, X, 1, XX, ABS( INCX ),
2158 $ 0, N - 1, RESET, TRANSL )
2160 X( N/2 ) = ZERO
2161 XX( 1 + ABS( INCX )*( N/2 - 1 ) ) = ZERO
2172 $ ABS( INCY ), 0, N - 1, RESET, TRANSL )
2175 YY( 1 + ABS( INCY )*( N/2 - 1 ) ) = ZERO
2179 ALPHA = ALF( IA )
2180 NULL = N.LE.0.OR.ALPHA.EQ.ZERO
2186 $ NMAX, AA, LDA, N - 1, N - 1, RESET,
2195 ALS = ALPHA
2214 $ ALPHA, INCX, INCY, LDA
2217 CALL ZHER2( UPLO, N, ALPHA, XX, INCX, YY, INCY,
2222 $ ALPHA, INCX, INCY
2225 CALL ZHPR2( UPLO, N, ALPHA, XX, INCX, YY, INCY,
2229 * Check if error-exit was taken incorrectly.
2241 ISAME( 3 ) = ALS.EQ.ALPHA
2275 Z( I, 1 ) = X( I )
2279 Z( I, 1 ) = X( N - I + 1 )
2288 Z( I, 2 ) = Y( N - I + 1 )
2293 W( 1 ) = ALPHA*DCONJG( Z( J, 2 ) )
2294 W( 2 ) = DCONJG( ALPHA )*DCONJG( Z( J, 1 ) )
2300 LJ = N - J + 1
2351 WRITE( NOUT, FMT = 9993 )NC, SNAME, UPLO, N, ALPHA, INCX,
2354 WRITE( NOUT, FMT = 9994 )NC, SNAME, UPLO, N, ALPHA, INCX, INCY
2362 9998 FORMAT( ' ******* FATAL ERROR - PARAMETER NUMBER ', I2, ' WAS CH',
2366 $ ' - SUSPECT *******' )
2369 9994 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',(', F4.1, ',',
2370 $ F4.1, '), X,', I2, ', Y,', I2, ', AP) ',
2372 9993 FORMAT( 1X, I6, ': ', A6, '(''', A1, ''',', I3, ',(', F4.1, ',',
2373 $ F4.1, '), X,', I2, ', Y,', I2, ', A,', I3, ') ',
2375 9992 FORMAT( ' ******* FATAL ERROR - ERROR-EXIT TAKEN ON VALID CALL *',
2384 * Requires a special version of the error-handling routine XERBLA.
2385 * ALPHA, RALPHA, BETA, A, X and Y should not need to be defined.
2389 * -- Written on 10-August-1987.
2400 COMPLEX*16 ALPHA, BETA local
2403 COMPLEX*16 A( 1, 1 ), X( 1 ), Y( 1 ) local
2415 * it is called, and is then tested and re-set by CHKXER.
2421 CALL ZGEMV( '/', 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2424 CALL ZGEMV( 'N', -1, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2427 CALL ZGEMV( 'N', 0, -1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2430 CALL ZGEMV( 'N', 2, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2433 CALL ZGEMV( 'N', 0, 0, ALPHA, A, 1, X, 0, BETA, Y, 1 )
2436 CALL ZGEMV( 'N', 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 0 )
2440 CALL ZGBMV( '/', 0, 0, 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2443 CALL ZGBMV( 'N', -1, 0, 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2446 CALL ZGBMV( 'N', 0, -1, 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2449 CALL ZGBMV( 'N', 0, 0, -1, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2452 CALL ZGBMV( 'N', 2, 0, 0, -1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2455 CALL ZGBMV( 'N', 0, 0, 1, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2458 CALL ZGBMV( 'N', 0, 0, 0, 0, ALPHA, A, 1, X, 0, BETA, Y, 1 )
2461 CALL ZGBMV( 'N', 0, 0, 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 0 )
2465 CALL ZHEMV( '/', 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2468 CALL ZHEMV( 'U', -1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2471 CALL ZHEMV( 'U', 2, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2474 CALL ZHEMV( 'U', 0, ALPHA, A, 1, X, 0, BETA, Y, 1 )
2477 CALL ZHEMV( 'U', 0, ALPHA, A, 1, X, 1, BETA, Y, 0 )
2481 CALL ZHBMV( '/', 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2484 CALL ZHBMV( 'U', -1, 0, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2487 CALL ZHBMV( 'U', 0, -1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2490 CALL ZHBMV( 'U', 0, 1, ALPHA, A, 1, X, 1, BETA, Y, 1 )
2493 CALL ZHBMV( 'U', 0, 0, ALPHA, A, 1, X, 0, BETA, Y, 1 )
2496 CALL ZHBMV( 'U', 0, 0, ALPHA, A, 1, X, 1, BETA, Y, 0 )
2500 CALL ZHPMV( '/', 0, ALPHA, A, X, 1, BETA, Y, 1 )
2503 CALL ZHPMV( 'U', -1, ALPHA, A, X, 1, BETA, Y, 1 )
2506 CALL ZHPMV( 'U', 0, ALPHA, A, X, 0, BETA, Y, 1 )
2509 CALL ZHPMV( 'U', 0, ALPHA, A, X, 1, BETA, Y, 0 )
2513 CALL ZTRMV( '/', 'N', 'N', 0, A, 1, X, 1 )
2516 CALL ZTRMV( 'U', '/', 'N', 0, A, 1, X, 1 )
2519 CALL ZTRMV( 'U', 'N', '/', 0, A, 1, X, 1 )
2522 CALL ZTRMV( 'U', 'N', 'N', -1, A, 1, X, 1 )
2525 CALL ZTRMV( 'U', 'N', 'N', 2, A, 1, X, 1 )
2528 CALL ZTRMV( 'U', 'N', 'N', 0, A, 1, X, 0 )
2532 CALL ZTBMV( '/', 'N', 'N', 0, 0, A, 1, X, 1 )
2535 CALL ZTBMV( 'U', '/', 'N', 0, 0, A, 1, X, 1 )
2538 CALL ZTBMV( 'U', 'N', '/', 0, 0, A, 1, X, 1 )
2541 CALL ZTBMV( 'U', 'N', 'N', -1, 0, A, 1, X, 1 )
2544 CALL ZTBMV( 'U', 'N', 'N', 0, -1, A, 1, X, 1 )
2547 CALL ZTBMV( 'U', 'N', 'N', 0, 1, A, 1, X, 1 )
2550 CALL ZTBMV( 'U', 'N', 'N', 0, 0, A, 1, X, 0 )
2554 CALL ZTPMV( '/', 'N', 'N', 0, A, X, 1 )
2557 CALL ZTPMV( 'U', '/', 'N', 0, A, X, 1 )
2560 CALL ZTPMV( 'U', 'N', '/', 0, A, X, 1 )
2563 CALL ZTPMV( 'U', 'N', 'N', -1, A, X, 1 )
2566 CALL ZTPMV( 'U', 'N', 'N', 0, A, X, 0 )
2570 CALL ZTRSV( '/', 'N', 'N', 0, A, 1, X, 1 )
2573 CALL ZTRSV( 'U', '/', 'N', 0, A, 1, X, 1 )
2576 CALL ZTRSV( 'U', 'N', '/', 0, A, 1, X, 1 )
2579 CALL ZTRSV( 'U', 'N', 'N', -1, A, 1, X, 1 )
2582 CALL ZTRSV( 'U', 'N', 'N', 2, A, 1, X, 1 )
2585 CALL ZTRSV( 'U', 'N', 'N', 0, A, 1, X, 0 )
2589 CALL ZTBSV( '/', 'N', 'N', 0, 0, A, 1, X, 1 )
2592 CALL ZTBSV( 'U', '/', 'N', 0, 0, A, 1, X, 1 )
2595 CALL ZTBSV( 'U', 'N', '/', 0, 0, A, 1, X, 1 )
2598 CALL ZTBSV( 'U', 'N', 'N', -1, 0, A, 1, X, 1 )
2601 CALL ZTBSV( 'U', 'N', 'N', 0, -1, A, 1, X, 1 )
2604 CALL ZTBSV( 'U', 'N', 'N', 0, 1, A, 1, X, 1 )
2607 CALL ZTBSV( 'U', 'N', 'N', 0, 0, A, 1, X, 0 )
2611 CALL ZTPSV( '/', 'N', 'N', 0, A, X, 1 )
2614 CALL ZTPSV( 'U', '/', 'N', 0, A, X, 1 )
2617 CALL ZTPSV( 'U', 'N', '/', 0, A, X, 1 )
2620 CALL ZTPSV( 'U', 'N', 'N', -1, A, X, 1 )
2623 CALL ZTPSV( 'U', 'N', 'N', 0, A, X, 0 )
2627 CALL ZGERC( -1, 0, ALPHA, X, 1, Y, 1, A, 1 )
2630 CALL ZGERC( 0, -1, ALPHA, X, 1, Y, 1, A, 1 )
2633 CALL ZGERC( 0, 0, ALPHA, X, 0, Y, 1, A, 1 )
2636 CALL ZGERC( 0, 0, ALPHA, X, 1, Y, 0, A, 1 )
2639 CALL ZGERC( 2, 0, ALPHA, X, 1, Y, 1, A, 1 )
2643 CALL ZGERU( -1, 0, ALPHA, X, 1, Y, 1, A, 1 )
2646 CALL ZGERU( 0, -1, ALPHA, X, 1, Y, 1, A, 1 )
2649 CALL ZGERU( 0, 0, ALPHA, X, 0, Y, 1, A, 1 )
2652 CALL ZGERU( 0, 0, ALPHA, X, 1, Y, 0, A, 1 )
2655 CALL ZGERU( 2, 0, ALPHA, X, 1, Y, 1, A, 1 )
2659 CALL ZHER( '/', 0, RALPHA, X, 1, A, 1 )
2662 CALL ZHER( 'U', -1, RALPHA, X, 1, A, 1 )
2665 CALL ZHER( 'U', 0, RALPHA, X, 0, A, 1 )
2668 CALL ZHER( 'U', 2, RALPHA, X, 1, A, 1 )
2672 CALL ZHPR( '/', 0, RALPHA, X, 1, A )
2675 CALL ZHPR( 'U', -1, RALPHA, X, 1, A )
2678 CALL ZHPR( 'U', 0, RALPHA, X, 0, A )
2682 CALL ZHER2( '/', 0, ALPHA, X, 1, Y, 1, A, 1 )
2685 CALL ZHER2( 'U', -1, ALPHA, X, 1, Y, 1, A, 1 )
2688 CALL ZHER2( 'U', 0, ALPHA, X, 0, Y, 1, A, 1 )
2691 CALL ZHER2( 'U', 0, ALPHA, X, 1, Y, 0, A, 1 )
2694 CALL ZHER2( 'U', 2, ALPHA, X, 1, Y, 1, A, 1 )
2698 CALL ZHPR2( '/', 0, ALPHA, X, 1, Y, 1, A )
2701 CALL ZHPR2( 'U', -1, ALPHA, X, 1, Y, 1, A )
2704 CALL ZHPR2( 'U', 0, ALPHA, X, 0, Y, 1, A )
2707 CALL ZHPR2( 'U', 0, ALPHA, X, 1, Y, 0, A )
2717 9999 FORMAT( ' ', A6, ' PASSED THE TESTS OF ERROR-EXITS' )
2718 9998 FORMAT( ' ******* ', A6, ' FAILED THE TESTS OF ERROR-EXITS *****',
2736 * -- Written on 10-August-1987.
2745 PARAMETER ( ROGUE = ( -1.0D10, 1.0D10 ) )
2749 PARAMETER ( RROGUE = -1.0D10 )
2780 IF( ( I.LE.J.AND.J - I.LE.KU ).OR.
2781 $ ( I.GE.J.AND.I - J.LE.KL ) )THEN
2808 AA( I + ( J - 1 )*LDA ) = A( I, J )
2811 AA( I + ( J - 1 )*LDA ) = ROGUE
2816 DO 60 I1 = 1, KU + 1 - J
2817 AA( I1 + ( J - 1 )*LDA ) = ROGUE
2819 DO 70 I2 = I1, MIN( KL + KU + 1, KU + 1 + M - J )
2820 AA( I2 + ( J - 1 )*LDA ) = A( I2 + J - KU - 1, J )
2823 AA( I3 + ( J - 1 )*LDA ) = ROGUE
2831 IEND = J - 1
2843 DO 100 I = 1, IBEG - 1
2844 AA( I + ( J - 1 )*LDA ) = ROGUE
2847 AA( I + ( J - 1 )*LDA ) = A( I, J )
2850 AA( I + ( J - 1 )*LDA ) = ROGUE
2853 JJ = J + ( J - 1 )*LDA
2861 IBEG = MAX( 1, KL + 2 - J )
2874 IEND = MIN( KL + 1, 1 + M - J )
2876 DO 140 I = 1, IBEG - 1
2877 AA( I + ( J - 1 )*LDA ) = ROGUE
2880 AA( I + ( J - 1 )*LDA ) = A( I + J - KK, J )
2883 AA( I + ( J - 1 )*LDA ) = ROGUE
2886 JJ = KK + ( J - 1 )*LDA
2917 SUBROUTINE ZMVCH( TRANS, M, N, ALPHA, A, NMAX, X, INCX, BETA, Y, argument
2924 * -- Written on 10-August-1987.
2934 COMPLEX*16 ALPHA, BETA local
2940 COMPLEX*16 A( NMAX, * ), X( * ), Y( * ), YT( * ), YY( * ) local
2965 INCXL = -1
2972 INCYL = -1
2978 * Compute expected result in YT using data in A, X and Y.
2988 YT( IY ) = YT( IY ) + A( J, I )*X( JX )
2989 G( IY ) = G( IY ) + ABS1( A( J, I ) )*ABS1( X( JX ) )
2994 YT( IY ) = YT( IY ) + DCONJG( A( J, I ) )*X( JX )
2995 G( IY ) = G( IY ) + ABS1( A( J, I ) )*ABS1( X( JX ) )
3000 YT( IY ) = YT( IY ) + A( I, J )*X( JX )
3001 G( IY ) = G( IY ) + ABS1( A( I, J ) )*ABS1( X( JX ) )
3005 YT( IY ) = ALPHA*YT( IY ) + BETA*Y( IY )
3006 G( IY ) = ABS1( ALPHA )*G( IY ) + ABS1( BETA )*ABS1( Y( IY ) )
3014 ERRI = ABS( YT( I ) - YY( 1 + ( I - 1 )*ABS( INCY ) ) )/EPS
3031 $ YY( 1 + ( I - 1 )*ABS( INCY ) )
3034 $ YY( 1 + ( I - 1 )*ABS( INCY ) ), YT( I )
3041 9999 FORMAT( ' ******* FATAL ERROR - COMPUTED RESULT IS LESS THAN HAL',
3044 9998 FORMAT( 1X, I7, 2( ' (', G15.6, ',', G15.6, ')' ) )
3055 * -- Written on 10-August-1987.
3087 * -- Written on 10-August-1987.
3118 DO 30 I = 1, IBEG - 1
3141 * distributed between -0.5 and 0.5.
3145 * -- Written on 10-August-1987.
3178 I = I - 1000*( I/1000 )
3179 J = J - 1000*( J/1000 )
3184 ZBEG = DCMPLX( ( I - 500 )/1001.0D0, ( J - 500 )/1001.0D0 )
3190 DOUBLE PRECISION FUNCTION DDIFF( X, Y )
3194 * -- Written on 10-August-1987.
3198 DOUBLE PRECISION X, Y local
3200 DDIFF = X - Y
3212 * -- Written on 10-August-1987.
3247 * -- Written on 10-August-1987.