create_rbody_clause.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/.
!||====================================================================
!||    create_rbody_clause    ../starter/source/model/sets/create_rbody_clause.f
!||--- called by ------------------------------------------------------
!||    hm_set                 ../starter/source/model/sets/hm_set.F
!||--- calls      -----------------------------------------------------
!||    create_rbody_box       ../starter/source/model/sets/create_rbody_box.F
!||    create_rbody_list      ../starter/source/model/sets/create_rbody_clause.F
!||    create_rbody_list_g    ../starter/source/model/sets/create_rbody_clause.F
!||--- uses       -----------------------------------------------------
!||    hm_option_read_mod     ../starter/share/modules1/hm_option_read_mod.F
!||    message_mod            ../starter/share/message_module/message_mod.F
!||    submodel_mod           ../starter/share/modules1/submodel_mod.F
!||====================================================================
      SUBROUTINE create_rbody_clause(
     .                        CLAUSE   ,IRBODYM  ,JCLAUSE ,OPT_G ,IS_AVAILABLE ,
     .                        LSUBMODEL,OPT_B   ,IBOX    ,X     ,SKEW         ,
     .                        SET_TITLE,KEYSET  ,RBY_MSN )
C-----------------------------------------------
C   ROUTINE DESCRIPTION :
C   ===================
C   Treat the RBODY Clause, read RBODYs from HM_READER & fill clause
C   Calls CREATE_RBODY_LIST (simple list)
C   Calls CREATE_RBODY_LIST_G (RBODY_G : All nodes from a MIN to MAX with increment)
C------------------------------------------------------------------
C   DUMMY ARGUMENTS DESCRIPTION:
C   ===================
C
C     NAME          DESCRIPTION                         
C
C     CLAUSE        (SET structure) Clause to be treated
C     ITABM1        MAP Table UID -> LocalID
C     JCLAUSE       parameter with HM_READER (current clause read)
C     Opt_G         Opt_G operator 1 if PART_G is set, 0 else
C     IS_AVAILABLE  Bool / Result of HM_interface
C     LSUBMODEL     SUBMODEL Structure.
C============================================================================
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE setdef_mod
      USE submodel_mod
      USE message_mod
      USE hm_option_read_mod
      USE optiondef_mod
      USE names_and_titles_mod , ONLY : nchartitle, ncharfield
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com04_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JCLAUSE,OPT_G,OPT_B
      LOGICAL :: IS_AVAILABLE
      INTEGER, INTENT(IN), DIMENSION(NRBODY,2) :: IRBODYM
      INTEGER, INTENT(IN), DIMENSION(2,NRBODY) :: RBY_MSN
      my_real x(3,*),skew(lskew,*)
      CHARACTER(LEN=NCHARFIELD) :: KEYSET
      CHARACTER(LEN=NCHARTITLE) :: SET_TITLE
C-----------------------------------------------
      TYPE (SET_) ::  CLAUSE
      TYPE(submodel_data),INTENT(IN):: LSUBMODEL(*)
      TYPE (BOX_)  , DIMENSION(NBBOX)  :: IBOX
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
!
        IF ( opt_g == 1 ) THEN
 
             CALL create_rbody_list_g(clause,  irbodym  ,jclause ,is_available ,lsubmodel)
 
        ELSEIF ( opt_g == 0 .AND. opt_b == 0 ) THEN
 
             CALL create_rbody_list(clause,  irbodym  ,jclause ,is_available ,lsubmodel)
 
        ELSEIF ( opt_b == 1 ) THEN
 
             CALL create_rbody_box(clause  ,irbodym  ,jclause ,is_available ,lsubmodel,
     .                             ibox    ,x        ,skew    ,set_title    ,keyset   ,
     .                             rby_msn )
        ENDIF
C-----------------------------------------------
      END
!||====================================================================
!||    create_rbody_list           ../starter/source/model/sets/create_rbody_clause.F
!||--- called by ------------------------------------------------------
!||    create_rbody_clause         ../starter/source/model/sets/create_rbody_clause.F
!||--- calls      -----------------------------------------------------
!||    ancmsg                      ../starter/source/output/message/message.F
!||    hm_get_int_array_2indexes   ../starter/source/devtools/hm_reader/hm_get_int_array_2indexes.F
!||    hm_get_int_array_index      ../starter/source/devtools/hm_reader/hm_get_int_array_index.F
!||    set_usrtos                  ../starter/source/model/sets/ipartm1.F
!||--- uses       -----------------------------------------------------
!||    hm_option_read_mod          ../starter/share/modules1/hm_option_read_mod.F
!||    message_mod                 ../starter/share/message_module/message_mod.F
!||    submodel_mod                ../starter/share/modules1/submodel_mod.F
!||====================================================================
      SUBROUTINE create_rbody_list(
     .               CLAUSE,  IRBODYM  ,JCLAUSE ,IS_AVAILABLE ,LSUBMODEL)
C-----------------------------------------------
C   ROUTINE DESCRIPTION :
C   ===================
C   Create RBODY Clause from LIST
C------------------------------------------------------------------
C   DUMMY ARGUMENTS DESCRIPTION:
C   ===================
C
C     NAME          DESCRIPTION                         
C
C     CLAUSE        (SET structure) Clause to be treated
C     IRBODYM       MAP Table UID -> LocalID
C     JCLAUSE       parameter with HM_READER (current clause read)
C     IS_AVAILABLE  Bool / Result of HM_interface
C     LSUBMODEL     SUBMODEL Structure.
C============================================================================
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE setdef_mod
      USE submodel_mod
      USE message_mod
      USE hm_option_read_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com04_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JCLAUSE
      LOGICAL :: IS_AVAILABLE
      INTEGER, INTENT(IN), DIMENSION(NRBODY,2) :: IRBODYM
!
      TYPE (SET_) ::  CLAUSE
      TYPE(submodel_data),INTENT(IN):: LSUBMODEL(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,IDS,NINDX,LIST_SIZE,IDS_MAX,NODSYS,RBYM
      INTEGER IWORK(70000)
!
      INTEGER, ALLOCATABLE, DIMENSION(:) :: RBYM_READ_TMP,INDEX,SORT
C
      INTEGER SET_USRTOS
      EXTERNAL SET_USRTOS
C=======================================================================
 
      CALL hm_get_int_array_index('idsmax' ,ids_max ,jclause,is_available,lsubmodel)
 
      ALLOCATE(rbym_read_tmp(ids_max))
      rbym_read_tmp(1:ids_max) = 0
      ALLOCATE(index(2*ids_max))
      ALLOCATE(sort(ids_max))
 
      index = 0
 
      nindx = 0
      list_size = 0
 
 
      ! Read & convert Nodes
      ! ---------------------
      DO i=1,ids_max
         CALL hm_get_int_array_2indexes('ids',ids,jclause,i,is_available,lsubmodel)
 
 
 
         rbym = set_usrtos(ids,irbodym,nrbody)
         IF (rbym == 0) THEN
           ! Rbody was not found. Issue a Warning & Skip.
           CALL ancmsg(msgid=1902,anmode=aninfo,
     .                             msgtype=msgwarning,
     .                             i1 = clause%SET_ID,
     .                             i2=ids,
     .                             c1=trim(clause%TITLE),
     .                             c2='RBODY')
         ELSE
 
           rbym=irbodym(rbym,2)
!
           nindx=nindx+1    !   nb of CLAUSE Rbody
           rbym_read_tmp(nindx) = rbym
         ENDIF
 
      ENDDO ! DO K=1,IDS_MAX
 
 
 
 
 
      ! Sort the Readed RBODYs and remove eventual duplicates
      ! ---------------------------------------------------- 
 
      DO i=1,nindx
        index(i) = i
        sort(i) = rbym_read_tmp(i)
      ENDDO
      CALL my_orders(0,iwork,sort,index,nindx,1)
 
      DO i=1,nindx
         rbym_read_tmp(i) = sort(index(i))
      ENDDO
 
      CALL remove_duplicates(rbym_read_tmp,nindx,list_size)
 
      ! Copy in final SET
      ! ------------------
      clause%NB_RBODY = list_size  
      ALLOCATE( clause%RBODY( list_size ) )
 
      DO i=1,list_size
        clause%RBODY(i) = rbym_read_tmp(i)
      ENDDO
C-------------------------
       DEALLOCATE(rbym_read_tmp)
       DEALLOCATE(sort)
       DEALLOCATE(index)
C-------------------------
      RETURN
      END
 
!||====================================================================
!||    create_rbody_list_g         ../starter/source/model/sets/create_rbody_clause.F
!||--- called by ------------------------------------------------------
!||    create_rbody_clause         ../starter/source/model/sets/create_rbody_clause.F
!||--- calls      -----------------------------------------------------
!||    hm_get_int_array_2indexes   ../starter/source/devtools/hm_reader/hm_get_int_array_2indexes.F
!||    hm_get_int_array_index      ../starter/source/devtools/hm_reader/hm_get_int_array_index.F
!||    set_usrtos_nearest          ../starter/source/model/sets/ipartm1.F
!||--- uses       -----------------------------------------------------
!||    hm_option_read_mod          ../starter/share/modules1/hm_option_read_mod.F
!||    message_mod                 ../starter/share/message_module/message_mod.F
!||    submodel_mod                ../starter/share/modules1/submodel_mod.F
!||====================================================================
      SUBROUTINE create_rbody_list_g(
     .               CLAUSE,  IRBODYM  ,JCLAUSE ,IS_AVAILABLE ,LSUBMODEL)
C-----------------------------------------------
C   ROUTINE DESCRIPTION :
C   ===================
C   Create RBODY Clause from Generation All parts from Min to Max with Increment (Opt_G)
C--------------------------------------------------------------------------------------
C   DUMMY ARGUMENTS DESCRIPTION:
C   ===================
C
C     NAME          DESCRIPTION                         
C
C     CLAUSE        (SET structure) Clause to be treated
C     IRBODYM       MAP Table UID -> LocalID
C     JCLAUSE       parameter with HM_READER (current clause read)
C     IS_AVAILABLE  Bool / Result of HM_interface
C     LSUBMODEL     SUBMODEL Structure.
C============================================================================
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE setdef_mod
      USE submodel_mod
      USE message_mod
      USE hm_option_read_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com04_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JCLAUSE
      LOGICAL :: IS_AVAILABLE
      INTEGER, INTENT(IN), DIMENSION(NRBODY,2) :: IRBODYM
!
      TYPE (SET_) ::  CLAUSE
      TYPE(submodel_data),INTENT(IN):: LSUBMODEL(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,IDS,LIST_SIZE,IDS_MAX,GENE_MAX,K,R,R1
      INTEGER START_GENE,END_GENE,INCR_GENE,RSTART,RSTOP,STACK,STACK_ONE,NB_RESULT
!-
      INTEGER, ALLOCATABLE, DIMENSION(:) :: RBYM_READ_TMP,
     .         RBYM_READ_ONE,RESULT
C
      INTEGER  SET_USRTOS_NEAREST
      EXTERNAL SET_USRTOS_NEAREST
C=======================================================================
      CALL hm_get_int_array_index('genemax' ,gene_max ,jclause,is_available,lsubmodel)
 
      ALLOCATE(rbym_read_tmp(nrbody))
      ALLOCATE(rbym_read_one(nrbody))
 
 
      IF (gene_max > 1) THEN
          ALLOCATE(result(nrbody))
      ENDIF
 
      stack=0
 
      DO k=1,gene_max
          CALL hm_get_int_array_2indexes('start'  ,start_gene,jclause,k,is_available,lsubmodel)
          CALL hm_get_int_array_2indexes('end'    ,end_gene  ,jclause,k,is_available,lsubmodel)
          CALL hm_get_int_array_2indexes('by'     ,incr_gene ,jclause,k,is_available,lsubmodel)
 
          ! set value by default for increment to 1
          IF (incr_gene == 0) incr_gene = 1
 
          rstart = set_usrtos_nearest(start_gene,irbodym,nrbody,1)
          rstop  = set_usrtos_nearest(end_gene,irbodym,nrbody,2)
 
          stack_one=0
          DO r=rstart, rstop
             r1 = irbodym(r,1)
             IF ( mod( r1-start_gene , incr_gene) == 0)THEN
                stack_one = stack_one+1
                rbym_read_one(stack_one) = irbodym(r,2)
             ENDIF
          ENDDO
 
          IF (stack==0) THEN
             rbym_read_tmp(1:stack_one) = rbym_read_one(1:stack_one)
             stack = stack_one
          ELSE
             ! This code will not go if  GENE_MAX == 1 / Result does not need to be allocated
             CALL union_2_sorted_sets( rbym_read_tmp, stack ,
     *                                 rbym_read_one, stack_one ,
     *                                 result,        nb_result )
 
            rbym_read_tmp(1:nb_result) = result(1:nb_result)
            stack = nb_result
          ENDIF
      ENDDO
 
      clause%NB_RBODY = stack
      ALLOCATE(clause%RBODY(stack))
      clause%RBODY(1:stack) = rbym_read_tmp(1:stack)
C-------------------------
      DEALLOCATE (rbym_read_tmp)
      DEALLOCATE (rbym_read_one)
      IF (ALLOCATED(result)) DEALLOCATE (result)
C-------------------------
      END