geom_hex_dihedral_angle.f90 Source File

!
!     (c) 2019 Guide Star Engineering, LLC
!     This Software was developed for the US Nuclear Regulatory Commission (US NRC)
!     under contract "Multi-Dimensional Physics Implementation into Fuel Analysis under
!     Steady-state and Transients (FAST)", contract # NRC-HQ-60-17-C-0007
!
!    NEMO - Numerical Engine (for) Multiphysics Operators
! Copyright (c) 2007, Stefano Toninel
!                     Gian Marco Bianchi  University of Bologna
!              David P. Schmidt    University of Massachusetts - Amherst
!              Salvatore Filippone University of Rome Tor Vergata
! All rights reserved.
!
! Redistribution and use in source and binary forms, with or without modification,
! are permitted provided that the following conditions are met:
!
!     1. Redistributions of source code must retain the above copyright notice,
!        this list of conditions and the following disclaimer.
!     2. Redistributions in binary form must reproduce the above copyright notice,
!        this list of conditions and the following disclaimer in the documentation
!        and/or other materials provided with the distribution.
!     3. Neither the name of the NEMO project nor the names of its contributors
!        may be used to endorse or promote products derived from this software
!        without specific prior written permission.
!
! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND
! ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
! WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
! DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR
! ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
! (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
! LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON
! ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
! (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
! SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
!
!---------------------------------------------------------------------------------
!
! $Id: geom_hex_dihedral_angle.f90 3093 2008-04-22 14:51:09Z sfilippo $
!
! Description:
!    Calculates the smallest dihedral angle (in radians).
!    Note that dihedral angle is equal to the angle between normal vectors.
!    WARNING! Only hexes are currently supported.
!
SUBMODULE(tools_mesh_basics) geom_hex_dihedral_angle_implementation
    USE class_vector, ONLY : vector, OPERATOR(.dot.)
    IMPLICIT NONE

    CONTAINS

        MODULE PROCEDURE geom_hex_dihedral_angle

            USE class_psblas, ONLY : psb_dpk_
            !USE class_vector, ONLY : vector, OPERATOR(.dot.)
            USE tools_math,   ONLY : Pi
            IMPLICIT NONE
            !
            INTEGER :: edge,if1, if2
            REAL(psb_dpk_) :: arg   !! argument of the acos function
            REAL(psb_dpk_) :: angle !! between two faces
            TYPE(vector)     :: unitnorm(6)

            unitnorm(:) = af(:)%unit()

            largest  = 0.d0 !! Initialize to smallest possible value
            smallest = pi   !! Initialize to largest possible value

            DO edge = 1,12  !! Loop over all twelve edges in a hex cell

                if1 = adjacent(edge,1)
                if2 = adjacent(edge,2)

                IF ( .NOT. ( (if1 == 0) .OR. (if2 == 0) ) ) THEN
                    ! cosine = dot product of two vectors/divided by the product of their magnitude

                    ! area contains the normal scaled so the mag. = the face area
                    arg = (unitnorm(if1) .dot. unitnorm(if2))

                    angle = ACOS(arg)

                    largest  = MAX(largest,angle)
                    smallest = MIN(smallest,angle)
                ENDIF

            ENDDO

        END PROCEDURE geom_hex_dihedral_angle

END SUBMODULE geom_hex_dihedral_angle_implementation