File: DIAMAN.FT of Disk: Disks/MyPDP/m8-2-rka1-rkb1

(Source file text) 

---

C	TEST PROGRAM FOR DIAGON SUBROUTINE
	COMMON /FITPAR/ JIT,FISTOP,AMARQI,UPMARQ,DNMARQ,VARI,
     $  UPVAR,DNVAR,VARMIN,MLOOP,LOOPLW,ISTART,ISTOP,ITEST
	COMMON /DIAG/ A(24,24),Z(24,24),D(24),E(24)
	COMMON /SIZE/ N,NP,K1,K3
	COMMON /PARAM/ DKHI,SHFACT,CHISQ,ITORQ,IDILA,DK
     $  ,HUP,PS,TEMP,AKL,AKV,CAP,EICH
	COMMON /VECT/ C(24),D1(6,4),Q(6,4),DIAGEL(24)
	COMMON /DATIN/ F(256)
	COMMON /CONST/ PI
C
	LOGICAL ITEST
C
	ITEST=.FALSE.
	WRITE (4,100)
100	FORMAT(1H1,'DIAGONALIZATION OF REAL SYMMETRIC MATRIX',/)
	WRITE(4,101)
101	FORMAT(1X,'GIVE SIZE OF MATRIX: ',$)
	READ(4,102)N
102	FORMAT(1I2,/)
	WRITE(4,103)
103	FORMAT(1X,'ENTER MATRIX ELEMENTS',/)
	DO 10 J=1,N
	DO 20 K=J,N
	WRITE(4,104)J,K
104	FORMAT(1H ,'A(',I2,',',I2,')= ',$)
	READ(4,105) A(J,K)
20	CONTINUE
	WRITE(4,106)
10	CONTINUE
105	FORMAT(1G14.6)
106	FORMAT(1H )
	CALL DIAGON
	WRITE (4,107) JIT
107	FORMAT(1X,'# OF ITERATIONS=',I2)
	WRITE (4,108)
108	FORMAT (1H0,'EIGENVALUES:',//)
	WRITE(4,106)
	WRITE (4,109) (D(J),J=1,N)
109	FORMAT (10F13.5)
	WRITE (4,110)
110	FORMAT(1X,//,1X,'EIGENVECTORS:',//)
	DO 40 J=1,N
	WRITE(4,106)
	WRITE (4,111) (Z(J,K),K=1,N)
111	FORMAT (10F13.5)
40	CONTINUE
	END

---