[u/mrichter/AliRoot.git] / MINICERN / mathlib / gen / f / eisch3.F

*
* $Id$
*
* $Log$
* 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
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(4N+1),WORK(5N+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(4N+1),WORK(5N+1),WORK(6N+1),WORK(7N+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