s10jacob1.F File Reference

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

Go to the source code of this file.

Functions/Subroutines

subroutine

s10jacob1 (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, nx1, nx2, nx3, nx4, nx5, nx6, nx7, nx8, nx9, nx10, vol, nnega, index, voldp)

Function/Subroutine Documentation

s10jacob1()

subroutine s10jacob1	(		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,
			nx1,
			nx2,
			nx3,
			nx4,
			nx5,
			nx6,
			nx7,
			nx8,
			nx9,
			nx10,
			vol,
		integer	nnega,
		integer, dimension(*)	index,
		double precision, dimension(mvsiz)	voldp )

Definition at line 28 of file s10jacob1.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 NNEGA,INDEX(*)
C     REAL
       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(mvsiz)
 
       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
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J
      DOUBLE PRECISION 
     .   D,A4MB4,A4,B4,BB,AB,B2,
     .   DXDR,DXDS,DXDT,DYDR,DYDS,DYDT,DZDR,DZDS,DZDT
     
      DOUBLE PRECISION 
     .   AA,A4M1(MVSIZ),B4M1(MVSIZ)
      
      DOUBLE PRECISION 
     .   DRDX(MVSIZ), DSDX(MVSIZ), DTDX(MVSIZ),
     .   DRDY(MVSIZ), DSDY(MVSIZ), DTDY(MVSIZ),
     .   DRDZ(MVSIZ), DSDZ(MVSIZ), DTDZ(MVSIZ),
     .   DET(MVSIZ)
C-----------------------------------------------
c      AA = (TWO*ALPH - ONE)*ALPH
c      BB = (TWO*BETA - ONE)*BETA
c      B2 = FOUR*BETA*BETA
c      AB = FOUR*ALPH*BETA      
C-----------------------------------------------
#include "vectorize.inc"
      DO j=1,nnega
       i = index(j)
       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(i)=dyds*dzdt-dzds*dydt
       dsdx(i)=dydt*dzdr-dzdt*dydr
       dtdx(i)=dydr*dzds-dzdr*dyds
C
       drdy(i)=dzds*dxdt-dxds*dzdt
       dsdy(i)=dzdt*dxdr-dxdt*dzdr
       dtdy(i)=dzdr*dxds-dxdr*dzds
C
       drdz(i)=dxds*dydt-dyds*dxdt
       dsdz(i)=dxdt*dydr-dydt*dxdr
       dtdz(i)=dxdr*dyds-dydr*dxds
C
       det(i) = dxdr * drdx(i)
     .        + dydr * drdy(i)
     .        + dzdr * drdz(i)
       voldp(i) = w * det(i) * one_over_6
       vol(i) = voldp(i)
C
      ENDDO
C      
#include "vectorize.inc"
      DO j=1,nnega
       i = index(j)
        d = one/max(em30,det(i))
C
        a4 = four * alph
        b4 = four * beta
        a4m1(i)  = a4 - one
        b4m1(i)  = b4 - one
        b4m1(i) = d * b4m1(i)
        a4m1(i) = d * a4m1(i)
        b4 = d * b4
        a4 = d * a4
        a4mb4  = a4 - b4
C
        px1(i) = b4m1(i) * drdx(i)
        py1(i) = b4m1(i) * drdy(i)
        pz1(i) = b4m1(i) * drdz(i)
C
        px2(i) = b4m1(i) * dsdx(i)
        py2(i) = b4m1(i) * dsdy(i)
        pz2(i) = b4m1(i) * dsdz(i)
C
        px3(i) = b4m1(i) * dtdx(i)
        py3(i) = b4m1(i) * dtdy(i)
        pz3(i) = b4m1(i) * dtdz(i)
C
        px4(i) =-a4m1(i) * (drdx(i)+dsdx(i)+dtdx(i))
        py4(i) =-a4m1(i) * (drdy(i)+dsdy(i)+dtdy(i))
        pz4(i) =-a4m1(i) * (drdz(i)+dsdz(i)+dtdz(i))
C
        px5(i) = b4 * (drdx(i)+dsdx(i))
        py5(i) = b4 * (drdy(i)+dsdy(i))
        pz5(i) = b4 * (drdz(i)+dsdz(i))
C
        px6(i) = b4 * (dsdx(i)+dtdx(i))
        py6(i) = b4 * (dsdy(i)+dtdy(i))
        pz6(i) = b4 * (dsdz(i)+dtdz(i))
C
        px7(i) = b4 * (dtdx(i)+drdx(i))
        py7(i) = b4 * (dtdy(i)+drdy(i))
        pz7(i) = b4 * (dtdz(i)+drdz(i))
C
        px8(i) = a4mb4  * drdx(i) - px6(i)
        py8(i) = a4mb4  * drdy(i) - py6(i)
        pz8(i) = a4mb4  * drdz(i) - pz6(i)
C
        px9(i) = a4mb4  * dsdx(i) - px7(i)
        py9(i) = a4mb4  * dsdy(i) - py7(i)
        pz9(i) = a4mb4  * dsdz(i) - pz7(i)
C
        px10(i)= a4mb4  * dtdx(i) - px5(i)
        py10(i)= a4mb4  * dtdy(i) - py5(i)
        pz10(i)= a4mb4  * dtdz(i) - pz5(i)
C
      ENDDO
C
      RETURN