1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
|
// Copyright (c) 2006-2008 Fernando Luis Cacciola Carballal. All rights reserved.
//
// This file is part of CGAL (www.cgal.org)
//
// $URL: https://github.com/CGAL/cgal/blob/v6.1/Straight_skeleton_2/include/CGAL/certified_numeric_predicates.h $
// $Id: include/CGAL/certified_numeric_predicates.h b26b07a1242 $
// SPDX-License-Identifier: GPL-3.0-or-later OR LicenseRef-Commercial
//
// Author(s) : Fernando Cacciola <fernando_cacciola@ciudad.com.ar>
//
#ifndef CGAL_CERTIFIED_NUMERIC_PREDICATES_H
#define CGAL_CERTIFIED_NUMERIC_PREDICATES_H
#include <CGAL/license/Straight_skeleton_2.h>
#include <CGAL/number_utils.h>
#include <CGAL/Interval_arithmetic.h>
#include <CGAL/Uncertain.h>
namespace CGAL {
inline Uncertain<bool> logical_or ( Uncertain<bool> a, Uncertain<bool> b ) { return a | b ; }
inline Uncertain<bool> logical_and( Uncertain<bool> a, Uncertain<bool> b ) { return a & b ; }
inline Uncertain<bool> logical_or ( Uncertain<bool> a, Uncertain<bool> b, Uncertain<bool> c ) { return a | b | c ; }
inline Uncertain<bool> logical_and( Uncertain<bool> a, Uncertain<bool> b, Uncertain<bool> c ) { return a & b & c ; }
template <class NT>
inline Uncertain<bool> certified_is_zero(const NT& x)
{
return CGAL_NTS is_valid(x) ? make_uncertain(CGAL_NTS is_zero(x)) : Uncertain<bool>::indeterminate() ;
}
template <class NT>
inline Uncertain<bool> certified_is_not_zero(const NT& x)
{
return CGAL_NTS is_valid(x) ? make_uncertain(!CGAL_NTS is_zero(x)) : Uncertain<bool>::indeterminate() ;
}
template <class NT>
inline Uncertain<bool> certified_is_one(const NT& x)
{
return CGAL_NTS is_valid(x) ? make_uncertain(CGAL_NTS is_one(x) ) : Uncertain<bool>::indeterminate() ;
}
template <class NT>
inline Uncertain<bool> certified_is_negative(const NT& x)
{
return CGAL_NTS is_valid(x) ? make_uncertain(CGAL_NTS is_negative(x) ) : Uncertain<bool>::indeterminate() ;
}
template <class NT>
inline Uncertain<bool> certified_is_positive(const NT& x)
{
return CGAL_NTS is_valid(x) ? make_uncertain(CGAL_NTS is_positive(x) ) : Uncertain<bool>::indeterminate() ;
}
template <class NT>
inline Uncertain<Sign> certified_sign(const NT& x)
{
return CGAL_NTS is_valid(x) ? make_uncertain(CGAL_NTS sign(x)) : Uncertain<Sign>::indeterminate() ;
}
template <class NT1, class NT2>
inline Uncertain<Comparison_result> certified_compare(const NT1& n1, const NT2& n2)
{
return CGAL_NTS is_valid(n1) && CGAL_NTS is_valid(n2) ? make_uncertain(CGAL_NTS compare(n1,n2))
: Uncertain<Comparison_result>::indeterminate() ;
}
inline Uncertain<bool> certified_is_smaller( Uncertain<Comparison_result> c )
{
return c == SMALLER;
}
inline Uncertain<bool> certified_is_equal( Uncertain<Comparison_result> c )
{
return c == EQUAL;
}
inline Uncertain<bool> certified_is_larger( Uncertain<Comparison_result> c )
{
return c == LARGER;
}
inline Uncertain<bool> certified_is_smaller_or_equal( Uncertain<Comparison_result> c )
{
return logical_or( c == SMALLER , c == EQUAL ) ;
}
inline Uncertain<bool> certified_is_larger_or_equal( Uncertain<Comparison_result> c )
{
return logical_or( c == LARGER , c == EQUAL ) ;
}
template <class NT1, class NT2>
inline Uncertain<bool> certified_is_smaller(const NT1& n1, const NT2& n2)
{
return certified_is_smaller(certified_compare(n1,n2));
}
template <class NT1, class NT2>
inline Uncertain<bool> certified_is_equal(const NT1& n1, const NT2& n2)
{
return certified_is_equal(certified_compare(n1,n2));
}
template <class NT1, class NT2>
inline Uncertain<bool> certified_is_larger(const NT1& n1, const NT2& n2)
{
return certified_is_larger(certified_compare(n1,n2));
}
template <class NT1, class NT2>
inline Uncertain<bool> certified_is_smaller_or_equal(const NT1& n1, const NT2& n2)
{
return certified_is_smaller_or_equal(certified_compare(n1,n2)) ;
}
template <class NT1, class NT2>
inline Uncertain<bool> certified_is_larger_or_equal(const NT1& n1, const NT2& n2)
{
return certified_is_larger_or_equal(certified_compare(n1,n2)) ;
}
template <class NT>
inline Uncertain<Sign> certified_sign_of_determinant2x2( const NT& a00
, const NT& a01
, const NT& a10
, const NT& a11
)
{
return certified_compare(a00*a11, a10*a01) ;
}
} // end namespace CGAL
#endif // CGAL_CERTIFIED_NUMERIC_PREDICATES_H
|
