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
