From 9cd8f0e8e14cc2a2d3127467afa2d636369fb96c Mon Sep 17 00:00:00 2001 From: Laurent Rineau Date: Thu, 15 Jun 2023 11:32:41 +0200 Subject: [PATCH] fix case 12 --- .../CGAL/Triangulation_segment_traverser_3.h | 13 ++++++++++--- 1 file changed, 10 insertions(+), 3 deletions(-) diff --git a/Triangulation_3/include/CGAL/Triangulation_segment_traverser_3.h b/Triangulation_3/include/CGAL/Triangulation_segment_traverser_3.h index 02a29f092a4..0860eaeee56 100644 --- a/Triangulation_3/include/CGAL/Triangulation_segment_traverser_3.h +++ b/Triangulation_3/include/CGAL/Triangulation_segment_traverser_3.h @@ -870,9 +870,13 @@ public: } break; } - case Locate_type::FACET: - _curr_simplex = chprev; //query goes through the cell + case Locate_type::FACET: { + if(facet_has_edge(Facet(chprev, liprev), get_edge())) + _curr_simplex = Facet(chprev, liprev); //edge-facet-outside + else + _curr_simplex = chprev; //query goes through the cell break; + } case Locate_type::CELL: _curr_simplex = ch; //edge-cell-outside break; @@ -915,7 +919,10 @@ public: } case Locate_type::FACET: - _curr_simplex = Cell_handle(_cell_iterator);//query goes through the cell + if(facet_has_edge(Facet(ch, li_exit), get_edge())) + _curr_simplex = Facet(ch, li_exit); //edge-facet-outside + else + _curr_simplex = Cell_handle(_cell_iterator);//query goes through the cell break; default: