unitVectors_p Subroutine

  1. korc_collisions.f90
  2. korc_collisions
  3. unitVectors_p

private subroutine unitVectors_p(b_unit_X, b_unit_Y, b_unit_Z, b1_X, b1_Y, b1_Z, b2_X, b2_Y, b2_Z, b3_X, b3_Y, b3_Z)

Arguments

Type IntentOptional AttributesName
real(kind=rp), intent(in), DIMENSION(8):: b_unit_X
real(kind=rp), intent(in), DIMENSION(8):: b_unit_Y
real(kind=rp), intent(in), DIMENSION(8):: b_unit_Z
real(kind=rp), intent(out), DIMENSION(8):: b1_X
real(kind=rp), intent(out), DIMENSION(8):: b1_Y
real(kind=rp), intent(out), DIMENSION(8):: b1_Z
real(kind=rp), intent(out), DIMENSION(8):: b2_X
real(kind=rp), intent(out), DIMENSION(8):: b2_Y
real(kind=rp), intent(out), DIMENSION(8):: b2_Z
real(kind=rp), intent(out), DIMENSION(8):: b3_X
real(kind=rp), intent(out), DIMENSION(8):: b3_Y
real(kind=rp), intent(out), DIMENSION(8):: b3_Z

Contents

Source Code

unitVectors_p

Source Code

  subroutine unitVectors_p(b_unit_X,b_unit_Y,b_unit_Z,b1_X,b1_Y,b1_Z, &
            b2_X,b2_Y,b2_Z,b3_X,b3_Y,b3_Z)
    REAL(rp), DIMENSION(8), INTENT(IN) 	:: b_unit_X,b_unit_Y,b_unit_Z
    REAL(rp), DIMENSION(8), INTENT(OUT) :: b1_X,b1_Y,b1_Z
    REAL(rp), DIMENSION(8), INTENT(OUT) :: b2_X,b2_Y,b2_Z
    REAL(rp), DIMENSION(8), INTENT(OUT) :: b3_X,b3_Y,b3_Z
    REAL(rp), DIMENSION(8) :: b2mag,b3mag
    integer(ip) :: cc

    !$OMP SIMD 
!    !$OMP& aligned(b1_X,b1_Y,b1_Z,b_unit_X,b_unit_Y,b_unit_Z, &
!    !$OMP& b2_X,b2_Y,b2_Z,b2mag,b3_X,b3_Y,b3_Z,b3mag)
    do cc=1_idef,8_idef
       b1_X(cc) = b_unit_X(cc)
       b1_Y(cc) = b_unit_Y(cc)
       b1_Z(cc) = b_unit_Z(cc)

       b2_X(cc) = b1_Y(cc)
       b2_Y(cc) = -b1_X(cc)
       b2_Z(cc) = 0._rp

       b2mag(cc)=sqrt(b2_X(cc)*b2_X(cc)+b2_Y(cc)*b2_Y(cc)+b2_Z(cc)*b2_Z(cc))

       b2_X(cc) = b2_X(cc)/b2mag(cc)
       b2_Y(cc) = b2_Y(cc)/b2mag(cc)
       b2_Z(cc) = b2_Z(cc)/b2mag(cc)

       b3_X(cc)=b1_Y(cc)*b2_Z(cc)-b1_Z(cc)*b2_Y(cc)
       b3_Y(cc)=b1_Z(cc)*b2_X(cc)-b1_X(cc)*b2_Z(cc)
       b3_Z(cc)=b1_X(cc)*b2_Y(cc)-b1_Y(cc)*b2_X(cc)

       b3mag(cc)=sqrt(b3_X(cc)*b3_X(cc)+b3_Y(cc)*b3_Y(cc)+b3_Z(cc)*b3_Z(cc))

       b3_X(cc) = b3_X(cc)/b3mag(cc)
       b3_Y(cc) = b3_Y(cc)/b3mag(cc)
       b3_Z(cc) = b3_Z(cc)/b3mag(cc)
    end do
    !$OMP END SIMD
    
  end subroutine unitVectors_p