From 32fb8ce9dd769156bad2c400cf4a5e5aade085b8 Mon Sep 17 00:00:00 2001 From: Simon Giraudot Date: Thu, 18 Aug 2016 11:29:02 +0200 Subject: [PATCH] Bugfix: loose bbox length was badly computed --- .../CGAL/Optimal_transportation_reconstruction_2.h | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/Optimal_transportation_reconstruction_2/include/CGAL/Optimal_transportation_reconstruction_2.h b/Optimal_transportation_reconstruction_2/include/CGAL/Optimal_transportation_reconstruction_2.h index a8386a6a6cb..aeb498f4c2a 100644 --- a/Optimal_transportation_reconstruction_2/include/CGAL/Optimal_transportation_reconstruction_2.h +++ b/Optimal_transportation_reconstruction_2/include/CGAL/Optimal_transportation_reconstruction_2.h @@ -441,17 +441,17 @@ public: if (m_verbose > 0) std::cerr << "insert loose bbox..."; - double xl = (bbox.xmax()-bbox.xmin()/2); - double yl = (bbox.ymax()-bbox.ymin()/2); + double dl = (std::max)((bbox.xmax()-bbox.xmin()) / 2., + (bbox.ymax()-bbox.ymin()) / 2.); timer.start(); int nb = static_cast(m_dt.number_of_vertices()); typename Traits::Construct_point_2 point_2 = m_traits.construct_point_2_object(); - insert_point(point_2(bbox.xmin()-xl, bbox.ymin()-yl), true, nb++); - insert_point(point_2(bbox.xmin()-xl, bbox.ymax()+yl), true, nb++); - insert_point(point_2(bbox.xmax()+xl, bbox.ymax()+yl), true, nb++); - insert_point(point_2(bbox.xmax()+xl, bbox.ymin()-yl), true, nb++); + insert_point(point_2(bbox.xmin()-dl, bbox.ymin()-dl), true, nb++); + insert_point(point_2(bbox.xmin()-dl, bbox.ymax()+dl), true, nb++); + insert_point(point_2(bbox.xmax()+dl, bbox.ymax()+dl), true, nb++); + insert_point(point_2(bbox.xmax()+dl, bbox.ymin()-dl), true, nb++); if (m_verbose > 0) std::cerr << "done (" << nb << " vertices, "