admnorm.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/.
!||====================================================================
!||    admnorm4      ../engine/source/model/remesh/admnorm.F
!||--- called by ------------------------------------------------------
!||    admdiv        ../engine/source/model/remesh/admdiv.F
!||--- uses       -----------------------------------------------------
!||    element_mod   ../common_source/modules/elements/element_mod.F90
!||    remesh_mod    ../engine/share/modules/remesh_mod.F
!||====================================================================
      SUBROUTINE admnorm4(N,IXC,X)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE remesh_mod
      use element_mod , only : nixc
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER N, IXC(NIXC,*)
      my_real
     .        x(3,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER SON,LEVEL,IP,NLEV,N1,N2,N3,N4
      my_real
     .        x13,y13,z13,x24,y24,z24,nx,ny,nz,aaa
C-----------------------------------------------
      n1=ixc(2,n)
      n2=ixc(3,n)
      n3=ixc(4,n)
      n4=ixc(5,n)
 
      tagnod(n1)=1
      tagnod(n2)=1
      tagnod(n3)=1
      tagnod(n4)=1
 
      x13 = x(1,n3) - x(1,n1)
      y13 = x(2,n3) - x(2,n1)
      z13 = x(3,n3) - x(3,n1)
 
      x24 = x(1,n4) - x(1,n2)
      y24 = x(2,n4) - x(2,n2)
      z24 = x(3,n4) - x(3,n2)
 
      nx = y13*z24 - z13*y24
      ny = z13*x24 - x13*z24
      nz = x13*y24 - y13*x24
 
      aaa=one/max(em30,sqrt(nx*nx+ny*ny+nz*nz))
      nx = nx * aaa
      ny = ny * aaa
      nz = nz * aaa
 
      nodnorm(1,n1)=nodnorm(1,n1)+nx
      nodnorm(2,n1)=nodnorm(2,n1)+ny
      nodnorm(3,n1)=nodnorm(3,n1)+nz
C  
      nodnorm(1,n2)=nodnorm(1,n2)+nx
      nodnorm(2,n2)=nodnorm(2,n2)+ny
      nodnorm(3,n2)=nodnorm(3,n2)+nz
C  
      nodnorm(1,n3)=nodnorm(1,n3)+nx
      nodnorm(2,n3)=nodnorm(2,n3)+ny
      nodnorm(3,n3)=nodnorm(3,n3)+nz
C  
      nodnorm(1,n4)=nodnorm(1,n4)+nx
      nodnorm(2,n4)=nodnorm(2,n4)+ny
      nodnorm(3,n4)=nodnorm(3,n4)+nz
 
      RETURN
      END     
 
!||====================================================================
!||    admnorm3      ../engine/source/model/remesh/admnorm.F
!||--- called by ------------------------------------------------------
!||    admdiv        ../engine/source/model/remesh/admdiv.f
!||--- uses       -----------------------------------------------------
!||    element_mod   ../common_source/modules/elements/element_mod.F90
!||    remesh_mod    ../engine/share/modules/remesh_mod.F
!||====================================================================
      SUBROUTINE admnorm3(N,IXTG,X)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE remesh_mod
      use element_mod , only : nixtg
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER N, IXTG(NIXTG,*)
      my_real
     .        x(3,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER N1,N2,N3
      my_real
     .       x12,y12,z12,x13,y13,z13,nx,ny,nz,aaa
C-----------------------------------------------
      n1=ixtg(2,n)
      n2=ixtg(3,n)
      n3=ixtg(4,n)
 
      tagnod(n1)=1
      tagnod(n2)=1
      tagnod(n3)=1
 
      x13 = x(1,n3) - x(1,n1)
      y13 = x(2,n3) - x(2,n1)
      z13 = x(3,n3) - x(3,n1)
 
      x12 = x(1,n2) - x(1,n1)
      y12 = x(2,n2) - x(2,n1)
      z12 = x(3,n2) - x(3,n1)
 
      nx = y12*z13 - z12*y13
      ny = z12*x13 - x12*z13
      nz = x12*y13 - y12*x13
 
      aaa=one/max(em30,sqrt(nx*nx+ny*ny+nz*nz))
      nx = nx * aaa
      ny = ny * aaa
      nz = nz * aaa
 
      nodnorm(1,n1)=nodnorm(1,n1)+nx
      nodnorm(2,n1)=nodnorm(2,n1)+ny
      nodnorm(3,n1)=nodnorm(3,n1)+nz
C  
      nodnorm(1,n2)=nodnorm(1,n2)+nx
      nodnorm(2,n2)=nodnorm(2,n2)+ny
      nodnorm(3,n2)=nodnorm(3,n2)+nz
C  
      nodnorm(1,n3)=nodnorm(1,n3)+nx
      nodnorm(2,n3)=nodnorm(2,n3)+ny
      nodnorm(3,n3)=nodnorm(3,n3)+nz
 
      RETURN
      END