]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 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 |