s10jacobnel.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	s10jacobnel (alph, beta, w, x1b, x2b, x3b, x4a, x5b, x6b, x7b, x8b, x9b, x10b, x8a, x9a, x10a, y1b, y2b, y3b, y4a, y5b, y6b, y7b, y8b, y9b, y10b, y8a, y9a, y10a, z1b, z2b, z3b, z4a, z5b, z6b, z7b, z8b, z9b, z10b, z8a, z9a, z10a, 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, vol, voldp, nel)

Function/Subroutine Documentation

s10jacobnel()

subroutine s10jacobnel	(		alph,
			beta,
			w,
		double precision, dimension(mvsiz)	x1b,
		double precision, dimension(mvsiz)	x2b,
		double precision, dimension(mvsiz)	x3b,
		double precision, dimension(mvsiz)	x4a,
		double precision, dimension(mvsiz)	x5b,
		double precision, dimension(mvsiz)	x6b,
		double precision, dimension(mvsiz)	x7b,
		double precision, dimension(mvsiz)	x8b,
		double precision, dimension(mvsiz)	x9b,
		double precision, dimension(mvsiz)	x10b,
		double precision, dimension(mvsiz)	x8a,
		double precision, dimension(mvsiz)	x9a,
		double precision, dimension(mvsiz)	x10a,
		double precision, dimension(mvsiz)	y1b,
		double precision, dimension(mvsiz)	y2b,
		double precision, dimension(mvsiz)	y3b,
		double precision, dimension(mvsiz)	y4a,
		double precision, dimension(mvsiz)	y5b,
		double precision, dimension(mvsiz)	y6b,
		double precision, dimension(mvsiz)	y7b,
		double precision, dimension(mvsiz)	y8b,
		double precision, dimension(mvsiz)	y9b,
		double precision, dimension(mvsiz)	y10b,
		double precision, dimension(mvsiz)	y8a,
		double precision, dimension(mvsiz)	y9a,
		double precision, dimension(mvsiz)	y10a,
		double precision, dimension(mvsiz)	z1b,
		double precision, dimension(mvsiz)	z2b,
		double precision, dimension(mvsiz)	z3b,
		double precision, dimension(mvsiz)	z4a,
		double precision, dimension(mvsiz)	z5b,
		double precision, dimension(mvsiz)	z6b,
		double precision, dimension(mvsiz)	z7b,
		double precision, dimension(mvsiz)	z8b,
		double precision, dimension(mvsiz)	z9b,
		double precision, dimension(mvsiz)	z10b,
		double precision, dimension(mvsiz)	z8a,
		double precision, dimension(mvsiz)	z9a,
		double precision, dimension(mvsiz)	z10a,
			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,
			vol,
		double precision, dimension(*)	voldp,
		integer	nel )

Definition at line 28 of file s10jacobnel.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
       double precision
     .   x1b(mvsiz),x2b(mvsiz),x3b(mvsiz),x4a(mvsiz),x5b(mvsiz),
     .   x6b(mvsiz),x7b(mvsiz),x8b(mvsiz),x9b(mvsiz),x10b(mvsiz),
     .   x8a(mvsiz),x9a(mvsiz),x10a(mvsiz),
     .   y1b(mvsiz),y2b(mvsiz),y3b(mvsiz),y4a(mvsiz),y5b(mvsiz),
     .   y6b(mvsiz),y7b(mvsiz),y8b(mvsiz),y9b(mvsiz),y10b(mvsiz),
     .   y8a(mvsiz),y9a(mvsiz),y10a(mvsiz),
     .   z1b(mvsiz),z2b(mvsiz),z3b(mvsiz),z4a(mvsiz),z5b(mvsiz),
     .   z6b(mvsiz),z7b(mvsiz),z8b(mvsiz),z9b(mvsiz),z10b(mvsiz),
     .   z8a(mvsiz),z9a(mvsiz),z10a(mvsiz),voldp(*)
 
       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),
     .   vol(mvsiz),alph,beta,w
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      DOUBLE PRECISION 
     .   D,A4MB4,A4,B4,BB,AB,B2,A4_0,B4_0,
     .   DXDR,DXDS,DXDT,DYDR,DYDS,DYDT,DZDR,DZDS,DZDT
     
      DOUBLE PRECISION AA,A4M1,B4M1
      
      DOUBLE PRECISION 
     .   DRDX, DSDX, DTDX,
     .   DRDY, DSDY, DTDY,
     .   DRDZ, DSDZ, DTDZ,
     .   DET,WUNSIX
C-----------------------------------------------
c      AA = (TWO*ALPH - ONE)*ALPH
c      BB = (TWO*BETA - ONE)*BETA
c      B2 = FOUR*BETA*BETA
c      AB = FOUR*ALPH*BETA      
cC-----------------------------------------------
C
       a4_0 = four * alph
       b4_0 = four * beta
       wunsix = w*one_over_6
#include   "nofusion.inc"
      DO i=1,nel
       aa = x5b(i) + x6b(i) + x7b(i)
     .    - x4a(i) - x8b(i) - x9b(i) - x10b(i)
       dxdr = x1b(i) +x8a(i) - x6b(i) + aa
       dxds = x2b(i) +x9a(i) - x7b(i) + aa
       dxdt = x3b(i) +x10a(i)- x5b(i) + aa
C
       aa = y5b(i) + y6b(i) + y7b(i)
     .    - y4a(i) - y8b(i) - y9b(i) - y10b(i)
       dydr = y1b(i) +y8a(i) - y6b(i) + aa
       dyds = y2b(i) +y9a(i) - y7b(i) + aa
       dydt = y3b(i) +y10a(i)- y5b(i) + aa
C
       aa = z5b(i) + z6b(i) + z7b(i)
     .    - z4a(i) - z8b(i) - z9b(i) - z10b(i)
       dzdr = z1b(i) +z8a(i) - z6b(i) + aa
       dzds = z2b(i) +z9a(i) - z7b(i) + aa
       dzdt = z3b(i) +z10a(i)- z5b(i) + aa
C
       drdx=dyds*dzdt-dzds*dydt
       dsdx=dydt*dzdr-dzdt*dydr
       dtdx=dydr*dzds-dzdr*dyds
C
       drdy=dzds*dxdt-dxds*dzdt
       dsdy=dzdt*dxdr-dxdt*dzdr
       dtdy=dzdr*dxds-dxdr*dzds
C
       drdz=dxds*dydt-dyds*dxdt
       dsdz=dxdt*dydr-dydt*dxdr
       dtdz=dxdr*dyds-dydr*dxds
C
       det = dxdr * drdx + dydr * drdy + dzdr * drdz
 
       d = one/max(em30,det)
C
       a4m1  = d *(a4_0 - one)
       b4m1  =  d *(b4_0 - one)
!
       b4 = d * b4_0
       a4 = d * a4_0
       a4mb4  = a4 - b4
       voldp(i) = wunsix * det 
       vol(i) = voldp(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