> Is there an easy way of working out whether a particular vertex is on the
> boundary of the mesh? You can check this for a particular cell with
>
> cell->at_boundary()
>
> and with a face
> cell->face(face)->at_boundary()
>
> but I can't seem to find out how to do this for a particular vertex.
There is no such function that retrieves this information from a place already
pre-computed. However, it should be relatively straightforward to write this:
std::vector<bool> vertex_is_at_boundary (tria.n_vertices, false);
for (cell=...)
for (face=...)
if (cell->face(face)->at_boundary())
for (vertex=...)
vertex_is_at_boundary[cell->face(face)->vertex_index(v)] = true;
You can then poke into this array next time you want to figure out whether a
vertex is at the boundary or not.
Best
W.
-------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
