aglo24.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	aglo24 (nel, sig, r11, r12, r13, r21, r22, r23, r31, r32, r33)

Function/Subroutine Documentation

aglo24()

subroutine aglo24	(	integer	nel,
			sig,
		intent(in)	r11,
		intent(in)	r12,
		intent(in)	r13,
		intent(in)	r21,
		intent(in)	r22,
		intent(in)	r23,
		intent(in)	r31,
		intent(in)	r32,
		intent(in)	r33 )

Definition at line 29 of file aglo24.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 NEL
      my_real, DIMENSION(NEL), INTENT(IN) :: r11,r12,r13,r21,r22,r23,r31,r32,r33
      my_real sig(nel,6)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real, DIMENSION(NEL) :: a11,a12,a13,a21,a22,a23,a31,a32,a33
C=======================================================================
      DO i=1,nel
        a11(i)=sig(i,1)*r11(i)+sig(i,4)*r12(i)+sig(i,6)*r13(i)
        a12(i)=sig(i,1)*r21(i)+sig(i,4)*r22(i)+sig(i,6)*r23(i)
        a13(i)=sig(i,1)*r31(i)+sig(i,4)*r32(i)+sig(i,6)*r33(i)
        a21(i)=sig(i,4)*r11(i)+sig(i,2)*r12(i)+sig(i,5)*r13(i)
        a22(i)=sig(i,4)*r21(i)+sig(i,2)*r22(i)+sig(i,5)*r23(i)
        a23(i)=sig(i,4)*r31(i)+sig(i,2)*r32(i)+sig(i,5)*r33(i)
        a31(i)=sig(i,6)*r11(i)+sig(i,5)*r12(i)+sig(i,3)*r13(i)
        a32(i)=sig(i,6)*r21(i)+sig(i,5)*r22(i)+sig(i,3)*r23(i)
        a33(i)=sig(i,6)*r31(i)+sig(i,5)*r32(i)+sig(i,3)*r33(i)
      ENDDO
c
      DO i=1,nel
        sig(i,1)=r11(i)*a11(i)+r12(i)*a21(i)+r13(i)*a31(i)
        sig(i,2)=r21(i)*a12(i)+r22(i)*a22(i)+r23(i)*a32(i)
        sig(i,3)=r31(i)*a13(i)+r32(i)*a23(i)+r33(i)*a33(i)
        sig(i,4)=r11(i)*a12(i)+r12(i)*a22(i)+r13(i)*a32(i)
        sig(i,5)=r21(i)*a13(i)+r22(i)*a23(i)+r23(i)*a33(i)
        sig(i,6)=r11(i)*a13(i)+r12(i)*a23(i)+r13(i)*a33(i)
      ENDDO
c-----------
      RETURN