Skip to Content.
Sympa Menu

cgal-discuss - Re: [cgal-discuss] incident_facets function of a Delaunay triangulation

Subject: CGAL users discussion list

List archive

Re: [cgal-discuss] incident_facets function of a Delaunay triangulation


Chronological Thread 
  • From: Manuel Caroli <>
  • To:
  • Subject: Re: [cgal-discuss] incident_facets function of a Delaunay triangulation
  • Date: Tue, 24 Mar 2009 19:34:29 +0100
Hi,

Andreas was asking for a "minimalistic example", i.e. an example that does compile and produces the error (see also http://www.cgal.org/bug_report.html).
Your example does not compile because
1) I don't have "stdafx..h"
2) I don't seem to have System::String (my fault?)

Anyway when uncommenting these two things it compiles fine and runs through without any error. So I'm sorry I cannot reproduce your problem.

Maybe you could try with CGAL 3.4 and if the problem still persists file a bug report as described in the link above.

best

Manuel


Query Cgal wrote:
Sorry that I missed them out before.

the following is the whole piece of the code. Thanks

#include "stdafx..h"

#include <fstream>
#include <vector>
#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
#include <CGAL/Delaunay_triangulation_3.h>
#include <CGAL/Triangulation_vertex_base_3.h>
#include <CGAL/Cartesian.h>
#include <iostream>

using namespace std;


typedef CGAL::Cartesian<double> Rep_class;
typedef CGAL::Vector_3<Rep_class> Vector;
typedef CGAL::Point_3<Rep_class> Point_3;

typedef CGAL::Object CGAL_object;



template < class GT, class Vb = CGAL::Triangulation_vertex_base_3<GT> >
class Vertex_base_for_ssr
: public Vb
{
public:
typedef typename Vb::Vertex_handle Vertex_handle;
typedef typename Vb::Cell_handle Cell_handle;
typedef typename Vb::Point Point;

template < class TDS2 >
struct Rebind_TDS {
typedef typename Vb::template Rebind_TDS<TDS2>::Other Vb2;
typedef Vertex_base_for_ssr<GT, Vb2> Other;
};


Vector surface_normal;
vector<Vertex_base_for_ssr<GT, Vb> *> neighbor_vertex_handle;
double fifth_nearest_dist;


int surface_label;
bool is_safe;



Vertex_base_for_ssr() : Vb()
{
surface_label = -1;
is_safe = false;
}

Vertex_base_for_ssr(const Point& p)
: Vb(p)
{
coordinate[0] = p.x();
coordinate[1] = p.y();
coordinate[2] = p.z();
surface_label = -1;
is_safe = false;
}


Vertex_base_for_ssr(const Point& p, Cell_handle c)
: Vb(p, c)
{
surface_label = -1;
is_safe = false;
}

void set_surface_normal(double x_i, double y_i, double z_i)
{
surface_normal = Vector(x_i, y_i, z_i);
}

Vector get_surface_normal()
{
return surface_normal;
}

Point get_point()
{
return Vb::point();
}

double coordinate[3];

void set_coordiante()
{
coordinate[0] = point().x();
coordinate[1] = point().y();
coordinate[2] = point().z();
}
};





typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
typedef CGAL::Triangulation_data_structure_3<Vertex_base_for_ssr<K> > Tds;
typedef CGAL::Delaunay_triangulation_3<K, Tds> Delaunay_triangulation;

typedef Delaunay_triangulation::Vertex_handle Vertex_handle;
typedef Delaunay_triangulation::Vertex Vertex;
typedef Delaunay_triangulation::Point Point;
typedef Delaunay_triangulation::Finite_vertices_iterator Finite_vertices_iterator;



typedef Delaunay_triangulation::Facet;


int main(array<System::String ^> ^args)
{
Delaunay_triangulation Dt;

Dt.insert(Point(0,0,0));
Dt.insert(Point(10,0,0));
Dt.insert(Point(5,10,0));
Dt.insert(Point(8,8,0));

vector<Vertex_handle> a;

Finite_vertices_iterator start = Dt.finite_vertices_begin();
Finite_vertices_iterator end = Dt.finite_vertices_end();
Finite_vertices_iterator curr;



list<Delaunay_triangulation::Facet> f;

Dt.incident_facets(start, back_inserter(f));

}

--- *2009年3月24日 星期二，Andreas Fabri /<>/* 寫道﹕


寄件人: Andreas Fabri
<>
主題: Re: [cgal-discuss] incident_facets function of a Delaunay
triangulation
收件人:

日期: 2009 3 24 星期二 下午 6:08


There are no #include statements




Query Cgal wrote:
> Hi. Thanks for your reply.
>
> I got the exception when I try the following piece of code.
Thanks for your help
>
>
>
> typedef CGAL::Cartesian<double> Rep_class;
> typedef CGAL::Vector_3<Rep_class> Vector;
> typedef CGAL::Point_3<Rep_class> Point_3;
>
> template < class GT, class Vb =
CGAL::Triangulation_vertex_base_3<GT> >
> class Vertex_base_for_ssr
> : public Vb
> {
> public:
> typedef typename Vb::Vertex_handle Vertex_handle;
> typedef typename Vb::Cell_handle Cell_handle;
> typedef typename Vb::Point Point;
>
> template < class TDS2 >
> struct Rebind_TDS {
> typedef typename Vb::template Rebind_TDS<TDS2>::Other Vb2;
> typedef Vertex_base_for_ssr<GT, Vb2> Other;
> };
>
>
> Vector surface_normal;
> vector<Vertex_base_for_ssr<GT, Vb> *> neighbor_vertex_handle;
> double fifth_nearest_dist;
>
>
> int surface_label;
> bool is_safe;
>
>
>
> Vertex_base_for_ssr() : Vb()
> {
> surface_label = -1;
> is_safe = false;
> }
>
> Vertex_base_for_ssr(const Point& p)
> : Vb(p)
> {
> coordinate[0] = p.x();
> coordinate[1] = p.y();
> coordinate[2] = p.z();
> surface_label = -1;
> is_safe = false;
> }
>
>
> Vertex_base_for_ssr(const Point& p, Cell_handle c)
> : Vb(p, c)
> {
> surface_label = -1;
> is_safe = false;
> }
>
> void set_surface_normal(double x_i, double y_i, double z_i)
> {
> surface_normal = Vector(x_i, y_i, z_i);
> }
>
> Vector get_surface_normal()
> {
> return surface_normal;
> }
>
> Point get_point()
> {
> return Vb::point();
> }
>
> double coordinate[3];
>
> void set_coordiante()
> {
> coordinate[0] = point().x();
> coordinate[1] = point().y();
> coordinate[2] = point().z();
> }
> };
>
> typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
> typedef
CGAL::Triangulation_data_structure_3<Vertex_base_for_ssr<K> > Tds;
> typedef CGAL::Delaunay_triangulation_3<K, Tds>
Delaunay_triangulation;
>
> int main(array<System::String ^> ^args)
> {
> Delaunay_triangulation Dt;
>
> Dt.insert(Point(0,0,0));
> Dt.insert(Point(10,0,0));
> Dt.insert(Point(5,10,0));
> Dt.insert(Point(8,8,0));
>
> vector<Vertex_handle> a;
>
> Finite_vertices_iterator start = Dt.finite_vertices_begin();
> Finite_vertices_iterator end = Dt.finite_vertices_end();
> Finite_vertices_iterator curr;
>
>
>
> list<Delaunay_triangulation::Facet> f;
>
> Dt.incident_facets(start, back_inserter(f));
> }
>
> --- *2009年3月24日 星期二，Andreas Fabri

/<

</mc/compose?>>/*
寫道﹕
>
>
> 寄件人: Andreas Fabri
<

</mc/compose?>>
> 主題: Re: [cgal-discuss] incident_facets function of a Delaunay
> triangulation
> 收件人:


</mc/compose?>
> 日期: 2009 3 24 星期二 下午 5:47
>
>
> Hi Query CGAL
>
> Can you please write a minimalistic example program that produces
> the error,
> so that we can try to reproduce it
>
>
> andreas
>
>
>
> Query Cgal wrote:
> > Excuse me, may I ask something about a strange problem I have
> asked beofre?
> >
> > The problem is that when I use the incident_facets()
function of
> a 3D Delaunay Triangulation, I get
> > System.NullReferenceException
> >
> >
> > and I found out that the exception happens when it is
executing
> line 662 of triangulation_data_structure_3.h
> >
> > I am using Microsoft Visual C++ 2005 and CGAL-3.3.1
> >
> > May I ask what may be the possible reason for the problem
then?
> >
> > Thank you very much
> >
> >
> > --- *2009年3月23日 星期一，Query Cgal

/<

</mc/compose?>
>
</mc/compose?

</mc/compose?>>>/*
寫道﹕
> >
> >
> > 寄件人: Query Cgal
<

</mc/compose?>
>
</mc/compose?

</mc/compose?>>>
> > 主題: Re: [cgal-discuss] incident_facets function of a
Delaunay
> > triangulation
> > 收件人:
..fr

</mc/compose?..fr>
>
</mc/compose?..fr

</mc/compose?..fr>>
> > 日期: 2009 3 23 星期一 下午 12:35
> >
> > Thank you for your reply.
> >
> > I am using C++.
> > The complier that I use is Visual Studio C++ 2005.
> > And the CGAL version that I use is 3.3.1
> >
> > Thanks
> >
> > --- *2009年3月23日 星期一，Andreas Fabri
> >
/<

</mc/compose?>
>
</mc/compose?

</mc/compose?>>>/*
寫道﹕
> >
> >
> > 寄件人: Andreas Fabri

<

</mc/compose?>
>
</mc/compose?

</mc/compose?>>>
> > 主題: Re: [cgal-discuss] incident_facets function of a
> Delaunay
> > triangulation
> > 收件人:
..fr

</mc/compose?..fr>
>
</mc/compose?...fr

</mc/compose?..fr>>
> > 日期: 2009 3 23 星期一 上午 9:44
> >
> > Hello,
> >
> > What programming language do you use? I ask because
> >
> > int main(array<System::String ^> ^args)
> >
> > and
> >
> > System.NullReferenceException
> >
> > don't look like C++
> >
> > andreas
> >
> >
> > naresh wrote:
> > > Hi
> > > Can you change kernel for Vertex_base_for_ssr
class as
> > Exact_predicates_inexact_constructions_kernel ?
> > > I donot anything wrong in code it should work
try to match
> > similar kernel .
> > > > ----- Original Message -----
> > > *From:* Query Cgal

<mailto:cgalquery@yahoo..com.hk

</mc/compose?to=cgalquery@yahoo..com.hk>
>
</mc/compose?to=cgalquery@yahoo..com.hk

</mc/compose?to=cgalquery@yahoo..com.hk>>>
> > > *To:*


</mc/compose?>
>
</mc/compose?

</mc/compose?>>
> > >
<mailto:

</mc/compose?>
>
</mc/compose?

</mc/compose?>>>
> > > *Sent:* Monday, March 23, 2009 10:13 AM
> > > *Subject:* Re: [cgal-discuss] incident_facets
> function of
> > a Delaunay
> > > triangulation
> > >
> > > Hi, Sorry that I have sent an empty reply
before
> because
> > I have hot
> > > the wrong button.
> > >
> > > The line that I got the exception is the line
> > >
> > > Dt.incident_facets(start, back_inserter(f));
> > >
> > > and after I have chnaged the vector to a
list, the
> > exception still
> > > occur..
> > >
> > > Thanks
> > >
> > > --- *2009年3月23日 星期一，Query Cgal
> >
/<

</mc/compose?>
>
</mc/compose?..hk

</mc/compose?..hk>>
> > >
<mailto:

</mc/compose?to=cgalquery@yahoo..com.hk>
>
</mc/compose?

</mc/compose?>>>>/*
寫道﹕
> > >
> > >
> > > 寄件人: Query Cgal

<cgalquery@yahoo..com.hk

</mc/compose?to=cgalquery@yahoo..com.hk>
>
</mc/compose?to=cgalquery@yahoo..com.hk

</mc/compose?to=cgalquery@yahoo..com.hk>>
> > >
<mailto:

</mc/compose?>
>
</mc/compose?

</mc/compose?>>>>
> > > 主題: Re: [cgal-discuss] incident_facets
> function of
> > a Delaunay
> > > triangulation
> > > 收件人:


</mc/compose?>
>
</mc/compose?

</mc/compose?>>
> > > 日期: 2009 3 23 星期一 上午 4:41
> > >
> > >
> > >
> > > --- *2009年3月23日 星期一，naresh
> >
/<alexdowson@hotmail...com

</mc/compose?to=alexdowson@hotmail...com>
>
</mc/compose?to=alexdowson@hotmail...com

</mc/compose?to=alexdowson@hotmail...com>>>/*
> > > 寫道﹕
> > >
> > >
> > > 寄件人: naresh

<

</mc/compose?>
>
</mc/compose?

</mc/compose?>>>
> > > 主題: Re: [cgal-discuss]
incident_facets
> function
> > of a
> > > Delaunay triangulation
> > > 收件人:



</mc/compose?>
>
</mc/compose?

</mc/compose?>>
> > > 日期: 2009 3 23 星期一 上午 4:33
> > >
> > > 
> > > Hi
> > > in Which line you are
getting
> > System.NullReferenceException
> > > of code ?
> > > Hmmm.... You can try this
> instead of this
> > > vector<Delaunay_triangulation::Facet> f; with
> > > list<Delaunay_triangulation::Facet> f;
> > > try list instead of
vector.
> > > > ----- Original Message -----
> > > *From:* Query Cgal
> > > *To:*



</mc/compose?>
>
</mc/compose?

</mc/compose?>>
> > > *Sent:* Monday, March 23, 2009
9:47 AM
> > > *Subject:* Re: [cgal-discuss]
> incident_facets
> > function
> > > of a Delaunay triangulation
> > >
> > > Hi.
> > >
> > > Thank you for your reply.
> > >
> > > Finally, I have tried with the
> following code
> > to use the
> > > incident_facets function.
> > >
> > > int main(array<System::String
^> ^args)
> > > {
> > > Delaunay_triangulation Dt;
> > >
> > > Dt.insert(Point(0,0,0));
> > > Dt.insert(Point(10,0,0));
> > > Dt.insert(Point(5,10,0));
> > > Dt.insert(Point(8,8,0));
> > > Dt...insert(Point(-3, -20, 0));
> > > Dt.insert(Point(-5,10,0));
> > >
> > >
> > > Finite_vertices_iterator
start =
> > > Dt.finite_vertices_begin();
> > >
> > > vector<Delaunay_triangulation::Facet> f;
> > >
> > > Dt.incident_facets(start,
> back_inserter(f));
> > >
> > > }
> > >
> > > However, when it has
> > "System.NullReferenceException'
> > > May I ask is it due to the
nature of the
> > points that I
> > > have inserted into the
triangulation?
> > >
> > > I am using CGAL-3.3.1 and
Visual C++ 2005
> > > And I am using the
triangulation with the
> > following
> > > declaration.
> > > Thank you very much for your help.
> > >
> > >
> > >
> > > typedef CGAL::Cartesian<double>
Rep_class;
> > > typedef
CGAL::Vector_3<Rep_class> Vector;
> > > typedef
CGAL::Point_3<Rep_class> Point_3;
> > >
> > > typedef CGAL::Object CGAL_object;
> > >
> > >
> > >
> > > template < class GT, class Vb =
> > > CGAL::Triangulation_vertex_base_3<GT> >
> > > class Vertex_base_for_ssr
> > > : public Vb
> > > {
> > > public:
> > > typedef typename
Vb::Vertex_handle Vertex_handle;
> > > typedef typename
Vb::Cell_handle Cell_handle;
> > > typedef typename Vb::Point Point;
> > >
> > > template < class TDS2 >
> > > struct Rebind_TDS {
> > > typedef typename Vb::template
> > > Rebind_TDS<TDS2>::Other Vb2;
> > > typedef Vertex_base_for_ssr<GT,
> > > Vb2> Other;
> > > };
> > >
> > >
> > > Vector surface_normal;
> > > vector<Vertex_base_for_ssr<GT, Vb> *>
> > > neighbor_vertex_handle;
> > > double fifth_nearest_dist;
> > >
> > >
> > > int surface_label;
> > > bool is_safe;
> > >
> > >
> > >
> > > Vertex_base_for_ssr() : Vb()
> > > {
> > > surface_label = -1;
> > > is_safe = false;
> > > }
> > >
> > > Vertex_base_for_ssr(const
Point& p)
> > > : Vb(p)
> > > {
> > > coordinate[0] = p.x();
> > > coordinate[1] = p.y();
> > > coordinate[2] = p.z();
> > > surface_label = -1;
> > > is_safe = false;
> > > }
> > >
> > >
> > > Vertex_base_for_ssr(const
Point& p,
> > Cell_handle c)
> > > : Vb(p, c)
> > > {
> > > surface_label = -1;
> > > is_safe = false;
> > > }
> > >
> > > void
set_surface_normal(double x_i,
> double
> > y_i, double
> > > z_i)
> > > {
> > > surface_normal = Vector(x_i,
> y_i, z_i);
> > > }
> > >
> > > Vector get_surface_normal()
> > > {
> > > return surface_normal;
> > > }
> > >
> > > Point get_point()
> > > {
> > > return Vb::point();
> > > }
> > >
> > > double coordinate[3];
> > >
> > > void set_coordiante()
> > > {
> > > coordinate[0] = point()..x();
> > > coordinate[1] = point().y();
> > > coordinate[2] = point().z();
> > > }
> > > };
> > >
> > > typedef
> > > CGAL::Exact_predicates_inexact_constructions_kernel K;
> > > typedef
> > > CGAL::Triangulation_data_structure_3<Vertex_base_for_ssr<K>
> > > > Tds;
> > > typedef
> CGAL::Delaunay_triangulation_3<K, Tds>
> > > Delaunay_triangulation;
> > >
> > >
> > >
> > >
> > >
> > > --- 2009年3月21日 星期六，naresh
> > >
<

</mc/compose?>
>
</mc/compose?

</mc/compose?>>>
寫道﹕
> > >
> > >
> > > 寄件人: naresh
>
<

</mc/compose?>

</mc/compose?

</mc/compose?>>>
> > > 主題: Re: [cgal-discuss]
> incident_facets
> > function of
> > > a Delaunay triangulation
> > > 收件人:
>


</mc/compose?>
>
</mc/compose?

</mc/compose?>>
> > > 日期: 2009 3 21 星期六 下午
1:47
> > >
> > > 
> > > Hi
> > > >>In
fact, I am
> > dealing with a 3D triangulation
> > > (sorry that I forgot to mention
that
> it is 3D
> > in the
> > > >>previous post.)
> > > >>and therefore, I cannot use
> Facet since
> > Facet is
> > > not represented in 3D
triangulation.
> > > >>So, what should I do to use
> > incident_facets() of a
> > > 3D triangulation.
> > > I donot
know which
> > triangulation are you using. for
> > > simple triangulation here there
you can
> > define Facet
> > > easily .. here is code..
> > > try to modifiy with ur code
> > > typedef
> > > CGAL::Exact_predicates_inexact_constructions_kernel K;
> > >
> > > typedef
CGAL::Triangulation_3<K>
> > Triangulation;
> > >
> > > typedef
Triangulation::Facet Facet;
> > >
> > > If there is still confusion
post some
> > code . But i
> > > will not able reply today
because saturday
> > now and i am
> > > off. :) I can reply on monday..
> > >
> > > > Hope this helps
> > >
> > > Alex
> > >
> > > > ----- Original Message -----
> > > From: Query Cgal
> > > To:
>
..fr

</mc/compose?..fr>
>
</mc/compose?..fr

</mc/compose?..fr>>
> > > Sent: Saturday, March
21, 2009
> 6:46 PM
> > > Subject: Re: [cgal-discuss]
> > incident_facets
> > > function of a Delaunay
triangulation
> > >
> > > THank you for your reply.
> > >
> > > In fact, I am dealing
with a 3D
> > triangulation
> > > (sorry that I forgot to mention
that
> it is 3D
> > in the
> > > previous post.)
> > > and therefore, I cannot
use Facet
> > since Facet is
> > > not represented in 3D
triangulation.
> > > So, what should I do to use
> > incident_facets() of
> > > a 3D triangulation.
> > >
> > > Thank you
> > >
> > >
> > >
> > > --- 2009年3月21日 星期
六，naresh
> > >
<

</mc/compose?>
>
</mc/compose?

</mc/compose?>>>
寫道﹕
> > >
> > >
> > > 寄件人: naresh
> >
<

</mc/compose?>
>
</mc/compose?

</mc/compose?>>>
> > > 主題: Re:
[cgal-discuss]
> > incident_facets
> > > function of a Delaunay
triangulation
> > > 收件人:
> >


</mc/compose?>
>
</mc/compose?

</mc/compose?>>
> > > 日期: 2009 3 21 星
期六 下
> 午 1:07
> > >
> > > 
> > > Hi
> > > I
> > think you should try like this
> > > typedef Trianluation_2::Facet Facet; //
> > > something here that you had
used in ur
> > programme may
> > > different i had showed only example
> > > vector<Facet> f;
> > > Dt.incident_facets(start,
> > back_inserter(f));
> > > > -----
Original Message -----
> > > From: Query Cgal
> > > To:
> >


</mc/compose?>
>
</mc/compose?

</mc/compose?>>
> > > Sent: Saturday,
March 21,
> > 2009 5:45 PM
> > > Subject:
[cgal-discuss]
> > incident_facets
> > > function of a Delaunay
triangulation
> > >
> > > Hi, I have tried to
> use the
> > function below
> > >
> > > template <class
> OutputIterator>
> > > OutputIterator t.incident_facets (
> > > Vertex_handle v, OutputIterator
facets)
> > >
> > > Copies all Facets incident to v to
> > > the output iterator facets.
Returns the
> > resulting output
> > > iterator.
> > > Precondition: t...dimension() >1 , v
> > > Vertex_handle(), t.is_vertex(v)..
> > > When I tried to used it in the following
> > > way but it didn't work,
> > >
> > > vector<Finite_facets_iterator> f;
> > > Dt.incident_facets(start,
> > back_inserter(f));
> > >
> > > May I ask in what
> other way
> > should I use
> > > the Output iterator then ?
> > >
> > > Thank you
> > >
> > >
> > >
> > >
> > >
> > >
> > >
> > > Yahoo!香港提供
網上安全
> 攻略，
> > 教你如何防範
> > > 黑客!了解更多
> > >
> > >
> > > Yahoo!香港提供網上安全
攻略，教
> 你如何
> > 防範黑客!了
> > > 解更多
> > >
> > >
> > > ------------------------------------------------------------------------
> > > Yahoo!香港提供網上安全攻略，教
你如何防
> 範黑客!
> > *了解更多*
> > > <http://hk.promo.yahoo.com/security/
> > <http://hk.promo..yahoo.com/security/>>
> > >
> > >
> > > ------------------------------------------------------------------------
> > > Yahoo!香港提供網上安全攻略，教你如何防
範黑客!*
> 了解更多*
> > > <http://hk.promo.yahoo.com/security/
> > <http://hk..promo.yahoo.com/security/>>
> > >
> > >
> > > ------------------------------------------------------------------------
> > > Yahoo!香港提供網上安全攻略，教你如何防範黑
客!*了解
> 更多*
> > > <http://hk.promo.yahoo.com/security/>
> >
> > -- You are currently subscribed to cgal-discuss.
> > To unsubscribe or access the archives, go to
> > https://lists-sop.inria.fr/wws/info/cgal-discuss
> >
> >
> > ------------------------------------------------------------------------
> > Yahoo!香港提供網上安全攻略，教你如何防範黑客!*了解更多*
> > <http://hk.promo.yahoo.com/security/>
> >
> >
> >
> ------------------------------------------------------------------------
> > Yahoo!香港提供網上安全攻略，教你如何防範黑客!*了解更多*
> <http://hk..promo.yahoo.com/security/>
>
> -- You are currently subscribed to cgal-discuss.
> To unsubscribe or access the archives, go to
> https://lists-sop.inria.fr/wws/info/cgal-discuss
>
>
>
------------------------------------------------------------------------
> Yahoo!香港提供網上安全攻略，教你如何防範黑客!*了解更多*
<http://hk.promo.yahoo.com/security/>

-- You are currently subscribed to cgal-discuss.
To unsubscribe or access the archives, go to
https://lists-sop.inria.fr/wws/info/cgal-discuss


------------------------------------------------------------------------
Yahoo!香港提供網上安全攻略，教你如何防範黑客!*了解更多* <http://hk.promo.yahoo.com/security/>


Archive powered by MHonArc 2.6.16.

Top of Page