setlenwa.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/.
!||====================================================================
!||    setlenwa   ../starter/source/restart/ddsplit/setlenwa.F
!||--- called by ------------------------------------------------------
!||    ddsplit    ../starter/source/restart/ddsplit/ddsplit.F
!||    lectur     ../starter/source/starter/lectur.F
!||====================================================================
      SUBROUTINE setlenwa(
     1      LENWA0 ,NTHWA0  ,NAIRWA0   ,NUMELS0  ,NUMELQ0,
     2      NUMELC0,NUMELTG0,NUMELT0   ,NUMELP0  ,NUMELR0,
     3      NUMNOD0,NMNT0   ,L_MUL_LAG1,L_MUL_LAG,MAXNX0 ,
     4      LWASPH0,NUMSPH0   ,LWASPIO  , NRCVVOIS0 ,
     5      LWAMP_L, LWANMP_L ,ITHERM   )
C Attention les dummy peuvent etre soit locaux soit globaux suivant appel
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      "com09_c.inc"
#include      "sphcom.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(INOUT) :: LWAMP_L
      INTEGER, INTENT(INOUT) :: LWANMP_L
      INTEGER, INTENT(IN)    :: ITHERM
      INTEGER LENWA0,NTHWA0,NAIRWA0,
     .        NUMELS0,NUMELQ0,NUMELC0,NUMELTG0,NUMELT0,NUMELP0,
     .        numelr0,numnod0,nmnt0,l_mul_lag1,l_mul_lag,maxnx0,
     .        lwasph0,numsph0,lwaspio,nrcvvois0
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NFLUWA, NBHOLWA, NAIRWA, NINTWA, NXELWA, NSPHWA,
     .        NASSNMP, NASSMP, NIDELWA, NWASPH, MSPR,
     .        IUN,LENVOISWA
C-----------------------------------------------
      IUN = 1
C---------------------------------------------------------------------
C SPH
      nsphwa = max(2*kvoisph,4*kvoisph,
     .             16*numsph0,lwaspio)   ! 19*NUMSPH0 enleve (tri en auto)
C work array WASPH, prive au SPH.
      lwasph0=17*numsph0+min(iun,nsphio)*3*numsph0
     .                            +min(iun,nsphsol)*numsph0
C---------------------------------------------------------------------
C ALE / CFD
      nfluwa  = max(iale,itherm,ieuler) * max(1,nmult) * 
     .          max(7 * (numels0+numelq0),3 * numnod0)
      nbholwa = max(iale,itherm,ieuler) * nmult * numelq0
      lenvoiswa = 3 * ( numnod0 + nrcvvois0 ) 
C AirBag 
      nairwa  = max(nairwa0, 4 * (numelc0 + numeltg0 + ibagsurf))   ! en toute rigueur NAIRWA0 et IBAGSURF locaux
C Interface
      nintwa  = nmnt0+l_mul_lag1                                    ! en toute rigueur NMNT local
C
C Assemblage SPMD : tout en dynamique
      nassnmp = 0
      nassmp  = 0
C shooting nodes
      nidelwa = 3*numnod0
C lenwa non parallele
      lwanmp_l   = max(nthwa0,nassnmp,l_mul_lag)
C---------------------------------------------------------------------
      nxelwa   = 22*maxnx0                                         ! en toute rigueur MAXNX local
C lenwa parallele
      lwamp_l    = max(nairwa,nintwa,nfluwa+nbholwa,nassmp,nxelwa,nsphwa,
     +               nidelwa,lenvoiswa)
C
C MSPR utilise dans le starter
      mspr    = max(4*numelr0,nparg*ngroup)
C
      lenwa0  = max(lwamp_l,lwanmp_l,mspr)
C
      RETURN
      END
C
!||====================================================================
!||    setmulti   ../starter/source/restart/ddsplit/setlenwa.F
!||--- called by ------------------------------------------------------
!||    lectur     ../starter/source/starter/lectur.F
!||====================================================================
      SUBROUTINE setmulti(IPARI)
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      "spmd_c.inc"
#include      "param_c.inc"
#include      "sphcom.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IPARI(NPARI,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, MULTIMP, ICONT, ITY
C-----------------------------------------------
C
C MULTIMAX : parametre servant a allouer les buffer de comm interface engine (SPMD)
C
      multimax = 0
      icont = 0
      DO i = 1, ninter
        multimp = ipari(23,i)
          IF(multimp<=4) THEN
            multimax = max(multimax,4*multimp)
          ELSEIF(multimp<=12)THEN
            multimax= max(multimax,3*multimp)
          ELSEIF(multimp>12)THEN
            multimax= max(multimax,2*multimp)
          END IF
        ity = ipari(7,i)
        IF(ity==7.OR.ity==10.OR.ity==11.OR. ity==22)icont=1
      END DO
      IF(numsph/=0.AND.multimax==0)multimax=1
      IF(numelx/=0.AND.multimax==0)multimax=4
      IF(ninter>0.AND.icont==0)multimax=100
 
C
      RETURN
      END