CGAL users discussion list

[Zhonggui Chen <>]

Hi Laurent,

Thanks for your reply!

Actually I want to initialize the C3t3 as the input for the CGAL::perturb_mesh_3() function. 
I tried your method, but it doesn't work. The output .mesh file containts only the points I inserted, no facets or tetrahehra. Any other information I should provide to the C3t3?

 Here are my codes:
------------------------------------------
    Polyhedron polyhedron;
    Build_mesh<HalfedgeDS> builder(boundary);
    polyhedron.delegate(builder);
    Mesh_domain domain(polyhedron);

    C3t3 c3t3;
    for (int i=0; i<tetmesh->nb_vertices(); i++)
    {
      // insert points
 c3t3.triangulation().insert(Point(tetmesh->vertices_[i].point_[0],tetmesh->vertices_[i].point_[1],tetmesh->vertices_[i].point_[2]));
     }

    CGAL::perturb_mesh_3(c3t3, domain, time_limit=300, sliver_bound=10);

    std::ofstream medit_file("out_1.mesh");
    c3t3.output_to_medit(medit_file);
------------------------------------------------------------------

--zhong

On Sun, Apr 11, 2010 at 8:50 PM, Laurent Rineau (GeometryFactory) <> wrote:

Le Dimanche 11 Avril 2010 14:33:24, Zhonggui Chen a écrit :

> Hi,
>
> I have some initial points in the domain bounded by a surface mesh. I'd
> like to use the C3t3 to optimize the locations of the points to get a
> better 3D meshing of the domain. How to insert the initial points into
> C3t3?
>
> Any suggestion would be appreciated!!!

Use :
 c3t3.triangulation().insert(p)

--
Laurent Rineau, PhD
Release Manager of the CGAL Project       http://www.cgal.org/
R&D Engineer at GeometryFactory           http://www.geometryfactory.com/

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