s8zjac_ic.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	s8zjac_ic (xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8, aj1, aj2, aj3, aj4, aj5, aj6, aj7, aj8, aj9, hx, hy, hz, pxc1, pxc2, pxc3, pxc4, pyc1, pyc2, pyc3, pyc4, pzc1, pzc2, pzc3, pzc4, jac_i, index, nel)

Function/Subroutine Documentation

s8zjac_ic()

subroutine s8zjac_ic	(	double precision, dimension(mvsiz)	xd1,
		double precision, dimension(mvsiz)	xd2,
		double precision, dimension(mvsiz)	xd3,
		double precision, dimension(mvsiz)	xd4,
		double precision, dimension(mvsiz)	xd5,
		double precision, dimension(mvsiz)	xd6,
		double precision, dimension(mvsiz)	xd7,
		double precision, dimension(mvsiz)	xd8,
		double precision, dimension(mvsiz)	yd1,
		double precision, dimension(mvsiz)	yd2,
		double precision, dimension(mvsiz)	yd3,
		double precision, dimension(mvsiz)	yd4,
		double precision, dimension(mvsiz)	yd5,
		double precision, dimension(mvsiz)	yd6,
		double precision, dimension(mvsiz)	yd7,
		double precision, dimension(mvsiz)	yd8,
		double precision, dimension(mvsiz)	zd1,
		double precision, dimension(mvsiz)	zd2,
		double precision, dimension(mvsiz)	zd3,
		double precision, dimension(mvsiz)	zd4,
		double precision, dimension(mvsiz)	zd5,
		double precision, dimension(mvsiz)	zd6,
		double precision, dimension(mvsiz)	zd7,
		double precision, dimension(mvsiz)	zd8,
			aj1,
			aj2,
			aj3,
			aj4,
			aj5,
			aj6,
			aj7,
			aj8,
			aj9,
			hx,
			hy,
			hz,
			pxc1,
			pxc2,
			pxc3,
			pxc4,
			pyc1,
			pyc2,
			pyc3,
			pyc4,
			pzc1,
			pzc2,
			pzc3,
			pzc4,
			jac_i,
		integer, dimension(*)	index,
		integer	nel )

Definition at line 28 of file s8zjac_ic.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,INDEX(*)
      my_real
     .   hx(4,*), hy(4,*), hz(4,*),  
     .   aj1(*),aj2(*),aj3(*),
     .   aj4(*),aj5(*),aj6(*),
     .   aj7(*),aj8(*),aj9(*),jac_i(10,*),
     .   pxc1(*), pxc2(*), pxc3(*), pxc4(*), 
     .   pyc1(*), pyc2(*), pyc3(*), pyc4(*), 
     .   pzc1(*), pzc2(*), pzc3(*), pzc4(*) 
      double precision
     .   xd1(mvsiz), xd2(mvsiz), xd3(mvsiz), xd4(mvsiz), 
     .   xd5(mvsiz), xd6(mvsiz), xd7(mvsiz), xd8(mvsiz), 
     .   yd1(mvsiz), yd2(mvsiz), yd3(mvsiz), yd4(mvsiz),
     .   yd5(mvsiz), yd6(mvsiz), yd7(mvsiz), yd8(mvsiz),
     .   zd1(mvsiz), zd2(mvsiz), zd3(mvsiz), zd4(mvsiz),
     .   zd5(mvsiz), zd6(mvsiz), zd7(mvsiz), zd8(mvsiz)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,  J
      
      my_real
     .   det(mvsiz) ,dett(mvsiz)
C                                                                     12
      my_real
     .   x17(mvsiz) , x28(mvsiz) , x35(mvsiz) , x46(mvsiz),
     .   y17(mvsiz) , y28(mvsiz) , y35(mvsiz) , y46(mvsiz),
     .   z17(mvsiz) , z28(mvsiz) , z35(mvsiz) , z46(mvsiz),
     .   jac_59_68(mvsiz), jac_67_49(mvsiz), jac_48_57(mvsiz),
     .   jac_38_29(mvsiz), jac_19_37(mvsiz), jac_27_18(mvsiz),
     .   jac_26_35(mvsiz), jac_34_16(mvsiz), jac_15_24(mvsiz),
     .   aji1, aji2, aji3,aji4, aji5, aji6,aji7, aji8, aji9,
     .   aj12, aj45, aj78,aj12p, aj45p, aj78p,
     .   a17(mvsiz) , a28(mvsiz) ,
     .   b17(mvsiz) , b28(mvsiz) ,
     .   c17(mvsiz) , c28(mvsiz) ,jac_1(10,mvsiz)
C-----------------------------------------------
      DO i=1,nel
      x17(i)=xd7(i)-xd1(i)
      x28(i)=xd8(i)-xd2(i)
      x35(i)=xd5(i)-xd3(i)
      x46(i)=xd6(i)-xd4(i)
      y17(i)=yd7(i)-yd1(i)
      y28(i)=yd8(i)-yd2(i)
      y35(i)=yd5(i)-yd3(i)
      y46(i)=yd6(i)-yd4(i)
      z17(i)=zd7(i)-zd1(i)
      z28(i)=zd8(i)-zd2(i)
      z35(i)=zd5(i)-zd3(i)
      z46(i)=zd6(i)-zd4(i)
      END DO
C
      DO i=1,nel
      aj4(i)=x17(i)+x28(i)-x35(i)-x46(i)
      aj5(i)=y17(i)+y28(i)-y35(i)-y46(i)
      aj6(i)=z17(i)+z28(i)-z35(i)-z46(i)
      a17(i)=x17(i)+x46(i)
      a28(i)=x28(i)+x35(i)
      b17(i)=y17(i)+y46(i)
      b28(i)=y28(i)+y35(i)
      c17(i)=z17(i)+z46(i)
      c28(i)=z28(i)+z35(i)
C
      aj7(i)=a17(i)+a28(i)
      aj8(i)=b17(i)+b28(i)
      aj9(i)=c17(i)+c28(i)
      aj1(i)=a17(i)-a28(i)
      aj2(i)=b17(i)-b28(i)
      aj3(i)=c17(i)-c28(i)
      END DO
C
C     JACOBIAN
C
      DO i=1,nel
      jac_59_68(i)=aj5(i)*aj9(i)-aj6(i)*aj8(i)
      jac_67_49(i)=aj6(i)*aj7(i)-aj4(i)*aj9(i)
      jac_38_29(i)=(-aj2(i)*aj9(i)+aj3(i)*aj8(i))
      jac_19_37(i)=( aj1(i)*aj9(i)-aj3(i)*aj7(i))
      jac_27_18(i)=(-aj1(i)*aj8(i)+aj2(i)*aj7(i))
      jac_26_35(i)=( aj2(i)*aj6(i)-aj3(i)*aj5(i))
      jac_34_16(i)=(-aj1(i)*aj6(i)+aj3(i)*aj4(i))
      jac_15_24(i)=( aj1(i)*aj5(i)-aj2(i)*aj4(i))
      jac_48_57(i)=aj4(i)*aj8(i)-aj5(i)*aj7(i)
C
      det(i)=one_over_64*(aj1(i)*jac_59_68(i)+aj2(i)*jac_67_49(i)+aj3(i)*jac_48_57(i))
      dett(i)=one_over_64/det(i)
      ENDDO
C
       DO i=1,nel
         aji1=dett(i)*jac_59_68(i)
         aji4=dett(i)*jac_67_49(i)
         aji7=dett(i)*jac_48_57(i)
         aji2=dett(i)*jac_38_29(i)
         aji5=dett(i)*jac_19_37(i)
         aji8=dett(i)*jac_27_18(i)
         aji3=dett(i)*jac_26_35(i)
         aji6=dett(i)*jac_34_16(i)
         aji9=dett(i)*jac_15_24(i)
         aj12=aji1-aji2
         aj45=aji4-aji5
         aj78=aji7-aji8
         pxc2(i)= aj12-aji3
         pyc2(i)= aj45-aji6
         pzc2(i)= aj78-aji9
         pxc4(i)=-aj12-aji3
         pyc4(i)=-aj45-aji6
         pzc4(i)=-aj78-aji9
         aj12p=aji1+aji2
         aj45p=aji4+aji5
         aj78p=aji7+aji8
         pxc1(i)=-aj12p-aji3
         pyc1(i)=-aj45p-aji6
         pzc1(i)=-aj78p-aji9
         pxc3(i)=aj12p-aji3
         pyc3(i)=aj45p-aji6
         pzc3(i)=aj78p-aji9
         jac_1(1,i)=aji1
         jac_1(4,i)=aji4
         jac_1(7,i)=aji7
         jac_1(2,i)=aji2
         jac_1(5,i)=aji5
         jac_1(8,i)=aji8
         jac_1(3,i)=aji3
         jac_1(6,i)=aji6
         jac_1(9,i)=aji9
         jac_1(10,i)=det(i)
       ENDDO
C mode 1
C 1 1 -1 -1 -1 -1 1 1
       DO i=1,nel
         hx(1,i)=(xd1(i)+xd2(i)-xd3(i)-xd4(i)-xd5(i)-xd6(i)+xd7(i)+xd8(i))
         hy(1,i)=(yd1(i)+yd2(i)-yd3(i)-yd4(i)-yd5(i)-yd6(i)+yd7(i)+yd8(i))
         hz(1,i)=(zd1(i)+zd2(i)-zd3(i)-zd4(i)-zd5(i)-zd6(i)+zd7(i)+zd8(i))
       ENDDO
C mode 2
C 1 -1 -1 1 -1 1 1 -1
       DO i=1,nel
         hx(2,i)=(xd1(i)-xd2(i)-xd3(i)+xd4(i)-xd5(i)+xd6(i)+xd7(i)-xd8(i))
         hy(2,i)=(yd1(i)-yd2(i)-yd3(i)+yd4(i)-yd5(i)+yd6(i)+yd7(i)-yd8(i))
         hz(2,i)=(zd1(i)-zd2(i)-zd3(i)+zd4(i)-zd5(i)+zd6(i)+zd7(i)-zd8(i))
       ENDDO
C mode 3
C 1 -1 1 -1 1 -1 1 -1
       DO i=1,nel
         hx(3,i)=(xd1(i)-xd2(i)+xd3(i)-xd4(i)+xd5(i)-xd6(i)+xd7(i)-xd8(i))
         hy(3,i)=(yd1(i)-yd2(i)+yd3(i)-yd4(i)+yd5(i)-yd6(i)+yd7(i)-yd8(i))
         hz(3,i)=(zd1(i)-zd2(i)+zd3(i)-zd4(i)+zd5(i)-zd6(i)+zd7(i)-zd8(i))
       ENDDO
C mode 4
C -1 1 -1 1 1 -1 1 -1
       DO i=1,nel
         hx(4,i)=(-xd1(i)+xd2(i)-xd3(i)+xd4(i)+xd5(i)-xd6(i)+xd7(i)-xd8(i))
         hy(4,i)=(-yd1(i)+yd2(i)-yd3(i)+yd4(i)+yd5(i)-yd6(i)+yd7(i)-yd8(i))
         hz(4,i)=(-zd1(i)+zd2(i)-zd3(i)+zd4(i)+zd5(i)-zd6(i)+zd7(i)-zd8(i))
       ENDDO
#include "vectorize.inc"
        DO j=1,nel
          i = index(j)
          jac_i(1:10,i) = jac_1(1:10,j)
        ENDDO      
C       
      RETURN
C