sysfus.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/.
C     REAL
#include "my_real.inc"
!||====================================================================
!||    sysfus        ../engine/source/system/sysfus.F
!||--- called by ------------------------------------------------------
!||    leccut        ../engine/source/tools/sect/leccut.F
!||--- calls      -----------------------------------------------------
!||    ancmsg        ../engine/source/output/message/message.F
!||--- uses       -----------------------------------------------------
!||    message_mod   ../engine/share/message_module/message_mod.F
!||====================================================================
       my_real  FUNCTION sysfus(IU,ITABM1,NUMNOD,MESS)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE message_mod
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 iu, numnod
      CHARACTER mess*40
      INTEGER itabm1(*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "warn_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER jinf, jsup, j
C-----------------------------------------------
c
      jinf=1
      jsup=numnod
C  045      J=NUMNOD/2
      j=min(1,numnod/2)
   10 IF(jsup<=jinf.AND.(iu-itabm1(j))/=0) THEN
        CALL ancmsg(msgid=187,anmode=aninfo,
     .            i1=iu,c1=mess)
        ierr=ierr+1
        sysfus=1.1
        RETURN
      ENDIF
      IF((iu-itabm1(j))==0)THEN
C     >CAS IU=TABM FIN DE LA RECHERCHE
         sysfus=itabm1(j+numnod)+0.1
         RETURN
      ELSE IF (iu-itabm1(j)<0) THEN
C     >CAS IU<TABM
         jsup=j-1
      ELSE
C     >CAS IU>TABM
         jinf=j+1
      ENDIF
      j=(jsup+jinf)/2
      GO TO 10
      END
C
C     REAL
!||====================================================================
!||    sysfus2       ../engine/source/system/sysfus.F
!||--- called by ------------------------------------------------------
!||    ale51_spmd2   ../engine/source/ale/ale51/ale51_spmd2.F
!||    ale51_spmd3   ../engine/source/ale/ale51/ale51_spmd3.F
!||    forcefingeo   ../engine/source/loads/general/forcefingeo.f
!||    fr_rlale      ../engine/source/mpi/kinematic_conditions/fr_rlink1.F
!||    fr_rlink1     ../engine/source/mpi/kinematic_conditions/fr_rlink1.F
!||    genh3d        ../engine/source/output/h3d/h3d_results/genh3d.F
!||    lcbcsf        ../engine/source/constraints/general/bcs/lcbcsf.f
!||    lecnoise      ../engine/source/general_controls/computation/lecnoise.F
!||====================================================================
      INTEGER FUNCTION sysfus2(IU,ITABM1,NUMNOD)
C cette fonction retourne le no interne correspondant a un no user ou 0
C si le noeud n'existe pas
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 iu, numnod
      INTEGER itabm1(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER jinf, jsup, j
C-----------------------------------------------
C
      ! Check exit parameters
 
      ! 1st NUMNOD=0
      IF (numnod==0) THEN
        sysfus2=0
        RETURN
      END IF
 
      ! 2nd NodeID is lower than smallest NodeID.
      IF ( iu-itabm1(1)<0 ) THEN
        sysfus2=0
        RETURN
      ENDIF
 
 
      jinf=1
      jsup=numnod
      j=min(1,numnod/2)
   10 IF(jsup<=jinf.AND.(iu-itabm1(j))/=0) THEN
C     >CAS NOEUD non trouve
        sysfus2=0
        RETURN
      ENDIF
      IF((iu-itabm1(j))==0)THEN
C     >CAS IU=TABM FIN DE LA RECHERCHE
         sysfus2=itabm1(j+numnod)
         RETURN
      ELSE IF (iu-itabm1(j)<0) THEN
C     >CAS IU<TABM
         jsup=j-1
      ELSE
C     >CAS IU>TABM
         jinf=j+1
      ENDIF
      j=(jsup+jinf)/2
      GO TO 10
      END