songsenand
/
cgal
mirror of https://github.com/CGAL/cgal
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

fix for an infinite edge conflict 

Here is the problematic input:

$ cat ~/Dropbox/cgal/demo/Segment_Delaunay_graph_Linf_2/September/29sept2seg1pt.cin
s -200 -100 200 100
p 200 -50
s 220 90 -130 -90
This commit is contained in:
Panagiotis Cheilaris 2013-01-23 00:04:01 +02:00
parent 1f2a84d650
commit 17b7ec7ffd
1 changed files with 73 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

73
include/CGAL/Segment_Delaunay_graph_Linf_2/Infinite_edge_interior_conflict_C2.h
View File

@ -24,6 +24,7 @@ public:
typedef typename K::Site_2 Site_2;
typedef typename K::Point_2 Point_2;
typedef typename K::Direction_2 Direction_2;
typedef typename K::RT RT;
typedef typename K::Boolean Boolean;
@ -138,6 +139,78 @@ public:
not (same_points(q, t.source_site()) or
same_points(q, t.target_site()) ) );
// here q is point
if ( not (t.segment().is_horizontal() or
t.segment().is_vertical() ) ) {
Line_2 lt = compute_supporting_line(t);
// Linf-project point q to line lt
Point_2 projq =
compute_linf_projection_nonhom(lt, q.point());
Line_2 lq = compute_line_from_to(projq, q.point());
Point_2 srep;
if (s.is_point()) {
srep = s.point();
} else {
// s is segment
CGAL_assertion(
not (s.segment().is_horizontal() or
s.segment().is_vertical() ) );
Direction_2 d (s.segment());
Line_2 ls = compute_supporting_line(s);
if (CGAL::sign(d.dx()) == CGAL::sign(d.dy())) {
srep = compute_horizontal_projection(
ls, q.point());
} else {
srep = compute_vertical_projection(
ls, q.point());
}
}
Oriented_side oss =
oriented_side_of_line(lq, srep);
Point_2 rrep;
if (r.is_point()) {
rrep = r.point();
} else {
// r is segment
CGAL_assertion(
not (r.segment().is_horizontal() or
r.segment().is_vertical() ) );
Direction_2 d (r.segment());
Line_2 lr = compute_supporting_line(r);
if (CGAL::sign(d.dx()) == CGAL::sign(d.dy())) {
rrep = compute_vertical_projection(
lr, q.point());
} else {
rrep = compute_horizontal_projection(
lr, q.point());
}
}
Oriented_side osr =
oriented_side_of_line(lq, rrep);
if ((oss == ON_NEGATIVE_SIDE) and
(osr == ON_POSITIVE_SIDE) ) {
CGAL_SDG_DEBUG(
std::cout
<< "debug infinite-edge-int-cf tocheck (q,s,r,t,sgn)= "
<< q << ' ' << s << ' ' << r << ' ' << t
<< ' ' << sgn << " returns "
<< false << std::endl;);
return false;
}
} // case where t is neither horizontal nor vertical
//CGAL_assertion(false);
// philaris: tocheck more