facnor.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/.
!||====================================================================
!||    facnor    ../engine/source/output/anim/generate/facnor.F
!||--- called by ------------------------------------------------------
!||    parsorc   ../engine/source/output/anim/generate/parsorc.F
!||====================================================================
      SUBROUTINE facnor(X,D,II,XNORM,CDG,INVERT)
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-----------------------------------------------
C     REAL
      my_real
     .   x(3,*),d(3,*),xnorm(3,*),cdg(*)
      INTEGER II(4),INVERT
C     REAL
      my_real
     .   x0(4),y0(4),z0(4),nx,ny,nz,xc,yc,zc,s
      INTEGER I,I1,I2,I3,I4
C-----------------------------------------------
C       FACE NORMAL
C-----------------------------------------------
      DO 100 i = 1,4
            x0(i) = x(1,ii(i))-d(1,ii(i))
            y0(i) = x(2,ii(i))-d(2,ii(i))
            z0(i) = x(3,ii(i))-d(3,ii(i))
 100  CONTINUE
C
      nx = (y0(3) - y0(1)) * (z0(4) - z0(2)) - 
     .     (z0(3) - z0(1)) * (y0(4) - y0(2))
      ny = (z0(3) - z0(1)) * (x0(4) - x0(2)) - 
     .     (x0(3) - x0(1)) * (z0(4) - z0(2))
      nz = (x0(3) - x0(1)) * (y0(4) - y0(2)) -
     .     (y0(3) - y0(1)) * (x0(4) - x0(2))
C
      xc = .25*(x0(1) + x0(2) + x0(3) + x0(4)) 
      yc = .25*(y0(1) + y0(2) + y0(3) + y0(4)) 
      zc = .25*(z0(1) + z0(2) + z0(3) + z0(4)) 
      xc = xc - cdg(1)
      yc = yc - cdg(2)
      zc = zc - cdg(3)
      s  = nx * xc + ny * yc + nz * zc
      invert = 1
c      IF(S<ZERO)THEN
c             INVERT = -1
c             I1    = II(2)
c             I2    = II(1)
c             I3    = II(4)
c             I4    = II(3)
c             II(1) = I1
c             II(2) = I2
c             II(3) = I3
c             II(4) = I4
c      ENDIF
C
      DO 200 i = 1,4
            x0(i) = x(1,ii(i))
            y0(i) = x(2,ii(i))
            z0(i) = x(3,ii(i))
 200  CONTINUE
C
      nx = (y0(3) - y0(1)) * (z0(4) - z0(2)) - 
     .     (z0(3) - z0(1)) * (y0(4) - y0(2))
      ny = (z0(3) - z0(1)) * (x0(4) - x0(2)) - 
     .     (x0(3) - x0(1)) * (z0(4) - z0(2))
      nz = (x0(3) - x0(1)) * (y0(4) - y0(2)) -
     .     (y0(3) - y0(1)) * (x0(4) - x0(2))
C
      DO 300 i = 1,4
          xnorm(1,ii(i)) = xnorm(1,ii(i)) + nx
          xnorm(2,ii(i)) = xnorm(2,ii(i)) + ny
          xnorm(3,ii(i)) = xnorm(3,ii(i)) + nz
 300  CONTINUE
C
C-----------------------------------------------
      RETURN
      END
C