icinv3.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/.
C
!||====================================================================
!||    icinv3        ../starter/source/interfaces/inter3d1/icinv3.F
!||--- called by ------------------------------------------------------
!||    i12chk3       ../starter/source/interfaces/inter3d1/i12chk3.F
!||--- uses       -----------------------------------------------------
!||====================================================================
      SUBROUTINE icinv3(IXS,IXC,IXTG,IADD,INVC)
      use element_mod , only :nixs,nixc,nixtg
C
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"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IXS(NIXS,*), IXC(NIXC,*), IXTG(NIXTG,*),
     .        IADD(*), INVC(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C---------------------`-------------------------
      INTEGER I, J, N
C-----------------------------------------------
      iadd(1) = 1
      iadd(2) = 1
C
      DO i=3,numnod+1
          iadd(i)=0
      ENDDO
C
      DO j=2,5
        DO i=1,numelc
          n = ixc(j,i) + 2
          iadd(n)=iadd(n)+1
        ENDDO
      ENDDO
C
      DO j=2,4
        DO i=1,numeltg
          n = ixtg(j,i) + 2
          iadd(n)=iadd(n)+1
        ENDDO
      ENDDO
C
      DO j=2,9
        DO i=1,numels
          n = ixs(j,i) + 2
          iadd(n)=iadd(n)+1
        ENDDO
      ENDDO
C
      DO i=3,numnod+1
          iadd(i)=iadd(i)+iadd(i-1)
      ENDDO
C
      DO i=1,numelc
 
          n = ixc(2,i) + 1
          invc(iadd(n)) = i
          iadd(n) = iadd(n) + 1
 
          n = ixc(3,i) + 1
          invc(iadd(n)) = i
          iadd(n) = iadd(n) + 1
 
          n = ixc(4,i) + 1
          invc(iadd(n)) = i
          iadd(n) = iadd(n) + 1
 
          n = ixc(5,i) + 1
          invc(iadd(n)) = i
          iadd(n) = iadd(n) + 1
 
      ENDDO
C
      DO i=1,numeltg
          n = ixtg(2,i) + 1
          invc(iadd(n)) = i+numelc
          iadd(n) = iadd(n) + 1
 
          n = ixtg(3,i) + 1
          invc(iadd(n)) = i+numelc
          iadd(n) = iadd(n) + 1
 
          n = ixtg(4,i) + 1
          invc(iadd(n)) = i+numelc
          iadd(n) = iadd(n) + 1
 
      ENDDO
C
      DO i=1,numels
 
          n = ixs(2,i) + 1
          invc(iadd(n)) = i+numelc+numeltg
          iadd(n) = iadd(n) + 1
 
          n = ixs(3,i) + 1
          invc(iadd(n)) = i+numelc+numeltg
          iadd(n) = iadd(n) + 1
 
          n = ixs(4,i) + 1
          invc(iadd(n)) = i+numelc+numeltg
          iadd(n) = iadd(n) + 1
 
          n = ixs(5,i) + 1
          invc(iadd(n)) = i+numelc+numeltg
          iadd(n) = iadd(n) + 1
 
          n = ixs(6,i) + 1
          invc(iadd(n)) = i+numelc+numeltg
          iadd(n) = iadd(n) + 1
 
          n = ixs(7,i) + 1
          invc(iadd(n)) = i+numelc+numeltg
          iadd(n) = iadd(n) + 1
 
          n = ixs(8,i) + 1
          invc(iadd(n)) = i+numelc+numeltg
          iadd(n) = iadd(n) + 1
 
          n = ixs(9,i) + 1
          invc(iadd(n)) = i+numelc+numeltg
          iadd(n) = iadd(n) + 1
      ENDDO
C-------------------------------------
C     APRES DO 70 
C     IADD(1:NUMNOD) <=> IADD(2:NUMNOD+1)
C------------------------------------- 
      RETURN
      END