lit
    JINRLIB

tEQN, tINV, tEQINV - ,

F010

: G.A.Erskine :

popaa tEQN ( t=D C) pa apo pa:

  A * X = B , (1)

oopo pca coo cc N coc
pa c pa ac:

F010

popaa tINV apao ap A.
popaa tEQINV pa cc (1), a ap A, o ao , tEQN.
c pc a c op ap A o p ocoao coo a (1) c o ap oo A, o c pa ac B, o c cooa popa F011.
pa a t aa ao popa aa :

t = D DOUBLE PRECISION
t = C COMPLEX*16

:

: - SUBROUTINE
: - tEQN tINV tEQINV
(t=D,C)
: - F010PR
: - tFACT tFEQN tFINV (F011)

: t = D C

CALL tEQN   (N,A,IDIMN,IR,IFAIL,K,B)

CALL tINV   (N,A,IDIMN,IR,IFAIL)

CALL tEQINV (N,A,IDIMN,IR,IFAIL,K,B)

apap:

N - (INTEGER) opo apao ap A;
A - ( cooc c t) p acc, pa apoc oopoo a IDIMN;
IDIMN - (INTEGER) pa papoc acca A ( acca B, c K>1 );
IR - (INTEGER) acc, o pa p N o;
IFAIL - (INTEGER) a o IFAIL=-1, c A cp, IFAIL=0;
K - (INTEGER) coo ap B X;
B - ( cooc c t) p acc, pa papoc oopoo IDIMN.
B o oop, c K=1.

, apa A cpa, IFAIL=0,
acc A B :

tEQN - X oac acc B, aa A copac;
tINV - ap A-1 acc A;
tEQINV - X acc B, ap cc A.

Ec apa A - cpa, o IFAIL=-1.
o ca copo acca A pcao, copo B cac .

:

coc poa aopa pcaoo .

:

c N<1, IDIMN<1, K<1, aac co o o,
popa ppaac .

:

poo, o apa A papoc 10*10 apa papoc
10*3 aoc, cooco, acca A popa, copa opaop:

DIMENSION IR(25)
DOUBLE PRECISION A(25,30), B(25,10)

Toa, A * X = B 10
, 10*3 ,
A,
A 100,
c o :

CALL DEQINV(10,A,25,IR,IFAIL,3,B)
IF ( IFAIL. NE. 0 ) GO TO 100


home up e-mail