lit
    JINRLIB

EISCH1, EISCH2, EISCH3, EISCH4 - .

F222

: E.Edberg :

p popa EISCH1,EISCH2,EISCH3,EISCH4 c c oop coc a coc op oco poo ap, co a popa EISPACK(F220).

EISCH1 c c coc a coc op.
EISCH2 c c coc a.
EISCH3 c coc a pa [RLB,RUB] cooc coc op.
EISCH4 c coc a pa [RLB,RUB].

:

: - SUBROUTINE
: - EISCH1, EISCH2, EISCH3, EISCH4
: - HTRIDI,TQL1,TQL2,BISECT,TINVIT,
HTRIBK (EISPACK,F220)

:

CALL EISCH1(NM,N,AR,AI,WR,ZR,ZI,IERR,WORK)
CALL EISCH2(NM,N,AR,AI,WR,IERR,WORK)
CALL EISCH3(NM,N,AR,AI,RLB,RUB,MM,M,WR,ZR,ZI,IERR,WORK,IWORK)
CALL EISCH4(NM,N,AR,AI,RLB,RUB,MM,M,IERR,WORK,IWORK), :

NM - acco AR,AI,ZR,ZI paop DIMENSION a popa;
N - opo ap, N ≤ NM;
AR,AI - c p acc papoc N*N, opa cooco c c coo ap;
WR - c oop acc, copa c coc a.
apoc WR oa N, EISCH3 EISCH4 - MM;
ZR,ZI - c p acc, copa oa cooco c ac c coc opo.
Koco coo oo N, EISCH3 - MM;
IERR - a pa, caa coo o a. c a o IERR pao 0, o c oop;
WORK - c oop pao acc,
apoc oopoo 2*N, 3*N, 10*N, 5*N+MM cooco opopaa;
RLB,RUB - aac pa [RLB,RUB] c oc a o pa;
MM - acao co oc coc a pa [RLB,RUB];
M - co a coc a.
c M > MM, o co IERR pao 0;
IWORK - oop pao acc papoc N.

:

ooo po poa ap AR a ao aoa ap AI 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 EISCH1(NM,N,AR,AI,WR,ZR,ZI,IERR,WORK)
       . . .
       CALL EISCH2(NM,N,BR,BI,WR,IERR,WORK)
       . . .
       CALL EISCH3(NM,N,AR,AI,RLB,RUB,MM,M,WR,ZR,ZI,IERR,WORK,IWORK)
       . . .
       CALL EISCH4(NM,N,AR,AI,RLB,RUB,MM,M,IERR,WORK,IWORK)
:
              COMPLEX MATRIX
 
    3.000    .000     1.000    .000      .000    .000      .000   2.000
    1.000    .000     3.000    .000      .000  -2.000      .000    .000
     .000    .000      .000   2.000     1.000    .000     1.000    .000
     .000  -2.000      .000    .000     1.000    .000     1.000    .000
 
              REAL EIGENVALUES
 
    -.8284271247E+00     .8681126901E-15     .4000000000E+01     .4828427125E+01
 
              COMPLEX EIGENVECTORS
 
     .0000000000E+00    -.2705980501E+00     .0000000000E+00    -.2705980501E+00
    -.6532814824E+00     .0000000000E+00     .6532814824E+00     .0000000000E+00
 
     .0000000000E+00     .5000000000E+00     .0000000000E+00    -.5000000000E+00
    -.5000000000E+00     .0000000000E+00    -.5000000000E+00     .0000000000E+00
 
     .0000000000E+00    -.5000000000E+00     .0000000000E+00     .5000000000E+00
    -.5000000000E+00     .0000000000E+00    -.5000000000E+00     .0000000000E+00
 
     .0000000000E+00     .6532814824E+00     .0000000000E+00     .6532814824E+00
    -.2705980501E+00     .0000000000E+00     .2705980501E+00     .0000000000E+00


home up e-mail