Hello Stephane,

yes, you are setting the boundary of the Triangulation to CylinderBoundary right at the beginning of the program. Then new vertices will be put on the "real" boundary of the cylinder, if you refine the triangulation. There is no need to use the CylinderBoundary in every refine step. You just have to asure that the CylinderBoundary object lives as long as the triangulation. Probably you want to have a look on tutorial step-1. There the same procedure is shown for a circle.

Best Regards,
Markus



Am 24.09.10 09:28, schrieb veys:
Hi Markus :)
Yes I know that, but what I want is to adapt the mesh and force it to
have the shape of cylinder, or by using triangulation.refine_global (x)
the shape will not be a "perfect" cylinder beacause boundaries of mesh
will not be adapted isn't it ?
That's why I use CylinderBoundary<dim>  (am I clear ?)
Bests
to refine a grid it suffices to use

triangulation.refine_global (x);

where x is an integer.

_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii

Reply via email to