Infinite loop due to zero-area triangle but positive oriented according to ExactPredicates.jl

[DanielVandH]

i  = 265 
rng = StableRNG(i)
points = [(0.0, 0.0), (1.0, 0.0), (1.0, 1.0), (0.0, 1.0), (0.5, 0.5), (0.2, 0.8), (0.1, 0.785)]
tri = triangulate(points; rng)

hits an infinite loop, since

julia> p
(0.5, 0.5)

julia> q
(0.2, 0.8)

julia> r
(0.0, 1.0)

julia> orient(p,q,r)
1

julia> DelaunayTriangulation.triangle_area(p,q,r)
0.0