intti2.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/.
!||====================================================================
!||    intti2          ../engine/source/interfaces/interf/intti2.F
!||--- called by ------------------------------------------------------
!||    resol           ../engine/source/engine/resol.F
!||--- calls      -----------------------------------------------------
!||    intti12v        ../engine/source/interfaces/interf/intti12.F
!||    intti2v         ../engine/source/interfaces/interf/intti2v.F
!||--- uses       -----------------------------------------------------
!||    intbufdef_mod   ../common_source/modules/interfaces/intbufdef_mod.F90
!||====================================================================
      SUBROUTINE intti2(IPARI,X     ,V      ,
     .                  A    ,VR    ,AR    ,KHIE  ,MS     ,
     .                  IN   ,WEIGHT,WA    ,SKEW  ,
     .                  INTBUF_TAB )
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE intbufdef_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"
#include      "param_c.inc"
#include      "task_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IPARI(NPARI,*), WEIGHT(*),KHIE
C     REAL
      my_real
     .   x(*), v(*), a(*),ar(*),vr(*), ms(*),in(*),
     .   wa(*),skew(*)
 
      TYPE(intbuf_struct_) INTBUF_TAB(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER N, NTY, JI, JB, NMN, NINT,ILAGM
      INTEGER :: ILEV,II,J
C-----------------------------------------------
C wa are no longer used!!!
C=======================================================================
      DO n=ninter,1,-1
        nty   = ipari(7,n)                                             
        ilagm = 0                                                     
        ilagm = ipari(33,n)                                           
        IF (ilagm/=0) cycle                                   
        IF (npari>=26)THEN                                          
          IF(ipari(26,n)/=khie) cycle                           
        ENDIF                                                         
        ji  =ipari(1,n)                                                 
        jb  =ipari(2,n)                                                 
        nmn =ipari(6,n)                                             
        nint=n                                                        
        IF (nty == 2) THEN                                            
          IF(ipari(5,n)>0) THEN !       check if NSN=IPARI(5,N)>0
        !       Optimization :
        !       If IPARI(5,N)=NSN=0, some variables are loaded in INTTI2V
        !       if the number of TYPE2 interface is important (>3000) and if the number of
        !       spmd domain is quite important (NSPMD>50), the initialisation time is important :
        !       for_array_copy_in and other initialisations represent up to 5% of total CPUTIME
        !       and break the scalability of the code
           CALL intti2v(                                               
     .          ipari(1,n),x         ,v        , 
     .          a         ,vr        ,ar        ,ms        ,in       , 
     .          weight    ,intbuf_tab(n) )
          ELSE
        !       WARNING : if NSN==0 and ILEV==2 or 4, ones needs to scale the mass 
                ilev = ipari(20,n)
                IF(ilev==2.OR.ilev==4) THEN
                        DO ii=1,nmn
                                j=intbuf_tab(n)%MSR(ii)
                                ms(j)=intbuf_tab(n)%NMAS(ii)
                        ENDDO
                ENDIF
          ENDIF
        ELSEIF (nty == 12) THEN                                       
          IF(ispmd==0)THEN                            
            CALL intti12v(                                            
     .           ipari(1,n),intbuf_tab(n)      ,x         ,v        , 
     .            a        ,ms       ,weight   ,wa        ,skew     ) 
          ENDIF                                                       
        ENDIF                                                         
      ENDDO
C-----------
      RETURN
      END