LAPACK 3.3.0

zdrvst.f

Go to the documentation of this file.
00001       SUBROUTINE ZDRVST( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
00002      $                   NOUNIT, A, LDA, D1, D2, D3, WA1, WA2, WA3, U,
00003      $                   LDU, V, TAU, Z, WORK, LWORK, RWORK, LRWORK,
00004      $                   IWORK, LIWORK, RESULT, INFO )
00005 *
00006 *  -- LAPACK test routine (version 3.1) --
00007 *     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd..
00008 *     November 2006
00009 *
00010 *     .. Scalar Arguments ..
00011       INTEGER            INFO, LDA, LDU, LIWORK, LRWORK, LWORK, NOUNIT,
00012      $                   NSIZES, NTYPES
00013       DOUBLE PRECISION   THRESH
00014 *     ..
00015 *     .. Array Arguments ..
00016       LOGICAL            DOTYPE( * )
00017       INTEGER            ISEED( 4 ), IWORK( * ), NN( * )
00018       DOUBLE PRECISION   D1( * ), D2( * ), D3( * ), RESULT( * ),
00019      $                   RWORK( * ), WA1( * ), WA2( * ), WA3( * )
00020       COMPLEX*16         A( LDA, * ), TAU( * ), U( LDU, * ),
00021      $                   V( LDU, * ), WORK( * ), Z( LDU, * )
00022 *     ..
00023 *
00024 *  Purpose
00025 *  =======
00026 *
00027 *       ZDRVST  checks the Hermitian eigenvalue problem drivers.
00028 *
00029 *               ZHEEVD computes all eigenvalues and, optionally,
00030 *               eigenvectors of a complex Hermitian matrix,
00031 *               using a divide-and-conquer algorithm.
00032 *
00033 *               ZHEEVX computes selected eigenvalues and, optionally,
00034 *               eigenvectors of a complex Hermitian matrix.
00035 *
00036 *               ZHEEVR computes selected eigenvalues and, optionally,
00037 *               eigenvectors of a complex Hermitian matrix
00038 *               using the Relatively Robust Representation where it can.
00039 *
00040 *               ZHPEVD computes all eigenvalues and, optionally,
00041 *               eigenvectors of a complex Hermitian matrix in packed
00042 *               storage, using a divide-and-conquer algorithm.
00043 *
00044 *               ZHPEVX computes selected eigenvalues and, optionally,
00045 *               eigenvectors of a complex Hermitian matrix in packed
00046 *               storage.
00047 *
00048 *               ZHBEVD computes all eigenvalues and, optionally,
00049 *               eigenvectors of a complex Hermitian band matrix,
00050 *               using a divide-and-conquer algorithm.
00051 *
00052 *               ZHBEVX computes selected eigenvalues and, optionally,
00053 *               eigenvectors of a complex Hermitian band matrix.
00054 *
00055 *               ZHEEV computes all eigenvalues and, optionally,
00056 *               eigenvectors of a complex Hermitian matrix.
00057 *
00058 *               ZHPEV computes all eigenvalues and, optionally,
00059 *               eigenvectors of a complex Hermitian matrix in packed
00060 *               storage.
00061 *
00062 *               ZHBEV computes all eigenvalues and, optionally,
00063 *               eigenvectors of a complex Hermitian band matrix.
00064 *
00065 *       When ZDRVST is called, a number of matrix "sizes" ("n's") and a
00066 *       number of matrix "types" are specified.  For each size ("n")
00067 *       and each type of matrix, one matrix will be generated and used
00068 *       to test the appropriate drivers.  For each matrix and each
00069 *       driver routine called, the following tests will be performed:
00070 *
00071 *       (1)     | A - Z D Z' | / ( |A| n ulp )
00072 *
00073 *       (2)     | I - Z Z' | / ( n ulp )
00074 *
00075 *       (3)     | D1 - D2 | / ( |D1| ulp )
00076 *
00077 *       where Z is the matrix of eigenvectors returned when the
00078 *       eigenvector option is given and D1 and D2 are the eigenvalues
00079 *       returned with and without the eigenvector option.
00080 *
00081 *       The "sizes" are specified by an array NN(1:NSIZES); the value of
00082 *       each element NN(j) specifies one size.
00083 *       The "types" are specified by a logical array DOTYPE( 1:NTYPES );
00084 *       if DOTYPE(j) is .TRUE., then matrix type "j" will be generated.
00085 *       Currently, the list of possible types is:
00086 *
00087 *       (1)  The zero matrix.
00088 *       (2)  The identity matrix.
00089 *
00090 *       (3)  A diagonal matrix with evenly spaced entries
00091 *            1, ..., ULP  and random signs.
00092 *            (ULP = (first number larger than 1) - 1 )
00093 *       (4)  A diagonal matrix with geometrically spaced entries
00094 *            1, ..., ULP  and random signs.
00095 *       (5)  A diagonal matrix with "clustered" entries 1, ULP, ..., ULP
00096 *            and random signs.
00097 *
00098 *       (6)  Same as (4), but multiplied by SQRT( overflow threshold )
00099 *       (7)  Same as (4), but multiplied by SQRT( underflow threshold )
00100 *
00101 *       (8)  A matrix of the form  U* D U, where U is unitary and
00102 *            D has evenly spaced entries 1, ..., ULP with random signs
00103 *            on the diagonal.
00104 *
00105 *       (9)  A matrix of the form  U* D U, where U is unitary and
00106 *            D has geometrically spaced entries 1, ..., ULP with random
00107 *            signs on the diagonal.
00108 *
00109 *       (10) A matrix of the form  U* D U, where U is unitary and
00110 *            D has "clustered" entries 1, ULP,..., ULP with random
00111 *            signs on the diagonal.
00112 *
00113 *       (11) Same as (8), but multiplied by SQRT( overflow threshold )
00114 *       (12) Same as (8), but multiplied by SQRT( underflow threshold )
00115 *
00116 *       (13) Symmetric matrix with random entries chosen from (-1,1).
00117 *       (14) Same as (13), but multiplied by SQRT( overflow threshold )
00118 *       (15) Same as (13), but multiplied by SQRT( underflow threshold )
00119 *       (16) A band matrix with half bandwidth randomly chosen between
00120 *            0 and N-1, with evenly spaced eigenvalues 1, ..., ULP
00121 *            with random signs.
00122 *       (17) Same as (16), but multiplied by SQRT( overflow threshold )
00123 *       (18) Same as (16), but multiplied by SQRT( underflow threshold )
00124 *
00125 *  Arguments
00126 *  =========
00127 *
00128 *  NSIZES  INTEGER
00129 *          The number of sizes of matrices to use.  If it is zero,
00130 *          ZDRVST does nothing.  It must be at least zero.
00131 *          Not modified.
00132 *
00133 *  NN      INTEGER array, dimension (NSIZES)
00134 *          An array containing the sizes to be used for the matrices.
00135 *          Zero values will be skipped.  The values must be at least
00136 *          zero.
00137 *          Not modified.
00138 *
00139 *  NTYPES  INTEGER
00140 *          The number of elements in DOTYPE.   If it is zero, ZDRVST
00141 *          does nothing.  It must be at least zero.  If it is MAXTYP+1
00142 *          and NSIZES is 1, then an additional type, MAXTYP+1 is
00143 *          defined, which is to use whatever matrix is in A.  This
00144 *          is only useful if DOTYPE(1:MAXTYP) is .FALSE. and
00145 *          DOTYPE(MAXTYP+1) is .TRUE. .
00146 *          Not modified.
00147 *
00148 *  DOTYPE  LOGICAL array, dimension (NTYPES)
00149 *          If DOTYPE(j) is .TRUE., then for each size in NN a
00150 *          matrix of that size and of type j will be generated.
00151 *          If NTYPES is smaller than the maximum number of types
00152 *          defined (PARAMETER MAXTYP), then types NTYPES+1 through
00153 *          MAXTYP will not be generated.  If NTYPES is larger
00154 *          than MAXTYP, DOTYPE(MAXTYP+1) through DOTYPE(NTYPES)
00155 *          will be ignored.
00156 *          Not modified.
00157 *
00158 *  ISEED   INTEGER array, dimension (4)
00159 *          On entry ISEED specifies the seed of the random number
00160 *          generator. The array elements should be between 0 and 4095;
00161 *          if not they will be reduced mod 4096.  Also, ISEED(4) must
00162 *          be odd.  The random number generator uses a linear
00163 *          congruential sequence limited to small integers, and so
00164 *          should produce machine independent random numbers. The
00165 *          values of ISEED are changed on exit, and can be used in the
00166 *          next call to ZDRVST to continue the same random number
00167 *          sequence.
00168 *          Modified.
00169 *
00170 *  THRESH  DOUBLE PRECISION
00171 *          A test will count as "failed" if the "error", computed as
00172 *          described above, exceeds THRESH.  Note that the error
00173 *          is scaled to be O(1), so THRESH should be a reasonably
00174 *          small multiple of 1, e.g., 10 or 100.  In particular,
00175 *          it should not depend on the precision (single vs. double)
00176 *          or the size of the matrix.  It must be at least zero.
00177 *          Not modified.
00178 *
00179 *  NOUNIT  INTEGER
00180 *          The FORTRAN unit number for printing out error messages
00181 *          (e.g., if a routine returns IINFO not equal to 0.)
00182 *          Not modified.
00183 *
00184 *  A       COMPLEX*16 array, dimension (LDA , max(NN))
00185 *          Used to hold the matrix whose eigenvalues are to be
00186 *          computed.  On exit, A contains the last matrix actually
00187 *          used.
00188 *          Modified.
00189 *
00190 *  LDA     INTEGER
00191 *          The leading dimension of A.  It must be at
00192 *          least 1 and at least max( NN ).
00193 *          Not modified.
00194 *
00195 *  D1      DOUBLE PRECISION array, dimension (max(NN))
00196 *          The eigenvalues of A, as computed by ZSTEQR simlutaneously
00197 *          with Z.  On exit, the eigenvalues in D1 correspond with the
00198 *          matrix in A.
00199 *          Modified.
00200 *
00201 *  D2      DOUBLE PRECISION array, dimension (max(NN))
00202 *          The eigenvalues of A, as computed by ZSTEQR if Z is not
00203 *          computed.  On exit, the eigenvalues in D2 correspond with
00204 *          the matrix in A.
00205 *          Modified.
00206 *
00207 *  D3      DOUBLE PRECISION array, dimension (max(NN))
00208 *          The eigenvalues of A, as computed by DSTERF.  On exit, the
00209 *          eigenvalues in D3 correspond with the matrix in A.
00210 *          Modified.
00211 *
00212 *  WA1     DOUBLE PRECISION array, dimension
00213 *
00214 *  WA2     DOUBLE PRECISION array, dimension
00215 *
00216 *  WA3     DOUBLE PRECISION array, dimension
00217 *
00218 *  U       COMPLEX*16 array, dimension (LDU, max(NN))
00219 *          The unitary matrix computed by ZHETRD + ZUNGC3.
00220 *          Modified.
00221 *
00222 *  LDU     INTEGER
00223 *          The leading dimension of U, Z, and V.  It must be at
00224 *          least 1 and at least max( NN ).
00225 *          Not modified.
00226 *
00227 *  V       COMPLEX*16 array, dimension (LDU, max(NN))
00228 *          The Housholder vectors computed by ZHETRD in reducing A to
00229 *          tridiagonal form.
00230 *          Modified.
00231 *
00232 *  TAU     COMPLEX*16 array, dimension (max(NN))
00233 *          The Householder factors computed by ZHETRD in reducing A
00234 *          to tridiagonal form.
00235 *          Modified.
00236 *
00237 *  Z       COMPLEX*16 array, dimension (LDU, max(NN))
00238 *          The unitary matrix of eigenvectors computed by ZHEEVD,
00239 *          ZHEEVX, ZHPEVD, CHPEVX, ZHBEVD, and CHBEVX.
00240 *          Modified.
00241 *
00242 *  WORK  - COMPLEX*16 array of dimension ( LWORK )
00243 *           Workspace.
00244 *           Modified.
00245 *
00246 *  LWORK - INTEGER
00247 *           The number of entries in WORK.  This must be at least
00248 *           2*max( NN(j), 2 )**2.
00249 *           Not modified.
00250 *
00251 *  RWORK   DOUBLE PRECISION array, dimension (3*max(NN))
00252 *           Workspace.
00253 *           Modified.
00254 *
00255 *  LRWORK - INTEGER
00256 *           The number of entries in RWORK.
00257 *
00258 *  IWORK   INTEGER array, dimension (6*max(NN))
00259 *          Workspace.
00260 *          Modified.
00261 *
00262 *  LIWORK - INTEGER
00263 *           The number of entries in IWORK.
00264 *
00265 *  RESULT  DOUBLE PRECISION array, dimension (??)
00266 *          The values computed by the tests described above.
00267 *          The values are currently limited to 1/ulp, to avoid
00268 *          overflow.
00269 *          Modified.
00270 *
00271 *  INFO    INTEGER
00272 *          If 0, then everything ran OK.
00273 *           -1: NSIZES < 0
00274 *           -2: Some NN(j) < 0
00275 *           -3: NTYPES < 0
00276 *           -5: THRESH < 0
00277 *           -9: LDA < 1 or LDA < NMAX, where NMAX is max( NN(j) ).
00278 *          -16: LDU < 1 or LDU < NMAX.
00279 *          -21: LWORK too small.
00280 *          If  DLATMR, SLATMS, ZHETRD, DORGC3, ZSTEQR, DSTERF,
00281 *              or DORMC2 returns an error code, the
00282 *              absolute value of it is returned.
00283 *          Modified.
00284 *
00285 *-----------------------------------------------------------------------
00286 *
00287 *       Some Local Variables and Parameters:
00288 *       ---- ----- --------- --- ----------
00289 *       ZERO, ONE       Real 0 and 1.
00290 *       MAXTYP          The number of types defined.
00291 *       NTEST           The number of tests performed, or which can
00292 *                       be performed so far, for the current matrix.
00293 *       NTESTT          The total number of tests performed so far.
00294 *       NMAX            Largest value in NN.
00295 *       NMATS           The number of matrices generated so far.
00296 *       NERRS           The number of tests which have exceeded THRESH
00297 *                       so far (computed by DLAFTS).
00298 *       COND, IMODE     Values to be passed to the matrix generators.
00299 *       ANORM           Norm of A; passed to matrix generators.
00300 *
00301 *       OVFL, UNFL      Overflow and underflow thresholds.
00302 *       ULP, ULPINV     Finest relative precision and its inverse.
00303 *       RTOVFL, RTUNFL  Square roots of the previous 2 values.
00304 *               The following four arrays decode JTYPE:
00305 *       KTYPE(j)        The general type (1-10) for type "j".
00306 *       KMODE(j)        The MODE value to be passed to the matrix
00307 *                       generator for type "j".
00308 *       KMAGN(j)        The order of magnitude ( O(1),
00309 *                       O(overflow^(1/2) ), O(underflow^(1/2) )
00310 *
00311 *  =====================================================================
00312 *
00313 *
00314 *     .. Parameters ..
00315       DOUBLE PRECISION   ZERO, ONE, TWO, TEN
00316       PARAMETER          ( ZERO = 0.0D+0, ONE = 1.0D+0, TWO = 2.0D+0,
00317      $                   TEN = 10.0D+0 )
00318       DOUBLE PRECISION   HALF
00319       PARAMETER          ( HALF = ONE / TWO )
00320       COMPLEX*16         CZERO, CONE
00321       PARAMETER          ( CZERO = ( 0.0D+0, 0.0D+0 ),
00322      $                   CONE = ( 1.0D+0, 0.0D+0 ) )
00323       INTEGER            MAXTYP
00324       PARAMETER          ( MAXTYP = 18 )
00325 *     ..
00326 *     .. Local Scalars ..
00327       LOGICAL            BADNN
00328       CHARACTER          UPLO
00329       INTEGER            I, IDIAG, IHBW, IINFO, IL, IMODE, INDWRK, INDX,
00330      $                   IROW, ITEMP, ITYPE, IU, IUPLO, J, J1, J2, JCOL,
00331      $                   JSIZE, JTYPE, KD, LGN, LIWEDC, LRWEDC, LWEDC,
00332      $                   M, M2, M3, MTYPES, N, NERRS, NMATS, NMAX,
00333      $                   NTEST, NTESTT
00334       DOUBLE PRECISION   ABSTOL, ANINV, ANORM, COND, OVFL, RTOVFL,
00335      $                   RTUNFL, TEMP1, TEMP2, TEMP3, ULP, ULPINV, UNFL,
00336      $                   VL, VU
00337 *     ..
00338 *     .. Local Arrays ..
00339       INTEGER            IDUMMA( 1 ), IOLDSD( 4 ), ISEED2( 4 ),
00340      $                   ISEED3( 4 ), KMAGN( MAXTYP ), KMODE( MAXTYP ),
00341      $                   KTYPE( MAXTYP )
00342 *     ..
00343 *     .. External Functions ..
00344       DOUBLE PRECISION   DLAMCH, DLARND, DSXT1
00345       EXTERNAL           DLAMCH, DLARND, DSXT1
00346 *     ..
00347 *     .. External Subroutines ..
00348       EXTERNAL           ALASVM, DLABAD, DLAFTS, XERBLA, ZHBEV, ZHBEVD,
00349      $                   ZHBEVX, ZHEEV, ZHEEVD, ZHEEVR, ZHEEVX, ZHET21,
00350      $                   ZHET22, ZHPEV, ZHPEVD, ZHPEVX, ZLACPY, ZLASET,
00351      $                   ZLATMR, ZLATMS
00352 *     ..
00353 *     .. Intrinsic Functions ..
00354       INTRINSIC          ABS, DBLE, INT, LOG, MAX, MIN, SQRT
00355 *     ..
00356 *     .. Data statements ..
00357       DATA               KTYPE / 1, 2, 5*4, 5*5, 3*8, 3*9 /
00358       DATA               KMAGN / 2*1, 1, 1, 1, 2, 3, 1, 1, 1, 2, 3, 1,
00359      $                   2, 3, 1, 2, 3 /
00360       DATA               KMODE / 2*0, 4, 3, 1, 4, 4, 4, 3, 1, 4, 4, 0,
00361      $                   0, 0, 4, 4, 4 /
00362 *     ..
00363 *     .. Executable Statements ..
00364 *
00365 *     1)      Check for errors
00366 *
00367       NTESTT = 0
00368       INFO = 0
00369 *
00370       BADNN = .FALSE.
00371       NMAX = 1
00372       DO 10 J = 1, NSIZES
00373          NMAX = MAX( NMAX, NN( J ) )
00374          IF( NN( J ).LT.0 )
00375      $      BADNN = .TRUE.
00376    10 CONTINUE
00377 *
00378 *     Check for errors
00379 *
00380       IF( NSIZES.LT.0 ) THEN
00381          INFO = -1
00382       ELSE IF( BADNN ) THEN
00383          INFO = -2
00384       ELSE IF( NTYPES.LT.0 ) THEN
00385          INFO = -3
00386       ELSE IF( LDA.LT.NMAX ) THEN
00387          INFO = -9
00388       ELSE IF( LDU.LT.NMAX ) THEN
00389          INFO = -16
00390       ELSE IF( 2*MAX( 2, NMAX )**2.GT.LWORK ) THEN
00391          INFO = -22
00392       END IF
00393 *
00394       IF( INFO.NE.0 ) THEN
00395          CALL XERBLA( 'ZDRVST', -INFO )
00396          RETURN
00397       END IF
00398 *
00399 *     Quick return if nothing to do
00400 *
00401       IF( NSIZES.EQ.0 .OR. NTYPES.EQ.0 )
00402      $   RETURN
00403 *
00404 *     More Important constants
00405 *
00406       UNFL = DLAMCH( 'Safe minimum' )
00407       OVFL = DLAMCH( 'Overflow' )
00408       CALL DLABAD( UNFL, OVFL )
00409       ULP = DLAMCH( 'Epsilon' )*DLAMCH( 'Base' )
00410       ULPINV = ONE / ULP
00411       RTUNFL = SQRT( UNFL )
00412       RTOVFL = SQRT( OVFL )
00413 *
00414 *     Loop over sizes, types
00415 *
00416       DO 20 I = 1, 4
00417          ISEED2( I ) = ISEED( I )
00418          ISEED3( I ) = ISEED( I )
00419    20 CONTINUE
00420 *
00421       NERRS = 0
00422       NMATS = 0
00423 *
00424       DO 1220 JSIZE = 1, NSIZES
00425          N = NN( JSIZE )
00426          IF( N.GT.0 ) THEN
00427             LGN = INT( LOG( DBLE( N ) ) / LOG( TWO ) )
00428             IF( 2**LGN.LT.N )
00429      $         LGN = LGN + 1
00430             IF( 2**LGN.LT.N )
00431      $         LGN = LGN + 1
00432             LWEDC = MAX( 2*N+N*N, 2*N*N )
00433             LRWEDC = 1 + 4*N + 2*N*LGN + 3*N**2
00434             LIWEDC = 3 + 5*N
00435          ELSE
00436             LWEDC = 2
00437             LRWEDC = 8
00438             LIWEDC = 8
00439          END IF
00440          ANINV = ONE / DBLE( MAX( 1, N ) )
00441 *
00442          IF( NSIZES.NE.1 ) THEN
00443             MTYPES = MIN( MAXTYP, NTYPES )
00444          ELSE
00445             MTYPES = MIN( MAXTYP+1, NTYPES )
00446          END IF
00447 *
00448          DO 1210 JTYPE = 1, MTYPES
00449             IF( .NOT.DOTYPE( JTYPE ) )
00450      $         GO TO 1210
00451             NMATS = NMATS + 1
00452             NTEST = 0
00453 *
00454             DO 30 J = 1, 4
00455                IOLDSD( J ) = ISEED( J )
00456    30       CONTINUE
00457 *
00458 *           2)      Compute "A"
00459 *
00460 *                   Control parameters:
00461 *
00462 *               KMAGN  KMODE        KTYPE
00463 *           =1  O(1)   clustered 1  zero
00464 *           =2  large  clustered 2  identity
00465 *           =3  small  exponential  (none)
00466 *           =4         arithmetic   diagonal, (w/ eigenvalues)
00467 *           =5         random log   Hermitian, w/ eigenvalues
00468 *           =6         random       (none)
00469 *           =7                      random diagonal
00470 *           =8                      random Hermitian
00471 *           =9                      band Hermitian, w/ eigenvalues
00472 *
00473             IF( MTYPES.GT.MAXTYP )
00474      $         GO TO 110
00475 *
00476             ITYPE = KTYPE( JTYPE )
00477             IMODE = KMODE( JTYPE )
00478 *
00479 *           Compute norm
00480 *
00481             GO TO ( 40, 50, 60 )KMAGN( JTYPE )
00482 *
00483    40       CONTINUE
00484             ANORM = ONE
00485             GO TO 70
00486 *
00487    50       CONTINUE
00488             ANORM = ( RTOVFL*ULP )*ANINV
00489             GO TO 70
00490 *
00491    60       CONTINUE
00492             ANORM = RTUNFL*N*ULPINV
00493             GO TO 70
00494 *
00495    70       CONTINUE
00496 *
00497             CALL ZLASET( 'Full', LDA, N, CZERO, CZERO, A, LDA )
00498             IINFO = 0
00499             COND = ULPINV
00500 *
00501 *           Special Matrices -- Identity & Jordan block
00502 *
00503 *                   Zero
00504 *
00505             IF( ITYPE.EQ.1 ) THEN
00506                IINFO = 0
00507 *
00508             ELSE IF( ITYPE.EQ.2 ) THEN
00509 *
00510 *              Identity
00511 *
00512                DO 80 JCOL = 1, N
00513                   A( JCOL, JCOL ) = ANORM
00514    80          CONTINUE
00515 *
00516             ELSE IF( ITYPE.EQ.4 ) THEN
00517 *
00518 *              Diagonal Matrix, [Eigen]values Specified
00519 *
00520                CALL ZLATMS( N, N, 'S', ISEED, 'H', RWORK, IMODE, COND,
00521      $                      ANORM, 0, 0, 'N', A, LDA, WORK, IINFO )
00522 *
00523             ELSE IF( ITYPE.EQ.5 ) THEN
00524 *
00525 *              Hermitian, eigenvalues specified
00526 *
00527                CALL ZLATMS( N, N, 'S', ISEED, 'H', RWORK, IMODE, COND,
00528      $                      ANORM, N, N, 'N', A, LDA, WORK, IINFO )
00529 *
00530             ELSE IF( ITYPE.EQ.7 ) THEN
00531 *
00532 *              Diagonal, random eigenvalues
00533 *
00534                CALL ZLATMR( N, N, 'S', ISEED, 'H', WORK, 6, ONE, CONE,
00535      $                      'T', 'N', WORK( N+1 ), 1, ONE,
00536      $                      WORK( 2*N+1 ), 1, ONE, 'N', IDUMMA, 0, 0,
00537      $                      ZERO, ANORM, 'NO', A, LDA, IWORK, IINFO )
00538 *
00539             ELSE IF( ITYPE.EQ.8 ) THEN
00540 *
00541 *              Hermitian, random eigenvalues
00542 *
00543                CALL ZLATMR( N, N, 'S', ISEED, 'H', WORK, 6, ONE, CONE,
00544      $                      'T', 'N', WORK( N+1 ), 1, ONE,
00545      $                      WORK( 2*N+1 ), 1, ONE, 'N', IDUMMA, N, N,
00546      $                      ZERO, ANORM, 'NO', A, LDA, IWORK, IINFO )
00547 *
00548             ELSE IF( ITYPE.EQ.9 ) THEN
00549 *
00550 *              Hermitian banded, eigenvalues specified
00551 *
00552                IHBW = INT( ( N-1 )*DLARND( 1, ISEED3 ) )
00553                CALL ZLATMS( N, N, 'S', ISEED, 'H', RWORK, IMODE, COND,
00554      $                      ANORM, IHBW, IHBW, 'Z', U, LDU, WORK,
00555      $                      IINFO )
00556 *
00557 *              Store as dense matrix for most routines.
00558 *
00559                CALL ZLASET( 'Full', LDA, N, CZERO, CZERO, A, LDA )
00560                DO 100 IDIAG = -IHBW, IHBW
00561                   IROW = IHBW - IDIAG + 1
00562                   J1 = MAX( 1, IDIAG+1 )
00563                   J2 = MIN( N, N+IDIAG )
00564                   DO 90 J = J1, J2
00565                      I = J - IDIAG
00566                      A( I, J ) = U( IROW, J )
00567    90             CONTINUE
00568   100          CONTINUE
00569             ELSE
00570                IINFO = 1
00571             END IF
00572 *
00573             IF( IINFO.NE.0 ) THEN
00574                WRITE( NOUNIT, FMT = 9999 )'Generator', IINFO, N, JTYPE,
00575      $            IOLDSD
00576                INFO = ABS( IINFO )
00577                RETURN
00578             END IF
00579 *
00580   110       CONTINUE
00581 *
00582             ABSTOL = UNFL + UNFL
00583             IF( N.LE.1 ) THEN
00584                IL = 1
00585                IU = N
00586             ELSE
00587                IL = 1 + INT( ( N-1 )*DLARND( 1, ISEED2 ) )
00588                IU = 1 + INT( ( N-1 )*DLARND( 1, ISEED2 ) )
00589                IF( IL.GT.IU ) THEN
00590                   ITEMP = IL
00591                   IL = IU
00592                   IU = ITEMP
00593                END IF
00594             END IF
00595 *
00596 *           Perform tests storing upper or lower triangular
00597 *           part of matrix.
00598 *
00599             DO 1200 IUPLO = 0, 1
00600                IF( IUPLO.EQ.0 ) THEN
00601                   UPLO = 'L'
00602                ELSE
00603                   UPLO = 'U'
00604                END IF
00605 *
00606 *              Call ZHEEVD and CHEEVX.
00607 *
00608                CALL ZLACPY( ' ', N, N, A, LDA, V, LDU )
00609 *
00610                NTEST = NTEST + 1
00611                CALL ZHEEVD( 'V', UPLO, N, A, LDU, D1, WORK, LWEDC,
00612      $                      RWORK, LRWEDC, IWORK, LIWEDC, IINFO )
00613                IF( IINFO.NE.0 ) THEN
00614                   WRITE( NOUNIT, FMT = 9999 )'ZHEEVD(V,' // UPLO //
00615      $               ')', IINFO, N, JTYPE, IOLDSD
00616                   INFO = ABS( IINFO )
00617                   IF( IINFO.LT.0 ) THEN
00618                      RETURN
00619                   ELSE
00620                      RESULT( NTEST ) = ULPINV
00621                      RESULT( NTEST+1 ) = ULPINV
00622                      RESULT( NTEST+2 ) = ULPINV
00623                      GO TO 130
00624                   END IF
00625                END IF
00626 *
00627 *              Do tests 1 and 2.
00628 *
00629                CALL ZHET21( 1, UPLO, N, 0, V, LDU, D1, D2, A, LDU, Z,
00630      $                      LDU, TAU, WORK, RWORK, RESULT( NTEST ) )
00631 *
00632                CALL ZLACPY( ' ', N, N, V, LDU, A, LDA )
00633 *
00634                NTEST = NTEST + 2
00635                CALL ZHEEVD( 'N', UPLO, N, A, LDU, D3, WORK, LWEDC,
00636      $                      RWORK, LRWEDC, IWORK, LIWEDC, IINFO )
00637                IF( IINFO.NE.0 ) THEN
00638                   WRITE( NOUNIT, FMT = 9999 )'ZHEEVD(N,' // UPLO //
00639      $               ')', IINFO, N, JTYPE, IOLDSD
00640                   INFO = ABS( IINFO )
00641                   IF( IINFO.LT.0 ) THEN
00642                      RETURN
00643                   ELSE
00644                      RESULT( NTEST ) = ULPINV
00645                      GO TO 130
00646                   END IF
00647                END IF
00648 *
00649 *              Do test 3.
00650 *
00651                TEMP1 = ZERO
00652                TEMP2 = ZERO
00653                DO 120 J = 1, N
00654                   TEMP1 = MAX( TEMP1, ABS( D1( J ) ), ABS( D3( J ) ) )
00655                   TEMP2 = MAX( TEMP2, ABS( D1( J )-D3( J ) ) )
00656   120          CONTINUE
00657                RESULT( NTEST ) = TEMP2 / MAX( UNFL,
00658      $                           ULP*MAX( TEMP1, TEMP2 ) )
00659 *
00660   130          CONTINUE
00661                CALL ZLACPY( ' ', N, N, V, LDU, A, LDA )
00662 *
00663                NTEST = NTEST + 1
00664 *
00665                IF( N.GT.0 ) THEN
00666                   TEMP3 = MAX( ABS( D1( 1 ) ), ABS( D1( N ) ) )
00667                   IF( IL.NE.1 ) THEN
00668                      VL = D1( IL ) - MAX( HALF*( D1( IL )-D1( IL-1 ) ),
00669      $                    TEN*ULP*TEMP3, TEN*RTUNFL )
00670                   ELSE IF( N.GT.0 ) THEN
00671                      VL = D1( 1 ) - MAX( HALF*( D1( N )-D1( 1 ) ),
00672      $                    TEN*ULP*TEMP3, TEN*RTUNFL )
00673                   END IF
00674                   IF( IU.NE.N ) THEN
00675                      VU = D1( IU ) + MAX( HALF*( D1( IU+1 )-D1( IU ) ),
00676      $                    TEN*ULP*TEMP3, TEN*RTUNFL )
00677                   ELSE IF( N.GT.0 ) THEN
00678                      VU = D1( N ) + MAX( HALF*( D1( N )-D1( 1 ) ),
00679      $                    TEN*ULP*TEMP3, TEN*RTUNFL )
00680                   END IF
00681                ELSE
00682                   TEMP3 = ZERO
00683                   VL = ZERO
00684                   VU = ONE
00685                END IF
00686 *
00687                CALL ZHEEVX( 'V', 'A', UPLO, N, A, LDU, VL, VU, IL, IU,
00688      $                      ABSTOL, M, WA1, Z, LDU, WORK, LWORK, RWORK,
00689      $                      IWORK, IWORK( 5*N+1 ), IINFO )
00690                IF( IINFO.NE.0 ) THEN
00691                   WRITE( NOUNIT, FMT = 9999 )'ZHEEVX(V,A,' // UPLO //
00692      $               ')', IINFO, N, JTYPE, IOLDSD
00693                   INFO = ABS( IINFO )
00694                   IF( IINFO.LT.0 ) THEN
00695                      RETURN
00696                   ELSE
00697                      RESULT( NTEST ) = ULPINV
00698                      RESULT( NTEST+1 ) = ULPINV
00699                      RESULT( NTEST+2 ) = ULPINV
00700                      GO TO 150
00701                   END IF
00702                END IF
00703 *
00704 *              Do tests 4 and 5.
00705 *
00706                CALL ZLACPY( ' ', N, N, V, LDU, A, LDA )
00707 *
00708                CALL ZHET21( 1, UPLO, N, 0, A, LDU, WA1, D2, Z, LDU, V,
00709      $                      LDU, TAU, WORK, RWORK, RESULT( NTEST ) )
00710 *
00711                NTEST = NTEST + 2
00712                CALL ZHEEVX( 'N', 'A', UPLO, N, A, LDU, VL, VU, IL, IU,
00713      $                      ABSTOL, M2, WA2, Z, LDU, WORK, LWORK, RWORK,
00714      $                      IWORK, IWORK( 5*N+1 ), IINFO )
00715                IF( IINFO.NE.0 ) THEN
00716                   WRITE( NOUNIT, FMT = 9999 )'ZHEEVX(N,A,' // UPLO //
00717      $               ')', IINFO, N, JTYPE, IOLDSD
00718                   INFO = ABS( IINFO )
00719                   IF( IINFO.LT.0 ) THEN
00720                      RETURN
00721                   ELSE
00722                      RESULT( NTEST ) = ULPINV
00723                      GO TO 150
00724                   END IF
00725                END IF
00726 *
00727 *              Do test 6.
00728 *
00729                TEMP1 = ZERO
00730                TEMP2 = ZERO
00731                DO 140 J = 1, N
00732                   TEMP1 = MAX( TEMP1, ABS( WA1( J ) ), ABS( WA2( J ) ) )
00733                   TEMP2 = MAX( TEMP2, ABS( WA1( J )-WA2( J ) ) )
00734   140          CONTINUE
00735                RESULT( NTEST ) = TEMP2 / MAX( UNFL,
00736      $                           ULP*MAX( TEMP1, TEMP2 ) )
00737 *
00738   150          CONTINUE
00739                CALL ZLACPY( ' ', N, N, V, LDU, A, LDA )
00740 *
00741                NTEST = NTEST + 1
00742 *
00743                CALL ZHEEVX( 'V', 'I', UPLO, N, A, LDU, VL, VU, IL, IU,
00744      $                      ABSTOL, M2, WA2, Z, LDU, WORK, LWORK, RWORK,
00745      $                      IWORK, IWORK( 5*N+1 ), IINFO )
00746                IF( IINFO.NE.0 ) THEN
00747                   WRITE( NOUNIT, FMT = 9999 )'ZHEEVX(V,I,' // UPLO //
00748      $               ')', IINFO, N, JTYPE, IOLDSD
00749                   INFO = ABS( IINFO )
00750                   IF( IINFO.LT.0 ) THEN
00751                      RETURN
00752                   ELSE
00753                      RESULT( NTEST ) = ULPINV
00754                      GO TO 160
00755                   END IF
00756                END IF
00757 *
00758 *              Do tests 7 and 8.
00759 *
00760                CALL ZLACPY( ' ', N, N, V, LDU, A, LDA )
00761 *
00762                CALL ZHET22( 1, UPLO, N, M2, 0, A, LDU, WA2, D2, Z, LDU,
00763      $                      V, LDU, TAU, WORK, RWORK, RESULT( NTEST ) )
00764 *
00765                NTEST = NTEST + 2
00766 *
00767                CALL ZHEEVX( 'N', 'I', UPLO, N, A, LDU, VL, VU, IL, IU,
00768      $                      ABSTOL, M3, WA3, Z, LDU, WORK, LWORK, RWORK,
00769      $                      IWORK, IWORK( 5*N+1 ), IINFO )
00770                IF( IINFO.NE.0 ) THEN
00771                   WRITE( NOUNIT, FMT = 9999 )'ZHEEVX(N,I,' // UPLO //
00772      $               ')', IINFO, N, JTYPE, IOLDSD
00773                   INFO = ABS( IINFO )
00774                   IF( IINFO.LT.0 ) THEN
00775                      RETURN
00776                   ELSE
00777                      RESULT( NTEST ) = ULPINV
00778                      GO TO 160
00779                   END IF
00780                END IF
00781 *
00782 *              Do test 9.
00783 *
00784                TEMP1 = DSXT1( 1, WA2, M2, WA3, M3, ABSTOL, ULP, UNFL )
00785                TEMP2 = DSXT1( 1, WA3, M3, WA2, M2, ABSTOL, ULP, UNFL )
00786                IF( N.GT.0 ) THEN
00787                   TEMP3 = MAX( ABS( WA1( 1 ) ), ABS( WA1( N ) ) )
00788                ELSE
00789                   TEMP3 = ZERO
00790                END IF
00791                RESULT( NTEST ) = ( TEMP1+TEMP2 ) /
00792      $                           MAX( UNFL, TEMP3*ULP )
00793 *
00794   160          CONTINUE
00795                CALL ZLACPY( ' ', N, N, V, LDU, A, LDA )
00796 *
00797                NTEST = NTEST + 1
00798 *
00799                CALL ZHEEVX( 'V', 'V', UPLO, N, A, LDU, VL, VU, IL, IU,
00800      $                      ABSTOL, M2, WA2, Z, LDU, WORK, LWORK, RWORK,
00801      $                      IWORK, IWORK( 5*N+1 ), IINFO )
00802                IF( IINFO.NE.0 ) THEN
00803                   WRITE( NOUNIT, FMT = 9999 )'ZHEEVX(V,V,' // UPLO //
00804      $               ')', IINFO, N, JTYPE, IOLDSD
00805                   INFO = ABS( IINFO )
00806                   IF( IINFO.LT.0 ) THEN
00807                      RETURN
00808                   ELSE
00809                      RESULT( NTEST ) = ULPINV
00810                      GO TO 170
00811                   END IF
00812                END IF
00813 *
00814 *              Do tests 10 and 11.
00815 *
00816                CALL ZLACPY( ' ', N, N, V, LDU, A, LDA )
00817 *
00818                CALL ZHET22( 1, UPLO, N, M2, 0, A, LDU, WA2, D2, Z, LDU,
00819      $                      V, LDU, TAU, WORK, RWORK, RESULT( NTEST ) )
00820 *
00821                NTEST = NTEST + 2
00822 *
00823                CALL ZHEEVX( 'N', 'V', UPLO, N, A, LDU, VL, VU, IL, IU,
00824      $                      ABSTOL, M3, WA3, Z, LDU, WORK, LWORK, RWORK,
00825      $                      IWORK, IWORK( 5*N+1 ), IINFO )
00826                IF( IINFO.NE.0 ) THEN
00827                   WRITE( NOUNIT, FMT = 9999 )'ZHEEVX(N,V,' // UPLO //
00828      $               ')', IINFO, N, JTYPE, IOLDSD
00829                   INFO = ABS( IINFO )
00830                   IF( IINFO.LT.0 ) THEN
00831                      RETURN
00832                   ELSE
00833                      RESULT( NTEST ) = ULPINV
00834                      GO TO 170
00835                   END IF
00836                END IF
00837 *
00838                IF( M3.EQ.0 .AND. N.GT.0 ) THEN
00839                   RESULT( NTEST ) = ULPINV
00840                   GO TO 170
00841                END IF
00842 *
00843 *              Do test 12.
00844 *
00845                TEMP1 = DSXT1( 1, WA2, M2, WA3, M3, ABSTOL, ULP, UNFL )
00846                TEMP2 = DSXT1( 1, WA3, M3, WA2, M2, ABSTOL, ULP, UNFL )
00847                IF( N.GT.0 ) THEN
00848                   TEMP3 = MAX( ABS( WA1( 1 ) ), ABS( WA1( N ) ) )
00849                ELSE
00850                   TEMP3 = ZERO
00851                END IF
00852                RESULT( NTEST ) = ( TEMP1+TEMP2 ) /
00853      $                           MAX( UNFL, TEMP3*ULP )
00854 *
00855   170          CONTINUE
00856 *
00857 *              Call ZHPEVD and CHPEVX.
00858 *
00859                CALL ZLACPY( ' ', N, N, V, LDU, A, LDA )
00860 *
00861 *              Load array WORK with the upper or lower triangular
00862 *              part of the matrix in packed form.
00863 *
00864                IF( IUPLO.EQ.1 ) THEN
00865                   INDX = 1
00866                   DO 190 J = 1, N
00867                      DO 180 I = 1, J
00868                         WORK( INDX ) = A( I, J )
00869                         INDX = INDX + 1
00870   180                CONTINUE
00871   190             CONTINUE
00872                ELSE
00873                   INDX = 1
00874                   DO 210 J = 1, N
00875                      DO 200 I = J, N
00876                         WORK( INDX ) = A( I, J )
00877                         INDX = INDX + 1
00878   200                CONTINUE
00879   210             CONTINUE
00880                END IF
00881 *
00882                NTEST = NTEST + 1
00883                INDWRK = N*( N+1 ) / 2 + 1
00884                CALL ZHPEVD( 'V', UPLO, N, WORK, D1, Z, LDU,
00885      $                      WORK( INDWRK ), LWEDC, RWORK, LRWEDC, IWORK,
00886      $                      LIWEDC, IINFO )
00887                IF( IINFO.NE.0 ) THEN
00888                   WRITE( NOUNIT, FMT = 9999 )'ZHPEVD(V,' // UPLO //
00889      $               ')', IINFO, N, JTYPE, IOLDSD
00890                   INFO = ABS( IINFO )
00891                   IF( IINFO.LT.0 ) THEN
00892                      RETURN
00893                   ELSE
00894                      RESULT( NTEST ) = ULPINV
00895                      RESULT( NTEST+1 ) = ULPINV
00896                      RESULT( NTEST+2 ) = ULPINV
00897                      GO TO 270
00898                   END IF
00899                END IF
00900 *
00901 *              Do tests 13 and 14.
00902 *
00903                CALL ZHET21( 1, UPLO, N, 0, A, LDA, D1, D2, Z, LDU, V,
00904      $                      LDU, TAU, WORK, RWORK, RESULT( NTEST ) )
00905 *
00906                IF( IUPLO.EQ.1 ) THEN
00907                   INDX = 1
00908                   DO 230 J = 1, N
00909                      DO 220 I = 1, J
00910                         WORK( INDX ) = A( I, J )
00911                         INDX = INDX + 1
00912   220                CONTINUE
00913   230             CONTINUE
00914                ELSE
00915                   INDX = 1
00916                   DO 250 J = 1, N
00917                      DO 240 I = J, N
00918                         WORK( INDX ) = A( I, J )
00919                         INDX = INDX + 1
00920   240                CONTINUE
00921   250             CONTINUE
00922                END IF
00923 *
00924                NTEST = NTEST + 2
00925                INDWRK = N*( N+1 ) / 2 + 1
00926                CALL ZHPEVD( 'N', UPLO, N, WORK, D3, Z, LDU,
00927      $                      WORK( INDWRK ), LWEDC, RWORK, LRWEDC, IWORK,
00928      $                      LIWEDC, IINFO )
00929                IF( IINFO.NE.0 ) THEN
00930                   WRITE( NOUNIT, FMT = 9999 )'ZHPEVD(N,' // UPLO //
00931      $               ')', IINFO, N, JTYPE, IOLDSD
00932                   INFO = ABS( IINFO )
00933                   IF( IINFO.LT.0 ) THEN
00934                      RETURN
00935                   ELSE
00936                      RESULT( NTEST ) = ULPINV
00937                      GO TO 270
00938                   END IF
00939                END IF
00940 *
00941 *              Do test 15.
00942 *
00943                TEMP1 = ZERO
00944                TEMP2 = ZERO
00945                DO 260 J = 1, N
00946                   TEMP1 = MAX( TEMP1, ABS( D1( J ) ), ABS( D3( J ) ) )
00947                   TEMP2 = MAX( TEMP2, ABS( D1( J )-D3( J ) ) )
00948   260          CONTINUE
00949                RESULT( NTEST ) = TEMP2 / MAX( UNFL,
00950      $                           ULP*MAX( TEMP1, TEMP2 ) )
00951 *
00952 *              Load array WORK with the upper or lower triangular part
00953 *              of the matrix in packed form.
00954 *
00955   270          CONTINUE
00956                IF( IUPLO.EQ.1 ) THEN
00957                   INDX = 1
00958                   DO 290 J = 1, N
00959                      DO 280 I = 1, J
00960                         WORK( INDX ) = A( I, J )
00961                         INDX = INDX + 1
00962   280                CONTINUE
00963   290             CONTINUE
00964                ELSE
00965                   INDX = 1
00966                   DO 310 J = 1, N
00967                      DO 300 I = J, N
00968                         WORK( INDX ) = A( I, J )
00969                         INDX = INDX + 1
00970   300                CONTINUE
00971   310             CONTINUE
00972                END IF
00973 *
00974                NTEST = NTEST + 1
00975 *
00976                IF( N.GT.0 ) THEN
00977                   TEMP3 = MAX( ABS( D1( 1 ) ), ABS( D1( N ) ) )
00978                   IF( IL.NE.1 ) THEN
00979                      VL = D1( IL ) - MAX( HALF*( D1( IL )-D1( IL-1 ) ),
00980      $                    TEN*ULP*TEMP3, TEN*RTUNFL )
00981                   ELSE IF( N.GT.0 ) THEN
00982                      VL = D1( 1 ) - MAX( HALF*( D1( N )-D1( 1 ) ),
00983      $                    TEN*ULP*TEMP3, TEN*RTUNFL )
00984                   END IF
00985                   IF( IU.NE.N ) THEN
00986                      VU = D1( IU ) + MAX( HALF*( D1( IU+1 )-D1( IU ) ),
00987      $                    TEN*ULP*TEMP3, TEN*RTUNFL )
00988                   ELSE IF( N.GT.0 ) THEN
00989                      VU = D1( N ) + MAX( HALF*( D1( N )-D1( 1 ) ),
00990      $                    TEN*ULP*TEMP3, TEN*RTUNFL )
00991                   END IF
00992                ELSE
00993                   TEMP3 = ZERO
00994                   VL = ZERO
00995                   VU = ONE
00996                END IF
00997 *
00998                CALL ZHPEVX( 'V', 'A', UPLO, N, WORK, VL, VU, IL, IU,
00999      $                      ABSTOL, M, WA1, Z, LDU, V, RWORK, IWORK,
01000      $                      IWORK( 5*N+1 ), IINFO )
01001                IF( IINFO.NE.0 ) THEN
01002                   WRITE( NOUNIT, FMT = 9999 )'ZHPEVX(V,A,' // UPLO //
01003      $               ')', IINFO, N, JTYPE, IOLDSD
01004                   INFO = ABS( IINFO )
01005                   IF( IINFO.LT.0 ) THEN
01006                      RETURN
01007                   ELSE
01008                      RESULT( NTEST ) = ULPINV
01009                      RESULT( NTEST+1 ) = ULPINV
01010                      RESULT( NTEST+2 ) = ULPINV
01011                      GO TO 370
01012                   END IF
01013                END IF
01014 *
01015 *              Do tests 16 and 17.
01016 *
01017                CALL ZHET21( 1, UPLO, N, 0, A, LDU, WA1, D2, Z, LDU, V,
01018      $                      LDU, TAU, WORK, RWORK, RESULT( NTEST ) )
01019 *
01020                NTEST = NTEST + 2
01021 *
01022                IF( IUPLO.EQ.1 ) THEN
01023                   INDX = 1
01024                   DO 330 J = 1, N
01025                      DO 320 I = 1, J
01026                         WORK( INDX ) = A( I, J )
01027                         INDX = INDX + 1
01028   320                CONTINUE
01029   330             CONTINUE
01030                ELSE
01031                   INDX = 1
01032                   DO 350 J = 1, N
01033                      DO 340 I = J, N
01034                         WORK( INDX ) = A( I, J )
01035                         INDX = INDX + 1
01036   340                CONTINUE
01037   350             CONTINUE
01038                END IF
01039 *
01040                CALL ZHPEVX( 'N', 'A', UPLO, N, WORK, VL, VU, IL, IU,
01041      $                      ABSTOL, M2, WA2, Z, LDU, V, RWORK, IWORK,
01042      $                      IWORK( 5*N+1 ), IINFO )
01043                IF( IINFO.NE.0 ) THEN
01044                   WRITE( NOUNIT, FMT = 9999 )'ZHPEVX(N,A,' // UPLO //
01045      $               ')', IINFO, N, JTYPE, IOLDSD
01046                   INFO = ABS( IINFO )
01047                   IF( IINFO.LT.0 ) THEN
01048                      RETURN
01049                   ELSE
01050                      RESULT( NTEST ) = ULPINV
01051                      GO TO 370
01052                   END IF
01053                END IF
01054 *
01055 *              Do test 18.
01056 *
01057                TEMP1 = ZERO
01058                TEMP2 = ZERO
01059                DO 360 J = 1, N
01060                   TEMP1 = MAX( TEMP1, ABS( WA1( J ) ), ABS( WA2( J ) ) )
01061                   TEMP2 = MAX( TEMP2, ABS( WA1( J )-WA2( J ) ) )
01062   360          CONTINUE
01063                RESULT( NTEST ) = TEMP2 / MAX( UNFL,
01064      $                           ULP*MAX( TEMP1, TEMP2 ) )
01065 *
01066   370          CONTINUE
01067                NTEST = NTEST + 1
01068                IF( IUPLO.EQ.1 ) THEN
01069                   INDX = 1
01070                   DO 390 J = 1, N
01071                      DO 380 I = 1, J
01072                         WORK( INDX ) = A( I, J )
01073                         INDX = INDX + 1
01074   380                CONTINUE
01075   390             CONTINUE
01076                ELSE
01077                   INDX = 1
01078                   DO 410 J = 1, N
01079                      DO 400 I = J, N
01080                         WORK( INDX ) = A( I, J )
01081                         INDX = INDX + 1
01082   400                CONTINUE
01083   410             CONTINUE
01084                END IF
01085 *
01086                CALL ZHPEVX( 'V', 'I', UPLO, N, WORK, VL, VU, IL, IU,
01087      $                      ABSTOL, M2, WA2, Z, LDU, V, RWORK, IWORK,
01088      $                      IWORK( 5*N+1 ), IINFO )
01089                IF( IINFO.NE.0 ) THEN
01090                   WRITE( NOUNIT, FMT = 9999 )'ZHPEVX(V,I,' // UPLO //
01091      $               ')', IINFO, N, JTYPE, IOLDSD
01092                   INFO = ABS( IINFO )
01093                   IF( IINFO.LT.0 ) THEN
01094                      RETURN
01095                   ELSE
01096                      RESULT( NTEST ) = ULPINV
01097                      RESULT( NTEST+1 ) = ULPINV
01098                      RESULT( NTEST+2 ) = ULPINV
01099                      GO TO 460
01100                   END IF
01101                END IF
01102 *
01103 *              Do tests 19 and 20.
01104 *
01105                CALL ZHET22( 1, UPLO, N, M2, 0, A, LDU, WA2, D2, Z, LDU,
01106      $                      V, LDU, TAU, WORK, RWORK, RESULT( NTEST ) )
01107 *
01108                NTEST = NTEST + 2
01109 *
01110                IF( IUPLO.EQ.1 ) THEN
01111                   INDX = 1
01112                   DO 430 J = 1, N
01113                      DO 420 I = 1, J
01114                         WORK( INDX ) = A( I, J )
01115                         INDX = INDX + 1
01116   420                CONTINUE
01117   430             CONTINUE
01118                ELSE
01119                   INDX = 1
01120                   DO 450 J = 1, N
01121                      DO 440 I = J, N
01122                         WORK( INDX ) = A( I, J )
01123                         INDX = INDX + 1
01124   440                CONTINUE
01125   450             CONTINUE
01126                END IF
01127 *
01128                CALL ZHPEVX( 'N', 'I', UPLO, N, WORK, VL, VU, IL, IU,
01129      $                      ABSTOL, M3, WA3, Z, LDU, V, RWORK, IWORK,
01130      $                      IWORK( 5*N+1 ), IINFO )
01131                IF( IINFO.NE.0 ) THEN
01132                   WRITE( NOUNIT, FMT = 9999 )'ZHPEVX(N,I,' // UPLO //
01133      $               ')', IINFO, N, JTYPE, IOLDSD
01134                   INFO = ABS( IINFO )
01135                   IF( IINFO.LT.0 ) THEN
01136                      RETURN
01137                   ELSE
01138                      RESULT( NTEST ) = ULPINV
01139                      GO TO 460
01140                   END IF
01141                END IF
01142 *
01143 *              Do test 21.
01144 *
01145                TEMP1 = DSXT1( 1, WA2, M2, WA3, M3, ABSTOL, ULP, UNFL )
01146                TEMP2 = DSXT1( 1, WA3, M3, WA2, M2, ABSTOL, ULP, UNFL )
01147                IF( N.GT.0 ) THEN
01148                   TEMP3 = MAX( ABS( WA1( 1 ) ), ABS( WA1( N ) ) )
01149                ELSE
01150                   TEMP3 = ZERO
01151                END IF
01152                RESULT( NTEST ) = ( TEMP1+TEMP2 ) /
01153      $                           MAX( UNFL, TEMP3*ULP )
01154 *
01155   460          CONTINUE
01156                NTEST = NTEST + 1
01157                IF( IUPLO.EQ.1 ) THEN
01158                   INDX = 1
01159                   DO 480 J = 1, N
01160                      DO 470 I = 1, J
01161                         WORK( INDX ) = A( I, J )
01162                         INDX = INDX + 1
01163   470                CONTINUE
01164   480             CONTINUE
01165                ELSE
01166                   INDX = 1
01167                   DO 500 J = 1, N
01168                      DO 490 I = J, N
01169                         WORK( INDX ) = A( I, J )
01170                         INDX = INDX + 1
01171   490                CONTINUE
01172   500             CONTINUE
01173                END IF
01174 *
01175                CALL ZHPEVX( 'V', 'V', UPLO, N, WORK, VL, VU, IL, IU,
01176      $                      ABSTOL, M2, WA2, Z, LDU, V, RWORK, IWORK,
01177      $                      IWORK( 5*N+1 ), IINFO )
01178                IF( IINFO.NE.0 ) THEN
01179                   WRITE( NOUNIT, FMT = 9999 )'ZHPEVX(V,V,' // UPLO //
01180      $               ')', IINFO, N, JTYPE, IOLDSD
01181                   INFO = ABS( IINFO )
01182                   IF( IINFO.LT.0 ) THEN
01183                      RETURN
01184                   ELSE
01185                      RESULT( NTEST ) = ULPINV
01186                      RESULT( NTEST+1 ) = ULPINV
01187                      RESULT( NTEST+2 ) = ULPINV
01188                      GO TO 550
01189                   END IF
01190                END IF
01191 *
01192 *              Do tests 22 and 23.
01193 *
01194                CALL ZHET22( 1, UPLO, N, M2, 0, A, LDU, WA2, D2, Z, LDU,
01195      $                      V, LDU, TAU, WORK, RWORK, RESULT( NTEST ) )
01196 *
01197                NTEST = NTEST + 2
01198 *
01199                IF( IUPLO.EQ.1 ) THEN
01200                   INDX = 1
01201                   DO 520 J = 1, N
01202                      DO 510 I = 1, J
01203                         WORK( INDX ) = A( I, J )
01204                         INDX = INDX + 1
01205   510                CONTINUE
01206   520             CONTINUE
01207                ELSE
01208                   INDX = 1
01209                   DO 540 J = 1, N
01210                      DO 530 I = J, N
01211                         WORK( INDX ) = A( I, J )
01212                         INDX = INDX + 1
01213   530                CONTINUE
01214   540             CONTINUE
01215                END IF
01216 *
01217                CALL ZHPEVX( 'N', 'V', UPLO, N, WORK, VL, VU, IL, IU,
01218      $                      ABSTOL, M3, WA3, Z, LDU, V, RWORK, IWORK,
01219      $                      IWORK( 5*N+1 ), IINFO )
01220                IF( IINFO.NE.0 ) THEN
01221                   WRITE( NOUNIT, FMT = 9999 )'ZHPEVX(N,V,' // UPLO //
01222      $               ')', IINFO, N, JTYPE, IOLDSD
01223                   INFO = ABS( IINFO )
01224                   IF( IINFO.LT.0 ) THEN
01225                      RETURN
01226                   ELSE
01227                      RESULT( NTEST ) = ULPINV
01228                      GO TO 550
01229                   END IF
01230                END IF
01231 *
01232                IF( M3.EQ.0 .AND. N.GT.0 ) THEN
01233                   RESULT( NTEST ) = ULPINV
01234                   GO TO 550
01235                END IF
01236 *
01237 *              Do test 24.
01238 *
01239                TEMP1 = DSXT1( 1, WA2, M2, WA3, M3, ABSTOL, ULP, UNFL )
01240                TEMP2 = DSXT1( 1, WA3, M3, WA2, M2, ABSTOL, ULP, UNFL )
01241                IF( N.GT.0 ) THEN
01242                   TEMP3 = MAX( ABS( WA1( 1 ) ), ABS( WA1( N ) ) )
01243                ELSE
01244                   TEMP3 = ZERO
01245                END IF
01246                RESULT( NTEST ) = ( TEMP1+TEMP2 ) /
01247      $                           MAX( UNFL, TEMP3*ULP )
01248 *
01249   550          CONTINUE
01250 *
01251 *              Call ZHBEVD and CHBEVX.
01252 *
01253                IF( JTYPE.LE.7 ) THEN
01254                   KD = 0
01255                ELSE IF( JTYPE.GE.8 .AND. JTYPE.LE.15 ) THEN
01256                   KD = MAX( N-1, 0 )
01257                ELSE
01258                   KD = IHBW
01259                END IF
01260 *
01261 *              Load array V with the upper or lower triangular part
01262 *              of the matrix in band form.
01263 *
01264                IF( IUPLO.EQ.1 ) THEN
01265                   DO 570 J = 1, N
01266                      DO 560 I = MAX( 1, J-KD ), J
01267                         V( KD+1+I-J, J ) = A( I, J )
01268   560                CONTINUE
01269   570             CONTINUE
01270                ELSE
01271                   DO 590 J = 1, N
01272                      DO 580 I = J, MIN( N, J+KD )
01273                         V( 1+I-J, J ) = A( I, J )
01274   580                CONTINUE
01275   590             CONTINUE
01276                END IF
01277 *
01278                NTEST = NTEST + 1
01279                CALL ZHBEVD( 'V', UPLO, N, KD, V, LDU, D1, Z, LDU, WORK,
01280      $                      LWEDC, RWORK, LRWEDC, IWORK, LIWEDC, IINFO )
01281                IF( IINFO.NE.0 ) THEN
01282                   WRITE( NOUNIT, FMT = 9998 )'ZHBEVD(V,' // UPLO //
01283      $               ')', IINFO, N, KD, JTYPE, IOLDSD
01284                   INFO = ABS( IINFO )
01285                   IF( IINFO.LT.0 ) THEN
01286                      RETURN
01287                   ELSE
01288                      RESULT( NTEST ) = ULPINV
01289                      RESULT( NTEST+1 ) = ULPINV
01290                      RESULT( NTEST+2 ) = ULPINV
01291                      GO TO 650
01292                   END IF
01293                END IF
01294 *
01295 *              Do tests 25 and 26.
01296 *
01297                CALL ZHET21( 1, UPLO, N, 0, A, LDA, D1, D2, Z, LDU, V,
01298      $                      LDU, TAU, WORK, RWORK, RESULT( NTEST ) )
01299 *
01300                IF( IUPLO.EQ.1 ) THEN
01301                   DO 610 J = 1, N
01302                      DO 600 I = MAX( 1, J-KD ), J
01303                         V( KD+1+I-J, J ) = A( I, J )
01304   600                CONTINUE
01305   610             CONTINUE
01306                ELSE
01307                   DO 630 J = 1, N
01308                      DO 620 I = J, MIN( N, J+KD )
01309                         V( 1+I-J, J ) = A( I, J )
01310   620                CONTINUE
01311   630             CONTINUE
01312                END IF
01313 *
01314                NTEST = NTEST + 2
01315                CALL ZHBEVD( 'N', UPLO, N, KD, V, LDU, D3, Z, LDU, WORK,
01316      $                      LWEDC, RWORK, LRWEDC, IWORK, LIWEDC, IINFO )
01317                IF( IINFO.NE.0 ) THEN
01318                   WRITE( NOUNIT, FMT = 9998 )'ZHBEVD(N,' // UPLO //
01319      $               ')', IINFO, N, KD, JTYPE, IOLDSD
01320                   INFO = ABS( IINFO )
01321                   IF( IINFO.LT.0 ) THEN
01322                      RETURN
01323                   ELSE
01324                      RESULT( NTEST ) = ULPINV
01325                      GO TO 650
01326                   END IF
01327                END IF
01328 *
01329 *              Do test 27.
01330 *
01331                TEMP1 = ZERO
01332                TEMP2 = ZERO
01333                DO 640 J = 1, N
01334                   TEMP1 = MAX( TEMP1, ABS( D1( J ) ), ABS( D3( J ) ) )
01335                   TEMP2 = MAX( TEMP2, ABS( D1( J )-D3( J ) ) )
01336   640          CONTINUE
01337                RESULT( NTEST ) = TEMP2 / MAX( UNFL,
01338      $                           ULP*MAX( TEMP1, TEMP2 ) )
01339 *
01340 *              Load array V with the upper or lower triangular part
01341 *              of the matrix in band form.
01342 *
01343   650          CONTINUE
01344                IF( IUPLO.EQ.1 ) THEN
01345                   DO 670 J = 1, N
01346                      DO 660 I = MAX( 1, J-KD ), J
01347                         V( KD+1+I-J, J ) = A( I, J )
01348   660                CONTINUE
01349   670             CONTINUE
01350                ELSE
01351                   DO 690 J = 1, N
01352                      DO 680 I = J, MIN( N, J+KD )
01353                         V( 1+I-J, J ) = A( I, J )
01354   680                CONTINUE
01355   690             CONTINUE
01356                END IF
01357 *
01358                NTEST = NTEST + 1
01359                CALL ZHBEVX( 'V', 'A', UPLO, N, KD, V, LDU, U, LDU, VL,
01360      $                      VU, IL, IU, ABSTOL, M, WA1, Z, LDU, WORK,
01361      $                      RWORK, IWORK, IWORK( 5*N+1 ), IINFO )
01362                IF( IINFO.NE.0 ) THEN
01363                   WRITE( NOUNIT, FMT = 9999 )'ZHBEVX(V,A,' // UPLO //
01364      $               ')', IINFO, N, KD, JTYPE, IOLDSD
01365                   INFO = ABS( IINFO )
01366                   IF( IINFO.LT.0 ) THEN
01367                      RETURN
01368                   ELSE
01369                      RESULT( NTEST ) = ULPINV
01370                      RESULT( NTEST+1 ) = ULPINV
01371                      RESULT( NTEST+2 ) = ULPINV
01372                      GO TO 750
01373                   END IF
01374                END IF
01375 *
01376 *              Do tests 28 and 29.
01377 *
01378                CALL ZHET21( 1, UPLO, N, 0, A, LDU, WA1, D2, Z, LDU, V,
01379      $                      LDU, TAU, WORK, RWORK, RESULT( NTEST ) )
01380 *
01381                NTEST = NTEST + 2
01382 *
01383                IF( IUPLO.EQ.1 ) THEN
01384                   DO 710 J = 1, N
01385                      DO 700 I = MAX( 1, J-KD ), J
01386                         V( KD+1+I-J, J ) = A( I, J )
01387   700                CONTINUE
01388   710             CONTINUE
01389                ELSE
01390                   DO 730 J = 1, N
01391                      DO 720 I = J, MIN( N, J+KD )
01392                         V( 1+I-J, J ) = A( I, J )
01393   720                CONTINUE
01394   730             CONTINUE
01395                END IF
01396 *
01397                CALL ZHBEVX( 'N', 'A', UPLO, N, KD, V, LDU, U, LDU, VL,
01398      $                      VU, IL, IU, ABSTOL, M2, WA2, Z, LDU, WORK,
01399      $                      RWORK, IWORK, IWORK( 5*N+1 ), IINFO )
01400                IF( IINFO.NE.0 ) THEN
01401                   WRITE( NOUNIT, FMT = 9998 )'ZHBEVX(N,A,' // UPLO //
01402      $               ')', IINFO, N, KD, JTYPE, IOLDSD
01403                   INFO = ABS( IINFO )
01404                   IF( IINFO.LT.0 ) THEN
01405                      RETURN
01406                   ELSE
01407                      RESULT( NTEST ) = ULPINV
01408                      GO TO 750
01409                   END IF
01410                END IF
01411 *
01412 *              Do test 30.
01413 *
01414                TEMP1 = ZERO
01415                TEMP2 = ZERO
01416                DO 740 J = 1, N
01417                   TEMP1 = MAX( TEMP1, ABS( WA1( J ) ), ABS( WA2( J ) ) )
01418                   TEMP2 = MAX( TEMP2, ABS( WA1( J )-WA2( J ) ) )
01419   740          CONTINUE
01420                RESULT( NTEST ) = TEMP2 / MAX( UNFL,
01421      $                           ULP*MAX( TEMP1, TEMP2 ) )
01422 *
01423 *              Load array V with the upper or lower triangular part
01424 *              of the matrix in band form.
01425 *
01426   750          CONTINUE
01427                NTEST = NTEST + 1
01428                IF( IUPLO.EQ.1 ) THEN
01429                   DO 770 J = 1, N
01430                      DO 760 I = MAX( 1, J-KD ), J
01431                         V( KD+1+I-J, J ) = A( I, J )
01432   760                CONTINUE
01433   770             CONTINUE
01434                ELSE
01435                   DO 790 J = 1, N
01436                      DO 780 I = J, MIN( N, J+KD )
01437                         V( 1+I-J, J ) = A( I, J )
01438   780                CONTINUE
01439   790             CONTINUE
01440                END IF
01441 *
01442                CALL ZHBEVX( 'V', 'I', UPLO, N, KD, V, LDU, U, LDU, VL,
01443      $                      VU, IL, IU, ABSTOL, M2, WA2, Z, LDU, WORK,
01444      $                      RWORK, IWORK, IWORK( 5*N+1 ), IINFO )
01445                IF( IINFO.NE.0 ) THEN
01446                   WRITE( NOUNIT, FMT = 9998 )'ZHBEVX(V,I,' // UPLO //
01447      $               ')', IINFO, N, KD, JTYPE, IOLDSD
01448                   INFO = ABS( IINFO )
01449                   IF( IINFO.LT.0 ) THEN
01450                      RETURN
01451                   ELSE
01452                      RESULT( NTEST ) = ULPINV
01453                      RESULT( NTEST+1 ) = ULPINV
01454                      RESULT( NTEST+2 ) = ULPINV
01455                      GO TO 840
01456                   END IF
01457                END IF
01458 *
01459 *              Do tests 31 and 32.
01460 *
01461                CALL ZHET22( 1, UPLO, N, M2, 0, A, LDU, WA2, D2, Z, LDU,
01462      $                      V, LDU, TAU, WORK, RWORK, RESULT( NTEST ) )
01463 *
01464                NTEST = NTEST + 2
01465 *
01466                IF( IUPLO.EQ.1 ) THEN
01467                   DO 810 J = 1, N
01468                      DO 800 I = MAX( 1, J-KD ), J
01469                         V( KD+1+I-J, J ) = A( I, J )
01470   800                CONTINUE
01471   810             CONTINUE
01472                ELSE
01473                   DO 830 J = 1, N
01474                      DO 820 I = J, MIN( N, J+KD )
01475                         V( 1+I-J, J ) = A( I, J )
01476   820                CONTINUE
01477   830             CONTINUE
01478                END IF
01479                CALL ZHBEVX( 'N', 'I', UPLO, N, KD, V, LDU, U, LDU, VL,
01480      $                      VU, IL, IU, ABSTOL, M3, WA3, Z, LDU, WORK,
01481      $                      RWORK, IWORK, IWORK( 5*N+1 ), IINFO )
01482                IF( IINFO.NE.0 ) THEN
01483                   WRITE( NOUNIT, FMT = 9998 )'ZHBEVX(N,I,' // UPLO //
01484      $               ')', IINFO, N, KD, JTYPE, IOLDSD
01485                   INFO = ABS( IINFO )
01486                   IF( IINFO.LT.0 ) THEN
01487                      RETURN
01488                   ELSE
01489                      RESULT( NTEST ) = ULPINV
01490                      GO TO 840
01491                   END IF
01492                END IF
01493 *
01494 *              Do test 33.
01495 *
01496                TEMP1 = DSXT1( 1, WA2, M2, WA3, M3, ABSTOL, ULP, UNFL )
01497                TEMP2 = DSXT1( 1, WA3, M3, WA2, M2, ABSTOL, ULP, UNFL )
01498                IF( N.GT.0 ) THEN
01499                   TEMP3 = MAX( ABS( WA1( 1 ) ), ABS( WA1( N ) ) )
01500                ELSE
01501                   TEMP3 = ZERO
01502                END IF
01503                RESULT( NTEST ) = ( TEMP1+TEMP2 ) /
01504      $                           MAX( UNFL, TEMP3*ULP )
01505 *
01506 *              Load array V with the upper or lower triangular part
01507 *              of the matrix in band form.
01508 *
01509   840          CONTINUE
01510                NTEST = NTEST + 1
01511                IF( IUPLO.EQ.1 ) THEN
01512                   DO 860 J = 1, N
01513                      DO 850 I = MAX( 1, J-KD ), J
01514                         V( KD+1+I-J, J ) = A( I, J )
01515   850                CONTINUE
01516   860             CONTINUE
01517                ELSE
01518                   DO 880 J = 1, N
01519                      DO 870 I = J, MIN( N, J+KD )
01520                         V( 1+I-J, J ) = A( I, J )
01521   870                CONTINUE
01522   880             CONTINUE
01523                END IF
01524                CALL ZHBEVX( 'V', 'V', UPLO, N, KD, V, LDU, U, LDU, VL,
01525      $                      VU, IL, IU, ABSTOL, M2, WA2, Z, LDU, WORK,
01526      $                      RWORK, IWORK, IWORK( 5*N+1 ), IINFO )
01527                IF( IINFO.NE.0 ) THEN
01528                   WRITE( NOUNIT, FMT = 9998 )'ZHBEVX(V,V,' // UPLO //
01529      $               ')', IINFO, N, KD, JTYPE, IOLDSD
01530                   INFO = ABS( IINFO )
01531                   IF( IINFO.LT.0 ) THEN
01532                      RETURN
01533                   ELSE
01534                      RESULT( NTEST ) = ULPINV
01535                      RESULT( NTEST+1 ) = ULPINV
01536                      RESULT( NTEST+2 ) = ULPINV
01537                      GO TO 930
01538                   END IF
01539                END IF
01540 *
01541 *              Do tests 34 and 35.
01542 *
01543                CALL ZHET22( 1, UPLO, N, M2, 0, A, LDU, WA2, D2, Z, LDU,
01544      $                      V, LDU, TAU, WORK, RWORK, RESULT( NTEST ) )
01545 *
01546                NTEST = NTEST + 2
01547 *
01548                IF( IUPLO.EQ.1 ) THEN
01549                   DO 900 J = 1, N
01550                      DO 890 I = MAX( 1, J-KD ), J
01551                         V( KD+1+I-J, J ) = A( I, J )
01552   890                CONTINUE
01553   900             CONTINUE
01554                ELSE
01555                   DO 920 J = 1, N
01556                      DO 910 I = J, MIN( N, J+KD )
01557                         V( 1+I-J, J ) = A( I, J )
01558   910                CONTINUE
01559   920             CONTINUE
01560                END IF
01561                CALL ZHBEVX( 'N', 'V', UPLO, N, KD, V, LDU, U, LDU, VL,
01562      $                      VU, IL, IU, ABSTOL, M3, WA3, Z, LDU, WORK,
01563      $                      RWORK, IWORK, IWORK( 5*N+1 ), IINFO )
01564                IF( IINFO.NE.0 ) THEN
01565                   WRITE( NOUNIT, FMT = 9998 )'ZHBEVX(N,V,' // UPLO //
01566      $               ')', IINFO, N, KD, JTYPE, IOLDSD
01567                   INFO = ABS( IINFO )
01568                   IF( IINFO.LT.0 ) THEN
01569                      RETURN
01570                   ELSE
01571                      RESULT( NTEST ) = ULPINV
01572                      GO TO 930
01573                   END IF
01574                END IF
01575 *
01576                IF( M3.EQ.0 .AND. N.GT.0 ) THEN
01577                   RESULT( NTEST ) = ULPINV
01578                   GO TO 930
01579                END IF
01580 *
01581 *              Do test 36.
01582 *
01583                TEMP1 = DSXT1( 1, WA2, M2, WA3, M3, ABSTOL, ULP, UNFL )
01584                TEMP2 = DSXT1( 1, WA3, M3, WA2, M2, ABSTOL, ULP, UNFL )
01585                IF( N.GT.0 ) THEN
01586                   TEMP3 = MAX( ABS( WA1( 1 ) ), ABS( WA1( N ) ) )
01587                ELSE
01588                   TEMP3 = ZERO
01589                END IF
01590                RESULT( NTEST ) = ( TEMP1+TEMP2 ) /
01591      $                           MAX( UNFL, TEMP3*ULP )
01592 *
01593   930          CONTINUE
01594 *
01595 *              Call ZHEEV
01596 *
01597                CALL ZLACPY( ' ', N, N, A, LDA, V, LDU )
01598 *
01599                NTEST = NTEST + 1
01600                CALL ZHEEV( 'V', UPLO, N, A, LDU, D1, WORK, LWORK, RWORK,
01601      $                     IINFO )
01602                IF( IINFO.NE.0 ) THEN
01603                   WRITE( NOUNIT, FMT = 9999 )'ZHEEV(V,' // UPLO // ')',
01604      $               IINFO, N, JTYPE, IOLDSD
01605                   INFO = ABS( IINFO )
01606                   IF( IINFO.LT.0 ) THEN
01607                      RETURN
01608                   ELSE
01609                      RESULT( NTEST ) = ULPINV
01610                      RESULT( NTEST+1 ) = ULPINV
01611                      RESULT( NTEST+2 ) = ULPINV
01612                      GO TO 950
01613                   END IF
01614                END IF
01615 *
01616 *              Do tests 37 and 38
01617 *
01618                CALL ZHET21( 1, UPLO, N, 0, V, LDU, D1, D2, A, LDU, Z,
01619      $                      LDU, TAU, WORK, RWORK, RESULT( NTEST ) )
01620 *
01621                CALL ZLACPY( ' ', N, N, V, LDU, A, LDA )
01622 *
01623                NTEST = NTEST + 2
01624                CALL ZHEEV( 'N', UPLO, N, A, LDU, D3, WORK, LWORK, RWORK,
01625      $                     IINFO )
01626                IF( IINFO.NE.0 ) THEN
01627                   WRITE( NOUNIT, FMT = 9999 )'ZHEEV(N,' // UPLO // ')',
01628      $               IINFO, N, JTYPE, IOLDSD
01629                   INFO = ABS( IINFO )
01630                   IF( IINFO.LT.0 ) THEN
01631                      RETURN
01632                   ELSE
01633                      RESULT( NTEST ) = ULPINV
01634                      GO TO 950
01635                   END IF
01636                END IF
01637 *
01638 *              Do test 39
01639 *
01640                TEMP1 = ZERO
01641                TEMP2 = ZERO
01642                DO 940 J = 1, N
01643                   TEMP1 = MAX( TEMP1, ABS( D1( J ) ), ABS( D3( J ) ) )
01644                   TEMP2 = MAX( TEMP2, ABS( D1( J )-D3( J ) ) )
01645   940          CONTINUE
01646                RESULT( NTEST ) = TEMP2 / MAX( UNFL,
01647      $                           ULP*MAX( TEMP1, TEMP2 ) )
01648 *
01649   950          CONTINUE
01650 *
01651                CALL ZLACPY( ' ', N, N, V, LDU, A, LDA )
01652 *
01653 *              Call ZHPEV
01654 *
01655 *              Load array WORK with the upper or lower triangular
01656 *              part of the matrix in packed form.
01657 *
01658                IF( IUPLO.EQ.1 ) THEN
01659                   INDX = 1
01660                   DO 970 J = 1, N
01661                      DO 960 I = 1, J
01662                         WORK( INDX ) = A( I, J )
01663                         INDX = INDX + 1
01664   960                CONTINUE
01665   970             CONTINUE
01666                ELSE
01667                   INDX = 1
01668                   DO 990 J = 1, N
01669                      DO 980 I = J, N
01670                         WORK( INDX ) = A( I, J )
01671                         INDX = INDX + 1
01672   980                CONTINUE
01673   990             CONTINUE
01674                END IF
01675 *
01676                NTEST = NTEST + 1
01677                INDWRK = N*( N+1 ) / 2 + 1
01678                CALL ZHPEV( 'V', UPLO, N, WORK, D1, Z, LDU,
01679      $                     WORK( INDWRK ), RWORK, IINFO )
01680                IF( IINFO.NE.0 ) THEN
01681                   WRITE( NOUNIT, FMT = 9999 )'ZHPEV(V,' // UPLO // ')',
01682      $               IINFO, N, JTYPE, IOLDSD
01683                   INFO = ABS( IINFO )
01684                   IF( IINFO.LT.0 ) THEN
01685                      RETURN
01686                   ELSE
01687                      RESULT( NTEST ) = ULPINV
01688                      RESULT( NTEST+1 ) = ULPINV
01689                      RESULT( NTEST+2 ) = ULPINV
01690                      GO TO 1050
01691                   END IF
01692                END IF
01693 *
01694 *              Do tests 40 and 41.
01695 *
01696                CALL ZHET21( 1, UPLO, N, 0, A, LDA, D1, D2, Z, LDU, V,
01697      $                      LDU, TAU, WORK, RWORK, RESULT( NTEST ) )
01698 *
01699                IF( IUPLO.EQ.1 ) THEN
01700                   INDX = 1
01701                   DO 1010 J = 1, N
01702                      DO 1000 I = 1, J
01703                         WORK( INDX ) = A( I, J )
01704                         INDX = INDX + 1
01705  1000                CONTINUE
01706  1010             CONTINUE
01707                ELSE
01708                   INDX = 1
01709                   DO 1030 J = 1, N
01710                      DO 1020 I = J, N
01711                         WORK( INDX ) = A( I, J )
01712                         INDX = INDX + 1
01713  1020                CONTINUE
01714  1030             CONTINUE
01715                END IF
01716 *
01717                NTEST = NTEST + 2
01718                INDWRK = N*( N+1 ) / 2 + 1
01719                CALL ZHPEV( 'N', UPLO, N, WORK, D3, Z, LDU,
01720      $                     WORK( INDWRK ), RWORK, IINFO )
01721                IF( IINFO.NE.0 ) THEN
01722                   WRITE( NOUNIT, FMT = 9999 )'ZHPEV(N,' // UPLO // ')',
01723      $               IINFO, N, JTYPE, IOLDSD
01724                   INFO = ABS( IINFO )
01725                   IF( IINFO.LT.0 ) THEN
01726                      RETURN
01727                   ELSE
01728                      RESULT( NTEST ) = ULPINV
01729                      GO TO 1050
01730                   END IF
01731                END IF
01732 *
01733 *              Do test 42
01734 *
01735                TEMP1 = ZERO
01736                TEMP2 = ZERO
01737                DO 1040 J = 1, N
01738                   TEMP1 = MAX( TEMP1, ABS( D1( J ) ), ABS( D3( J ) ) )
01739                   TEMP2 = MAX( TEMP2, ABS( D1( J )-D3( J ) ) )
01740  1040          CONTINUE
01741                RESULT( NTEST ) = TEMP2 / MAX( UNFL,
01742      $                           ULP*MAX( TEMP1, TEMP2 ) )
01743 *
01744  1050          CONTINUE
01745 *
01746 *              Call ZHBEV
01747 *
01748                IF( JTYPE.LE.7 ) THEN
01749                   KD = 0
01750                ELSE IF( JTYPE.GE.8 .AND. JTYPE.LE.15 ) THEN
01751                   KD = MAX( N-1, 0 )
01752                ELSE
01753                   KD = IHBW
01754                END IF
01755 *
01756 *              Load array V with the upper or lower triangular part
01757 *              of the matrix in band form.
01758 *
01759                IF( IUPLO.EQ.1 ) THEN
01760                   DO 1070 J = 1, N
01761                      DO 1060 I = MAX( 1, J-KD ), J
01762                         V( KD+1+I-J, J ) = A( I, J )
01763  1060                CONTINUE
01764  1070             CONTINUE
01765                ELSE
01766                   DO 1090 J = 1, N
01767                      DO 1080 I = J, MIN( N, J+KD )
01768                         V( 1+I-J, J ) = A( I, J )
01769  1080                CONTINUE
01770  1090             CONTINUE
01771                END IF
01772 *
01773                NTEST = NTEST + 1
01774                CALL ZHBEV( 'V', UPLO, N, KD, V, LDU, D1, Z, LDU, WORK,
01775      $                     RWORK, IINFO )
01776                IF( IINFO.NE.0 ) THEN
01777                   WRITE( NOUNIT, FMT = 9998 )'ZHBEV(V,' // UPLO // ')',
01778      $               IINFO, N, KD, JTYPE, IOLDSD
01779                   INFO = ABS( IINFO )
01780                   IF( IINFO.LT.0 ) THEN
01781                      RETURN
01782                   ELSE
01783                      RESULT( NTEST ) = ULPINV
01784                      RESULT( NTEST+1 ) = ULPINV
01785                      RESULT( NTEST+2 ) = ULPINV
01786                      GO TO 1140
01787                   END IF
01788                END IF
01789 *
01790 *              Do tests 43 and 44.
01791 *
01792                CALL ZHET21( 1, UPLO, N, 0, A, LDA, D1, D2, Z, LDU, V,
01793      $                      LDU, TAU, WORK, RWORK, RESULT( NTEST ) )
01794 *
01795                IF( IUPLO.EQ.1 ) THEN
01796                   DO 1110 J = 1, N
01797                      DO 1100 I = MAX( 1, J-KD ), J
01798                         V( KD+1+I-J, J ) = A( I, J )
01799  1100                CONTINUE
01800  1110             CONTINUE
01801                ELSE
01802                   DO 1130 J = 1, N
01803                      DO 1120 I = J, MIN( N, J+KD )
01804                         V( 1+I-J, J ) = A( I, J )
01805  1120                CONTINUE
01806  1130             CONTINUE
01807                END IF
01808 *
01809                NTEST = NTEST + 2
01810                CALL ZHBEV( 'N', UPLO, N, KD, V, LDU, D3, Z, LDU, WORK,
01811      $                     RWORK, IINFO )
01812                IF( IINFO.NE.0 ) THEN
01813                   WRITE( NOUNIT, FMT = 9998 )'ZHBEV(N,' // UPLO // ')',
01814      $               IINFO, N, KD, JTYPE, IOLDSD
01815                   INFO = ABS( IINFO )
01816                   IF( IINFO.LT.0 ) THEN
01817                      RETURN
01818                   ELSE
01819                      RESULT( NTEST ) = ULPINV
01820                      GO TO 1140
01821                   END IF
01822                END IF
01823 *
01824  1140          CONTINUE
01825 *
01826 *              Do test 45.
01827 *
01828                TEMP1 = ZERO
01829                TEMP2 = ZERO
01830                DO 1150 J = 1, N
01831                   TEMP1 = MAX( TEMP1, ABS( D1( J ) ), ABS( D3( J ) ) )
01832                   TEMP2 = MAX( TEMP2, ABS( D1( J )-D3( J ) ) )
01833  1150          CONTINUE
01834                RESULT( NTEST ) = TEMP2 / MAX( UNFL,
01835      $                           ULP*MAX( TEMP1, TEMP2 ) )
01836 *
01837                CALL ZLACPY( ' ', N, N, A, LDA, V, LDU )
01838                NTEST = NTEST + 1
01839                CALL ZHEEVR( 'V', 'A', UPLO, N, A, LDU, VL, VU, IL, IU,
01840      $                      ABSTOL, M, WA1, Z, LDU, IWORK, WORK, LWORK,
01841      $                      RWORK, LRWORK, IWORK( 2*N+1 ), LIWORK-2*N,
01842      $                      IINFO )
01843                IF( IINFO.NE.0 ) THEN
01844                   WRITE( NOUNIT, FMT = 9999 )'ZHEEVR(V,A,' // UPLO //
01845      $               ')', IINFO, N, JTYPE, IOLDSD
01846                   INFO = ABS( IINFO )
01847                   IF( IINFO.LT.0 ) THEN
01848                      RETURN
01849                   ELSE
01850                      RESULT( NTEST ) = ULPINV
01851                      RESULT( NTEST+1 ) = ULPINV
01852                      RESULT( NTEST+2 ) = ULPINV
01853                      GO TO 1170
01854                   END IF
01855                END IF
01856 *
01857 *              Do tests 45 and 46 (or ... )
01858 *
01859                CALL ZLACPY( ' ', N, N, V, LDU, A, LDA )
01860 *
01861                CALL ZHET21( 1, UPLO, N, 0, A, LDU, WA1, D2, Z, LDU, V,
01862      $                      LDU, TAU, WORK, RWORK, RESULT( NTEST ) )
01863 *
01864                NTEST = NTEST + 2
01865                CALL ZHEEVR( 'N', 'A', UPLO, N, A, LDU, VL, VU, IL, IU,
01866      $                      ABSTOL, M2, WA2, Z, LDU, IWORK, WORK, LWORK,
01867      $                      RWORK, LRWORK, IWORK( 2*N+1 ), LIWORK-2*N,
01868      $                      IINFO )
01869                IF( IINFO.NE.0 ) THEN
01870                   WRITE( NOUNIT, FMT = 9999 )'ZHEEVR(N,A,' // UPLO //
01871      $               ')', IINFO, N, JTYPE, IOLDSD
01872                   INFO = ABS( IINFO )
01873                   IF( IINFO.LT.0 ) THEN
01874                      RETURN
01875                   ELSE
01876                      RESULT( NTEST ) = ULPINV
01877                      GO TO 1170
01878                   END IF
01879                END IF
01880 *
01881 *              Do test 47 (or ... )
01882 *
01883                TEMP1 = ZERO
01884                TEMP2 = ZERO
01885                DO 1160 J = 1, N
01886                   TEMP1 = MAX( TEMP1, ABS( WA1( J ) ), ABS( WA2( J ) ) )
01887                   TEMP2 = MAX( TEMP2, ABS( WA1( J )-WA2( J ) ) )
01888  1160          CONTINUE
01889                RESULT( NTEST ) = TEMP2 / MAX( UNFL,
01890      $                           ULP*MAX( TEMP1, TEMP2 ) )
01891 *
01892  1170          CONTINUE
01893 *
01894                NTEST = NTEST + 1
01895                CALL ZLACPY( ' ', N, N, V, LDU, A, LDA )
01896                CALL ZHEEVR( 'V', 'I', UPLO, N, A, LDU, VL, VU, IL, IU,
01897      $                      ABSTOL, M2, WA2, Z, LDU, IWORK, WORK, LWORK,
01898      $                      RWORK, LRWORK, IWORK( 2*N+1 ), LIWORK-2*N,
01899      $                      IINFO )
01900                IF( IINFO.NE.0 ) THEN
01901                   WRITE( NOUNIT, FMT = 9999 )'ZHEEVR(V,I,' // UPLO //
01902      $               ')', IINFO, N, JTYPE, IOLDSD
01903                   INFO = ABS( IINFO )
01904                   IF( IINFO.LT.0 ) THEN
01905                      RETURN
01906                   ELSE
01907                      RESULT( NTEST ) = ULPINV
01908                      RESULT( NTEST+1 ) = ULPINV
01909                      RESULT( NTEST+2 ) = ULPINV
01910                      GO TO 1180
01911                   END IF
01912                END IF
01913 *
01914 *              Do tests 48 and 49 (or +??)
01915 *
01916                CALL ZLACPY( ' ', N, N, V, LDU, A, LDA )
01917 *
01918                CALL ZHET22( 1, UPLO, N, M2, 0, A, LDU, WA2, D2, Z, LDU,
01919      $                      V, LDU, TAU, WORK, RWORK, RESULT( NTEST ) )
01920 *
01921                NTEST = NTEST + 2
01922                CALL ZLACPY( ' ', N, N, V, LDU, A, LDA )
01923                CALL ZHEEVR( 'N', 'I', UPLO, N, A, LDU, VL, VU, IL, IU,
01924      $                      ABSTOL, M3, WA3, Z, LDU, IWORK, WORK, LWORK,
01925      $                      RWORK, LRWORK, IWORK( 2*N+1 ), LIWORK-2*N,
01926      $                      IINFO )
01927                IF( IINFO.NE.0 ) THEN
01928                   WRITE( NOUNIT, FMT = 9999 )'ZHEEVR(N,I,' // UPLO //
01929      $               ')', IINFO, N, JTYPE, IOLDSD
01930                   INFO = ABS( IINFO )
01931                   IF( IINFO.LT.0 ) THEN
01932                      RETURN
01933                   ELSE
01934                      RESULT( NTEST ) = ULPINV
01935                      GO TO 1180
01936                   END IF
01937                END IF
01938 *
01939 *              Do test 50 (or +??)
01940 *
01941                TEMP1 = DSXT1( 1, WA2, M2, WA3, M3, ABSTOL, ULP, UNFL )
01942                TEMP2 = DSXT1( 1, WA3, M3, WA2, M2, ABSTOL, ULP, UNFL )
01943                RESULT( NTEST ) = ( TEMP1+TEMP2 ) /
01944      $                           MAX( UNFL, ULP*TEMP3 )
01945  1180          CONTINUE
01946 *
01947                NTEST = NTEST + 1
01948                CALL ZLACPY( ' ', N, N, V, LDU, A, LDA )
01949                CALL ZHEEVR( 'V', 'V', UPLO, N, A, LDU, VL, VU, IL, IU,
01950      $                      ABSTOL, M2, WA2, Z, LDU, IWORK, WORK, LWORK,
01951      $                      RWORK, LRWORK, IWORK( 2*N+1 ), LIWORK-2*N,
01952      $                      IINFO )
01953                IF( IINFO.NE.0 ) THEN
01954                   WRITE( NOUNIT, FMT = 9999 )'ZHEEVR(V,V,' // UPLO //
01955      $               ')', IINFO, N, JTYPE, IOLDSD
01956                   INFO = ABS( IINFO )
01957                   IF( IINFO.LT.0 ) THEN
01958                      RETURN
01959                   ELSE
01960                      RESULT( NTEST ) = ULPINV
01961                      RESULT( NTEST+1 ) = ULPINV
01962                      RESULT( NTEST+2 ) = ULPINV
01963                      GO TO 1190
01964                   END IF
01965                END IF
01966 *
01967 *              Do tests 51 and 52 (or +??)
01968 *
01969                CALL ZLACPY( ' ', N, N, V, LDU, A, LDA )
01970 *
01971                CALL ZHET22( 1, UPLO, N, M2, 0, A, LDU, WA2, D2, Z, LDU,
01972      $                      V, LDU, TAU, WORK, RWORK, RESULT( NTEST ) )
01973 *
01974                NTEST = NTEST + 2
01975                CALL ZLACPY( ' ', N, N, V, LDU, A, LDA )
01976                CALL ZHEEVR( 'N', 'V', UPLO, N, A, LDU, VL, VU, IL, IU,
01977      $                      ABSTOL, M3, WA3, Z, LDU, IWORK, WORK, LWORK,
01978      $                      RWORK, LRWORK, IWORK( 2*N+1 ), LIWORK-2*N,
01979      $                      IINFO )
01980                IF( IINFO.NE.0 ) THEN
01981                   WRITE( NOUNIT, FMT = 9999 )'ZHEEVR(N,V,' // UPLO //
01982      $               ')', IINFO, N, JTYPE, IOLDSD
01983                   INFO = ABS( IINFO )
01984                   IF( IINFO.LT.0 ) THEN
01985                      RETURN
01986                   ELSE
01987                      RESULT( NTEST ) = ULPINV
01988                      GO TO 1190
01989                   END IF
01990                END IF
01991 *
01992                IF( M3.EQ.0 .AND. N.GT.0 ) THEN
01993                   RESULT( NTEST ) = ULPINV
01994                   GO TO 1190
01995                END IF
01996 *
01997 *              Do test 52 (or +??)
01998 *
01999                TEMP1 = DSXT1( 1, WA2, M2, WA3, M3, ABSTOL, ULP, UNFL )
02000                TEMP2 = DSXT1( 1, WA3, M3, WA2, M2, ABSTOL, ULP, UNFL )
02001                IF( N.GT.0 ) THEN
02002                   TEMP3 = MAX( ABS( WA1( 1 ) ), ABS( WA1( N ) ) )
02003                ELSE
02004                   TEMP3 = ZERO
02005                END IF
02006                RESULT( NTEST ) = ( TEMP1+TEMP2 ) /
02007      $                           MAX( UNFL, TEMP3*ULP )
02008 *
02009                CALL ZLACPY( ' ', N, N, V, LDU, A, LDA )
02010 *
02011 *
02012 *
02013 *
02014 *              Load array V with the upper or lower triangular part
02015 *              of the matrix in band form.
02016 *
02017  1190          CONTINUE
02018 *
02019  1200       CONTINUE
02020 *
02021 *           End of Loop -- Check for RESULT(j) > THRESH
02022 *
02023             NTESTT = NTESTT + NTEST
02024             CALL DLAFTS( 'ZST', N, N, JTYPE, NTEST, RESULT, IOLDSD,
02025      $                   THRESH, NOUNIT, NERRS )
02026 *
02027  1210    CONTINUE
02028  1220 CONTINUE
02029 *
02030 *     Summary
02031 *
02032       CALL ALASVM( 'ZST', NOUNIT, NERRS, NTESTT, 0 )
02033 *
02034  9999 FORMAT( ' ZDRVST: ', A, ' returned INFO=', I6, / 9X, 'N=', I6,
02035      $      ', JTYPE=', I6, ', ISEED=(', 3( I5, ',' ), I5, ')' )
02036  9998 FORMAT( ' ZDRVST: ', A, ' returned INFO=', I6, / 9X, 'N=', I6,
02037      $      ', KD=', I6, ', JTYPE=', I6, ', ISEED=(', 3( I5, ',' ), I5,
02038      $      ')' )
02039 *
02040       RETURN
02041 *
02042 *     End of ZDRVST
02043 *
02044       END
 All Files Functions