mrotens.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	mrotens (lft, llt, e1, e2, e3, e4, dyz, e6, r11, r12, r13, r21, r22, r23, r31, r32, r33)
subroutine	mrotensns (lft, llt, dxx, dxy, dxz, dyx, dyy, dyz, dzx, dzy, dzz, r11, r12, r13, r21, r22, r23, r31, r32, r33)

Function/Subroutine Documentation

mrotens()

subroutine mrotens	(	integer	lft,
		integer	llt,
			e1,
			e2,
			e3,
			e4,
			dyz,
			e6,
			r11,
			r12,
			r13,
			r21,
			r22,
			r23,
			r31,
			r32,
			r33 )

Definition at line 32 of file mrotens.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER LFT,LLT
      my_real e1(*),e2(*),e3(*),e4(*),dyz(*),e6(*),
     .   r11(*), r12(*), r13(*), 
     .   r21(*), r22(*), r23(*),
     .   r31(*), r32(*), r33(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real
     .   a11(mvsiz), a12(mvsiz), a13(mvsiz), 
     .   a21(mvsiz), a22(mvsiz), a23(mvsiz),
     .   a31(mvsiz), a32(mvsiz), a33(mvsiz)
C-----------------------------------------------
      DO 200 i=lft,llt
      a11(i)=e1(i)*r11(i)+e4(i)*r21(i)+e6(i)*r31(i)
      a12(i)=e1(i)*r12(i)+e4(i)*r22(i)+e6(i)*r32(i)
      a13(i)=e1(i)*r13(i)+e4(i)*r23(i)+e6(i)*r33(i)
      a21(i)=e4(i)*r11(i)+e2(i)*r21(i)+dyz(i)*r31(i)
      a22(i)=e4(i)*r12(i)+e2(i)*r22(i)+dyz(i)*r32(i)
      a23(i)=e4(i)*r13(i)+e2(i)*r23(i)+dyz(i)*r33(i)
      a31(i)=e6(i)*r11(i)+dyz(i)*r21(i)+e3(i)*r31(i)
      a32(i)=e6(i)*r12(i)+dyz(i)*r22(i)+e3(i)*r32(i)
      a33(i)=e6(i)*r13(i)+dyz(i)*r23(i)+e3(i)*r33(i)
 200  CONTINUE
      DO 201 i=lft,llt
      e1(i)=r11(i)*a11(i)+r21(i)*a21(i)+r31(i)*a31(i)
      e2(i)=r12(i)*a12(i)+r22(i)*a22(i)+r32(i)*a32(i)
      e3(i)=r13(i)*a13(i)+r23(i)*a23(i)+r33(i)*a33(i)
      e4(i)=r11(i)*a12(i)+r21(i)*a22(i)+r31(i)*a32(i)
      dyz(i)=r12(i)*a13(i)+r22(i)*a23(i)+r32(i)*a33(i)
      e6(i)=r11(i)*a13(i)+r21(i)*a23(i)+r31(i)*a33(i)
 201  CONTINUE
C
      RETURN

mrotensns()

subroutine mrotensns	(	integer	lft,
		integer	llt,
			dxx,
			dxy,
			dxz,
			dyx,
			dyy,
			dyz,
			dzx,
			dzy,
			dzz,
			r11,
			r12,
			r13,
			r21,
			r22,
			r23,
			r31,
			r32,
			r33 )

Definition at line 89 of file mrotens.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-----------------------------------------------
      INTEGER LFT,LLT
      my_real
     .   dxy(*) ,dyx(*) ,dyz(*) ,dzy(*) ,dzx(*) ,dxz(*)
      my_real dxx(*),dyy(*),dzz(*),
     .   r11(*), r12(*), r13(*), 
     .   r21(*), r22(*), r23(*),
     .   r31(*), r32(*), r33(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real
     .   a11, a12, a13, 
     .   a21, a22, a23,
     .   a31, a32, a33
C-----------------------------------------------
      DO i=lft,llt
        a11 = dxx(i)*r11(i) + dxy(i)*r21(i) + dxz(i)*r31(i)
        a12 = dxx(i)*r12(i) + dxy(i)*r22(i) + dxz(i)*r32(i)
        a13 = dxx(i)*r13(i) + dxy(i)*r23(i) + dxz(i)*r33(i)
 
        a21 = dyx(i)*r11(i) + dyy(i)*r21(i) + dyz(i)*r31(i)
        a22 = dyx(i)*r12(i) + dyy(i)*r22(i) + dyz(i)*r32(i)
        a23 = dyx(i)*r13(i) + dyy(i)*r23(i) + dyz(i)*r33(i)
 
        a31 = dzx(i)*r11(i) + dzy(i)*r21(i) + dzz(i)*r31(i)
        a32 = dzx(i)*r12(i) + dzy(i)*r22(i) + dzz(i)*r32(i)
        a33 = dzx(i)*r13(i) + dzy(i)*r23(i) + dzz(i)*r33(i)
 
        dxx(i) = r11(i)*a11 + r21(i)*a21 + r31(i)*a31
        dxy(i) = r11(i)*a12 + r21(i)*a22 + r31(i)*a32
        dxz(i) = r11(i)*a13 + r21(i)*a23 + r31(i)*a33
 
        dyx(i) = r12(i)*a11 + r22(i)*a21 + r32(i)*a31
        dyy(i) = r12(i)*a12 + r22(i)*a22 + r32(i)*a32
        dyz(i) = r12(i)*a13 + r22(i)*a23 + r32(i)*a33
 
        dzx(i) = r13(i)*a11 + r23(i)*a21 + r33(i)*a31
        dzy(i) = r13(i)*a12 + r23(i)*a22 + r33(i)*a32
        dzz(i) = r13(i)*a13 + r23(i)*a23 + r33(i)*a33
      ENDDO
C
      RETURN