i12cor3.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
!||====================================================================
!||    i12cor3    ../engine/source/interfaces/inter3d/i12cor3.F
!||--- called by ------------------------------------------------------
!||    intti12a   ../engine/source/interfaces/interf/intti12.F
!||====================================================================
      SUBROUTINE i12cor3(
     1   XS,      XM,      IRECT,   MSR,
     2   IRTL,    ILEV,    X,       TETS,
     3   IX1,     IX2,     IX3,     IX4,
     4   X1,      X2,      X3,      X4,
     5   XI,      Y1,      Y2,      Y3,
     6   Y4,      YI,      Z1,      Z2,
     7   Z3,      Z4,      ZI,      XFACE,
     8   LFT,     LLT,     NFT)
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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(INOUT) :: LFT
      INTEGER, INTENT(INOUT) :: LLT
      INTEGER, INTENT(INOUT) :: NFT
      INTEGER IRECT(4,*), MSR(*), IRTL(*),ILEV
      my_real
     .   XS(3,*),XM(3,*),X(3,*),TETS(*)
      INTEGER, DIMENSION(MVSIZ), INTENT(OUT) ::  IX1,IX2,IX3,IX4
      my_real, DIMENSION(MVSIZ), INTENT(OUT) ::  X1,X2,X3,X4,XI
      my_real, DIMENSION(MVSIZ), INTENT(OUT) ::  y1,y2,y3,y4,yi
      my_real, DIMENSION(MVSIZ), INTENT(OUT) ::  z1,z2,z3,z4,zi
      my_real, DIMENSION(MVSIZ), INTENT(OUT) ::  xface
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, IL, IG, L, NN
      my_real
     .   deupi,tetmin,tetmax,rmin
C-----------------------------------------------
 
      deupi=two*pi
 
      DO i=lft,llt
        xface(i)=one
        il=i+nft
        l=irtl(il)
 
        xi(i)=xs(1,il)
        yi(i)=xs(2,il)
        zi(i)=xs(3,il)
 
        nn=irect(1,l)
        ix1(i)=msr(nn)
        x1(i)=xm(1,nn)
        y1(i)=xm(2,nn)
        z1(i)=xm(3,nn)
 
        nn=irect(2,l)
        ix2(i)=msr(nn)
        x2(i)=xm(1,nn)
        y2(i)=xm(2,nn)
        z2(i)=xm(3,nn)
 
        nn=irect(3,l)
        ix3(i)=msr(nn)
        x3(i)=xm(1,nn)
        y3(i)=xm(2,nn)
        z3(i)=xm(3,nn)
 
        nn=irect(4,l)
        ix4(i)=msr(nn)
        x4(i)=xm(1,nn)
        y4(i)=xm(2,nn)
        z4(i)=xm(3,nn)
      ENDDO
      
      IF(ilev==1)THEN
      !--------------------------------------------------!
      !   CORRECTION DES TET SI A CHEVAL AXE Y NEGATIF   !
      !--------------------------------------------------!
        DO i=lft,llt
          tetmin=min(z1(i),z2(i),z3(i),z4(i),zi(i))
          tetmax=max(z1(i),z2(i),z3(i),z4(i),zi(i))
          IF(tetmax-tetmin >= pi)THEN
            IF(z1(i)<zero)z1(i)=z1(i)+deupi
            IF(z2(i)<zero)z2(i)=z2(i)+deupi
            IF(z3(i)<zero)z3(i)=z3(i)+deupi
            IF(z4(i)<zero)z4(i)=z4(i)+deupi
            IF(zi(i)<zero)zi(i)=zi(i)+deupi
          ENDIF
        ENDDO
 
        !--------------------------------------------------!
        !   CAS D'UNE FACETTE CENTRALE                     !
        !--------------------------------------------------!
        DO i=lft,llt
          rmin=min(y1(i),y2(i),y3(i),y4(i))
          IF(yi(i)<rmin)THEN
            tetmin=min(z1(i),z2(i),z3(i),z4(i))
            tetmax=max(z1(i),z2(i),z3(i),z4(i))
            IF(tetmax-tetmin >= pi)THEN
              il=i+nft
              l=irtl(il)
              nn=msr(irect(1,l))
              ix1(i)=nn
              x1(i)=x(1,nn)
              y1(i)=x(2,nn)
              z1(i)=x(3,nn)
 
              nn=msr(irect(2,l))
              ix2(i)=nn
              x2(i)=x(1,nn)
              y2(i)=x(2,nn)
              z2(i)=x(3,nn)
 
              nn=msr(irect(3,l))
              ix3(i)=nn
              x3(i)=x(1,nn)
              y3(i)=x(2,nn)
              z3(i)=x(3,nn)
 
              nn=msr(irect(4,l))
              ix4(i)=nn
              x4(i)=x(1,nn)
              y4(i)=x(2,nn)
              z4(i)=x(3,nn)
              tets(i)=10001.
            ENDIF          
          ENDIF
        ENDDO !I=LFT,LLT
      ENDIF !(ILEV==1)
 
      RETURN
      END