i12loc3.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/.
!||====================================================================
!||    i12loc3    ../engine/source/interfaces/inter3d/i12loc3.F
!||--- called by ------------------------------------------------------
!||    intti12a   ../engine/source/interfaces/interf/intti12.F
!||====================================================================
      SUBROUTINE i12loc3(
     1   XS,      XM,      IRECT,   LMSR,
     2   MSR,     NSV,     ILOC,    NSEG,
     3   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   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,*), LMSR(*), MSR(*), NSV(*), ILOC(*), NSEG(*)
      my_real
     .   xs(3,*),xm(3,*)
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,  JL, JLNEW, LL2, LL1, LL, LG, J, K, M, N
      my_real
     .   GMS, CMS,
     .   dms, ems, fms
C-----------------------------------------------
      DO i=lft,llt
        il=i+nft
        jl=iloc(il)
        jlnew=jl
        ll2=nseg(jl+1)-1
        ll1=nseg(jl)
        gms=1.e30
        DO ll=ll1,ll2
          lg=lmsr(ll)
          j=irect(1,lg)
          k=irect(2,lg)
          m=irect(3,lg)
          n=irect(4,lg)
          cms=(xs(1,il)-xm(1,j))**2+(xs(2,il)-xm(2,j))**2
     &                             +(xs(3,il)-xm(3,j))**2
          dms=(xs(1,il)-xm(1,k))**2+(xs(2,il)-xm(2,k))**2
     &                             +(xs(3,il)-xm(3,k))**2
          ems=(xs(1,il)-xm(1,m))**2+(xs(2,il)-xm(2,m))**2
     &                             +(xs(3,il)-xm(3,m))**2
          fms=(xs(1,il)-xm(1,n))**2+(xs(2,il)-xm(2,n))**2
     &                             +(xs(3,il)-xm(3,n))**2
          IF(cms<gms) THEN
            gms=cms
            jlnew=j
          ENDIF
          IF(dms<gms) THEN
            gms=dms
            jlnew=k
          ENDIF
          IF(ems<gms) THEN
            gms=ems
            jlnew=m
          ENDIF
          IF(fms<gms) THEN
            gms=fms
            jlnew=n
          ENDIF
        END DO !LL=LL1,LL2
        iloc(il)=jlnew      
      END DO !I=LFT,LLT
 
      RETURN
      END