CGAL users discussion list

[Song Lin <>]

Hi,

As Weighted_alpha_shape_3 is a subset of Regular_triangulation_3, anyone
knows how to decide which cell(face,edge...) of Triangulation_3 is
included in Alpha_shape_3 and which is not?

My codes are mainly like these:(based on the example
"ex_weighted_alpha_shapes_3.cpp")

typedef CGAL::Alpha_shape_vertex_base_3<Gt>         Vb;
typedef CGAL::Alpha_shape_cell_base_3<Gt>           Fb;
typedef CGAL::Triangulation_data_structure_3<Vb,Fb> Tds;
typedef CGAL::Regular_triangulation_3<Gt,Tds>       Triangulation_3;
typedef CGAL::Alpha_shape_3<Triangulation_3>        Alpha_shape_3;

int main()
{
list<Weighted_point> lwp;

Alpha_shape_3  as(lwp.begin(), lwp.end(), 0, Alpha_shape_3::GENERAL);

list<Edge> edges;
as.get_alpha_shape_edges(back_inserter(edges),Alpha_shape_3::REGULAR);
as.get_alpha_shape_edges(back_inserter(edges),Alpha_shape_3::SINGULAR);

Alpha_shape_3::Finite_edges_iterator vit;
for (vit = as.finite_edges_begin(); vit != as.finite_edges_end();vit++)
{
        ##########(I want to decide which edge of Triangulation_3 is included
in Alpha_shape_3 but I don't know how. Seems that the two iterators are
not of same kind)##########
}

I'm a new user of CGAL. Sorry if noob. Maybe my idea was totaly wrong
and there were more simple solutions. Anyone can help me with this
problem? Thank you!


Lin

-- 
Center for Computational and Systems Biology, 
Institute of Biophysics, CAS.
15 Datun Road, Chaoyang District
Beijing, 100101
P.R.China