s8xref_imp.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	s8xref_imp (offg, xref, xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8, r, nel)

Function/Subroutine Documentation

s8xref_imp()

subroutine s8xref_imp	(		offg,
		double precision, dimension(nel,21)	xref,
		double precision, dimension(mvsiz)	xd1,
		double precision, dimension(mvsiz)	xd2,
		double precision, dimension(mvsiz)	xd3,
		double precision, dimension(mvsiz)	xd4,
		double precision, dimension(mvsiz)	xd5,
		double precision, dimension(mvsiz)	xd6,
		double precision, dimension(mvsiz)	xd7,
		double precision, dimension(mvsiz)	xd8,
		double precision, dimension(mvsiz)	yd1,
		double precision, dimension(mvsiz)	yd2,
		double precision, dimension(mvsiz)	yd3,
		double precision, dimension(mvsiz)	yd4,
		double precision, dimension(mvsiz)	yd5,
		double precision, dimension(mvsiz)	yd6,
		double precision, dimension(mvsiz)	yd7,
		double precision, dimension(mvsiz)	yd8,
		double precision, dimension(mvsiz)	zd1,
		double precision, dimension(mvsiz)	zd2,
		double precision, dimension(mvsiz)	zd3,
		double precision, dimension(mvsiz)	zd4,
		double precision, dimension(mvsiz)	zd5,
		double precision, dimension(mvsiz)	zd6,
		double precision, dimension(mvsiz)	zd7,
		double precision, dimension(mvsiz)	zd8,
		double precision, dimension(3,3,mvsiz)	r,
		integer	nel )

Definition at line 30 of file s8xref_imp.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 NEL
C     REAL
      my_real
     .   offg(*)
      
      DOUBLE PRECISION 
     .   XD1(MVSIZ), XD2(MVSIZ), XD3(MVSIZ), XD4(MVSIZ), 
     .   XD5(MVSIZ), XD6(MVSIZ), XD7(MVSIZ), XD8(MVSIZ), 
     .   YD1(MVSIZ), YD2(MVSIZ), YD3(MVSIZ), YD4(MVSIZ), 
     .   YD5(MVSIZ), YD6(MVSIZ), YD7(MVSIZ), YD8(MVSIZ), 
     .   ZD1(MVSIZ), ZD2(MVSIZ), ZD3(MVSIZ), ZD4(MVSIZ), 
     .   ZD5(MVSIZ), ZD6(MVSIZ), ZD7(MVSIZ), ZD8(MVSIZ),
     .   R(3,3,MVSIZ), XREF(NEL,21)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      double precision
     .    u,v,w,ul,vl,wl
C=======================================================================
        DO i=1,nel
          !IF(ABS(OFFG(I)) <= ONE )THEN           !sb voir comment faire avec OFFG
            u = xd2(i)-xd1(i)
            v = yd2(i)-yd1(i)
            w = zd2(i)-zd1(i)
            ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
            vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
            wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
            xref(i,1) = ul
            xref(i,2) = vl
            xref(i,3) = wl
 
            u = xd3(i)-xd1(i)
            v = yd3(i)-yd1(i)
            w = zd3(i)-zd1(i)
            ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
            vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
            wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
            xref(i,4) = ul
            xref(i,5) = vl
            xref(i,6) = wl
      
            u = xd4(i)-xd1(i)
            v = yd4(i)-yd1(i)
            w = zd4(i)-zd1(i)
            ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
            vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
            wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
            xref(i,7) = ul
            xref(i,8) = vl
            xref(i,9) = wl
      
            u = xd5(i)-xd1(i)
            v = yd5(i)-yd1(i)
            w = zd5(i)-zd1(i)
            ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
            vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
            wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
            xref(i,10) = ul
            xref(i,11) = vl
            xref(i,12) = wl
      
            u = xd6(i)-xd1(i)
            v = yd6(i)-yd1(i)
            w = zd6(i)-zd1(i)
            ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
            vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
            wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
            xref(i,13) = ul
            xref(i,14) = vl
            xref(i,15) = wl
      
            u = xd7(i)-xd1(i)
            v = yd7(i)-yd1(i)
            w = zd7(i)-zd1(i)
            ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
            vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
            wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
            xref(i,16) = ul
            xref(i,17) = vl
            xref(i,18) = wl
      
            u = xd8(i)-xd1(i)
            v = yd8(i)-yd1(i)
            w = zd8(i)-zd1(i)
            ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
            vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
            wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
            xref(i,19) = ul
            xref(i,20) = vl
            xref(i,21) = wl
        
        ENDDO
C-----------
      RETURN