s8edericto3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	s8edericto3 (off, det, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, jac_i, nel)

Function/Subroutine Documentation

s8edericto3()

subroutine s8edericto3	(		off,
			det,
			px1,
			px2,
			px3,
			px4,
			py1,
			py2,
			py3,
			py4,
			pz1,
			pz2,
			pz3,
			pz4,
			jac_i,
		integer, intent(in)	nel )

Definition at line 28 of file s8edericto3.F.

C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
C   C o m m o n   B l o c k 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) :: NEL
C     REAL   
      my_real
     .   off(*),det(*),  
     .   px1(*), px2(*), px3(*), px4(*),  
     .   py1(*), py2(*), py3(*), py4(*),  
     .   pz1(*), pz2(*), pz3(*), pz4(*),  
     .   jac_i(10,mvsiz)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J ,ICOR
C     REAL
C                                                                     12
      my_real
     .   aji1(mvsiz), aji2(mvsiz), aji3(mvsiz),
     .   aji4(mvsiz), aji5(mvsiz), aji6(mvsiz),
     .   aji7(mvsiz), aji8(mvsiz), aji9(mvsiz),
     .   aj12(mvsiz), aj45(mvsiz), aj78(mvsiz)
C-----------------------------------------------
C  
      DO i=1,nel
      aji1(i)=jac_i(1,i)
      aji4(i)=jac_i(4,i)
      aji7(i)=jac_i(7,i)
      aji2(i)=jac_i(2,i)
      aji5(i)=jac_i(5,i)
      aji8(i)=jac_i(8,i)
      aji3(i)=jac_i(3,i)
      aji6(i)=jac_i(6,i)
      aji9(i)=jac_i(9,i)
      det(i) =jac_i(10,i)
      END DO
C
      DO 160 i=1,nel
      aj12(i)=aji1(i)-aji2(i)
      aj45(i)=aji4(i)-aji5(i)
  160 aj78(i)=aji7(i)-aji8(i)
      DO 170 i=1,nel
      px2(i)= aj12(i)-aji3(i)
      py2(i)= aj45(i)-aji6(i)
      pz2(i)= aj78(i)-aji9(i)
      px4(i)=-aj12(i)-aji3(i)
      py4(i)=-aj45(i)-aji6(i)
  170 pz4(i)=-aj78(i)-aji9(i)
      DO 180 i=1,nel
      aj12(i)=aji1(i)+aji2(i)
      aj45(i)=aji4(i)+aji5(i)
  180 aj78(i)=aji7(i)+aji8(i)
      DO 190 i=1,nel
      px1(i)=-aj12(i)-aji3(i)
      py1(i)=-aj45(i)-aji6(i)
      pz1(i)=-aj78(i)-aji9(i)
      px3(i)=aj12(i)-aji3(i)
      py3(i)=aj45(i)-aji6(i)
  190 pz3(i)=aj78(i)-aji9(i)
C       
      RETURN