ig3daire.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	ige3daire (itel, n, xxi, yyi, zzi, wwi, idx, idy, idz, aire, nctrl, gaussx, gaussy, gaussz, kx, ky, kz, px, py, pz)

Function/Subroutine Documentation

ige3daire()

subroutine ige3daire	(	integer	itel,
		integer	n,
		dimension(*)	xxi,
		dimension(*)	yyi,
		dimension(*)	zzi,
		dimension(*)	wwi,
		integer	idx,
		integer	idy,
		integer	idz,
		dimension(*)	aire,
		integer	nctrl,
			gaussx,
			gaussy,
			gaussz,
		dimension(*)	kx,
		dimension(*)	ky,
		dimension(*)	kz,
		integer	px,
		integer	py,
		integer	pz )

Definition at line 32 of file ig3daire.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE message_mod
C--------------------------------------------------------------------------------------------------------
C Implicite Types
C--------------------------------------------------------------------------------------------------------
#include     "implicit_f.inc"
C--------------------------------------------------------------------------------------------------------
C Dummy Arguments 
C--------------------------------------------------------------------------------------------------------
      INTEGER IDX, IDY, IDZ, NCTRL, PX, PY, PZ, ITEL, N
      my_real
     .  gaussx, gaussy, gaussz
      my_real, 
     .  DIMENSION(*)  :: aire
      my_real,
     .  DIMENSION(*) :: xxi,yyi,zzi,wwi
      my_real,
     .  DIMENSION(*) :: kx, ky, kz
C--------------------------------------------------------------------------------------------------------
C Local Variables
C--------------------------------------------------------------------------------------------------------
      INTEGER NUMLOC, I, J, K, NA, NB, NC
      my_real,
     .  DIMENSION(NCTRL) :: r
      my_real,
     .  DIMENSION(NCTRL,3) :: drdxi
      my_real
     .  sumtot
      my_real, 
     .  DIMENSION(PX+1) :: fn, dndxi
      my_real, 
     .  DIMENSION(PY+1) :: fm, dmdxi
      my_real, 
     .  DIMENSION(PZ+1) :: fl, dldxi
      my_real, 
     .  DIMENSION(3) :: xi, sumxi
      my_real, 
     .  DIMENSION(3,3)  :: dxdxi, ajmat, dxidtildexi
C--------------------------------------------------------------------------------------------------------
 
C INITIALISATION DES DIFFERENTES MATRICES LOCALES
 
 
      dxdxi(:,:)=zero
      dxidtildexi(:,:)=zero
      ajmat(:,:)=zero
 
 
      xi(1) = ((kx(idx+1)-kx(idx))*gaussx + (kx(idx+1)+(kx(idx))))/two
      xi(2) = ((ky(idy+1)-ky(idy))*gaussy + (ky(idy+1)+(ky(idy))))/two
      xi(3) = ((kz(idz+1)-kz(idz))*gaussz + (kz(idz+1)+(kz(idz))))/two
 
C CALCULATE UNIVARIATE B-SPLINE FUNCTION AT XI POINT
 
      CALL dersbasisfuns(idx, px, xi(1), kx, fn, dndxi)
      CALL dersbasisfuns(idy, py, xi(2), ky, fm, dmdxi)
      CALL dersbasisfuns(idz, pz, xi(3), kz, fl, dldxi)
 
C BUILD NUMERATORS AND DENOMINATORS
 
      numloc=0
      sumtot=zero
      sumxi(1)=zero
      sumxi(2)=zero
      sumxi(3)=zero
 
C BOUCLE SUR TOUT LES POINTS DE CONTROLE 
 
      DO k = 0,pz
        DO j = 0,py
          DO i = 0,px
            numloc = numloc+1
            r(numloc)=fn(px+1-i)*fm(py+1-j)*fl(pz+1-k)*wwi(numloc)
            sumtot=sumtot+r(numloc)
            drdxi(numloc,1)=dndxi(px+1-i)*fm(py+1-j)*fl(pz+1-k)*
     .                       wwi(numloc)
            sumxi(1)=sumxi(1)+drdxi(numloc,1)
            drdxi(numloc,2)=fn(px+1-i)*dmdxi(py+1-j)*fl(pz+1-k)*
     .                       wwi(numloc)
            sumxi(2)=sumxi(2)+drdxi(numloc,2)
            drdxi(numloc,3)=fn(px+1-i)*fm(py+1-j)*dldxi(pz+1-k)*
     .                       wwi(numloc)
            sumxi(3)=sumxi(3)+drdxi(numloc,3)
          ENDDO
        ENDDO
      ENDDO
 
C DIVIDE BY DENOMINATOR TO COMPLETE DEFINITION OF FUNCTION AND DERIVATIVES
 
      DO numloc=1,nctrl
        r(numloc)=r(numloc)/sumtot
      ENDDO
 
      DO i=1,3
        DO numloc=1,nctrl
          drdxi(numloc,i)=(drdxi(numloc,i)-r(numloc)*sumxi(i))/sumtot
        ENDDO
      ENDDO
 
C GRADIENT OF MAPPING FROM PARAMETER SPACE TO PHYSICAL SPACE 
 
      DO numloc=1,nctrl
        DO nb=1,3
          dxdxi(1,nb)=dxdxi(1,nb)+xxi(numloc)*drdxi(numloc,nb)
          dxdxi(2,nb)=dxdxi(2,nb)+yyi(numloc)*drdxi(numloc,nb)
          dxdxi(3,nb)=dxdxi(3,nb)+zzi(numloc)*drdxi(numloc,nb)
        ENDDO                             
      ENDDO
 
C GRADIENT OF MAPPING FROM PARENT ELEMENT TO PARAMETER SPACE
 
      dxidtildexi(1,1)=(kx(idx+1)-kx(idx))/two
      dxidtildexi(2,2)=(ky(idy+1)-ky(idy))/two
      dxidtildexi(3,3)=(kz(idz+1)-kz(idz))/two
 
      DO na=1,3
        DO nb=1,3
          DO nc=1,3
            ajmat(na,nb)=ajmat(na,nb)+dxdxi(na,nc)*dxidtildexi(nc,nb)
          ENDDO
        ENDDO
      ENDDO
 
      aire(1) = sqrt(ajmat(1,1)*ajmat(1,1)+ajmat(2,1)*ajmat(2,1)
     .              +ajmat(3,1)*ajmat(3,1))
     .         *sqrt(ajmat(1,2)*ajmat(1,2)+ajmat(2,2)*ajmat(2,2)
     .              +ajmat(3,2)*ajmat(3,2)) 
      aire(2) = sqrt(ajmat(1,1)*ajmat(1,1)+ajmat(2,1)*ajmat(2,1)
     .              +ajmat(3,1)*ajmat(3,1))
     .         *sqrt(ajmat(1,3)*ajmat(1,3)+ajmat(2,3)*ajmat(2,3)
     .              +ajmat(3,3)*ajmat(3,3))
      aire(3) = sqrt(ajmat(1,2)*ajmat(1,2)+ajmat(2,2)*ajmat(2,2)
     .              +ajmat(3,2)*ajmat(3,2))
     .         *sqrt(ajmat(1,3)*ajmat(1,3)+ajmat(2,3)*ajmat(2,3)
     .              +ajmat(3,3)*ajmat(3,3)) 
    
      RETURN