int18_law151_nsv_shift.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/.
!||====================================================================
!||    int18_law151_nsv_shift_mod   ../common_source/interf/int18_law151_nsv_shift.F
!||--- called by ------------------------------------------------------
!||    inintr                       ../starter/source/interfaces/interf1/inintr.F
!||    intfop2                      ../engine/source/interfaces/interf/intfop2.F
!||    inttri                       ../engine/source/interfaces/intsort/inttri.F
!||====================================================================
      module int18_law151_nsv_shift_mod
      contains
! ======================================================================================================================
!                                                   procedures
! ======================================================================================================================
!! \brief index shifting for law151
!! \details law151 has collocated numerical scheme node identifier is 1:NUMELS starting from NUMNOD+1 to NUMNOD+NUMELS
!||====================================================================
!||    int18_law151_nsv_shift   ../common_source/interf/int18_law151_nsv_shift.F
!||--- called by ------------------------------------------------------
!||    inintr                   ../starter/source/interfaces/interf1/inintr.F
!||    intfop2                  ../engine/source/interfaces/interf/intfop2.F
!||    inttri                   ../engine/source/interfaces/intsort/inttri.F
!||--- uses       -----------------------------------------------------
!||    groupdef_mod             ../common_source/modules/groupdef_mod.F
!||    intbufdef_mod            ../common_source/modules/interfaces/intbufdef_mod.F90
!||    multi_fvm_mod            ../common_source/modules/ale/multi_fvm_mod.F90
!||====================================================================
        subroutine int18_law151_nsv_shift(mode,itask,nthread,multi_fvm,ipari,intbuf_tab,npari,ninter,numnod,opt_int_id)
!$COMMENT
!       INT18_LAW151_NSV_SHIFT description
!       shift the secondary node ID
!       for /INT18 + LAW151
!
!       INT18_LAW151_NSV_SHIFT organization :
!       - at the beginning of inttri/intfop2, shift with +numnod
!       - at the end of inttri/intfop2, shift with -numnod
!$ENDCOMMENT
! ----------------------------------------------------------------------------------------------------------------------
!                                                   Modules
! ----------------------------------------------------------------------------------------------------------------------
          use intbufdef_mod
          use multi_fvm_mod
          use groupdef_mod
! ----------------------------------------------------------------------------------------------------------------------
!                                                   Implicit None
! ----------------------------------------------------------------------------------------------------------------------
          implicit none
! ----------------------------------------------------------------------------------------------------------------------
!                                                   Arguments
! ----------------------------------------------------------------------------------------------------------------------
          integer,intent(in) :: npari,ninter,numnod,nthread
          character(len=*), intent(in) :: mode
          integer, intent(in) :: itask
          type(intbuf_struct_),dimension(ninter), intent(inout) :: intbuf_tab
          integer, dimension(npari,*), intent(in) :: ipari
          type(multi_fvm_struct) :: multi_fvm
          integer,intent(in),optional :: opt_int_id
! ----------------------------------------------------------------------------------------------------------------------
!                                                   Arguments
! ----------------------------------------------------------------------------------------------------------------------
          INTEGER :: II,N,NN
          INTEGER :: NSN
          INTEGER :: NSNF,NSNL
          INTEGER :: INTER_ID
C-----------------------------------------------
          inter_id = 0
          IF(PRESENT(opt_int_id))THEN
            inter_id = max(0,opt_int_id)
          ENDIF
          DO nn=1,multi_fvm%NUMBER_INT18
            n = multi_fvm%INT18_LIST(nn)
            IF(inter_id > 0)THEN
              IF(inter_id /= n) cycle
            ENDIF
            nsn = ipari(5,n)    ! number of secondary nodes
            nsnf = 1 + itask * nsn / nthread
            nsnl = (1 + itask) * nsn / nthread
            !   add/substract NUMNOD to NSV for phantom nodes
            !   if NSV > 0 --> local phantom node
            !   else if NSV < 0 --> reote node
            IF(mode=='+') THEN
              DO ii = nsnf,nsnl
                IF(intbuf_tab(n)%NSV(ii)>0) intbuf_tab(n)%NSV(ii) = intbuf_tab(n)%NSV(ii) + numnod
              ENDDO
            ELSE
              DO ii = nsnf,nsnl
                IF(intbuf_tab(n)%NSV(ii)>0) intbuf_tab(n)%NSV(ii) = intbuf_tab(n)%NSV(ii) - numnod
              ENDDO
            ENDIF
          ENDDO
 
          RETURN
        end subroutine int18_law151_nsv_shift
 
!===============================================================================
      end module int18_law151_nsv_shift_mod