lit
    JINRLIB

tMMLT, tMMLA, tMNLA, CMMLTC... -

F004

: H.Lipps :

popa c pp ap Z=XY Z=XY,
Y - o cop oc ap Y oo ap pa:

Z=XY+Z, Z=XY-Z, Z=-XY+Z, Z=-XY-Z

pa a t aa popa aa apo:

t = I INTEGER
t = D DOUBLE PRECISION
t = C COMPLEX*16

:

: - SUBROUTINE
: - tMMLA tMMLS tMMLT
tMNMA tMNMS (t=I,D,C) CMMLTC

: t = I, D, C

CALL tMMLT(M,N,K,X11,X12,X21,Y11,Y12,Y21,Z11,Z12,Z21,W) Z=XY
CALL tMMLA(M,N,K,X11,X12,X21,Y11,Y12,Y21,Z11,Z12,Z21) Z=XY+Z
CALL tMMLS(M,N,K,X11,X12,X21,Y11,Y12,Y21,Z11,Z12,Z21) Z=XY-Z
CALL tMNLA(M,N,K,X11,X21,X12,Y11,Y21,Y12,Z11,Z12,Z21) Z=-XY+Z
CALL tMNLS(M,N,K,X11,X12,X21,Y11,Y12,Y21,Z11,Z12,Z21) Z=-XY-Z

t = C oo:

CALL MMLT(M,N,K,X11,X12,X21,Y11,Y12,Y21,Z11,Z12,Z21,W) Z=XY

:

M,N,K - (INTEGER) poc a X(M,N), Y(N,K), Z(M,K);
X11,X12,X21 - ( cooc t) acca,
copa x11,x12,x21 ap X;
Y11,Y12,Y21 - ( cooc t) acca,
copa y11,y12,y21 ap Y;
Z11,Z12,Z21 - ( cooc t) acca,
copa z11,z12,z21 ap Z;
W - ( cooc t) ao acc, , c Z ppa X Y.
po ca o coa pa.

M<1, N<1 K<1 c opopa o opa.

:

  1. ap X, Y, Z ac pa coc c coa oppaa, papa o paoao pacoo cpo coo. acoc, ap o pac o cpoa. aa opopaa o ppoa c pacopoao ap. oo, o o o ooo, aa apa .
    app, a opopaa oppoa c ap A = (aij), c 11,12,21 ac a11, a12, a21 ppoa c pacopoao ap A', c 11,12,21 ac a11, a21, a12.
  2. pa apa Z o ppa X Y oo c ca:
    tMMLT: X = XY Y = Y'Y , acc W K o.
    Y = XY X = XX' p co, o acc W M o.
    CMMLTC: X = XY Y = Y'Y p co, o acc W K o.
    Y = XY, X = XX', X'= XX', Y = Y'Y , acc W .
  3. op oo o ap LOC.
    IBM/PC- 16- LOC . , F004, 64 .

:

poo, o p acc A,B,C,D,E, op acc W
ap V cac :

COMPLEX*16 A(9,9),B(9,9),C(9,9),D(9,9),E(9,9),V,W(99)

o apa A papoc 4*5, apa B apoc 5*7
apa c papoc 7*3 pac cooc paa oppaa acca c a aa, oa:

  1. c D = AB :

    CALL CMMLT(4,5,7,A,A(1,2),A(2,1),B,B(1,2),B(2,1),D,D(1,2),D(2,1),V)

    ao pp ap AB papoc 4*7 ocpoo acc W:

    CALL CMMLT(4,5,7,A,A(1,2),A(2,1),B,B(2,1),B(1,2),W,W(2),W(8),V)

    a, o Z11 ooo W(1), Z12 - W(2) Z21 - W(8).

    po ac ooa:
      A(1,1),A(1,2),A(2,1) p A,
      A(1,1),A(2,1),A(1,2) p A',
    - acco B,C,D,E.
    Toa c D = AB:

    CALL CMMLT(4,5,7,A,B,D,V)

  2. c D = A'B' = (AB)':

    CALL CMMLT(7,5,4,B',A',D,V)

    CALL CMMLT(4,5,7,A,B,D',V)

  3. c D = AA' E = A'A:

    CALL CMMLT(4,5,4,A,A',D,V)

    CALL CMMLT(5,4,5,A',A,E,V)

  4. a A a AB a AA':

    CALL CMMLT(4,5,7,A,B,A,W)

    CALL CMMLT(4,5,4,A,A',A,W)

    po opa p pao op W, copa 7 c o, opo p pao op, copa 4 oc a.

  5. c D = AB E = BC = (C'B')':

    CALL CMMLTC(4,5,7,A,B,D,V)

    CALL CMMLTC(3,7,5,C',B',E',V)



home up e-mail