spmd_cell_size_exchange_init.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/.
!||====================================================================
!||    spmd_cell_size_exchange_init   ../engine/source/mpi/interfaces/spmd_cell_size_exchange_init.F
!||--- called by ------------------------------------------------------
!||    spmd_cell_list_exchange        ../engine/source/mpi/interfaces/spmd_cell_list_exchange.F
!||--- uses       -----------------------------------------------------
!||    inter_sorting_mod              ../engine/share/modules/inter_sorting_mod.F
!||    inter_struct_mod               ../engine/share/modules/inter_struct_mod.F
!||    tri7box                        ../engine/share/modules/tri7box.F
!||====================================================================
        SUBROUTINE spmd_cell_size_exchange_init(IRCVFROM,ISENDTO,IPARI,NIN,INTER_STRUCT,SORT_COMM)
!$COMMENT
!       SPMD_CELL_SIZE_EXCHANGE_INIT description :
!       initialization of buffer size + allocation
!       SPMD_CELL_SIZE_EXCHANGE_INIT organization :
!$ENDCOMMENT
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
        USE inter_sorting_mod
        USE inter_struct_mod
        USE tri7box
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      "com01_c.inc"
#include      "com04_c.inc"
#include      "task_c.inc"
#include      "sms_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
        INTEGER, INTENT(in) :: NIN
        INTEGER, DIMENSION(NPARI,NINTER), INTENT(in) ::  IPARI
        INTEGER, DIMENSION(NINTER+1,NSPMD+1), INTENT(in) :: ISENDTO,IRCVFROM
        TYPE(inter_struct_type), DIMENSION(NINTER), INTENT(inout) :: INTER_STRUCT   !   structure for interface
        TYPE(sorting_comm_type), DIMENSION(NINTER), INTENT(inout) :: SORT_COMM   ! structure for interface sorting comm
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C----------------------------------------------- 
#ifdef MPI        
        INTEGER :: I,J,NOD,L,L2,KK,IJK,KJI
        INTEGER :: P,P_LOC
        INTEGER :: ADRESS,SHIFT_
        INTEGER :: ISIZ,RSIZ,IDEB,JDEB
        INTEGER :: NSN,NMN,IGAP,INTTH,INTFRIC,ITYP,ITIED
        INTEGER :: IFQ,INACTI
        INTEGER :: LOC_PROC
!   --------------------------------------------------------------------     
 
        loc_proc = ispmd + 1
        sort_comm(nin)%RSIZ = 0
        sort_comm(nin)%ISIZ = 0
 
        igap = ipari(21,nin)
        intth = ipari(47,nin)
        intfric = ipari(72,nin)
        ityp = ipari(7,nin)
        itied = ipari(85,nin)
        nmn = ipari(6,nin)
        nsn = ipari(5,nin)
        inacti = ipari(22,nin)
        ifq =ipari(31,nin)
 
        IF(inacti==5.OR.inacti==6.OR.inacti==7.OR.ifq>0.OR.itied/=0) THEN
            IF(.NOT.ALLOCATED(inter_struct(nin)%NSNFIOLD) ) THEN
                ALLOCATE(inter_struct(nin)%NSNFIOLD(nspmd))
                inter_struct(nin)%NSNFIOLD(1:nspmd) = 0
            ENDIF
            DO p = 1, nspmd
                inter_struct(nin)%NSNFIOLD(p) = nsnfi(nin)%P(p)
            END DO
        END IF
        sort_comm(nin)%NBSEND_NB = 0
        sort_comm(nin)%NBRECV_NB = 0
 
        IF(ircvfrom(nin,loc_proc)/=0.OR.isendto(nin,loc_proc)/=0) THEN
 
            IF(.NOT.ALLOCATED(sort_comm(nin)%REQUEST_NB_R)) THEN
                ALLOCATE(sort_comm(nin)%REQUEST_NB_R(nspmd) )
            ENDIF
            IF(.NOT.ALLOCATED(sort_comm(nin)%REQUEST_NB_S)) THEN
                ALLOCATE(sort_comm(nin)%REQUEST_NB_S(nspmd) )
            ENDIF
            IF(.NOT.ALLOCATED(sort_comm(nin)%NB)) THEN
                ALLOCATE(sort_comm(nin)%NB(nspmd) )
                sort_comm(nin)%NB(1:nspmd) = 0
            ENDIF
            IF(.NOT.ALLOCATED(sort_comm(nin)%INDEX_RCV)) THEN
                ALLOCATE( sort_comm(nin)%INDEX_RCV(sort_comm(nin)%PROC_NUMBER) )
            ENDIF
        
            rsiz = 8    
            isiz = 6
 
c specific cases 
c IGAP=1 or IGAP=2
            IF(igap==1 .OR. igap==2)THEN
                rsiz = rsiz + 1
c IGAP=3
            ELSEIF(igap==3)THEN
                rsiz = rsiz + 2
            ENDIF
 
C thermic      
            IF(intth > 0 ) THEN    
                rsiz = rsiz + 2
                      isiz = isiz + 1
            ENDIF
C Friction      
            IF(intfric > 0 ) THEN    
                      isiz = isiz + 1
            ENDIF
 
C -- IDTMINS==2      
            IF(idtmins == 2)THEN     
                isiz = isiz + 2
C -- IDTMINS_INT /= 0
            ELSEIF(idtmins_int/=0)THEN    
                isiz = isiz + 1
            END IF
 
            sort_comm(nin)%RSIZ = rsiz
            sort_comm(nin)%ISIZ = isiz
 
            IF(isendto(nin,loc_proc)>0) THEN
                IF(.NOT.ALLOCATED(sort_comm(nin)%SEND_NB)) ALLOCATE(sort_comm(nin)%SEND_NB(nspmd))
            ENDIF
            IF(ircvfrom(nin,loc_proc)>0) THEN
                IF(.NOT.ALLOCATED(sort_comm(nin)%RECV_NB)) ALLOCATE(sort_comm(nin)%RECV_NB(nspmd))
            ENDIF
        ENDIF
#endif
        RETURN
        END SUBROUTINE spmd_cell_size_exchange_init