lit
    JINRLIB

EISRG1, EISRG2 - .

F223

: E.Edberg :

popa EISRG1 EISRG2 c coc a oc op co ap oo a, co a popa EISPACK(F220).

EISRG1 c c coc a coc op.
EISRG2 c oo coc a.

:

: - SUBROUTINE
: - EISRG1, EISRG2
: - BALANC,ELMHES,ELTRAN,HQR,
BALBAK,HQR2 (EISPACK,F220)

:

CALL EISRG1(NM,N,AR,WR,WI,ZR,IERR,WORK,IWORK)
CALL EISRG2(NM,N,AR,WR,WI,IERR,IWORK), :

NM - acco AR ZR opaop DIMENSION a popa;
N - opo ap, N ≤ NM;
AR - c p acc papoc N*N, opa co ap;
WR,WI - c oop acc papoc N, copa cooco c ac c coc a.
ap cop oc coc a a ocoao;
ZR - c p acc papoc N*N, opa coa coc op.
Koc coc op, cooc
(I)-o (I+1)-o coc a o a a a
(I)-o coa + i*(I+1)-o coa ooco
(I)-o coa - i*(I+1)-o oa acca ZR;
IERR - a pa, caa coo o a. c IERR pao 0, o c oop;
WORK - c oop pao acc papoc N;
IWORK - oop pao acc papoc N.

:

coa apa 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(10),BR(10,10),BI(10,10),IWORK(10)
       . . .
     C   
       . . .
       CALL EISRG1(NM,N,AR,WR,WI,ZR,IERR,WORK,IWORK)
       . . .
       CALL EISRG2(NM,N,AR,WR,WI,IERR,IWORK)
:
              REAL MATRIX

      4.000    -5.000      .000     3.000
       .000     4.000    -3.000    -5.000
      5.000    -3.000     4.000      .000
      3.000      .000     5.000     4.000

              COMPLEX EIGENVALUES

    .1200000000E+02     .0000000000E+00     .1000000000E+01     .5000000000E+01
    .1000000000E+01    -.5000000000E+01     .2000000000E+01     .0000000000E+00

              REAL EIGENVECTORS

    .5600876295E+00    -.5600876295E+00     .5600876295E+00     .5600876295E+00

    .2182380549E+00     .8047631930E+00     .8047631930E+00    -.2182380549E+00

    .8047631930E+00    -.2182380549E+00    -.2182380549E+00    -.8047631930E+00

    .3697784622E+00     .3697784622E+00    -.3697784622E+00     .3697784622E+00


home up e-mail