Hi all, In two dimensional geometry I have the following rectangle [-1,1]x[-1,1]. I am going to specify a boundary conditions on the boundary edges, I would like to choose the part which located on y=1, and to take the part (y=1 && x \in [-0.5,0.5] ), so I need the set of points which have this condition.
I looped over the whole cells, then I looped over the whole faces (lines in 2d), then I put my condition, then the points which satisfy this condition will have the name (boundary_indicator(1)). I tried to write this in c++: for (unsigned int face=0; f<GeometryInfo<dim>::faces_per_cell; ++f) if (std::fabs(cell->face(f)->( Point< dim >::unit_vector(1)- (1) < 1e-12) && (-0.5< Point< dim >::unit_vector(0)<2>(0)<0.5))) cell->face(f)->set_boundary_indicator (1); but unfortunately it gave me error msg., please help me to fix this error... Best Regards, Ahmad
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
