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