Hi, I attach my dataset ( note that it is in Point3D format, but it shouldn't

matter as the z is 0 always and you can easily convert to Point2D).

And this is the gist of my code:

typedef CGAL::Triangulation_vertex_base_2<Kernel> VbT;

typedef CGAL::Triangulation_face_base_with_info_2<FaceInfo2,Kernel> Fbb;

typedef CGAL::Constrained_triangulation_face_base_2<Kernel,Fbb> FbT;

typedef CGAL::Triangulation_data_structure_2<VbT,FbT> TDS;

typedef CGAL::Exact_predicates_tag Itag;

typedef CGAL::Constrained_Delaunay_triangulation_2<Kernel, TDS, Itag> CDT;

typedef CGAL::Alpha_shape_vertex_base_2<Kernel> Vb;

typedef CGAL::Alpha_shape_face_base_2<Kernel> Fb;

typedef CGAL::Triangulation_data_structure_2<Vb,Fb> Tds;

typedef CGAL::Delaunay_triangulation_2<Kernel,Tds> Alpha_Triangulation_2;

typedef CGAL::Alpha_shape_2<Alpha_Triangulation_2> Alpha_shape_2;

Alpha_shape_2 *A =new Alpha_shape_2(first,last, 100000,Mode.Regular);

auto maxLevel = A->number_of_alphas()-1;

auto alphaLevel = A->get_nth_alpha(maxLevel);

int i=0;

for(fit=A->finite_faces_begin(); fit!=A->finite_faces_end(); ++fit)

{

A->set_alpha(alphaLevel);

auto classType2 = A->classify(fit); //i==504 is

EXTERIOR

i++;

}

From what I can tell, there is an element at i==504, which is EXTERIOR. (

there are 545 elements).

Do let me know whether you can reproduce the problem with the dataset and

the above pseudo code. If not I will prepare a fully working sample.

I am using CGAL 4.6.2. Not sure whether it matters or not

coordinate.csv

<

http://cgal-discuss.949826.n4.nabble.com/file/t212754/coordinate.csv>

--

Sent from:

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

You are currently subscribed to cgal-discuss.

To unsubscribe or access the archives, go to

https://sympa.inria.fr/sympa/info/cgal-discuss