1
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 |
SUBROUTINE ZERRTZ( PATH, NUNIT )
* * -- LAPACK test routine (version 3.3.1) -- * Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. * -- April 2011 -- * * .. Scalar Arguments .. CHARACTER*3 PATH INTEGER NUNIT * .. * * Purpose * ======= * * ZERRTZ tests the error exits for ZTZRQF and ZTZRZF. * * Arguments * ========= * * PATH (input) CHARACTER*3 * The LAPACK path name for the routines to be tested. * * NUNIT (input) INTEGER * The unit number for output. * * ===================================================================== * * .. Parameters .. INTEGER NMAX PARAMETER ( NMAX = 2 ) * .. * .. Local Scalars .. CHARACTER*2 C2 INTEGER INFO * .. * .. Local Arrays .. COMPLEX*16 A( NMAX, NMAX ), TAU( NMAX ), W( NMAX ) * .. * .. External Functions .. LOGICAL LSAMEN EXTERNAL LSAMEN * .. * .. External Subroutines .. EXTERNAL ALAESM, CHKXER, ZTZRQF, ZTZRZF * .. * .. Scalars in Common .. LOGICAL LERR, OK CHARACTER*32 SRNAMT INTEGER INFOT, NOUT * .. * .. Common blocks .. COMMON / INFOC / INFOT, NOUT, OK, LERR COMMON / SRNAMC / SRNAMT * .. * .. Intrinsic Functions .. INTRINSIC DCMPLX * .. * .. Executable Statements .. * NOUT = NUNIT C2 = PATH( 2: 3 ) A( 1, 1 ) = DCMPLX( 1.D+0, -1.D+0 ) A( 1, 2 ) = DCMPLX( 2.D+0, -2.D+0 ) A( 2, 2 ) = DCMPLX( 3.D+0, -3.D+0 ) A( 2, 1 ) = DCMPLX( 4.D+0, -4.D+0 ) W( 1 ) = DCMPLX( 0.D+0, 0.D+0 ) W( 2 ) = DCMPLX( 0.D+0, 0.D+0 ) OK = .TRUE. * * Test error exits for the trapezoidal routines. * WRITE( NOUT, FMT = * ) IF( LSAMEN( 2, C2, 'TZ' ) ) THEN * * ZTZRQF * SRNAMT = 'ZTZRQF' INFOT = 1 CALL ZTZRQF( -1, 0, A, 1, TAU, INFO ) CALL CHKXER( 'ZTZRQF', INFOT, NOUT, LERR, OK ) INFOT = 2 CALL ZTZRQF( 1, 0, A, 1, TAU, INFO ) CALL CHKXER( 'ZTZRQF', INFOT, NOUT, LERR, OK ) INFOT = 4 CALL ZTZRQF( 2, 2, A, 1, TAU, INFO ) CALL CHKXER( 'ZTZRQF', INFOT, NOUT, LERR, OK ) * * ZTZRZF * SRNAMT = 'ZTZRZF' INFOT = 1 CALL ZTZRZF( -1, 0, A, 1, TAU, W, 1, INFO ) CALL CHKXER( 'ZTZRZF', INFOT, NOUT, LERR, OK ) INFOT = 2 CALL ZTZRZF( 1, 0, A, 1, TAU, W, 1, INFO ) CALL CHKXER( 'ZTZRZF', INFOT, NOUT, LERR, OK ) INFOT = 4 CALL ZTZRZF( 2, 2, A, 1, TAU, W, 1, INFO ) CALL CHKXER( 'ZTZRZF', INFOT, NOUT, LERR, OK ) INFOT = 7 CALL ZTZRZF( 2, 2, A, 2, TAU, W, 0, INFO ) CALL CHKXER( 'ZTZRZF', INFOT, NOUT, LERR, OK ) INFOT = 7 CALL ZTZRZF( 2, 3, A, 2, TAU, W, 1, INFO ) CALL CHKXER( 'ZTZRZF', INFOT, NOUT, LERR, OK ) END IF * * Print a summary line. * CALL ALAESM( PATH, OK, NOUT ) * RETURN * * End of ZERRTZ * END |