comput_mesh_neighbour.F File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	comput_mesh_neighbour (dir, dir2, ipid, ipid2, id_meshsurf, tab_coinknot, l_tab_coinknot, tab_coinknot_test, l_tab_coinknot_test, knot_insere, knot_insere2, tab_meshsurfcut, l_tab_meshsurfcut, tab_meshsurfcut2, l_tab_meshsurfcut2, flag)

Function/Subroutine Documentation

comput_mesh_neighbour()

subroutine comput_mesh_neighbour	(	integer	dir,
		integer	dir2,
		integer	ipid,
		integer	ipid2,
		integer	id_meshsurf,
			tab_coinknot,
		integer	l_tab_coinknot,
			tab_coinknot_test,
		integer	l_tab_coinknot_test,
			knot_insere,
			knot_insere2,
		integer, dimension(*)	tab_meshsurfcut,
		integer	l_tab_meshsurfcut,
		integer, dimension(*)	tab_meshsurfcut2,
		integer	l_tab_meshsurfcut2,
		integer	flag )

Definition at line 28 of file comput_mesh_neighbour.F.

C----------------------------------------------------------
C     routine that will test and count the meshsurfs that
C     intersect with the meshsurfsurf being processed
C     NB : MESHSURF = SURFACE DE COUPE D'ELEMENTS DEFINISSANT
C          one refinement in one direction
C----------------------------------------------------------
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER TAB_MESHSURFCUT(*),TAB_MESHSURFCUT2(*)
      INTEGER DIR,DIR2,IPID,IPID2,L_TAB_MESHSURFCUT,L_TAB_MESHSURFCUT2,
     .        L_TAB_COINKNOT,L_TAB_COINKNOT_TEST,ID_MESHSURF,FLAG
      my_real tab_coinknot(2,l_tab_coinknot),
     .        tab_coinknot_test(2,l_tab_coinknot_test),
     .        knot_insere,knot_insere2
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER DIRTANG1, DIRTANG2
      my_real coinenglob(2,2), coinenglob_test(2,2)
C=======================================================================
C
C     in the parametric space
C
C     RAFFINEMENT X :      RAFFINEMENT Y :      RAFFINEMENT Z :   
C
C          Z                    X                    Y         
C          |                    |                    |
C          |_ _ Y               |_ _ Z               |_ _ X
C         /                    /                    /
C        X                    Y                    Z
C
C=======================================================================
C                                    |    
C    C1 ___________c2_ _ _ _ _p1 |               
C      |            |           /    |                 
C      |            |                |  C1 ___________C2
C      |            |         /      |    |            |_ _ _  _P1 _ _ _ _ P2   
C      |____________|_ _ _ _ _ _ P2  |    |            |      /           /
C     C4           C3      /         |    |            |
C                         C2    /    |    |____________|    /           /
C                        /|          |  C4             C3   
C                       / |   /      |                    /           /
C                      /  |          |                         
C                     /   | /        |                  /           /
C                  C1/    |C3        |                C2__________ C3
C                    |   /           |                /           /
C                    |  /            |               /           /
C                    | /  <= TESTEE  |              /           / <= TESTEE
C                    |/              |             /___________/
C                   C4               |           C1           C4
C                                    |
C=======================================================================
C----------------------------------------------------------------------
C   TRAITEMENT SIMPLE QUI PERMET D'ECARTER LES MESHSURF D'AUTRES PATCHS
C   or the parallel meshsurf
C----------------------------------------------------------------------
c
      IF(ipid2/=ipid) RETURN ! not the same patches
      IF(dir2==dir) RETURN ! Parallel Messiturf
c
C----------------------------------------------------------------------
C   construction of bounding boxes for the meshsurfs
C   PERMET DE SAVOIR SI IL Y A POTENTIELLEMENT INTERSECTION
C   (IL Y A ONE TEST SUR LES FONCTIONS EN AVAL, 
C   LA CERTITUDE D'INTERSECTION N'EST PAS NECESSAIRE)
C----------------------------------------------------------------------
c
      IF(dir==1) THEN
        dirtang1 = 2
        dirtang2 = 3
      ELSEIF(dir==2) THEN
        dirtang1 = 3
        dirtang2 = 1
      ELSEIF(dir==3) THEN
        dirtang1 = 1
        dirtang2 = 2
      ENDIF
c
      coinenglob(1,1) = minval(tab_coinknot(1,1:(l_tab_coinknot)))
      coinenglob(2,1) = minval(tab_coinknot(2,1:(l_tab_coinknot)))
      coinenglob(1,2) = maxval(tab_coinknot(1,1:(l_tab_coinknot)))
      coinenglob(2,2) = maxval(tab_coinknot(2,1:(l_tab_coinknot)))
c
      coinenglob_test(1,1) = minval(tab_coinknot_test(1,1:(l_tab_coinknot_test)))
      coinenglob_test(2,1) = minval(tab_coinknot_test(2,1:(l_tab_coinknot_test)))
      coinenglob_test(1,2) = maxval(tab_coinknot_test(1,1:(l_tab_coinknot_test)))
      coinenglob_test(2,2) = maxval(tab_coinknot_test(2,1:(l_tab_coinknot_test))) 
c
C----------------------------------------------------------------------
C   comparison of the intervals of the bounding boxes and the inserted knots,
C   PUIS VERIFICATION SUR LA TROISIEME DIRECTION COMPLEMENTAIRE AUX TWO DIRECTIONS
C   of the inserted knots from the superposition of the meshsurfs
C   TEST QUI EXCLUE : PLUS SIMPLE
C----------------------------------------------------------------------
c
      IF(dir2==dirtang1) THEN
        IF(knot_insere2<coinenglob(1,1).OR.knot_insere2>coinenglob(1,2)) RETURN
        IF(knot_insere<coinenglob_test(2,1).OR.knot_insere>coinenglob_test(2,2)) RETURN
        IF(coinenglob_test(1,2)<coinenglob(2,1)) RETURN
        IF(coinenglob_test(1,1)>coinenglob(2,2)) RETURN
      ELSEIF(dir2==dirtang2) THEN
        IF(knot_insere2<coinenglob(2,1).OR.knot_insere2>coinenglob(2,2)) RETURN
        IF(knot_insere<coinenglob_test(1,1).OR.knot_insere>coinenglob_test(1,2)) RETURN
        IF(coinenglob_test(2,2)<coinenglob(1,1)) RETURN
        IF(coinenglob_test(2,1)>coinenglob(1,2)) RETURN
      ENDIF
c
C----------------------------------------------------------------------
C   SI ON ARRIVE LA, IL Y A ONE CROISEMENT PROBABLE
C----------------------------------------------------------------------
c
      l_tab_meshsurfcut = l_tab_meshsurfcut + 1
      IF(flag==1) THEN
        tab_meshsurfcut(l_tab_meshsurfcut) = id_meshsurf
      ENDIF
c
C----------------------------------------------------------------------
C   we can take this meshsurf and all the meshsurfs it intersects
C   (TESTONS SANS POUR L'INSTANT CAR ON TESTE TOUTES LES MESHSURFS 
C   anyway, with respect to the meshsurf we are on
C----------------------------------------------------------------------
c
c      DO I=1,L_TAB_MESHSURFCUT2
c        IF(FLAG==0) THEN
c        L_TAB_MESHSURFCUT = L_TAB_MESHSURFCUT + 1
c        ENDIF
c        IF(FLAG==1) THEN
c          IOUT=0
c          DO J=1,L_TAB_MESHSURFCUT
c            IF(TAB_MESHSURFCUT(J)==TAB_MESHSURFCUT2(I)) IOUT=1
c          ENDDO
c          IF(IOUT==0) THEN
c            L_TAB_MESHSURFCUT = L_TAB_MESHSURFCUT + 1
c            TAB_MESHSURFCUT(L_TAB_MESHSURFCUT) = TAB_MESHSURFCUT2(I)
c          ENDIF
c        ENDIF
c      ENDDO
c
C----------------------------------------------------------------------
c
      RETURN