Subject: CGAL users discussion list
List archive
- From: Mariette Yvinec <>
- To:
- Subject: Re: [cgal-discuss] meshing 2d
- Date: Mon, 17 May 2010 14:38:35 +0200
A better way, may be to mark the cells
with some indexes of the subdomain they belong to.
Using a propagation algorithm through cell adjacency
it is easy to mark the cells.
Then knowing on which boundary
are the edges is straightforward.
The main advantage is to avoid using
geometric tests like is_colinear()
which might yield a false result if you use floating point computations
and is costly when using an EPIC Kernel.
Hansjoerg Seybold wrote:
hello,
thanks a lot. This would help to determine weather boundary or not.
For a fem mesh i further have to know on which boundary the
nodes are to apply the appropriate boundary conditions in the
calculations. Also for internal boundary conditions the is_in_domain()
is not sufficient.
I was thinking to check all vvertices with all edges of the PLGS with
the basic functions
are_ordered_along line() or is_colinear() because i could not find a
better
point location method which returns the edge of the PLGS for a given
vertex. Please let me know if you have a better suggestion for me.
best hj
On Mon, May 17, 2010 at 10:51 AM, Mariette
Yvinec <>
wrote:
You may identify the boundary edges as edges incident
to a single face that is in the doamain.
This can be tested using the is_in_domain() member function
of the facet.
Hansjoerg Seybold wrote:
Hello,
i could not find a function /flag in the mesher that identifies the
boundary
edges, meaning edges that are part of the initial PLGS.
I saw that for a 3D mesh a IS_BOUNDARY type exists. but for 2d meshes i
could not find an equivalence.
Do i have to check the boundary query by hand using geometric queries
for the vertex points of the mesh and the edges of the original PLGS?
thanks a lot
hansjoerg
On Thu, Apr 22, 2010 at 3:31 PM,
Hansjoerg
Seybold <>
wrote:
Thanks
for the reply.
the mesh_domain thus is a simple triangulation object.
One question remains: how do i identify if the newly inserted vertices
(edge) is lying on the original boundary
triangulation identify the boudnary of the mesh domain
thanks
On Thu, Apr 22, 2010 at 2:01 PM, Ramin
H <>
wrote:
Look
here:
http://www.cgal.org/Manual/last/doc_html/cgal_manual/Triangulation_2_ref/Class_Triangulation_2.html
typedef
Tds::Face_iterator |
All_faces_iterator; |
iterator
over
all faces. |
typedef
Tds::Edge_iterator |
All_edges_iterator; |
iterator
over
all edges |
typedef
Tds::Vertex_iterator |
All_vertices_iterator; |
iterator
over
all vertices |
-Ramin
On Thu, Apr 22, 2010 at 5:44 AM,
Hansjoerg Seybold <>
wrote:
hello,
I am trying to make a simple 2D mesher for a FEM code using the cgal
meshing capabilities.
I experimenteeeed with the meshing examples and demos. What i could
not find in the
documentation and the examples is
How can i iterate through the mesh, Faces, Vertices, Edges?
Is there a method to determine if a edge/vertex lies on a boundary
(meaning part of the CDT which is used as input
for the mesher?
Are there default io streams for the meshs?
It would be nice if anybody has some suggestions references for
documentation etc on this topic.
Thanks a lot
best hansjoerg
--
Mariette Yvinec
Geometrica project team
INRIA Sophia-Antipolis
--
Mariette Yvinec
Geometrica project team
INRIA Sophia-Antipolis
|
Archive powered by MHonArc 2.6.16.