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

endpoint of segment code completion 

case of predicate pqrt of type PPSP:
if new point t is on the same side of line pq
as the other endpoint of r, then NEGATIVE, else POSITIVE.

This fixes inputs of the following form (zppsp1.cin):
s -100 50 0 0
p 100 0
p 0 100

With the fix, point p 0 0 has as neighbors both
p 100 0 and p 0 100. Before, p 0 100 was not a neighbor.

Signed-off-by: Panagiotis Cheilaris <philaris@cs.ntua.gr>
This commit is contained in:
Panagiotis Cheilaris 2013-08-09 08:56:05 +03:00
parent e3a376ebab
commit f81ca24f1b
1 changed files with 24 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

24
Segment_Delaunay_graph_Linf_2/include/CGAL/Segment_Delaunay_graph_Linf_2/Voronoi_vertex_sqrt_field_new_C2.h
View File

@ -2122,6 +2122,30 @@ private:
CGAL_assertion(v_type == PPS);
// if new point t is on the same side of line pq
// as the other endpoint of r, then CONFLICT
Line_2 l = compute_line_from_to(p.point(), q.point());
Oriented_side ost = oriented_side_of_line(l, t.point());
CGAL_assertion(ost != ON_ORIENTED_BOUNDARY);
Point_2 other_of_r;
if (is_p_endp_of_r) {
other_of_r = (same_points(p, r.source_site()))?
(r.target_site().point()) : (r.source_site().point());
} else { // is_q_endp_of_r
other_of_r = (same_points(q, r.source_site()))?
(r.target_site().point()) : (r.source_site().point());
}
Oriented_side osother = oriented_side_of_line(l, other_of_r);
CGAL_assertion(osother != ON_ORIENTED_BOUNDARY);
if (osother == ost) {
return NEGATIVE;
} else {
return POSITIVE;
}
} // case r is non-hv and has endpoint p or q
CGAL_assertion(num_same_quadrant_as_t == 0);