Returns true if the two vectors are perpendicular.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=wp), | intent(in), | dimension(:) | :: | v1 | ||
real(kind=wp), | intent(in), | dimension(:) | :: | v2 |
pure function perpendicular(v1, v2) result(is_parallel) implicit none real(wp),dimension(:),intent(in) :: v1 real(wp),dimension(:),intent(in) :: v2 logical :: is_parallel real(wp),parameter :: tol = 10.0_wp * epsilon(1.0_wp) !! tolerance is_parallel = abs(dot_product(unit(v1), unit(v2))) <= tol end function perpendicular