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