> 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?!
You can ask cell->face(face)->vertex(0) cell->face(face)->vertex(1) to get at the two vertices of a face in 2d. Best W. -- ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
