Skip to Content.
Sympa Menu

cgal-discuss - Re: [cgal-discuss] delaunay 2d triangulation/meshing with fixed points

Subject: CGAL users discussion list

List archive

Re: [cgal-discuss] delaunay 2d triangulation/meshing with fixed points


Chronological Thread 
  • From: rytis <>
  • To:
  • Subject: Re: [cgal-discuss] delaunay 2d triangulation/meshing with fixed points
  • Date: Fri, 26 Feb 2016 05:21:01 +0100
  • Authentication-results: mail2-smtp-roc.national.inria.fr; spf=None ; spf=Pass ; spf=None
  • Ironport-phdr: 9a23:A/fDIB2Xpt8OVR6LsmDT+DRfVm0co7zxezQtwd8ZsegeKfad9pjvdHbS+e9qxAeQG96LtLQa0aGJ6+jJYi8p39WoiDg6aptCVhsI2409vjcLJ4q7M3D9N+PgdCcgHc5PBxdP9nC/NlVJSo6lPwWB6kO74TNaIBjjLw09fr2zQd6NyZrtnL/qs7ToICx2xxOFKYtoKxu3qQiD/uI3uqBFbpgL9x3Sv3FTcP5Xz247bXianhL7+9vitMU7q3cY6Loc8dVdW/D6Y7ggVu4fSy83Nng8osztrxjKCwWVoWANV30f1RtODQ+C5x7zWtL9szDxq/FmixWdaMb5RLRxVTW55LpwUzfpjj0GPng36jL5kMt12bhKvBO670h13pzGOpuNHPV7d6LZO9gdQDwSDY5qSyVdD9bkPMM0BO0bMLMAog==

Thank you. I'll try your suggestion.

On 24 Feb 2016 15:36, "Sebastien Loriot (GeometryFactory)" <> wrote:
Vertex handles inserted will remain valid in both scenarios
(refine_Delaunay_mesh_2 is using Delaunay_mesher_2 internally).
About the process, you can create your domain, insert your fixed points
and then only once the mesher.

Sebastien.

On 02/10/2016 10:52 AM, rytis wrote:
Hello,

I am a beginner with CGAL which I like very much so far. I think there
is an issue (of documentation) with constrained Delaunay triangulation.
The issue is below, but first, I want to describe the problem. If
someone could indicate how solve my indexation problem, that would be
nice too!

Problem:
I need to create a 2d mesh from "thin air" + boundaries (as if by
Delaunay_mesher_2) but I must insert a few points which remain constant
(that don't disappear or change their coordinates -- as if they were
"punctual" constraints).


Issue:
Is it correct that,
A) refine_Delaunay_mesh_2 takes vertices as constants, so that, in
particular, vertex_handles remain valid throughout the refinement
procedure, and that no new vertices are added, whereas
B) Delaunay_mesher_2 invalidates vertex_handles, by possibly adding
and/or removing some of them in the process of refinement?

So far this is what I tried to solve my problem:

1: insert boundary nodes and set boundary constraints,
2: run B) to generate a "new" triangulation
3: insert my "custom" vertices,
4: run A) (expecting that vertices remain constant)

I have to generate two types of output for nodes and links, respectively
ID,x,y
...

ID,ID
...

For that, I keep vertex handles in a map like this:
map<string,Vertex_handle> node_map;

this map is updated  after each step like so:
node_map[ID]=Vertex_handle(Vertex_iterator);

which allows the retrieval of the vertex coordinates by their ID.

I noticed that during step 4 some of the vertex handles were not matched
by an ID which makes me think that A) is incorrect (Help!) I checked
that those vertices inserted at step 3 were still present though (Help....)

Thanks,
Rytis


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





Archive powered by MHonArc 2.6.18.

Top of Page