inter_deallocate_wait.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/.
!||====================================================================
!||    inter_deallocate_wait   ../engine/source/interfaces/generic/inter_deallocate_wait.F
!||--- called by ------------------------------------------------------
!||    inttri                  ../engine/source/interfaces/intsort/inttri.f
!||--- calls      -----------------------------------------------------
!||    my_barrier              ../engine/source/system/machine.F
!||--- uses       -----------------------------------------------------
!||    h3d_mod                 ../engine/share/modules/h3d_mod.F
!||    intbufdef_mod           ../common_source/modules/interfaces/intbufdef_mod.F90
!||    inter_sorting_mod       ../engine/share/modules/inter_sorting_mod.F
!||    inter_struct_mod        ../engine/share/modules/inter_struct_mod.F
!||    metric_mod              ../common_source/modules/interfaces/metric_mod.F
!||    multi_fvm_mod           ../common_source/modules/ale/multi_fvm_mod.F90
!||    sensor_mod              ../common_source/modules/sensor_mod.F90
!||====================================================================
        SUBROUTINE inter_deallocate_wait( ITASK,NB_INTER_SORTED,LIST_INTER_SORTED,IPARI,
     1                          NSENSOR,IRECVFROM,SENSOR_TAB,INTER_STRUCT,SORT_COMM )
!$COMMENT
!       INTER_DEALLOCATE_WAIT description
!       wait the message "send the secondary nodes" & deallocation of array
!       INTER_DEALLOCATE_WAIT organization :
!$ENDCOMMENT
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE multi_fvm_mod
      USE h3d_mod
      USE metric_mod
      USE intbufdef_mod  
      USE inter_struct_mod
      USE inter_sorting_mod
      USE sensor_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      "com01_c.inc"
#include      "com04_c.inc"
#include      "com08_c.inc"
#include      "param_c.inc"
#include      "task_c.inc"
#include "spmd.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(in) :: ITASK    !   omp thread ID
      INTEGER, INTENT(in) :: NB_INTER_SORTED        !   number of interfaces that need to be sorted
      INTEGER, INTENT(in) :: NSENSOR
      INTEGER, DIMENSION(NB_INTER_SORTED), INTENT(in) :: LIST_INTER_SORTED   !   list of interfaces that need to be sorted
      INTEGER, DIMENSION(NPARI,NINTER), INTENT(in) :: IPARI    !   interface data
      INTEGER, DIMENSION(NINTER+1,NSPMD+1), INTENT(in) :: IRECVFROM 
      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
      TYPE (SENSOR_STR_) ,DIMENSION(NSENSOR) :: SENSOR_TAB
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
        INTEGER :: KK,N,P
        INTEGER :: NTY,INACTI,ISENS
        INTEGER :: NSN,NMN
        LOGICAL :: TYPE18
#ifdef MPI
        INTEGER :: IERROR,STATUS(MPI_STATUS_SIZE)
#endif
        my_real :: ts
!   ----------------------------------------
        IF(itask==0.AND.nb_inter_sorted>0) THEN
            DO kk=1,nb_inter_sorted
                n = list_inter_sorted(kk)
 
                nty = ipari(7,n)
                inacti = ipari(22,n)
                type18=.false.
                IF(nty==7.AND.inacti==7)type18=.true.
 
                isens = 0
                IF(nty == 7.OR.nty == 11.OR.nty == 24.OR.nty == 25) isens = ipari(64,n)  
                IF (isens > 0) THEN    ! IF INTERFACE IS ACTIVATED BY SENSOR 
                    ts = sensor_tab(isens)%TSTART
                ELSE
                    ts = tt
                ENDIF
 
                nsn = ipari(5,n)
                nmn = ipari(6,n)
                ! -----------------------------------------------------
                ! type 7
                IF((nty==7.AND.tt>=ts).AND.(.NOT.type18))THEN
                    ! wait the message only for proc with only secondary nodes 
                    IF(nsn>0.AND.nmn==0) THEN   !   local nsn>0 && local nmn == 0
                        DO p = 1, nspmd
                            IF(irecvfrom(n,p)/=0) THEN     !   nmn >0
                                IF(p/=ispmd+1) THEN
                                    IF(sort_comm(n)%NB(p)/=0) THEN
#ifdef MPI
                                        CALL mpi_wait(sort_comm(n)%REQ_SD2(p),status,ierror) 
                                        DEALLOCATE(sort_comm(n)%DATA_PROC(p)%RBUF) 
                                        CALL mpi_wait(sort_comm(n)%REQ_SD3(p),status,ierror) 
                                        DEALLOCATE(sort_comm(n)%DATA_PROC(p)%IBUF)  
                                        sort_comm(n)%NB(p) = 0 
#endif                     
                                    END IF
                                ENDIF
                            ENDIF
                        ENDDO
                    ENDIF
                ENDIF
                ! -----------------------------------------------------
                ! deallocation
                IF(ALLOCATED(sort_comm(n)%COARSE_GRID) )DEALLOCATE( sort_comm(n)%COARSE_GRID )
                IF(ALLOCATED(sort_comm(n)%MAIN_COARSE_GRID) )DEALLOCATE( sort_comm(n)%MAIN_COARSE_GRID )
                IF(ALLOCATED(sort_comm(n)%GLOBAL_COARSE_CELL) )DEALLOCATE( sort_comm(n)%GLOBAL_COARSE_CELL )
 
                IF(ALLOCATED(sort_comm(n)%SEND_SIZE_COARSE_CELL) ) DEALLOCATE( sort_comm(n)%SEND_SIZE_COARSE_CELL )
                IF(ALLOCATED(sort_comm(n)%RCV_SIZE_COARSE_CELL) ) DEALLOCATE( sort_comm(n)%RCV_SIZE_COARSE_CELL )
                IF(ALLOCATED(sort_comm(n)%SEND_DISPLS_COARSE_CELL) ) DEALLOCATE( sort_comm(n)%SEND_DISPLS_COARSE_CELL )
                IF(ALLOCATED(sort_comm(n)%RCV_DISPLS_COARSE_CELL) ) DEALLOCATE( sort_comm(n)%RCV_DISPLS_COARSE_CELL )
               
 
                IF(ALLOCATED(sort_comm(n)%KEEP_PROC) )DEALLOCATE( sort_comm(n)%KEEP_PROC )
                IF(ALLOCATED(sort_comm(n)%CELL_LIST) )DEALLOCATE( sort_comm(n)%CELL_LIST )
                IF(ALLOCATED(sort_comm(n)%NB_CELL_PROC) )DEALLOCATE( sort_comm(n)%NB_CELL_PROC )
                IF(ALLOCATED(sort_comm(n)%CELL) )DEALLOCATE( sort_comm(n)%CELL )
 
                IF(ALLOCATED(sort_comm(n)%SEND_NB_CELL) ) DEALLOCATE( sort_comm(n)%SEND_NB_CELL )
                IF(ALLOCATED(sort_comm(n)%RCV_NB_CELL) ) DEALLOCATE( sort_comm(n)%RCV_NB_CELL )
                IF(ALLOCATED(sort_comm(n)%SEND_DISPLS_NB_CELL) ) DEALLOCATE( sort_comm(n)%SEND_DISPLS_NB_CELL )
                IF(ALLOCATED(sort_comm(n)%RCV_DISPLS_NB_CELL) ) DEALLOCATE( sort_comm(n)%RCV_DISPLS_NB_CELL )
 
                IF(ALLOCATED(sort_comm(n)%SEND_SIZE_CELL) ) DEALLOCATE( sort_comm(n)%SEND_SIZE_CELL )
                IF(ALLOCATED(sort_comm(n)%RCV_SIZE_CELL) ) DEALLOCATE( sort_comm(n)%RCV_SIZE_CELL )
                IF(ALLOCATED(sort_comm(n)%SEND_DISPLS_CELL) ) DEALLOCATE( sort_comm(n)%SEND_DISPLS_CELL )
                IF(ALLOCATED(sort_comm(n)%RCV_DISPLS_CELL) ) DEALLOCATE( sort_comm(n)%RCV_DISPLS_CELL )
 
                IF(ALLOCATED(sort_comm(n)%REQUEST_NB_S) ) DEALLOCATE( sort_comm(n)%REQUEST_NB_S )
                IF(ALLOCATED(sort_comm(n)%REQUEST_NB_R) ) DEALLOCATE( sort_comm(n)%REQUEST_NB_R )
                IF(ALLOCATED(sort_comm(n)%NB) ) DEALLOCATE( sort_comm(n)%NB )
                IF(ALLOCATED(sort_comm(n)%INDEX_RCV) ) DEALLOCATE( sort_comm(n)%INDEX_RCV )
 
                IF(ALLOCATED(sort_comm(n)%SEND_NB) ) DEALLOCATE( sort_comm(n)%SEND_NB )
                IF(ALLOCATED(sort_comm(n)%RECV_NB) ) DEALLOCATE( sort_comm(n)%RECV_NB )
 
                IF(ALLOCATED(sort_comm(n)%REQUEST_CELL_SEND) ) DEALLOCATE( sort_comm(n)%REQUEST_CELL_SEND )
                IF(ALLOCATED(sort_comm(n)%REQUEST_CELL_RCV) ) DEALLOCATE( sort_comm(n)%REQUEST_CELL_RCV )
 
                IF(ALLOCATED(sort_comm(n)%IIX) )DEALLOCATE( sort_comm(n)%IIX )
                IF(ALLOCATED(sort_comm(n)%IIY) )DEALLOCATE( sort_comm(n)%IIY )
                IF(ALLOCATED(sort_comm(n)%IIZ) )DEALLOCATE( sort_comm(n)%IIZ )
                IF(ALLOCATED(sort_comm(n)%LAST_NOD) )DEALLOCATE( sort_comm(n)%LAST_NOD )
                IF(ALLOCATED(sort_comm(n)%NEXT_NOD) )DEALLOCATE( sort_comm(n)%NEXT_NOD )
                IF(ALLOCATED(sort_comm(n)%VOXEL) )DEALLOCATE( sort_comm(n)%VOXEL )
 
                IF(ALLOCATED(inter_struct(n)%CURV_MAX)) DEALLOCATE(inter_struct(n)%CURV_MAX)
                ! -----------------------------------------------------
            ENDDO
            IF(ALLOCATED(nb_local_cell) )DEALLOCATE( nb_local_cell )
            IF(ALLOCATED(cell_bool) )DEALLOCATE( cell_bool )
        ENDIF
        CALL my_barrier()
 
        RETURN
        END SUBROUTINE inter_deallocate_wait