強大距離比較謂詞

Question

我需要通過下面的代碼定義一個強大的斷言：

CompareResult compareDistance(Point a, Point b, Point c, Point d) { 
    if (distance(a, b) > distance(c, d)) 
     return Larger; 
    else if (distance(a, b) == distance(c, d)) 
     return Equal; 
    else 
     return Smaller; 
} 

由於浮點運算的限制，我們無法計算究竟distance（即使它的平方），所以如果我們只是直接實現這個代碼，謂詞將不會健壯。我試圖在CGAL庫中找到它，但不能。

有點接近謂語我需要的是compare_distance_to_point(Point p, Point q, Point r)predicate。它返回Smaller如果distance(p, q) < distance(p, r)，Equal如果distance(p, q) == distance(p, r)和Larger否則。首先想到的是通過(c - a)向量轉移c和d，所以我們可以稱之爲compare_distance_to_point(a, b, d + (c - a))，但這將再次違反的魯棒性。那麼，有沒有人有適應它的想法？

Answer 1

如果你用準確的謂詞如 Exact_predicates_inexact_constructions_kernel， 可以使用仿函數Compare_distance_3這是概念CompareDistance_3模型的內核。

Kernel::Compare_distance_3 cmp;

return cmp(a,b,c,d);