workaround a case when the intersection polyline is not split

the polyline goes from border to interior but the node is
of degree 2 and the polyline is not split

  ***************
  *            *
  *           *
  *          *
 -*---------*-*---
| *           *   |
| *           *   |
| *************   |
|                 |
|                 |
|                 |
|_________________|
This commit is contained in:
Sébastien Loriot 2019-03-01 18:06:54 +01:00
parent a83dac3826
commit 9040b59cb8
1 changed files with 3 additions and 2 deletions

View File

@ -1103,8 +1103,9 @@ public:
// triangle which is tangent at its 3 vertices
// \todo improve this part which is not robust with a kernel
// with inexact constructions.
Bounded_side position = inside_tm2(midpoint(get(vpm1, source(h, tm1)),
get(vpm1, target(h, tm1)) ));
Bounded_side position = inside_tm2(centroid(get(vpm1, source(h, tm1)),
get(vpm1, target(h, tm1)),
get(vpm1, target(next(h, tm1), tm1)) ));
CGAL_assertion( position != ON_BOUNDARY);
if ( position == in_tm2 )
is_patch_inside_tm2.set(patch_id);