sjacidp.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	sjacidp (xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8, jac_i, nel)

Function/Subroutine Documentation

sjacidp()

subroutine sjacidp	(	double precision, dimension(mvsiz), intent(in)	xd1,
		double precision, dimension(mvsiz), intent(in)	xd2,
		double precision, dimension(mvsiz), intent(in)	xd3,
		double precision, dimension(mvsiz), intent(in)	xd4,
		double precision, dimension(mvsiz), intent(in)	xd5,
		double precision, dimension(mvsiz), intent(in)	xd6,
		double precision, dimension(mvsiz), intent(in)	xd7,
		double precision, dimension(mvsiz), intent(in)	xd8,
		double precision, dimension(mvsiz), intent(in)	yd1,
		double precision, dimension(mvsiz), intent(in)	yd2,
		double precision, dimension(mvsiz), intent(in)	yd3,
		double precision, dimension(mvsiz), intent(in)	yd4,
		double precision, dimension(mvsiz), intent(in)	yd5,
		double precision, dimension(mvsiz), intent(in)	yd6,
		double precision, dimension(mvsiz), intent(in)	yd7,
		double precision, dimension(mvsiz), intent(in)	yd8,
		double precision, dimension(mvsiz), intent(in)	zd1,
		double precision, dimension(mvsiz), intent(in)	zd2,
		double precision, dimension(mvsiz), intent(in)	zd3,
		double precision, dimension(mvsiz), intent(in)	zd4,
		double precision, dimension(mvsiz), intent(in)	zd5,
		double precision, dimension(mvsiz), intent(in)	zd6,
		double precision, dimension(mvsiz), intent(in)	zd7,
		double precision, dimension(mvsiz), intent(in)	zd8,
		intent(out)	jac_i,
		integer, intent(in)	nel )

Definition at line 30 of file sjacidp.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE message_mod
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
      my_real , DIMENSION(10,NEL),INTENT(OUT) :: jac_i     
      DOUBLE PRECISION , DIMENSION(MVSIZ), INTENT(IN) ::
     .   XD1, XD2, XD3, XD4, XD5, XD6, XD7, XD8,
     .   YD1, YD2, YD3, YD4, YD5, YD6, YD7, YD8,
     .   ZD1, ZD2, ZD3, ZD4, ZD5, ZD6, ZD7, ZD8
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: I
C
      my_real
     .   jac1(mvsiz), jac2(mvsiz), jac3(mvsiz), 
     .   jac4(mvsiz), jac5(mvsiz), jac6(mvsiz),
     .   jac7(mvsiz), jac8(mvsiz), jac9(mvsiz),
     .   jac_59_68(mvsiz), jac_67_49(mvsiz), jac_48_57(mvsiz)
C
      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)
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)
C
        y17(i)=yd7(i)-yd1(i)
        y28(i)=yd8(i)-yd2(i)
        y35(i)=yd5(i)-yd3(i)
        y46(i)=yd6(i)-yd4(i)
C
        z17(i)=zd7(i)-zd1(i)
        z28(i)=zd8(i)-zd2(i)
        z35(i)=zd5(i)-zd3(i)
        z46(i)=zd6(i)-zd4(i)
      ENDDO
C
C Jacobian Matrix
      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)
C
        jac4(i)=x17(i)+x46(i)+x28(i)+x35(i)
        jac5(i)=y17(i)+y46(i)+y28(i)+y35(i)
        jac6(i)=z17(i)+z46(i)+z28(i)+z35(i)
C
        jac7(i)=x17(i)+x46(i)-x28(i)-x35(i)
        jac8(i)=y17(i)+y46(i)-y28(i)-y35(i)
        jac9(i)=z17(i)+z46(i)-z28(i)-z35(i)
      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))
       dett(i)=one_over_64/max(det(i),em20)
      ENDDO     
C
C Jacobian matrix inverse
      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