n2edge3l.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/.
!||====================================================================
!||    n2edge3l     ../engine/source/interfaces/int24/n2edge3l.f
!||--- called by ------------------------------------------------------
!||    intersecsh   ../engine/source/interfaces/int24/i24dst3.F
!||====================================================================
        SUBROUTINE n2edge3l(
     1                           XXI    ,YYI      ,ZZI    ,
     2                           XXJ    ,YYJ      ,ZZJ    ,
     3                           NX     ,NY       ,NZ     ,
     4                           XI     ,YI       ,ZI     ,BBB )
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-----------------------------------------------
C     REAL
      my_real , INTENT(IN)     :: 
     1                           xxi    ,yyi      ,zzi    ,
     2                           xxj    ,yyj      ,zzj    ,
     3                           nx     ,ny       ,nz     ,
     4                           xi     ,yi       ,zi     
      my_real , INTENT(OUT)    :: bbb
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      my_real
     .     xmij(3),nij(3),di(3),aa,norm
C-----------------------------------------------
             xmij(1)=xxj-xxi
             xmij(2)=yyj-yyi
             xmij(3)=zzj-zzi
             nij(1)= ny*xmij(3) - nz*xmij(2)
             nij(2)= nz*xmij(1) - nx*xmij(3)
             nij(3)= nx*xmij(2) - ny*xmij(1)
             aa = nij(1)*nij(1)+nij(2)*nij(2)+nij(3)*nij(3)
             norm=max(em20,sqrt(aa))
             di(1)=half*(xxi+xxj)-xi
             di(2)=half*(yyi+yyj)-yi
             di(3)=half*(zzi+zzj)-zi
             bbb = (di(1)*nij(1)+di(2)*nij(2)+di(3)*nij(3))/norm
             IF (abs(bbb)<em3) then
               xmij(1)=xxj-xi
               xmij(2)=yyj-yi
               xmij(3)=zzj-zi
               aa = di(1)*xmij(1)+di(2)*xmij(2)+di(3)*xmij(3)
               IF (aa>abs(bbb)) bbb = aa
             END IF
C
      RETURN
      END