s8jac_i.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	s8jac_i (x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, jac_i, nel)

Function/Subroutine Documentation

s8jac_i()

subroutine s8jac_i	(	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,
			jac_i,
		integer	nel )

Definition at line 30 of file s8jac_i.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE message_mod
C------------------------------------------------------
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
      my_real :: jac_i(10,*)
      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(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J     
      my_real
     .   jac1(mvsiz), jac2(mvsiz), jac3(mvsiz), 
     .   jac4(mvsiz), jac5(mvsiz), jac6(mvsiz),
     .   jac7(mvsiz), jac8(mvsiz), jac9(mvsiz)
      my_real
     .   det(mvsiz), dett(mvsiz), 
     .   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),
     .   a17_46(mvsiz),
     .   a28_35(mvsiz),
     .   b17_46(mvsiz),
     .   b28_35(mvsiz),
     .   c17_46(mvsiz),
     .   c28_35(mvsiz)
C=======================================================================
      DO i=1,nel
        x17(i)=x7(i)-x1(i)
        x28(i)=x8(i)-x2(i)
        x35(i)=x5(i)-x3(i)
        x46(i)=x6(i)-x4(i)
C
        y17(i)=y7(i)-y1(i)
        y28(i)=y8(i)-y2(i)
        y35(i)=y5(i)-y3(i)
        y46(i)=y6(i)-y4(i)
C
        z17(i)=z7(i)-z1(i)
        z28(i)=z8(i)-z2(i)
        z35(i)=z5(i)-z3(i)
        z46(i)=z6(i)-z4(i)
      ENDDO
 
      DO i=1,nel
        jac1(i)=x17(i)+x28(i)-x35(i)-x46(i)
        jac2(i)=y17(i)+y28(i)-y35(i)-y46(i)
        jac3(i)=z17(i)+z28(i)-z35(i)-z46(i)
        a17_46(i)=x17(i)+x46(i)
        a28_35(i)=x28(i)+x35(i)
        b17_46(i)=y17(i)+y46(i)
        b28_35(i)=y28(i)+y35(i)
        c17_46(i)=z17(i)+z46(i)
        c28_35(i)=z28(i)+z35(i)
      ENDDO
         
      DO i=1,nel
        jac4(i)=a17_46(i)+a28_35(i)
        jac5(i)=b17_46(i)+b28_35(i)
        jac6(i)=c17_46(i)+c28_35(i)
        jac7(i)=a17_46(i)-a28_35(i)
        jac8(i)=b17_46(i)-b28_35(i)
        jac9(i)=c17_46(i)-c28_35(i)
      ENDDO
 
      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))
       dett(i)=one_over_64/max(det(i),em20)
      ENDDO     
C
C Inverse Jacobian Matrix
      DO i=1,nel
        jac_i(1,i)=dett(i)*jac_59_68(i)
        jac_i(4,i)=dett(i)*jac_67_49(i)
        jac_i(7,i)=dett(i)*jac_48_57(i)
        jac_i(2,i)=dett(i)*(-jac2(i)*jac9(i)+jac3(i)*jac8(i))
        jac_i(5,i)=dett(i)*( jac1(i)*jac9(i)-jac3(i)*jac7(i))
        jac_i(8,i)=dett(i)*(-jac1(i)*jac8(i)+jac2(i)*jac7(i))
        jac_i(3,i)=dett(i)*( jac2(i)*jac6(i)-jac3(i)*jac5(i))
        jac_i(6,i)=dett(i)*(-jac1(i)*jac6(i)+jac3(i)*jac4(i))
        jac_i(9,i)=dett(i)*( jac1(i)*jac5(i)-jac2(i)*jac4(i))
      ENDDO
C
      DO i=1,nel
        jac_i(10,i) = det(i)
      ENDDO
C
      RETURN