CylindricalManifold set manually does not work correctly

[Chengjiang Yin]

Hi all,

I'm trying to make a cylinder with a hole inside it.

I use the subdivided_cylinder to generate a cylinder with 3 slice and remove the central cell by create_triangulation_with_removed_cells. However, after I setup the boundary and manifold, the mesh refinement does not seem to work correctly.

Here I show the mesh with two global refinements respectively. It looks like the manifold is now working for the first refinement (like flat manifold) and so the second refinement could not generate as a single cylinder. 

Attached is a demo to generate these two meshes. I'm not sure if this is a bug of CylindricalManifold or I do not use it in the right way.

Regards,

Chengjiang Yin

[Chengjiang Yin]

I think I just find the problem. According to "problem with CylindricalManifold<3> after removing cells from a triangulation", using set_all_boundary_ids instead of set_boundary_id to set all of its children (line and vertex) to given boundary id.

Besides, you should pay extra attention to the intersection line of the bottom and the side boundary, which should be cylindrical manifold rather than flat manifold.