Artifact 6fb920137e0d6a57c29fdbe2614e361e2866072ed41a27c69621c7d9cd3f3b2b:

• Executable file r37/packages/gentran/gentran.tst — part of check-in [f2fda60abd] at 2011-09-02 18:13:33 on branch master — Some historical releases purely for archival purposes

  git-svn-id: https://svn.code.sf.net/p/reduce-algebra/code/trunk/historical@1375 2bfe0521-f11c-4a00-b80e-6202646ff360 (user: arthurcnorman@users.sourceforge.net, size: 1372) [annotate] [blame] [check-ins using] [more...]

• Executable file r38/packages/gentran/gentran.tst — part of check-in [f2fda60abd] at 2011-09-02 18:13:33 on branch master — Some historical releases purely for archival purposes

  git-svn-id: https://svn.code.sf.net/p/reduce-algebra/code/trunk/historical@1375 2bfe0521-f11c-4a00-b80e-6202646ff360 (user: arthurcnorman@users.sourceforge.net, size: 1372) [annotate] [blame] [check-ins using]

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$