CGAL users discussion list

[gp <>]

I've noticed that CGAL::insert on an Arrangement_2< Arr_Bezier_curve_traits_2
> will occasionally violate an assertion. I can't tell exactly what the
erroneous scenarios have in common, but I have a rough idea that they
involve one curve terminating on top of another. I've put together a minimal
example. For me (Windows 8.1, MSVC2015 running in QtCreator, CGAL-4.14),
this fails in debug mode on the assertion of line 234 of
No_intersection_surface_sweep_2_impl.h (in release mode, my system just
crashes):

         CGAL_assertion((m_statusLine.size() == 0));

Can anyone reproduce this?


main.cpp:

#include <CGAL/Cartesian.h>
#include <CGAL/CORE_algebraic_number_traits.h>
#include <CGAL/Arr_Bezier_curve_traits_2.h>
#include <CGAL/Arrangement_2.h>

#include <sstream>

using NtTraits = CGAL::CORE_algebraic_number_traits;
using Rational = NtTraits::Rational;
using Algebraic = NtTraits::Algebraic;
using RationalKernel = CGAL::Cartesian< Rational >;
using AlgebraicKernel = CGAL::Cartesian< Algebraic >;
using ArrTraits = CGAL::Arr_Bezier_curve_traits_2< RationalKernel,
AlgebraicKernel, NtTraits >;
using BezierCurve2 = ArrTraits::Curve_2;
using Arrangement = CGAL::Arrangement_2< ArrTraits >;

int main()
{
    // Change the 10.1 to a 10 and there is no longer an assertion
violation.
    const std::string curvesData =
                "2            \
                2 -10 0 10.1 0 \
                2 10 10 10 0 ";

    std::istringstream stream( curvesData );
    std::vector< BezierCurve2 > curvesToInsert;
    size_t numCurves = 0;
    stream >> numCurves;
    for( size_t i = 0; i < numCurves; i++ ) {
        BezierCurve2 curve;
        stream >> curve;
        curvesToInsert.push_back( curve );
    }

    Arrangement arr;
    // Assertion violation here.
    CGAL::insert( arr, curvesToInsert.begin(), curvesToInsert.end() );
    return 0;
} 



--
Sent from: http://cgal-discuss.949826.n4.nabble.com/