suderi3.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/.
!||====================================================================
!||    suderi3   ../starter/source/elements/elbuf_init/suderi3.F
!||--- called by ------------------------------------------------------
!||    suinit3   ../starter/source/elements/elbuf_init/suinit3.F
!||====================================================================
      SUBROUTINE suderi3(NEL   ,VOL  ,
     .              X1   ,X2   ,X3   ,X4   ,X5   ,X6   ,X7   ,X8   ,
     .              Y1   ,Y2   ,Y3   ,Y4   ,Y5   ,Y6   ,Y7   ,Y8   ,
     .              Z1   ,Z2   ,Z3   ,Z4   ,Z5   ,Z6   ,Z7   ,Z8)
C------------------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NEL
      my_real
     .   VOL(*), 
     .   x1(*), x2(*), x3(*), x4(*), x5(*), x6(*),
     .   x7(*), x8(*), y1(*), y2(*), y3(*), y4(*), y5(*), y6(*), y7(*),
     .   y8(*), z1(*), z2(*), z3(*), z4(*), z5(*), z6(*), z7(*), z8(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J
      
      DOUBLE PRECISION
     .   xd1(mvsiz), xd2(mvsiz), xd3(mvsiz), xd4(mvsiz),
     .   xd5(mvsiz), xd6(mvsiz), xd7(mvsiz), xd8(mvsiz),
     .   yd1(mvsiz), yd2(mvsiz), yd3(mvsiz), yd4(mvsiz),
     .   yd5(mvsiz), yd6(mvsiz), yd7(mvsiz), yd8(mvsiz),
     .   zd1(mvsiz), zd2(mvsiz), zd3(mvsiz), zd4(mvsiz),
     .   zd5(mvsiz), zd6(mvsiz), zd7(mvsiz), zd8(mvsiz)
C     REAL
      my_real
     .   jac1(mvsiz), jac2(mvsiz), jac3(mvsiz), jac4(mvsiz), jac5(mvsiz), jac6(mvsiz),
     .   jac7(mvsiz), jac8(mvsiz) , jac9(mvsiz),
     .   a17(mvsiz) , a28(mvsiz) ,
     .   b17(mvsiz) , b28(mvsiz) ,
     .   c17(mvsiz) , c28(mvsiz) 
      my_real
     .   x17(mvsiz), x28(mvsiz), x35(mvsiz), x46(mvsiz),
     .   y17(mvsiz), y28(mvsiz), y35(mvsiz), y46(mvsiz),
     .   z17(mvsiz), z28(mvsiz), z35(mvsiz), z46(mvsiz),
     .   jac_59_68(mvsiz), jac_67_49(mvsiz), jac_48_57(mvsiz)
C=======================================================================
      DO i=1,nel
        xd1(i)=x1(i)
        xd2(i)=x2(i)
        xd3(i)=x3(i)
        xd4(i)=x4(i)
        xd5(i)=x5(i)
        xd6(i)=x6(i)
        xd7(i)=x7(i)
        xd8(i)=x8(i)
        yd1(i)=y1(i)
        yd2(i)=y2(i)
        yd3(i)=y3(i)
        yd4(i)=y4(i)
        yd5(i)=y5(i)
        yd6(i)=y6(i)
        yd7(i)=y7(i)
        yd8(i)=y8(i)
        zd1(i)=z1(i)
        zd2(i)=z2(i)
        zd3(i)=z3(i)
        zd4(i)=z4(i)
        zd5(i)=z5(i)
        zd6(i)=z6(i)
        zd7(i)=z7(i)
        zd8(i)=z8(i)
      ENDDO 
 
      DO i=1,nel
        x17(i)=xd7(i)-xd1(i)
        x28(i)=xd8(i)-xd2(i)
        x35(i)=xd5(i)-xd3(i)
        x46(i)=xd6(i)-xd4(i)
        y17(i)=yd7(i)-yd1(i)
        y28(i)=yd8(i)-yd2(i)
        y35(i)=yd5(i)-yd3(i)
        y46(i)=yd6(i)-yd4(i)
        z17(i)=zd7(i)-zd1(i)
        z28(i)=zd8(i)-zd2(i)
        z35(i)=zd5(i)-zd3(i)
        z46(i)=zd6(i)-zd4(i)
      ENDDO 
C
      DO i=1,nel               
        jac1(i)=x17(i)+x28(i)-x35(i)-x46(i)
        jac2(i)=y17(i)+y28(i)-y35(i)-y46(i)
        jac3(i)=z17(i)+z28(i)-z35(i)-z46(i)
        a17(i)=x17(i)+x46(i)
        a28(i)=x28(i)+x35(i)
        b17(i)=y17(i)+y46(i)
        b28(i)=y28(i)+y35(i)
        c17(i)=z17(i)+z46(i)
        c28(i)=z28(i)+z35(i)
      ENDDO 
         
      DO i=1,nel
        jac4(i)=a17(i)+a28(i)
        jac5(i)=b17(i)+b28(i)
        jac6(i)=c17(i)+c28(i)
        jac7(i)=a17(i)-a28(i)
        jac8(i)=b17(i)-b28(i)
        jac9(i)=c17(i)-c28(i)
      ENDDO 
C
      DO i=1,nel
        jac_59_68(i)=jac5(i)*jac9(i)-jac6(i)*jac8(i)
        jac_67_49(i)=jac6(i)*jac7(i)-jac4(i)*jac9(i)
        jac_48_57(i)=jac4(i)*jac8(i)-jac5(i)*jac7(i)
      ENDDO
C
      DO i=1,nel
        vol(i)=one_over_64*(jac1(i)*jac_59_68(i)+jac2(i)*jac_67_49(i)+jac3(i)*jac_48_57(i))
      ENDDO     
C-----------
      RETURN
      END