next up previous contents
Next: Graphics Up: Code generation Previous: Code generation

Generated FORTRAN program

Generated program

C
C --- Calculate Matrix Values ---
C
      M(1,1)=-(9*M30*P**2*SIN(Q3)**2)+J30Z*SIN(Q3)**2-(J30Y*SIN(Q3)**2)+
     . 18*M30*P**2*COS(Q2)*COS(Q3)+18*M30*P**2+M10*P**2+J30Y+J10Y
      M(1,2)=-(9*M30*P**2*SIN(Q3)**2)+J30Z*SIN(Q3)**2-(J30Y*SIN(Q3)**2)+
     . 9*M30*P**2*COS(Q2)*COS(Q3)+9*M30*P**2+J30Y
      M(1,3)=-(9*M30*P**2*SIN(Q2)*SIN(Q3))
      M(2,2)=-(9*M30*P**2*SIN(Q3)**2)+J30Z*SIN(Q3)**2-(J30Y*SIN(Q3)**2)+
     . 9*M30*P**2+J30Y
      M(2,3)=0
      M(3,3)=9*M30*P**2+J30X
C
C --- Assign Non-Zero Matrix Values to Temporary Variables ---
C
      T0=M(1,1)
      T1=M(1,2)
      T2=M(1,3)
      T3=M(2,2)
      T4=M(3,3)
C
C --- Calculate Inverse Matrix Values ---
C
      MINV(1,1)=(T3*T4)/((T0*T3-T1**2)*T4-(T2**2*T3))
      MINV(1,2)=-((T1*T4)/((T0*T3-T1**2)*T4-(T2**2*T3)))
      MINV(1,3)=-((T2*T3)/((T0*T3-T1**2)*T4-(T2**2*T3)))
      MINV(2,2)=(T0*T4-T2**2)/((T0*T3-T1**2)*T4-(T2**2*T3))
      MINV(2,3)=(T1*T2)/((T0*T3-T1**2)*T4-(T2**2*T3))
      MINV(3,3)=(T0*T3-T1**2)/((T0*T3-T1**2)*T4-(T2**2*T3))
C
C --- Copy Entries Across Main Diagonals ---
C
      DO 25001 I=1,3
          DO 25002 J=I+1,3
              M(J,I)=M(I,J)
              MINV(J,I)=MINV(I,J)
25002     CONTINUE
25001 CONTINUE



Richard Liska