Dear Timo, many hanks for your reply, but I have one question, as I understand: center means the center of the edge, and the condition here if ( cell->face(face)->at_boundary() && cell->face(face)->center()(2)>0.0 && abs(cell->face(face)->center()(1))<0.5 ) he looks for one point in each edge (the center of the edge) which satisfy this condition, but what I need is : the whole points in each edge which satisfy this condition, I mean if we can use some thing like Point insted of center?! Best Regards, Ahmad
--- Timo Heister <[email protected]> schrieb am Mo, 24.1.2011: Von: Timo Heister <[email protected]> Betreff: Re: [deal.II] boundary_indicator An: "Ahmad Al Nabulsi" <[email protected]> CC: "deal.ii" <[email protected]> Datum: Montag, 24. Januar, 2011 07:48 Uhr Hello Ahmad, > [snip] Use something like this: for (; cell!=endc; ++cell) for (unsigned int face=0; face<GeometryInfo<dim>::faces_per_cell; ++face) if ( cell->face(face)->at_boundary() && cell->face(face)->center()(2)>0.0 && abs(cell->face(face)->center()(1))<0.5 ) cell->face(face)->set_boundary_indicator (1); Best, Timo -- Timo Heister http://num.math.uni-goettingen.de/~heister
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
