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
|
#include "../base/SRC_FIRST.hpp"
#include "robust_orientation.hpp"
#include "../base/macros.hpp"
extern "C"
{
#include "../3party/robust/predicates.c"
}
namespace m2 { namespace robust
{
bool Init()
{
exactinit();
return true;
}
double OrientedS(PointD const & p1, PointD const & p2, PointD const & p)
{
static bool res = Init();
ASSERT_EQUAL ( res, true, () );
UNUSED_VALUE(res);
double a[] = { p1.x, p1.y };
double b[] = { p2.x, p2.y };
double c[] = { p.x, p.y };
return orient2d(a, b, c);
}
bool SegmentsIntersect(PointD const & a, PointD const & b,
PointD const & c, PointD const & d)
{
return
max(a.x, b.x) >= min(c.x, d.x) &&
min(a.x, b.x) <= max(c.x, d.x) &&
max(a.y, b.y) >= min(c.y, d.y) &&
min(a.y, b.y) <= max(c.y, d.y) &&
OrientedS(a, b, c) * OrientedS(a, b, d) <= 0.0 &&
OrientedS(c, d, a) * OrientedS(c, d, b) <= 0.0;
}
}
}
|
