s4rota3.F File Reference

#include "implicit_f.inc"
#include "vect01_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	s4rota3 (r11, r12, r13, r21, r22, r23, r31, r32, r33, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4)
subroutine	s4rotat3 (r11, r12, r13, r21, r22, r23, r31, r32, r33, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4)

Function/Subroutine Documentation

s4rota3()

subroutine s4rota3	(		r11,
			r12,
			r13,
			r21,
			r22,
			r23,
			r31,
			r32,
			r33,
		double precision, dimension(*)	x1,
		double precision, dimension(*)	x2,
		double precision, dimension(*)	x3,
		double precision, dimension(*)	x4,
		double precision, dimension(*)	y1,
		double precision, dimension(*)	y2,
		double precision, dimension(*)	y3,
		double precision, dimension(*)	y4,
		double precision, dimension(*)	z1,
		double precision, dimension(*)	z2,
		double precision, dimension(*)	z3,
		double precision, dimension(*)	z4 )

Definition at line 28 of file s4rota3.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
C     REAL
      double precision
     .   x1(*), x2(*), x3(*), x4(*),
     .   y1(*), y2(*), y3(*), y4(*),
     .   z1(*), z2(*), z3(*), z4(*)     
      my_real
     .   r11(*), r12(*), r13(*),
     .   r21(*), r22(*), r23(*),
     .   r31(*), r32(*), r33(*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "vect01_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C     REAL
      my_real
     .   x,y,z
C-----------------------------------------------
C     X=RX' <=> X'=t(R)X chgt de base.
      DO i=lft,llt
        x=r11(i)*x1(i)+r21(i)*y1(i)+r31(i)*z1(i)
        y=r12(i)*x1(i)+r22(i)*y1(i)+r32(i)*z1(i)
        z=r13(i)*x1(i)+r23(i)*y1(i)+r33(i)*z1(i)
        x1(i)=x
        y1(i)=y
        z1(i)=z
        x=r11(i)*x2(i)+r21(i)*y2(i)+r31(i)*z2(i)
        y=r12(i)*x2(i)+r22(i)*y2(i)+r32(i)*z2(i)
        z=r13(i)*x2(i)+r23(i)*y2(i)+r33(i)*z2(i)
        x2(i)=x
        y2(i)=y
        z2(i)=z
        x=r11(i)*x3(i)+r21(i)*y3(i)+r31(i)*z3(i)
        y=r12(i)*x3(i)+r22(i)*y3(i)+r32(i)*z3(i)
        z=r13(i)*x3(i)+r23(i)*y3(i)+r33(i)*z3(i)
        x3(i)=x
        y3(i)=y
        z3(i)=z
        x=r11(i)*x4(i)+r21(i)*y4(i)+r31(i)*z4(i)
        y=r12(i)*x4(i)+r22(i)*y4(i)+r32(i)*z4(i)
        z=r13(i)*x4(i)+r23(i)*y4(i)+r33(i)*z4(i)
        x4(i)=x
        y4(i)=y
        z4(i)=z
      ENDDO
C-----------
      RETURN

s4rotat3()

subroutine s4rotat3	(		r11,
			r12,
			r13,
			r21,
			r22,
			r23,
			r31,
			r32,
			r33,
			x1,
			x2,
			x3,
			x4,
			y1,
			y2,
			y3,
			y4,
			z1,
			z2,
			z3,
			z4 )

Definition at line 95 of file s4rota3.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
C     REAL
      my_real
     .   x1(*), x2(*), x3(*), x4(*),
     .   y1(*), y2(*), y3(*), y4(*),
     .   z1(*), z2(*), z3(*), z4(*),     
     .   r11(*), r12(*), r13(*),
     .   r21(*), r22(*), r23(*),
     .   r31(*), r32(*), r33(*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "vect01_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C     REAL
      my_real
     .   x,y,z
C-----------------------------------------------
C     X=RX' <=> X'=t(R)X chgt de base.
      DO i=lft,llt
        x=r11(i)*x1(i)+r21(i)*y1(i)+r31(i)*z1(i)
        y=r12(i)*x1(i)+r22(i)*y1(i)+r32(i)*z1(i)
        z=r13(i)*x1(i)+r23(i)*y1(i)+r33(i)*z1(i)
        x1(i)=x
        y1(i)=y
        z1(i)=z
        x=r11(i)*x2(i)+r21(i)*y2(i)+r31(i)*z2(i)
        y=r12(i)*x2(i)+r22(i)*y2(i)+r32(i)*z2(i)
        z=r13(i)*x2(i)+r23(i)*y2(i)+r33(i)*z2(i)
        x2(i)=x
        y2(i)=y
        z2(i)=z
        x=r11(i)*x3(i)+r21(i)*y3(i)+r31(i)*z3(i)
        y=r12(i)*x3(i)+r22(i)*y3(i)+r32(i)*z3(i)
        z=r13(i)*x3(i)+r23(i)*y3(i)+r33(i)*z3(i)
        x3(i)=x
        y3(i)=y
        z3(i)=z
        x=r11(i)*x4(i)+r21(i)*y4(i)+r31(i)*z4(i)
        y=r12(i)*x4(i)+r22(i)*y4(i)+r32(i)*z4(i)
        z=r13(i)*x4(i)+r23(i)*y4(i)+r33(i)*z4(i)
        x4(i)=x
        y4(i)=y
        z4(i)=z
      ENDDO
C-----------
      RETURN