CGAL users discussion list

[rdzhao <>]

Hi everyone, 

I working working on generating tetrahedral mesh from input polyhedral
surface mesh. 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
in the final generated tet mesh.
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 post 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 certain H atom).

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.

Suppose that the executable file is called &quot;Meshing&quot;. The command
to run the code is:
./Mesh benzene.xyz benzene.off 0.2

.xyz file provides the atom locations.
.off file provides the model containing 12 balls (6 C atoms and 6 H atoms).
0.2 is a parameter to roughly control the element size.

The output contains a .off file for the surface of the tet meshing, where we
can find that there is missing components compared with the input .off file. 

Is there anyone that can 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/