reduce-historical/r34/xlog/gentran.log

Sat Jun 29 14:11:06 PDT 1991
REDUCE 3.4, 15-Jul-91 ...

1: 1: 
2: 2: 
3: 3: 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;


                      2  2              2               2
M(1,1) :=  - 9*SIN(Q3) *P *M30 - SIN(Q3) *J30Y + SIN(Q3) *J30Z

                                 2           2        2
           + 18*COS(Q3)*COS(Q2)*P *M30 + 18*P *M30 + P *M10 + J30Y

           + J10Y

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;


                                2  2              2
M(2,1) := M(1,2) :=  - 9*SIN(Q3) *P *M30 - SIN(Q3) *J30Y

          2                           2          2
 + SIN(Q3) *J30Z + 9*COS(Q3)*COS(Q2)*P *M30 + 9*P *M30 + J30Y

M(3,1) :=
M(1,3) := -9*SIN(Q3)*SIN(Q2)*M30*P**2;


                                          2
M(3,1) := M(1,3) :=  - 9*SIN(Q3)*SIN(Q2)*P *M30

M(2,2) := -SIN(Q3)**2*J30Y + SIN(Q3)**2*J30Z - 9*SIN(Q3)**2
           *M30*P**2 + J30Y + 9*M30*P**2;


                      2  2              2               2
M(2,2) :=  - 9*SIN(Q3) *P *M30 - SIN(Q3) *J30Y + SIN(Q3) *J30Z

                2
           + 9*P *M30 + J30Y

M(3,2) :=
M(2,3) := 0;


M(3,2) := M(2,3) := 0

M(3,3) := 9*M30*P**2 + J30X;


             2
M(3,3) := 9*P *M30 + J30X


GENTRANLANG!* := 'FORTRAN$


FORTLINELEN!* := 72$



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

C
C     *** COMPUTE VALUES FOR MATRIX M ***
C


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

      M(1,1)=-(9.0*SIN(REAL(Q3))**2*P**2*M30)-(SIN(REAL(Q3))**2*J30Y)+
     . SIN(REAL(Q3))**2*J30Z+18.0*COS(REAL(Q3))*COS(REAL(Q2))*P**2*M30+
     . 18.0*P**2*M30+P**2*M10+J30Y+J10Y
      M(1,2)=-(9.0*SIN(REAL(Q3))**2*P**2*M30)-(SIN(REAL(Q3))**2*J30Y)+
     . SIN(REAL(Q3))**2*J30Z+9.0*COS(REAL(Q3))*COS(REAL(Q2))*P**2*M30+
     . 9.0*P**2*M30+J30Y
      M(1,3)=-(9.0*SIN(REAL(Q3))*SIN(REAL(Q2))*P**2*M30)
      M(2,2)=-(9.0*SIN(REAL(Q3))**2*P**2*M30)-(SIN(REAL(Q3))**2*J30Y)+
     . SIN(REAL(Q3))**2*J30Z+9.0*P**2*M30+J30Y
      M(2,3)=0.0
      M(3,3)=9.0*P**2*M30+J30X


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

C
C     *** COMPUTE VALUES FOR INVERSE MATRIX ***
C


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))
        >>$

      T0=M(1,1)
      T1=M(1,2)
      T2=M(1,3)
      T3=M(2,2)
      T4=M(3,3)


COMMENT -- Contents of Matrix M: --$

M := M;


     [T0  T1  T2]
     [          ]
M := [T1  T3  0 ]
     [          ]
     [T2  0   T4]



MATRIX MXINV(3,3)$


MXINV := M**(-1)$



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

      MXINV(1,1)=(T3*T4)/(T0*T3*T4-(T1**2*T4)-(T2**2*T3))
      MXINV(1,2)=-(T1*T4)/(T0*T3*T4-(T1**2*T4)-(T2**2*T3))
      MXINV(1,3)=-(T2*T3)/(T0*T3*T4-(T1**2*T4)-(T2**2*T3))
      MXINV(2,2)=(T0*T4-T2**2)/(T0*T3*T4-(T1**2*T4)-(T2**2*T3))
      MXINV(2,3)=(T1*T2)/(T0*T3*T4-(T1**2*T4)-(T2**2*T3))
      MXINV(3,3)=(T0*T3-T1**2)/(T0*T3*T4-(T1**2*T4)-(T2**2*T3))


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

      DO 25001 J=1,3
          DO 25002 K=J+1,3
              M(K,J)=M(J,K)
              MXINV(K,J)=MXINV(J,K)
25002     CONTINUE
25001 CONTINUE


END$

4: 4: 
Quitting
Sat Jun 29 14:11:15 PDT 1991