lit
    JINRLIB

EISRS1, EISRS2, EISRS3, EISRS4 - .

F224

: E.Edberg :

p opopa EISRS1,EISRS2,EISRS3,EISRS4 c c oop coc a coc op co cpo ap, co a popa EISPACK(F220).

EISRS1 c c coc a coc op.
EISRS2 c c coc a.
EISRS3 c coc a pa [RLB,RUB] cooc coc op.
EISRS4 c coc a pa [RLB,RUB].

:

: - SUBROUTINE
: - EISRS1, EISRS2, EISRS3, EISRS4
: - TRED1,TRED2,TQL1,TQL2,BISECT,
TINVIT,TRBAK1 (EISPACK,F220)

:

CALL EISRS1(NM,N,AR,WR,ZR,IERR,WORK)
CALL EISRS2(NM,N,AR,WR,IERR,WORK)
CALL EISRS3(NM,N,AR,RLB,RUB,MM,M,WR,ZR,IERR,WORK,IWORK)
CALL EISRS4(NM,N,AR,RLB,RUB,MM,M,WR,IERR,WORK,IWORK), :

NM - acco AR ZR opaop DIMENSION a popa;
N - opo ap;
AR - ca cpa apa papoc N*N;
WR - c oop acc, copa oc a. apoc WR oa N, a EISRS3 EISRS4 - MM;
ZR - c p acc, copa oa coc op. ZR o N coo, a EISRS3 - MM;
IERR - a pa, caa coo o oa. c a o IERR ao 0, o c oop;
WORK - c oop pao acc, apoc oopoo N,N,8*N+MM 5*N+MM cooco popaa;
RLB,RUB - aac pa [RLB,RUB] c coc a o pa;
MM - ao acao co oc a pa [RLB,RUB];
M - co a coc a. c M > MM, o IERR pao 0;
IWORK - oop pao acc papoc N.

:

o p po ap AR copac ( EISRS1 copac oa apa).

:

  1. J.M. Bol, B.S.Garow, B.T.Smit, J.Ik, V.C.Klema, C.B.Moler.
    Matrix Eigensystem Routines - EISPACK Guide. Scond dition.
    Lcture Notes in Computer Science. Vol.6.
    Springr-Vrlag, Nw Yor. 1976.
  2. . EISPACK (F220).
  3. http://www.netlib.org/eispack
:
       IMPLICIT REAL*8 (A-H,O-Z)
       DIMENSION AR(10,10),AI(10,10),WR(10),WI(10),ZR(10,10),ZI(10,10),
      *WORK(100),BR(10,10),BI(10,10),CR(10,10),CI(10,10),
      *DR(10,10),DI(10,10),IWORK(10)
       . . .
     C   
       . . .
       CALL EISRS1(NM,N,AR,WR,ZR,IERR,WORK)
       . . .
       CALL EISRS2(NM,N,BR,WR,IERR,WORK)
       . . .
       CALL EISRS3(NM,N,AR,RLB,RUB,MM,M,WR,ZR,IERR,WORK,IWORK)
       . . .
       CALL EISRS4(NM,N,AR,RLB,RUB,MM,M,WR,IERR,WORK,IWORK)
:
              REAL MATRIX
 
       5.000     4.000     1.000     1.000
       4.000     5.000     1.000     1.000
       1.000     1.000     4.000     2.000
       1.000     1.000     2.000     4.000
 
              REAL EIGENVALUES
 
    .1000000000E+01     .2000000000E+01     .5000000000E+01     .1000000000E+02
 
              REAL EIGENVECTORS
 
    .7071067812E+00    -.7071067812E+00     .4532466518E-16     .0000000000E+00
 
    .1313155178E-15     .5075252013E-16    -.7071067812E+00     .7071067812E+00
  
    .3162277660E+00     .3162277660E+00    -.6324555320E+00    -.6324555320E+00
 
   -.6324555320E+00    -.6324555320E+00    -.3162277660E+00    -.3162277660E+00


home up e-mail