lit
    JINRLIB

tMADD, tMSET, tMRAN, tMUTL... -

F003

: H.Lipps :

popa ap opa c apa.
pa a t aa popa aa
SUBROUTINE FUNCTION ao FUNCTION.

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

:

: - SUBROUTINE FUNCTION
: - tMADD tMCPY tMBIL tMMNA tMMNS
tMMPA tMMPS tMMPY tMRAN tMSCL
tMSET tMSUB tMUTL tUMNA tUMNS
tUMPA tUMPS tUMPY tRSCL
(t=I,D,C)
: - D_RANDOM (V107)

: t = I, D, C

CALL tMSET(M,N,S,Z11,Z12,Z21) zij=s
CALL tMRAN(M,N,A,B,Z11,Z12,Z21) zij=c (..2)
CALL tMCPY(M,N,X11,X12,X21,Z11,Z12,Z21) zij=xij
CALL tMUTL(N,X11,X12,X21) xjk=xkj (j>k)(..3)
CALL tMADD(M,N,X11,X12,X21,Y11,Y12,Y21,Z11,Z12,Z21) zij=xij+yij
CALL tMSUB(M,N,X11,X12,X21,Y11,Y12,Y21,Z11,Z12,Z21) zij=xij-yij
CALL tRSCL(M,N,D1,D2,X11,X12,X21,Z11,Z12,Z21) zij=di*xij
CALL tMSCL(M,N,S,X11,X12,X21,Z11,Z12,Z21) zj=s*xij
CALL tMMPY(M,N,X11,X12,X21,Y1,Y2,Z1,Z2) zi=xi1*y1+...+xin*yn
CALL tMMPA(M,N,X11,X12,X21,Y1,Y2,Z1,Z2) zi=xi1*y1+...+xin*yn+zi
CALL tMMPS(M,N,X11,X12,X21,Y1,Y2,Z1,Z2) zi=xi1*y1+...+xin*yn-zi
CALL tMMNA(M,N,X11,X12,X21,Y1,Y2,Z1,Z2) zi=-xi1*y1-...-xin*yn+zi
CALL tMMNS(M,N,X11,X12,X21,Y1,Y2,Z1,Z2) zi=-xi1*y1-...-xin*yn-zi
CALL tUMPY(N,U11,U12,U22,Y1,Y2,Z1,Z2) zj=ujj*yj+...+ujn*yn
CALL tUMPA(N,U11,U12,U22,Y1,Y2,Z1,Z2) zj=ujj*yj+...+ujn*yn+zj
CALL tUMPS(N,U11,U12,U22,Y1,Y2,Z1,Z2) zj=ujj*yj+...+ujn*yn-zj
CALL tUMNA(N,U11,U12,U22,Y1,Y2,Z1,Z2) zj=-ujj*yj-...-ujn*yn+zj
CALL tUMNS(N,U11,U12,U22,Y1,Y2,Z1,Z2) zj=-ujj*yj-...-ujn*yn-zj
F = tMBIL(N,V1,V2,X11,X12,X21,Y1,Y2) f003

:

M, N - (INTEGER) apoc opo
(i=1,...,M; j,k=1,...,N);
S,A,B - ( cooc t) cap ;
X11,X12,X21 - ( cooc t) acca,
copa x11,x12,x21 ap (xij);
Y11,Y12,Y21 - ( cooc t) acca,
copa y11,y12,y21 ap (yij);
Y1,Y2 - ( cooc t) acca,
copa y1,y2 opa (yj);
V1,V2 - ( cooc t) acca,
copa v1,v2 opa (vk);
D1,D2 - ( cooc t) acca,
copa d1,d2 opa (di);
U11,U12,U22 - ( cooc t) acca,
copa u11,u12,u22 p poo ap (ujk);
Z11,Z12,Z21 - ( cooc t) acca,
copa z11,z12,z21 p ap (zij);
Z1,Z2 - ( cooc t) acca,
copa z1,z2 (zj).

M<1 N<1 c opopa a opa pa a-o c, c oaac pa .

:

  1. op (yj) p ac ao:
    Papac o paoao pacoo o.
    opopa op oo a opa Yj o ac apa Y1 Y2.
    aoo, ap (Xij) p ac pa cooc c oa oppaa;
    papac o paoao pacoo cpo coo.
    B acoc, ap o pac o cpoa.
    opopa op oo o ap (Xij) o ac X11, X12 21.
  2. tVRAN pcaa zij a caoo ca a t, paopo pacpoo pa (A,B).
    CMRAN ca a ac Zij paopo aco pacp paa (DBLE(A),DBLE(B)) (DIMAG(A),DIMAG(B)).
  3. tMULT op p po apao ap (Xjk) opa N o ap, coaa a opao p ap.
  4. op oo o ap LOC.
    IBM/PC- 16- LOC . , F003, 64 .
:
       IMPLICIT REAL*8 (A-H,O-Z)
       DIMENSION X(3,3),Y(3,3),Z(3,3)
       DATA X/9*7.D0/,Y/9*3.D0/
       CALL DMSUB(3,3,X(1,1),X(1,2),X(2,1),
      *Y(1,1),Y(1,2),Y(2,1),Z(1,1),Z(1,2),Z(2,1))
       PRINT 2,Z
     2 FORMAT(/2X,3F8.2)
       END
:
       4.00 4.00 4.00
       4.00 4.00 4.00
       4.00 4.00 4.00


home up e-mail