s4upd11t12.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	s4upd11t12 (offg, offg0, xd1, xd2, xd3, xd4, yd1, yd2, yd3, yd4, zd1, zd2, zd3, zd4, jac_1, sig, sigl, nel, jcvt)

Function/Subroutine Documentation

s4upd11t12()

subroutine s4upd11t12	(		offg,
			offg0,
		double precision, dimension(mvsiz)	xd1,
		double precision, dimension(mvsiz)	xd2,
		double precision, dimension(mvsiz)	xd3,
		double precision, dimension(mvsiz)	xd4,
		double precision, dimension(mvsiz)	yd1,
		double precision, dimension(mvsiz)	yd2,
		double precision, dimension(mvsiz)	yd3,
		double precision, dimension(mvsiz)	yd4,
		double precision, dimension(mvsiz)	zd1,
		double precision, dimension(mvsiz)	zd2,
		double precision, dimension(mvsiz)	zd3,
		double precision, dimension(mvsiz)	zd4,
			jac_1,
			sig,
			sigl,
		integer	nel,
		integer, intent(in)	jcvt )

Definition at line 30 of file s4upd11t12.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   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) :: JCVT
      INTEGER NEL
C     REAL
      my_real
     .   offg(*),offg0(*),jac_1(10,nel),sig(nel,6),sigl(nel,6)
      DOUBLE PRECISION 
     .   XD1(MVSIZ), XD2(MVSIZ), XD3(MVSIZ), XD4(MVSIZ),
     .   YD1(MVSIZ), YD2(MVSIZ), YD3(MVSIZ), YD4(MVSIZ),
     .   ZD1(MVSIZ), ZD2(MVSIZ), ZD3(MVSIZ), ZD4(MVSIZ)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I ,NCH,INDEX(MVSIZ),J    
      my_real
     .   jac_i(10,mvsiz)
      DOUBLE PRECISION 
     .   X1(MVSIZ), X2(MVSIZ), X3(MVSIZ), X4(MVSIZ),
     .   Y1(MVSIZ), Y2(MVSIZ), Y3(MVSIZ), Y4(MVSIZ),
     .   Z1(MVSIZ), Z2(MVSIZ), Z3(MVSIZ), Z4(MVSIZ)
C=======================================================================
C-----Always in global system and fixed then (not like Ismstr=11
C-----for JCVT>0 Fixed SAV is rotated w/ local system not like Ismstr=1
        nch=0
        DO i=1,nel
         IF(offg(i)/=offg0(i).AND.abs(offg(i)) > one ) THEN
           nch = nch +1
           index(nch) = i
         END IF
          END DO
          IF (nch==0) RETURN
        DO j=1,nch
          i= index(j)
          x1(j)=xd1(i)
          x2(j)=xd2(i)
          x3(j)=xd3(i)
          x4(j)=xd4(i)
          y1(j)=yd1(i)
          y2(j)=yd2(i)
          y3(j)=yd3(i)
          y4(j)=yd4(i)
          z1(j)=zd1(i)
          z2(j)=zd2(i)
          z3(j)=zd3(i)
          z4(j)=zd4(i)
        END DO
C
       CALL s4jac_i(
     .              x1   ,x2   ,x3   ,x4   ,
     .              y1   ,y2   ,y3   ,y4   ,
     .              z1   ,z2   ,z3   ,z4   ,
     .              jac_i ,nch  )
#include "vectorize.inc"
        DO j=1,nch
          i = index(j)
          jac_1(1:10,i) = jac_i(1:10,j)
          sigl(i,1:6) = sig(i,1:6)
        ENDDO      
C-----------
      RETURN