CGAL users discussion list

["Sebastien Loriot (GeometryFactory)" <>]

Could you provide a minimal example that we can run so that we can
reproduce the issue?

I don't see any issue in the snippet provided.

Thanks,

Sebastien.

On 9/7/20 8:33 AM, MariusZeinhofer wrote:
> Dear all,
>
> I am trying to mesh a 3d mesh using the create_implicit_mesh_domain method
> and I am running into problems with the sharp representations of 1d
> features.
>
> My 3d object is a set union of a cylinder with a segement of another
> cylinder of larger radius (see the image attached), and I want these edges
> to be sharp. I am working my way along the example 3.6.2 in the CGAL manual
> (Implicit Domain with 1D Features, the intersecting balls). The problem I'm
> getting is that my edges are "too sharp" represented which makes the mesh
> too big for my purposes. The mesh I get is extremly fine at the edges, see
> the following picture:
>
> <http://cgal-discuss.949826.n4.nabble.com/file/t376248/CylinderSmall.png>
>
> What confuses me is, that not all of the sharp edges are problematic and
> also if I am just meshing the smaller cylinder I'm not getting this issue at
> the top and bottom of the cylinder, even though I am using the same code.
>
> My mesh criteria are:
> Mesh_criteria criteria(edge_size = 0.29, facet_angle = 25, facet_size = 0.3,
> cell_radius_edge_ratio = 3.5, cell_size = .5, facet_distance = 0.4);
>
> The way I'm introducing the sharp features is:
>
> typedef std::vector<Point>        Polyline_3;
> typedef std::list<Polyline_3>       Polylines;
>
> for (int i = 0; i < 360; ++i)
>      {
>          Point p(15.0, 10 * std::cos(i * CGAL_PI / 180), 10 * std::sin(i *
> CGAL_PI / 360));
>          polyline.push_back(p);
>      }
>      polyline.push_back(polyline.front()); // close the line
>      polylines.push_back(polyline);
>
> and in a similar fashion for the other sharp features.
>
> Thank you very much for your help in advance!
> Cheers
> Marius
>
>
>
> --
> Sent from: http://cgal-discuss.949826.n4.nabble.com/