spmd_anim_ply_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_anim_ply_init   ../engine/source/mpi/anim/spmd_anim_ply_init.F
!||--- called by ------------------------------------------------------
!||    resol_init           ../engine/source/engine/resol_init.F
!||--- calls      -----------------------------------------------------
!||--- uses       -----------------------------------------------------
!||    plyxfem_mod          ../engine/share/modules/plyxfem_mod.F
!||    spmd_mod             ../engine/source/mpi/spmd_mod.F90
!||====================================================================
      SUBROUTINE spmd_anim_ply_init()
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE plyxfem_mod
      USE spmd_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------------------------
C   M e s s a g e   P a s s i n g
C-----------------------------------------------
#include "spmd.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "param_c.inc"
#include      "task_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
#ifdef MPI
      INTEGER, DIMENSION(:),ALLOCATABLE :: PLYNODTAG,SENDBUF,RECBUF
      INTEGER MAXSENDSZ,RECSZ,PLYSZ
      INTEGER I,J,P,ND, ITAG, IDCMAXP(NSPMD-1), REQ(NSPMD-1),NODCOUNT
      INTEGER MSGOFF,MSGOFF2,STAT(MPI_STATUS_SIZE,NSPMD-1), IERR
C-----------------------------------------------
      DATA msgoff/7054/
C-----------------------------------------------
C
C Allocate Global PLY Structure
C
      IF (ispmd==0)THEN
        ALLOCATE (plyiadnod(nplymax,nspmd))
        plyiadnod = 0
      ELSE
         ALLOCATE (plyiadnod(1,1))
      ENDIF
      plyiadnod = 0
 
 
C Need to rebuild a global PLYNOD Structure on processor 0 to have
C a conversion array from Global Node ID to Internal PLYNOD ID
 
      DO i=1,nplymax
        IF (ispmd==0)THEN
C PROC 0, RECEIVE NODES NUM for each ply
           DO p=2,nspmd
              itag=msgoff
              CALL mpi_recv(recsz,1, mpi_integer, it_spmd(p),
     .                      itag, spmd_comm_world, stat, ierr)
              plyiadnod(i,p) = recsz
           ENDDO
        ELSE
           itag=msgoff
           CALL mpi_send(plynod(i)%PLYNUMNODS,1,mpi_integer,it_spmd(1),
     .                   itag,spmd_comm_world,ierr)
        ENDIF
      ENDDO
 
C
#endif
      RETURN
      END