lit
    JINRLIB

EISCG1, EISCG2 - .

F221

: E.Edberg :

opopa EISCG1 EISCG2 c coc a coc op oco ap o a, co a popa EISPACK(F220).
EISCG1 c c coc a coc op.
EISCG2 c oo coc a.

:

: - SUBROUTINE
: - EISCG1, EISCG2
: - CBAL,COMHES,COMLR2,CBABK2,
COMLR (EISPACK,F220)

:

CALL EISCG1(NM,N,AR,AI,WR,WI,ZR,ZI,IERR,WORK,IWORK)
CALL EISCG2(NM,N,AR,AI,WR,WI,IERR,IWORK), :

NM - acco AR,AI,ZR,ZI opaop DIMENSION a popa;
N - opo ap, N ≤ NM;
AR,AI - c p acc papoc N*N,
opa cooco c ac coo ap;
WR,WI - c oop acc papoc N,
copa cooco c ac c coc a;
ZR,ZI - c p acc papoc N*N,
opa coa cooco c ac c coc opo;
IERR - a pa, caa oo o oa, oap o pao popa.
c a o IERR pao 0, o c oop;
WORK - c oop pao acc papoc N;
IWORK - oop pao acc papoc N.

:

co ap AR AI oc pao popa opac.

:

  1. J.M. Bol, B.S.Garow, B.T.Smit et.al.
    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(10),IWORK(10),BR(10,10),BI(10,10)
       . . .
    C    
       . . .
       CALL EISCG1(NM,N,AR,AI,WR,WI,ZR,ZI,IERR,WORK,IWORK)
       . . .
       CALL EISCG2(NM,N,AR,AI,WR,WI,IERR,IWORK)
:
              COMPLEX MATRIX
 
    5.000    9.000      5.000    5.000     -6.000   -6.000     -7.000   -7.000
    3.000    3.000      6.000   10.000     -5.000   -5.000     -6.000   -6.000
    2.000    2.000      3.000    3.000     -1.000    3.000     -5.000   -5.000
    1.000    1.000      2.000    2.000     -3.000   -3.000       .000    4.000
 
              COMPLEX EIGENVALUES
 
    .2000000000E+01     .6000000000E+01     .4000000000E+01     .8000000000E+01
    .3000000000E+01     .7000000000E+01     .1000000000E+01     .5000000000E+01
 
              COMPLEX EIGENVECTORS
 
    .5000000000E+00    -.4497138827E-16     .1000000000E+01     .0000000000E+00
    .5000000000E+00     .7845719978E-16     .5000000000E+00    -.5863178863E-16
 
   -.2000000000E+01    -.6657527432E-15    -.2000000000E+01    -.1110223025E-15
   -.2000000000E+01    -.9347308387E-16     .4440892099E-15     .6175242603E-16
 
    .2000000000E+01    -.9011494476E-15     .2000000000E+01    -.1110223025E-14
    .1887379142E-14     .1648445640E-15     .2000000000E+01    -.1116464299E-14
    
    .2000000000E+01     .1785516568E-15     .1000000000E+01     .7771561172E-15
    .1000000000E+01     .6005263602E-15     .1000000000E+01     .2559314017E-15


home up e-mail