srepiso12.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	srepiso12 (off, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, rx, ry, rz, sx, sy, sz, tx, ty, tz, nel)

Function/Subroutine Documentation

srepiso12()

subroutine srepiso12	(	dimension(mvsiz)	off,
		double precision, dimension(mvsiz), intent(in)	x1,
		double precision, dimension(mvsiz), intent(in)	x2,
		double precision, dimension(mvsiz), intent(in)	x3,
		double precision, dimension(mvsiz), intent(in)	x4,
		double precision, dimension(mvsiz), intent(in)	x5,
		double precision, dimension(mvsiz), intent(in)	x6,
		double precision, dimension(mvsiz), intent(in)	x7,
		double precision, dimension(mvsiz), intent(in)	x8,
		double precision, dimension(mvsiz), intent(in)	y1,
		double precision, dimension(mvsiz), intent(in)	y2,
		double precision, dimension(mvsiz), intent(in)	y3,
		double precision, dimension(mvsiz), intent(in)	y4,
		double precision, dimension(mvsiz), intent(in)	y5,
		double precision, dimension(mvsiz), intent(in)	y6,
		double precision, dimension(mvsiz), intent(in)	y7,
		double precision, dimension(mvsiz), intent(in)	y8,
		double precision, dimension(mvsiz), intent(in)	z1,
		double precision, dimension(mvsiz), intent(in)	z2,
		double precision, dimension(mvsiz), intent(in)	z3,
		double precision, dimension(mvsiz), intent(in)	z4,
		double precision, dimension(mvsiz), intent(in)	z5,
		double precision, dimension(mvsiz), intent(in)	z6,
		double precision, dimension(mvsiz), intent(in)	z7,
		double precision, dimension(mvsiz), intent(in)	z8,
		dimension(mvsiz)	rx,
		dimension(mvsiz)	ry,
		dimension(mvsiz)	rz,
		dimension(mvsiz)	sx,
		dimension(mvsiz)	sy,
		dimension(mvsiz)	sz,
		dimension(mvsiz)	tx,
		dimension(mvsiz)	ty,
		dimension(mvsiz)	tz,
		integer, intent(in)	nel )

Definition at line 28 of file srepiso12.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, INTENT(IN) :: NEL
C     REAL
      DOUBLE PRECISION, INTENT(IN),  DIMENSION(MVSIZ) ::
     .   X1, X2, X3, X4, X5, X6, X7, X8,
     .   Y1, Y2, Y3, Y4, Y5, Y6, Y7, Y8,  
     .   Z1, Z2, Z3, Z4, Z5, Z6, Z7, Z8
     
      my_real, INTENT(IN), DIMENSION(MVSIZ) ::
     .   off        
      my_real, INTENT(OUT), DIMENSION(MVSIZ) ::
     .   rx, ry, rz, sx, sy, sz, tx, ty, tz        
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C     REAL
      my_real
     .   x17 , x28 , x35 , x46, 
     .   y17 , y28 , y35 , y46, 
     .   z17 , z28 , z35 , z46,
     .   a17 , a28 , b17 , b28 , 
     .   c17 , c28 
C=======================================================================
      DO i=1,nel
        IF (off(i) <= one ) cycle 
        x17 = x7(i) - x1(i)
        x28 = x8(i) - x2(i)
        x35 = x5(i) - x3(i)
        x46 = x6(i) - x4(i)
        y17 = y7(i) - y1(i)
        y28 = y8(i) - y2(i)
        y35 = y5(i) - y3(i)
        y46 = y6(i) - y4(i)
        z17 = z7(i) - z1(i)
        z28 = z8(i) - z2(i)
        z35 = z5(i) - z3(i)
        z46 = z6(i) - z4(i)
 
        a17 = x17 + x46
        a28 = x28 + x35
        b17 = y17 + y46
        b28 = y28 + y35
        c17 = z17 + z46
        c28 = z28 + z35
C
        rx(i) = x17 + x28 - x35 - x46
        ry(i) = y17 + y28 - y35 - y46
        rz(i) = z17 + z28 - z35 - z46
        sx(i) = a17 + a28
        sy(i) = b17 + b28
        sz(i) = c17 + c28
        tx(i) = a17 - a28
        ty(i) = b17 - b28
        tz(i) = c17 - c28
      ENDDO
C-------------
      RETURN