File r34.1/xmpl/gentran.tst artifact 4d4adb65e1 part of check-in b5833487d7


MATRIX M(3,3)$
M(1,1) := 18*COS(Q3)*COS(Q2)*M30*P**2 - 9*SIN(Q3)**2*P**2*M30
          - SIN(Q3)**2*J30Y + SIN(Q3)**2*J30Z + P**2*M10
          + 18*P**2*M30 + J10Y + J30Y;
M(2,1) :=
M(1,2) := 9*COS(Q3)*COS(Q2)*M30*P**2 - SIN(Q3)**2*J30Y +
          SIN(Q3)**2*J30Z - 9*SIN(Q3)**2*M30*P**2 + J30Y +
          9*M30*P**2;
M(3,1) :=
M(1,3) := -9*SIN(Q3)*SIN(Q2)*M30*P**2;
M(2,2) := -SIN(Q3)**2*J30Y + SIN(Q3)**2*J30Z - 9*SIN(Q3)**2
           *M30*P**2 + J30Y + 9*M30*P**2;
M(3,2) :=
M(2,3) := 0;
M(3,3) := 9*M30*P**2 + J30X;

GENTRANLANG!* := 'FORTRAN$
FORTLINELEN!* := 72$

GENTRAN LITERAL "C", CR!*,
                "C", TAB!*, "*** COMPUTE VALUES FOR MATRIX M ***", CR!*,
                "C", CR!*$

FOR j:=1:3 DO
    FOR k:=j:3 DO
         GENTRAN M(j,k) ::=: M(j,k)$

GENTRAN LITERAL "C", CR!*,
                "C", TAB!*, "*** COMPUTE VALUES FOR INVERSE MATRIX ***",
                     CR!*,
                "C", CR!*$

SHARE var$
FOR j:=1:3 DO
    FOR k:=j:3 DO
        IF M(j,k) NEQ 0 THEN
        <<
            var := TEMPVAR NIL;
            MARKVAR var;
            M(j,k) := var;
            M(k,j) := var;
            GENTRAN
                EVAL(var) := M(EVAL(j),EVAL(k))
        >>$

COMMENT -- Contents of Matrix M: --$
M := M;

MATRIX MXINV(3,3)$
MXINV := M**(-1)$

FOR j:=1:3 DO
    FOR k:=j:3 DO
        GENTRAN MXINV(j,k) ::=: MXINV(j,k)$

GENTRAN
   for j:=1:3 do
       for k:=j+1:3 do
       <<
           m(k,j) := m(j,k);
           mxinv(k,j) := mxinv(j,k)
        >>$

END$


REDUCE Historical
REDUCE Sourceforge Project | Historical SVN Repository | GitHub Mirror | SourceHut Mirror | NotABug Mirror | Chisel Mirror | Chisel RSS ]