tetra4v.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/.
!||====================================================================
!||    tetra4v   ../engine/source/elements/solid/solide/tetra4v.F
!||--- called by ------------------------------------------------------
!||    deges4v   ../engine/source/elements/solid/solide/deges4v.F
!||====================================================================
      SUBROUTINE tetra4v(DET,
     .   X1, X2, X3, X4, 
     .   Y1, Y2, Y3, Y4,
     .   Z1, Z2, Z3, Z4)
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      my_real
     .   det,
     .   x1, x2, x3, x4, y1, y2, y3, y4,z1, z2, z3, z4
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
C                                                                     12
      my_real
     .   x41, y41, z41, x42, y42, z42, x43, y43, z43
      my_real
     .   rx, ry, rz, sx, sy, sz,tx, ty, tz,
     .   tsx, tsy,tsz
C-----------------------------------------------
       x43 = x4 - x3
       y43 = y4 - y3
       z43 = z4 - z3
       x41 = x4 - x1
       y41 = y4 - y1
       z41 = z4 - z1
       x42 = x4 - x2
       y42 = y4 - y2
       z42 = z4 - z2
C
       rx =  -x41
       ry =  -y41
       rz =  -z41
C
       sx =  -x42
       sy =  -y42
       sz =  -z42
C
       tx =  -x43
       ty =  -y43
       tz =  -z43
C
       tsx =  y43*z42 - y42*z43
       tsy =  z43*x42 - z42*x43
       tsz =  x43*y42 - x42*y43
C------V=r*(s^t)/6
       det = abs(rx*tsx + ry*tsy + rz*tsz)*one_over_6
C      
      RETURN
      END