lit
    JINRLIB

EISST1, EISST2, EISST3, EISST4 - .

F225

: E.Edberg :

popa c c oop coc a coc op co cpo paoao ap, co a popa EISPACK(F220).

EISST1 c c coc a coc op.
EISST2 c c coc a.
EISST3 c coc a pa [RLB,RUB] cooc coc op.
EISST4 c coc a pa [RLB,RUB].

:

: - SUBROUTINE
: - EISST1 EISST2 EISST3 EISST4
: - IMTQL1,IMTQL2,BISECT,TINVIT,
(EISPACK,F220)

:

CALL EISST1(NM,N,AR,WR,ZR,IERR,WORK)
CALL EISST2(NM,N,AR,WR,IERR,WORK)
CALL EISST3(NM,N,PR,RLB,RUB,MM,M,WR,ZR,IERR,WORK,IWORK)
CALL EISST4(NM,N,AR,RLB,RUB,MM,M,WR,IERR,WORK,IWORK), :

NM - acco AR ZR opaop DIMENSION a popa;
N - opo ap, N ≤ NM;
AR - c p acc papoc N*2.
Bopo co cop ao aoa,
a A(2,1),A(3,1),...,A(N,1) aoc oaoa ap;
WR - c oop acc, copa c coc a. apoc WR oa N, a EISST3, EISST4 - MM;
ZR - c p acc, copa c coc op. ZR o N coo, a EISST3 - MM;
IERR - a pa, caa coo o a. c a o IERR pao 0, o c oop;
WORK - pao acc papoc N,N,8*N+MM 5*N+MM cooco popaa;
RLB,RUB - aac pa [RLB,RUB] c oc 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.

:

coa apa AR oc pao popa copac.

:

  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 EISST1(NM,N,AR,WR,ZR,IERR,WORK)
       . . .
       CALL EISST2(NM,N,BR,WR,IERR,WORK)
       . . .
       CALL EISST3(NM,N,PR,RLB,RUB,MM,M,WR,ZR,IERR,WORK,IWORK)
       . . .
       CALL EISST4(NM,N,AR,RLB,RUB,MM,M,WR,IERR,WORK,IWORK)
:
      SYMMETRIC TRIDIAGONAL MATRIX, SUBDIAGONAL AND DIAGONAL
 
       4.000     6.000     6.000     4.000
      -4.000   -10.000   -12.000   -10.000    -4.000
 
              REAL EIGENVALUES
 
-.2000000000E+02 -.1200000000E+02 -.6000000000E+01 -.2000000000E+01 -.4387162225E-15
 
              REAL EIGENVECTORS
 
-.1195228609E+00  .4780914437E+00 -.7171371656E+00  .4780914437E+00 -.1195228609E+00
 
-.3162277660E+00  .6324555320E+00  .3295974604E-15 -.6324555320E+00  .3162277660E+00
 
 .5345224838E+00 -.2672612419E+00 -.5345224838E+00 -.2672612419E+00  .5345224838E+00
 
-.6324555320E+00 -.3162277660E+00  .0000000000E+00  .3162277660E+00  .6324555320E+00
 
 .4472135955E+00  .4472135955E+00  .4472135955E+00  .4472135955E+00  .4472135955E+00


home up e-mail