s8_is17jac_i.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	s8_is17jac_i (elbuf_tab, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, index, nch, nel)

Function/Subroutine Documentation

s8_is17jac_i()

subroutine s8_is17jac_i	(	type(elbuf_struct_), target	elbuf_tab,
		double precision, dimension(*)	x1,
		double precision, dimension(*)	x2,
		double precision, dimension(*)	x3,
		double precision, dimension(*)	x4,
		double precision, dimension(*)	x5,
		double precision, dimension(*)	x6,
		double precision, dimension(*)	x7,
		double precision, dimension(*)	x8,
		double precision, dimension(*)	y1,
		double precision, dimension(*)	y2,
		double precision, dimension(*)	y3,
		double precision, dimension(*)	y4,
		double precision, dimension(*)	y5,
		double precision, dimension(*)	y6,
		double precision, dimension(*)	y7,
		double precision, dimension(*)	y8,
		double precision, dimension(*)	z1,
		double precision, dimension(*)	z2,
		double precision, dimension(*)	z3,
		double precision, dimension(*)	z4,
		double precision, dimension(*)	z5,
		double precision, dimension(*)	z6,
		double precision, dimension(*)	z7,
		double precision, dimension(*)	z8,
		integer, dimension(*)	index,
		integer	nch,
		integer	nel )

Definition at line 33 of file s8_is17jac_i.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE elbufdef_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 INDEX(*),NCH   ,NEL
      double precision
     .   x1(*), x2(*), x3(*), x4(*), x5(*), x6(*),
     .   x7(*), x8(*), y1(*), y2(*), y3(*), y4(*), y5(*), y6(*), y7(*),
     .   y8(*), z1(*), z2(*), z3(*), z4(*), z5(*), z6(*), z7(*), z8(*)
      TYPE(ELBUF_STRUCT_), TARGET :: ELBUF_TAB
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER IP,IR, IS, IT,NPTR,NPTS,NPTT,
     . LLPIJ
C     REAL
      TYPE(L_BUFEL_) ,POINTER :: LBUF     
      TYPE(G_BUFEL_) ,POINTER :: GBUF     
      my_real
     .   ajc1(mvsiz) , ajc2(mvsiz) , ajc3(mvsiz) ,
     .   ajc4(mvsiz) , ajc5(mvsiz) , ajc6(mvsiz) ,
     .   ajc7(mvsiz) , ajc8(mvsiz) , ajc9(mvsiz) ,
     .   hx(4,mvsiz) , hy(4,mvsiz), hz(4,mvsiz),
     .   pxc1(mvsiz),pxc2(mvsiz),pxc3(mvsiz),pxc4(mvsiz),
     .   pyc1(mvsiz),pyc2(mvsiz),pyc3(mvsiz),pyc4(mvsiz),
     .   pzc1(mvsiz),pzc2(mvsiz),pzc3(mvsiz),pzc4(mvsiz)
      my_real
     . wt,zr,zs,zt,wi
C-----------------------------------------------
      my_real
     .  w_gauss(9,9),a_gauss(9,9)
      DATA w_gauss / 
c---
     1 2.d0               ,0.d0               ,0.d0               ,
     1 0.d0               ,0.d0               ,0.d0               ,
     1 0.d0               ,0.d0               ,0.d0               ,
     2 1.d0               ,1.d0               ,0.d0               ,
     2 0.d0               ,0.d0               ,0.d0               ,
     2 0.d0               ,0.d0               ,0.d0               ,
     3 0.555555555555556d0,0.888888888888889d0,0.555555555555556d0,
     3 0.d0               ,0.d0               ,0.d0               ,
     3 0.d0               ,0.d0               ,0.d0               ,
     4 0.347854845137454d0,0.652145154862546d0,0.652145154862546d0,
     4 0.347854845137454d0,0.d0               ,0.d0               ,
     4 0.d0               ,0.d0               ,0.d0               ,
     5 0.236926885056189d0,0.478628670499366d0,0.568888888888889d0,
     5 0.478628670499366d0,0.236926885056189d0,0.d0               ,
     5 0.d0               ,0.d0               ,0.d0               ,
     6 0.171324492379170d0,0.360761573048139d0,0.467913934572691d0,
     6 0.467913934572691d0,0.360761573048139d0,0.171324492379170d0,
     6 0.d0               ,0.d0               ,0.d0               ,
     7 0.129484966168870d0,0.279705391489277d0,0.381830050505119d0,
     7 0.417959183673469d0,0.381830050505119d0,0.279705391489277d0,
     7 0.129484966168870d0,0.d0               ,0.d0               ,
     8 0.101228536290376d0,0.222381034453374d0,0.313706645877887d0,
     8 0.362683783378362d0,0.362683783378362d0,0.313706645877887d0,
     8 0.222381034453374d0,0.101228536290376d0,0.d0               ,
     9 0.081274388361574d0,0.180648160694857d0,0.260610696402935d0,
     9 0.312347077040003d0,0.330239355001260d0,0.312347077040003d0,
     9 0.260610696402935d0,0.180648160694857d0,0.081274388361574d0/
c------------------------------------------------------------
      DATA a_gauss / 
     1 0.d0               ,0.d0               ,0.d0               ,
     1 0.d0               ,0.d0               ,0.d0               ,
     1 0.d0               ,0.d0               ,0.d0               ,
     2 -.577350269189625d0,0.577350269189625d0,0.d0               ,
     2 0.d0               ,0.d0               ,0.d0               ,
     2 0.d0               ,0.d0               ,0.d0               , 
     3 -.774596669241483d0,0.d0               ,0.774596669241483d0,
     3 0.d0               ,0.d0               ,0.d0               ,
     3 0.d0               ,0.d0               ,0.d0               ,
     4 -.861136311594053d0,-.339981043584856d0,0.339981043584856d0,
     4 0.861136311594053d0,0.d0               ,0.d0               ,
     4 0.d0               ,0.d0               ,0.d0               ,
     5 -.906179845938664d0,-.538469310105683d0,0.d0               ,
     5 0.538469310105683d0,0.906179845938664d0,0.d0               ,
     5 0.d0               ,0.d0               ,0.d0               ,
     6 -.932469514203152d0,-.661209386466265d0,-.238619186083197d0,
     6 0.238619186083197d0,0.661209386466265d0,0.932469514203152d0,
     6 0.d0               ,0.d0               ,0.d0               ,
     7 -.949107912342759d0,-.741531185599394d0,-.405845151377397d0,
     7 0.d0               ,0.405845151377397d0,0.741531185599394d0,
     7 0.949107912342759d0,0.d0               ,0.d0               ,
     8 -.960289856497536d0,-.796666477413627d0,-.525532409916329d0,
     8 -.183434642495650d0,0.183434642495650d0,0.525532409916329d0,
     8 0.796666477413627d0,0.960289856497536d0,0.d0               ,
     9 -.968160239507626d0,-.836031107326636d0,-.613371432700590d0,
     9 -.324253423403809d0,0.d0               ,0.324253423403809d0,
     9 0.613371432700590d0,0.836031107326636d0,0.968160239507626d0/
C
C-----------------------------------------------
      gbuf => elbuf_tab%GBUF
C   cas GBUF%JAC_I  for all case
        CALL s8zjac_ic(
     .           x1  ,x2  ,x3  ,x4  ,x5  ,x6  ,x7  ,x8   ,
     .           y1  ,y2  ,y3  ,y4  ,y5  ,y6  ,y7  ,y8   ,
     .           z1  ,z2  ,z3  ,z4  ,z5  ,z6  ,z7  ,z8   ,
     .           ajc1 ,ajc2 ,ajc3 ,
     .           ajc4 ,ajc5 ,ajc6 ,
     .           ajc7 ,ajc8 ,ajc9 ,
     .           hx, hy, hz,
     .           pxc1, pxc2, pxc3, pxc4, 
     .           pyc1, pyc2, pyc3, pyc4, 
     .           pzc1, pzc2, pzc3, pzc4,
     .           gbuf%JAC_I,index,nch )
        nptr = elbuf_tab%NPTR
        npts = elbuf_tab%NPTS
        nptt = elbuf_tab%NPTT
        llpij = elbuf_tab%BUFLY(1)%L_PIJ
C        
        DO ir=1,nptr
        DO is=1,npts
        DO it=1,nptt
C-----------
          lbuf => elbuf_tab%BUFLY(1)%LBUF(ir,is,it)
c
           zr = a_gauss(ir,nptr)
           zs = a_gauss(is,npts)
           zt = a_gauss(it,nptt)
           wt = w_gauss(it,nptt)
           ip = ir + ( (is-1) + (it-1)*npts )*nptr
           wi = w_gauss(ir,nptr)*w_gauss(is,npts)*wt
C
C   cas LBUF%L_PIJ=24 global system w/o assumed strain for Isolid=17,18 only
          CALL s8zjac_i3(
     .        zr,zs,zt,wi,
     .        hx, hy, hz,
     .        ajc1,ajc2,ajc3,
     .        ajc4,ajc5,ajc6,
     .        ajc7,ajc8,ajc9,lbuf%JAC_I,llpij,lbuf%PIJ,lbuf%SIG,lbuf%SIGL,
     .        index,nch,nel)
c
        ENDDO
        ENDDO
        ENDDO
C       
      RETURN
C