SplitQuadFace.F90

Go to the documentation of this file.

! *********************************************************************
! * Rocstar Simulation Suite                                          *
! * Copyright@2015, Illinois Rocstar LLC. All rights reserved.        *
! *                                                                   *
! * Illinois Rocstar LLC                                              *
! * Champaign, IL                                                     *
! * www.illinoisrocstar.com                                           *
! * sales@illinoisrocstar.com                                         *
! *                                                                   *
! * License: See LICENSE file in top level of distribution package or *
! * http://opensource.org/licenses/NCSA                               *
! *********************************************************************
! *********************************************************************
! * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,   *
! * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES   *
! * OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND          *
! * NONINFRINGEMENT.  IN NO EVENT SHALL THE CONTRIBUTORS OR           *
! * COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER       *
! * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,   *
! * Arising FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE    *
! * USE OR OTHER DEALINGS WITH THE SOFTWARE.                          *
! *********************************************************************
!******************************************************************************
!
! Purpose: Indicate which way quad face should be split into two triangular 
!   faces.
!
! Description: None.
!
! Input:
!   xyz         Coordinates of the four vertices
!
! Output:
!   splitFlag   Flag indicating whether face should be split along diagonal
!               13 or 24. If splitFlag == FACE_SPLIT_13, face should be split
!               along diagonal 13, giving the two triangular faces 123 and 134.
!               IF splitFlag == FACE_SPLIT_24, face should be split along
!               diagonal 24, giving the two triangular faces 124 and 234.
!
! Notes: None.
!
!******************************************************************************
!
! $Id: SplitQuadFace.F90,v 1.3 2008/12/06 08:44:10 mtcampbe Exp $
!
! Copyright: (c) 2002 by the University of Illinois
!
!******************************************************************************

SUBROUTINE splitquadface( global,xyz1,xyz2,xyz3,xyz4,splitFlag )

  USE moddatatypes
  USE modparameters
  USE moderror
  USE modglobal, ONLY     : t_global
  USE modinterfaces, ONLY : facevectortria
   
  IMPLICIT NONE

! ... parameters
  TYPE(t_global), POINTER :: global

  INTEGER, INTENT(OUT) :: splitflag

  REAL(RFREAL), DIMENSION(3), INTENT(IN) :: xyz1, xyz2, xyz3, xyz4

! ... local variables
  CHARACTER(CHRLEN) :: rcsidentstring

  REAL(RFREAL) :: fvecx1, fvecx2, fvecy1, fvecy2, fvecz1, fvecz2, mag1, mag2, & 
                  prod1, prod2
  REAL(RFREAL) :: xyznodes(3,3)

!******************************************************************************

  rcsidentstring = '$RCSfile: SplitQuadFace.F90,v $ $Revision: 1.3 $'

  CALL registerfunction( global,'SplitQuadFace',&
  'SplitQuadFace.F90' )

! Start -----------------------------------------------------------------------
 
! Split along diagonal 13: subdivision 123-134 
 
  xyznodes(1:3,1) = xyz1(1:3)
  xyznodes(1:3,2) = xyz2(1:3)
  xyznodes(1:3,3) = xyz3(1:3)    
 
  CALL facevectortria( xyznodes,fvecx1,fvecy1,fvecz1 )

  xyznodes(1:3,1) = xyz1(1:3)
  xyznodes(1:3,2) = xyz3(1:3)
  xyznodes(1:3,3) = xyz4(1:3)    

  CALL facevectortria( xyznodes,fvecx2,fvecy2,fvecz2 )

  mag1  = sqrt(fvecx1*fvecx1 + fvecy1*fvecy1 + fvecz1*fvecz1)
  mag2  = sqrt(fvecx2*fvecx2 + fvecy2*fvecy2 + fvecz2*fvecz2)
  prod1 = (fvecx1*fvecx2 + fvecy1*fvecy2 + fvecz1*fvecz2)/(mag1*mag2) 

  IF (prod1 < 0.0_rfreal) THEN 
    CALL errorstop( global,err_face_split,__line__ )
  ENDIF ! prod1

! Split along diagonal 24: subdivision 124-234 
 
  xyznodes(1:3,1) = xyz1(1:3)
  xyznodes(1:3,2) = xyz2(1:3)
  xyznodes(1:3,3) = xyz4(1:3)    
 
  CALL facevectortria( xyznodes,fvecx1,fvecy1,fvecz1 )

  xyznodes(1:3,1) = xyz2(1:3)
  xyznodes(1:3,2) = xyz3(1:3)
  xyznodes(1:3,3) = xyz4(1:3)    

  CALL facevectortria( xyznodes,fvecx2,fvecy2,fvecz2 )

  mag1  = sqrt(fvecx1*fvecx1 + fvecy1*fvecy1 + fvecz1*fvecz1)
  mag2  = sqrt(fvecx2*fvecx2 + fvecy2*fvecy2 + fvecz2*fvecz2)
  prod2 = (fvecx1*fvecx2 + fvecy1*fvecy2 + fvecz1*fvecz2)/(mag1*mag2)

  IF (prod2 < 0.0_rfreal) THEN 
    CALL errorstop( global,err_face_split,__line__ )
  ENDIF ! prod2

! Better subdivision is that which has smaller dotproduct

  IF (prod1 < prod2) THEN 
    splitflag = face_split_13
  ELSE 
    splitflag = face_split_24
  ENDIF ! prod1<prod2

! finalize --------------------------------------------------------------------

  CALL deregisterfunction( global )

END SUBROUTINE splitquadface

!******************************************************************************
!
! RCS Revision history:
!
! $Log: SplitQuadFace.F90,v $
! Revision 1.3  2008/12/06 08:44:10  mtcampbe
! Updated license.
!
! Revision 1.2  2008/11/19 22:17:23  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.1  2004/12/01 16:51:33  haselbac
! Initial revision after changing case
!
! Revision 1.3  2002/10/27 18:48:59  haselbac
! Removed tabs
!
! Revision 1.2  2002/09/05 17:40:20  jblazek
! Variable global moved into regions().
!
! Revision 1.1  2002/03/01 17:07:13  haselbac
! Initial revision
!
!******************************************************************************