MINICERN/mathlib/gen/f/eisch3.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/04/01 15:02:39  mclareni
   6 * Mathlib gen
   7 *
   8 *
   9 #include "gen/pilot.h"
  10       SUBROUTINE EISCH3(NM,N,AR,AI,RLB,RUB,MM,M,WR,ZR,ZI,IERR,WORK)
  11 C     SOME EIGENVALUES AND CORRESPONDING EIGENVECTORS OF A COMPLEX
  12 C     HERMITIAN MATRIX
  13 C
  14       DIMENSION AR(NM,NM),AI(NM,NM),WR(MM),ZR(NM,NM),ZI(NM,NM),WORK(1)
  15       CALL HTRIDI(NM,N,AR,AI,ZI,WORK,WORK(N+1),WORK(2*N+1))
  16       EPS1=0.0
  17       CALL BISECT(N,EPS1,ZI,WORK,WORK(N+1),RLB,RUB,MM,M,WR,WORK(9*N+1),
  18      X IERR,WORK(4*N+1),WORK(5*N+1))
  19       IF(IERR.NE.0) RETURN
  20       IF(M.EQ.0) RETURN
  21       CALL TINVIT(NM,N,ZI,WORK,WORK(N+1),M,WR,WORK(9*N+1),ZR,IERR,
  22      X WORK(4*N+1),WORK(5*N+1),WORK(6*N+1),WORK(7*N+1),WORK(8*N+1))
  23       IF(IERR.NE.0) RETURN
  24       CALL HTRIBK(NM,N,AR,AI,WORK(2*N+1),M,ZR,ZI)
  25       RETURN
  26       END