sderit3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	sderit3 (off, det, ngl, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, px1h1, px1h2, px1h3, px2h1, px2h2, px2h3, px3h1, px3h2, px3h3, px4h1, px4h2, px4h3, jac1, jac2, jac3, jac4, jac5, jac6, nel, jhbe)

Function/Subroutine Documentation

sderit3()

subroutine sderit3	(		off,
			det,
		integer, dimension(*), intent(in)	ngl,
		double precision, dimension(*)	x1,
		double precision, dimension(*)	x2,
		double precision, dimension(*)	x3,
		double precision, dimension(*)	x4,
		double precision, dimension(*)	x5,
		double precision, dimension(*)	x6,
		double precision, dimension(*)	x7,
		double precision, dimension(*)	x8,
		double precision, dimension(*)	y1,
		double precision, dimension(*)	y2,
		double precision, dimension(*)	y3,
		double precision, dimension(*)	y4,
		double precision, dimension(*)	y5,
		double precision, dimension(*)	y6,
		double precision, dimension(*)	y7,
		double precision, dimension(*)	y8,
		double precision, dimension(*)	z1,
		double precision, dimension(*)	z2,
		double precision, dimension(*)	z3,
		double precision, dimension(*)	z4,
		double precision, dimension(*)	z5,
		double precision, dimension(*)	z6,
		double precision, dimension(*)	z7,
		double precision, dimension(*)	z8,
			px1,
			px2,
			px3,
			px4,
			py1,
			py2,
			py3,
			py4,
			pz1,
			pz2,
			pz3,
			pz4,
			px1h1,
			px1h2,
			px1h3,
			px2h1,
			px2h2,
			px2h3,
			px3h1,
			px3h2,
			px3h3,
			px4h1,
			px4h2,
			px4h3,
			jac1,
			jac2,
			jac3,
			jac4,
			jac5,
			jac6,
		integer, intent(in)	nel,
		integer, intent(in)	jhbe )

Definition at line 30 of file sderit3.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.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-----------------------------------------------
      double precision
     .   x1(*), x2(*), x3(*), x4(*), x5(*), x6(*), x7(*), x8(*),
     .   y1(*), y2(*), y3(*), y4(*), y5(*), y6(*), y7(*), y8(*),  
     .   z1(*), z2(*), z3(*), z4(*), z5(*), z6(*), z7(*), z8(*)   
      INTEGER, INTENT(IN) :: NGL(*), NEL, JHBE 
      my_real
     .   off(*),det(*),
     .   px1(*), px2(*), px3(*), px4(*),  
     .   py1(*), py2(*), py3(*), py4(*),  
     .   pz1(*), pz2(*), pz3(*), pz4(*),  
     .   px1h1(*), px1h2(*), px1h3(*),  
     .   px2h1(*), px2h2(*), px2h3(*),  
     .   px3h1(*), px3h2(*), px3h3(*),  
     .   px4h1(*), px4h2(*), px4h3(*),  
     .   jac1(*), jac2(*), jac3(*),
     .   jac4(*), jac5(*), jac6(*)    
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J
      my_real
     .   dett ,  jac7(mvsiz), jac8(mvsiz) , jac9(mvsiz),
     .   jaci1, jaci2, jaci3,
     .   jaci4, jaci5, jaci6,
     .   jaci7, jaci8, jaci9,
     .   x17 , x28 , x35 , x46,
     .   y17 , y28 , y35 , y46,
     .   z17 , z28 , z35 , z46,
     .   jac_59_68(mvsiz), jac_67_49(mvsiz), jac_48_57(mvsiz),
     .   x_17_46, x_28_35,
     .   y_17_46, y_28_35,
     .   z_17_46, z_28_35,
     .   hx, hy, hz
C
C Jacobian matrix
      DO i=1,nel
        x17=x7(i)-x1(i)
        x28=x8(i)-x2(i)
        x35=x5(i)-x3(i)
        x46=x6(i)-x4(i)
C
        y17=y7(i)-y1(i)             
        y28=y8(i)-y2(i)
        y35=y5(i)-y3(i)
        y46=y6(i)-y4(i)
C
        z17=z7(i)-z1(i)
        z28=z8(i)-z2(i)
        z35=z5(i)-z3(i)
        z46=z6(i)-z4(i)
C            
        jac1(i)=x17+x28-x35-x46
        jac2(i)=y17+y28-y35-y46
        jac3(i)=z17+z28-z35-z46
C
        x_17_46=x17+x46
        x_28_35=x28+x35
        y_17_46=y17+y46
        y_28_35=y28+y35
        z_17_46=z17+z46
        z_28_35=z28+z35
C
        jac4(i)=x_17_46+x_28_35
        jac5(i)=y_17_46+y_28_35
        jac6(i)=z_17_46+z_28_35
C
        jac7(i)=x_17_46-x_28_35
        jac8(i)=y_17_46-y_28_35
        jac9(i)=z_17_46-z_28_35
      ENDDO
C                  
      DO i=1,nel      
        jac_59_68(i)=jac5(i)*jac9(i)-jac6(i)*jac8(i)
        jac_67_49(i)=jac6(i)*jac7(i)-jac4(i)*jac9(i)
        jac_48_57(i)=jac4(i)*jac8(i)-jac5(i)*jac7(i)
      ENDDO
C                  
      DO i=1,nel
        det(i)=one_over_64*(jac1(i)*jac_59_68(i)+jac2(i)*jac_67_49(i)+jac3(i)*jac_48_57(i))
      ENDDO   
C
C Check volume
      CALL schkjab3(
     1   off,     det,     ngl,     nel)
C
C Jacobian matrix inverse
      DO i=1,nel
        dett=one_over_64/det(i)
        jaci1=dett*jac_59_68(i)
        jaci4=dett*jac_67_49(i)
        jaci7=dett*jac_48_57(i)
        jaci2=dett*(jac3(i)*jac8(i)-jac2(i)*jac9(i))
        jaci5=dett*(jac1(i)*jac9(i)-jac3(i)*jac7(i))
        jaci8=dett*(jac2(i)*jac7(i)-jac1(i)*jac8(i))
        jaci3=dett*(jac2(i)*jac6(i)-jac3(i)*jac5(i))
        jaci6=dett*(jac3(i)*jac4(i)-jac1(i)*jac6(i))
        jaci9=dett*(jac1(i)*jac5(i)-jac2(i)*jac4(i))
C
        px1(i)=-jaci1-jaci2-jaci3
        py1(i)=-jaci4-jaci5-jaci6
        pz1(i)=-jaci7-jaci8-jaci9
C
        px2(i)=-jaci1-jaci2+jaci3
        py2(i)=-jaci4-jaci5+jaci6
        pz2(i)=-jaci7-jaci8+jaci9
C
        px3(i)= jaci1-jaci2+jaci3
        py3(i)= jaci4-jaci5+jaci6
        pz3(i)= jaci7-jaci8+jaci9
C
        px4(i)= jaci1-jaci2-jaci3
        py4(i)= jaci4-jaci5-jaci6
        pz4(i)= jaci7-jaci8-jaci9
      ENDDO
C
      IF(jhbe/=0)THEN
        DO i=1,nel
         hx=(x1(i)-x2(i)+x3(i)-x4(i)+x5(i)-x6(i)+x7(i)-x8(i))
         hy=(y1(i)-y2(i)+y3(i)-y4(i)+y5(i)-y6(i)+y7(i)-y8(i))
         hz=(z1(i)-z2(i)+z3(i)-z4(i)+z5(i)-z6(i)+z7(i)-z8(i))
         px1h1(i)=px1(i)*hx+ py1(i)*hy+pz1(i)*hz
         px2h1(i)=px2(i)*hx+ py2(i)*hy+pz2(i)*hz
         px3h1(i)=px3(i)*hx+ py3(i)*hy+pz3(i)*hz
         px4h1(i)=px4(i)*hx+ py4(i)*hy+pz4(i)*hz
        ENDDO
C
        DO i=1,nel
         hx=(x1(i)+x2(i)-x3(i)-x4(i)-x5(i)-x6(i)+x7(i)+x8(i))
         hy=(y1(i)+y2(i)-y3(i)-y4(i)-y5(i)-y6(i)+y7(i)+y8(i))
         hz=(z1(i)+z2(i)-z3(i)-z4(i)-z5(i)-z6(i)+z7(i)+z8(i))
         px1h2(i)=px1(i)*hx+ py1(i)*hy+pz1(i)*hz
         px2h2(i)=px2(i)*hx+ py2(i)*hy+pz2(i)*hz
         px3h2(i)=px3(i)*hx+ py3(i)*hy+pz3(i)*hz
         px4h2(i)=px4(i)*hx+ py4(i)*hy+pz4(i)*hz
        ENDDO
C
        DO i=1,nel
         hx=(x1(i)-x2(i)-x3(i)+x4(i)-x5(i)+x6(i)+x7(i)-x8(i))
         hy=(y1(i)-y2(i)-y3(i)+y4(i)-y5(i)+y6(i)+y7(i)-y8(i))
         hz=(z1(i)-z2(i)-z3(i)+z4(i)-z5(i)+z6(i)+z7(i)-z8(i))
         px1h3(i)=px1(i)*hx+ py1(i)*hy+pz1(i)*hz
         px2h3(i)=px2(i)*hx+ py2(i)*hy+pz2(i)*hz
         px3h3(i)=px3(i)*hx+ py3(i)*hy+pz3(i)*hz
         px4h3(i)=px4(i)*hx+ py4(i)*hy+pz4(i)*hz
        ENDDO
      ENDIF
      RETURN
C-----------