s10jacobt.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	s10jacobt (alph, beta, w, px1, px2, px3, px4, px5, px6, px7, px8, px9, px10, py1, py2, py3, py4, py5, py6, py7, py8, py9, py10, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, pz9, pz10, nx1, nx2, nx3, nx4, nx5, nx6, nx7, nx8, nx9, nx10, vol, jac_i, nel)

Function/Subroutine Documentation

s10jacobt()

subroutine s10jacobt	(		alph,
			beta,
			w,
			px1,
			px2,
			px3,
			px4,
			px5,
			px6,
			px7,
			px8,
			px9,
			px10,
			py1,
			py2,
			py3,
			py4,
			py5,
			py6,
			py7,
			py8,
			py9,
			py10,
			pz1,
			pz2,
			pz3,
			pz4,
			pz5,
			pz6,
			pz7,
			pz8,
			pz9,
			pz10,
			nx1,
			nx2,
			nx3,
			nx4,
			nx5,
			nx6,
			nx7,
			nx8,
			nx9,
			nx10,
			vol,
			jac_i,
		integer, intent(in)	nel )

Definition at line 26 of file s10jacobt.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, INTENT(IN) :: NEL
       my_real
     .   px1(mvsiz),px2(mvsiz),px3(mvsiz),px4(mvsiz),px5(mvsiz),
     .   px6(mvsiz),px7(mvsiz),px8(mvsiz),px9(mvsiz),px10(mvsiz),
     .   py1(mvsiz),py2(mvsiz),py3(mvsiz),py4(mvsiz),py5(mvsiz),
     .   py6(mvsiz),py7(mvsiz),py8(mvsiz),py9(mvsiz),py10(mvsiz),
     .   pz1(mvsiz),pz2(mvsiz),pz3(mvsiz),pz4(mvsiz),pz5(mvsiz),
     .   pz6(mvsiz),pz7(mvsiz),pz8(mvsiz),pz9(mvsiz),pz10(mvsiz),
     .   nx1(mvsiz),nx2(mvsiz),nx3(mvsiz),nx4(mvsiz),nx5(mvsiz),
     .   nx6(mvsiz),nx7(mvsiz),nx8(mvsiz),nx9(mvsiz),nx10(mvsiz),
     .   vol(mvsiz),alph,beta,w,jac_i(10,mvsiz)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real
     .   d,a4mb4,a4,b4,bb,ab,b2,
     .   dxdr,dxds,dxdt,dydr,dyds,dydt,dzdr,dzds,dzdt
     
      DOUBLE PRECISION 
     .   AA,A4M1,B4M1
      
      my_real
     .   drdx, dsdx, dtdx,
     .   drdy, dsdy, dtdy,
     .   drdz, dsdz, dtdz
C-----------------------------------------------
c      AA = (TWO*ALPH - ONE)*ALPH
c      BB = (TWO*BETA - ONE)*BETA
c      B2 = FOUR*BETA*BETA
c      AB = FOUR*ALPH*BETA      
C-----------------------------------------------
C
      a4 = four * alph
      b4 = four * beta
      a4m1  = a4 - one
      b4m1  = b4 - one
      a4mb4  = a4 - b4
      DO i=1,nel
       drdx=jac_i(1,i)
       dsdx=jac_i(4,i)
       dtdx=jac_i(7,i)
C
       drdy=jac_i(2,i)
       dsdy=jac_i(5,i)
       dtdy=jac_i(8,i)
C
       drdz=jac_i(3,i)
       dsdz=jac_i(6,i)
       dtdz=jac_i(9,i)
C
       vol(i) = jac_i(10,i)
C
        px1(i) = b4m1 * drdx
        py1(i) = b4m1 * drdy
        pz1(i) = b4m1 * drdz
C
        px2(i) = b4m1 * dsdx
        py2(i) = b4m1 * dsdy
        pz2(i) = b4m1 * dsdz
C
        px3(i) = b4m1 * dtdx
        py3(i) = b4m1 * dtdy
        pz3(i) = b4m1 * dtdz
C
        px4(i) =-a4m1 * (drdx+dsdx+dtdx)
        py4(i) =-a4m1 * (drdy+dsdy+dtdy)
        pz4(i) =-a4m1 * (drdz+dsdz+dtdz)
C
        px5(i) = b4 * (drdx+dsdx)
        py5(i) = b4 * (drdy+dsdy)
        pz5(i) = b4 * (drdz+dsdz)
C
        px6(i) = b4 * (dsdx+dtdx)
        py6(i) = b4 * (dsdy+dtdy)
        pz6(i) = b4 * (dsdz+dtdz)
C
        px7(i) = b4 * (dtdx+drdx)
        py7(i) = b4 * (dtdy+drdy)
        pz7(i) = b4 * (dtdz+drdz)
C
        px8(i) = a4mb4  * drdx - px6(i)
        py8(i) = a4mb4  * drdy - py6(i)
        pz8(i) = a4mb4  * drdz - pz6(i)
C
        px9(i) = a4mb4  * dsdx - px7(i)
        py9(i) = a4mb4  * dsdy - py7(i)
        pz9(i) = a4mb4  * dsdz - pz7(i)
C
        px10(i)= a4mb4  * dtdx - px5(i)
        py10(i)= a4mb4  * dtdy - py5(i)
        pz10(i)= a4mb4  * dtdz - pz5(i)
C
      ENDDO
C
      RETURN