ZLAEIN uses inverse iteration to find a right or left eigenvector
 corresponding to the eigenvalue W of a complex upper Hessenberg
 matrix H.

          RIGHTV is LOGICAL
          = .TRUE. : compute right eigenvector;
          = .FALSE.: compute left eigenvector.

          NOINIT is LOGICAL
          = .TRUE. : no initial vector supplied in V
          = .FALSE.: initial vector supplied in V.

          N is INTEGER
          The order of the matrix H.  N >= 0.

          H is COMPLEX*16 array, dimension (LDH,N)
          The upper Hessenberg matrix H.

          LDH is INTEGER
          The leading dimension of the array H.  LDH >= max(1,N).

          W is COMPLEX*16
          The eigenvalue of H whose corresponding right or left
          eigenvector is to be computed.

          V is COMPLEX*16 array, dimension (N)
          On entry, if NOINIT = .FALSE., V must contain a starting
          vector for inverse iteration; otherwise V need not be set.
          On exit, V contains the computed eigenvector, normalized so
          that the component of largest magnitude has magnitude 1; here
          the magnitude of a complex number (x,y) is taken to be
          |x| + |y|.

          B is COMPLEX*16 array, dimension (LDB,N)

          LDB is INTEGER
          The leading dimension of the array B.  LDB >= max(1,N).

          RWORK is DOUBLE PRECISION array, dimension (N)

          EPS3 is DOUBLE PRECISION
          A small machine-dependent value which is used to perturb
          close eigenvalues, and to replace zero pivots.

          SMLNUM is DOUBLE PRECISION
          A machine-dependent value close to the underflow threshold.

          INFO is INTEGER
          = 0:  successful exit
          = 1:  inverse iteration did not converge; V is set to the
                last iterate.