CGAL users discussion list

[rdzhao <>]

Hi Sebastien!

I'm facing another problem. 

I would like to create a tetrahedral mesh from a polyhedral surface with
make_mesh_3 and refine_mesh_3. The model is a benzene with 6 C atoms and 6 H
atoms. Those atoms are disconnected balls in the input polyhedral surface
mesh (benzene.off). 

However, by directly calling make_mesh_3, there will be missing components.
So I resort to another post:
http://cgal-discuss.949826.n4.nabble.com/make-mesh-3-with-features-initial-points-for-C3T3-td4661118.html
This pose tells that we need to set initial points as features to guide the
meshing process. So I modified the code with:

C3t3 c3t3 = CGAL::make_mesh_3<C3t3>(domain, criteria, features(domain),
no_lloyd(), no_odt(), no_perturb(), no_exude());
        
typedef C3t3::Triangulation::Point Weighted_point;
   for(int i=0; i<atoms.size(); i+=3){
        Weighted_point p(Weighted_point::Point(atoms[i], atoms[i+1], 
atoms[i+2]));
        Tr::Vertex_handle vh = c3t3.triangulation().insert(p);
        c3t3.set_dimension(vh, 3);
}
CGAL::refine_mesh_3(c3t3, domain, criteria);

The idea is adding atom positions as initial points. But the code still
misses one component.

A minimal code is attached here with model file (benzene.off) and atoms file
(benzene.xyz). It should be easy to directly compile and execute.
Can you help dig out the problem?
Thank you so much!
Rundong

&lt;nabble_a href=&quot;Archive.zip&quot;>Archive.zip 





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