szorth3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	storth3 (isorth, nel, g1x, g1y, g1z, g2x, g2y, g2z, g3x, g3y, g3z, gama)
subroutine	szordef3 (nel, dxx, dyy, dzz, d4, d5, d6, g1x, g1y, g1z, g2x, g2y, g2z, g3x, g3y, g3z)
subroutine	sroto3 (nel, sig, sign, g1x, g1y, g1z, g2x, g2y, g2z, g3x, g3y, g3z)

Function/Subroutine Documentation

sroto3()

subroutine sroto3	(	integer	nel,
			sig,
			sign,
			g1x,
			g1y,
			g1z,
			g2x,
			g2y,
			g2z,
			g3x,
			g3y,
			g3z )

Definition at line 134 of file szorth3.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
C
      my_real
     .   sig(6,*),sign(6,*), g1x(*),g1y(*),g1z(*),
     .   g2x(*),g2y(*),g2z(*),g3x(*),g3y(*),g3z(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C
      my_real
     .   sx,sy,sz
C-----------------------------------------------
      DO i=1,nel
         sx = sig(1,i)*g1x(i)+sig(4,i)*g1y(i)+sig(6,i)*g1z(i)
         sy = sig(4,i)*g1x(i)+sig(2,i)*g1y(i)+sig(5,i)*g1z(i)
         sz = sig(6,i)*g1x(i)+sig(5,i)*g1y(i)+sig(3,i)*g1z(i)
         sign(1,i) = sx*g1x(i)+sy*g1y(i)+sz*g1z(i)
         sign(4,i) = sx*g2x(i)+sy*g2y(i)+sz*g2z(i)
         sign(6,i) = sx*g3x(i)+sy*g3y(i)+sz*g3z(i)
         sx = sig(1,i)*g2x(i)+sig(4,i)*g2y(i)+sig(6,i)*g2z(i)
         sy = sig(4,i)*g2x(i)+sig(2,i)*g2y(i)+sig(5,i)*g2z(i)
         sz = sig(6,i)*g2x(i)+sig(5,i)*g2y(i)+sig(3,i)*g2z(i)
         sign(2,i) = sx*g2x(i)+sy*g2y(i)+sz*g2z(i)
         sign(5,i) = sx*g3x(i)+sy*g3y(i)+sz*g3z(i)
         sx = sig(1,i)*g3x(i)+sig(4,i)*g3y(i)+sig(6,i)*g3z(i)
         sy = sig(4,i)*g3x(i)+sig(2,i)*g3y(i)+sig(5,i)*g3z(i)
         sz = sig(6,i)*g3x(i)+sig(5,i)*g3y(i)+sig(3,i)*g3z(i)
         sign(3,i) = sx*g3x(i)+sy*g3y(i)+sz*g3z(i)
      ENDDO
C
      RETURN

storth3()

subroutine storth3	(	integer	isorth,
		integer	nel,
			g1x,
			g1y,
			g1z,
			g2x,
			g2y,
			g2z,
			g3x,
			g3y,
			g3z,
			gama )

Definition at line 28 of file szorth3.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 :: ISORTH, NEL
C
      my_real
     .   g1x(*),g1y(*),g1z(*),
     .   g2x(*),g2y(*),g2z(*),g3x(*),g3y(*),g3z(*),gama(nel,6)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C-----------------------------------------------
      IF (isorth == 1) THEN 
        DO i=1,nel
         g1x(i)=gama(i,1)
         g1y(i)=gama(i,2)
         g1z(i)=gama(i,3)
         g2x(i)=gama(i,4)
         g2y(i)=gama(i,5)
         g2z(i)=gama(i,6)
         g3x(i)=g1y(i)*g2z(i)-g1z(i)*g2y(i)
         g3y(i)=g1z(i)*g2x(i)-g1x(i)*g2z(i)
         g3z(i)=g1x(i)*g2y(i)-g1y(i)*g2x(i)
        END DO
      ELSE
        DO i=1,nel
         g1x(i)=one
         g1y(i)=zero
         g1z(i)=zero
         g2x(i)=zero
         g2y(i)=one
         g2z(i)=zero
         g3x(i)=zero
         g3y(i)=zero
         g3z(i)=one
        END DO
      ENDIF 
C
      RETURN

szordef3()

subroutine szordef3	(	integer	nel,
			dxx,
			dyy,
			dzz,
			d4,
			d5,
			d6,
			g1x,
			g1y,
			g1z,
			g2x,
			g2y,
			g2z,
			g3x,
			g3y,
			g3z )

Definition at line 82 of file szorth3.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 :: NEL
C
      my_real
     .   dxx(*),dyy(*),dzz(*),d4(*),d5(*),d6(*), g1x(*),g1y(*),g1z(*),
     .   g2x(*),g2y(*),g2z(*),g3x(*),g3y(*),g3z(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C
      my_real
     .   def(6,mvsiz),defn(6,mvsiz)
C-----------------------------------------------
      DO i=1,nel
        def(1,i) = dxx(i)
        def(2,i) = dyy(i)
        def(3,i) = dzz(i)
        def(4,i) = half*d4(i)
        def(5,i) = half*d5(i)
        def(6,i) = half*d6(i)
      ENDDO
      CALL sroto3(nel,def,defn,
     .            g1x,g1y,g1z,g2x,g2y,g2z,g3x,g3y,g3z)
      DO i=1,nel
        dxx(i) = defn(1,i)
        dyy(i) = defn(2,i)
        dzz(i) = defn(3,i)
        d4(i) = two*defn(4,i)
        d5(i) = two*defn(5,i)
        d6(i) = two*defn(6,i)
      ENDDO
C
      RETURN