* * $Id: eisch3.F,v 1.1.1.1 1996/04/01 15:02:39 mclareni Exp $ * * $Log: eisch3.F,v $ * Revision 1.1.1.1 1996/04/01 15:02:39 mclareni * Mathlib gen * * #include "gen/pilot.h" SUBROUTINE EISCH3(NM,N,AR,AI,RLB,RUB,MM,M,WR,ZR,ZI,IERR,WORK) C SOME EIGENVALUES AND CORRESPONDING EIGENVECTORS OF A COMPLEX C HERMITIAN MATRIX C DIMENSION AR(NM,NM),AI(NM,NM),WR(MM),ZR(NM,NM),ZI(NM,NM),WORK(1) CALL HTRIDI(NM,N,AR,AI,ZI,WORK,WORK(N+1),WORK(2*N+1)) EPS1=0.0 CALL BISECT(N,EPS1,ZI,WORK,WORK(N+1),RLB,RUB,MM,M,WR,WORK(9*N+1), X IERR,WORK(4*N+1),WORK(5*N+1)) IF(IERR.NE.0) RETURN IF(M.EQ.0) RETURN CALL TINVIT(NM,N,ZI,WORK,WORK(N+1),M,WR,WORK(9*N+1),ZR,IERR, X WORK(4*N+1),WORK(5*N+1),WORK(6*N+1),WORK(7*N+1),WORK(8*N+1)) IF(IERR.NE.0) RETURN CALL HTRIBK(NM,N,AR,AI,WORK(2*N+1),M,ZR,ZI) RETURN END