s8ejacit3.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	s8ejacit3 (voldp, ksi, eta, zeta, wi, ajp1, ajp2, ajp3, ajp4, ajp5, ajp6, ajp7, ajp8, ajp9, aji1, aji2, aji3, aji4, aji5, aji6, aji7, aji8, aji9, jac_i, nel, ismstr)

Function/Subroutine Documentation

s8ejacit3()

subroutine s8ejacit3	(	double precision, dimension(mvsiz)	voldp,
			ksi,
			eta,
			zeta,
			wi,
			ajp1,
			ajp2,
			ajp3,
			ajp4,
			ajp5,
			ajp6,
			ajp7,
			ajp8,
			ajp9,
			aji1,
			aji2,
			aji3,
			aji4,
			aji5,
			aji6,
			aji7,
			aji8,
			aji9,
			jac_i,
		integer, intent(in)	nel,
		integer, intent(in)	ismstr )

Definition at line 28 of file s8ejacit3.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   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      INTEGER, INTENT(IN) :: ISMSTR
C     REAL
      my_real
     .   ksi,eta,zeta,wi,
     .   ajp1(*),ajp2(*),ajp3(*),
     .   ajp4(*),ajp5(*),ajp6(*),
     .   ajp7(*),ajp8(*),ajp9(*),
     .   aji1(mvsiz), aji2(mvsiz), aji3(mvsiz),
     .   aji4(mvsiz), aji5(mvsiz), aji6(mvsiz),
     .   aji7(mvsiz), aji8(mvsiz), aji9(mvsiz),
     .   jac_i(10,mvsiz)
      DOUBLE PRECISION 
     .   VOLDP(MVSIZ)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C     REAL
C                                                                     12
      my_real
     .   det(mvsiz) ,dett(mvsiz) ,  
     .   aj1(mvsiz),aj2(mvsiz),aj3(mvsiz),
     .   aj4(mvsiz),aj5(mvsiz),aj6(mvsiz),
     .   aj7(mvsiz),aj8(mvsiz),aj9(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)
      DOUBLE PRECISION 
     .   DETDP
C-----------------------------------------------
      IF (ismstr==11) THEN
       aj1(1:nel)=ajp1(1:nel)
       aj2(1:nel)=ajp2(1:nel)
       aj3(1:nel)=ajp3(1:nel)
       aj4(1:nel)=ajp4(1:nel)
       aj5(1:nel)=ajp5(1:nel)
       aj6(1:nel)=ajp6(1:nel)
       aj7(1:nel)=ajp7(1:nel)
       aj8(1:nel)=ajp8(1:nel)
       aj9(1:nel)=ajp9(1:nel)
 
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)
      ENDDO
C
      DO i=1,nel
      detdp=one_over_512*(aj1(i)*jac_59_68(i)+aj2(i)*jac_67_49(i)+aj3(i)*jac_48_57(i))
      det(i) = detdp
      voldp(i)= wi*detdp
      dett(i)=one_over_512/det(i)
      ENDDO
C
      DO i=1,nel
       aji1(i)=dett(i)*jac_59_68(i)
       aji4(i)=dett(i)*jac_67_49(i)
       aji7(i)=dett(i)*jac_48_57(i)
       aji2(i)=dett(i)*jac_38_29(i)
       aji5(i)=dett(i)*jac_19_37(i)
       aji8(i)=dett(i)*jac_27_18(i)
       aji3(i)=dett(i)*jac_26_35(i)
       aji6(i)=dett(i)*jac_34_16(i)
       aji9(i)=dett(i)*jac_15_24(i)
      ENDDO
C---------don't forget this for Ismstr11 and R,S,T
c      DO I=1,NEL
c       IF (IDEG(I)==0) DELTAX(I)=MIN(DELTAX(I),HUNDRED28*DET(I)*SMAX(I))
c      ENDDO
C------ ISMSTR =10
      ELSE
      DO i=1,nel
       aji1(i)=jac_i(1,i)
       aji4(i)=jac_i(4,i)
       aji7(i)=jac_i(7,i)
       aji2(i)=jac_i(2,i)
       aji5(i)=jac_i(5,i)
       aji8(i)=jac_i(8,i)
       aji3(i)=jac_i(3,i)
       aji6(i)=jac_i(6,i)
       aji9(i)=jac_i(9,i)
       voldp(i) =jac_i(10,i)
      END DO
      END IF !(ISMSTR==11)
C
      RETURN