// should be smaller than the smallest cell
// diameter.
- typename Triangulation<dim>::face_iterator face = tria.begin_face(),
- endface = tria.end_face();
+ typename Triangulation<dim,spacedim>::face_iterator face = tria.begin_face(),
+ endface = tria.end_face();
for (; face!=endface; ++face)
{
- if (face->boundary_indicator() == 0)
+ if (face->boundary_id() == 0)
{
- const Point<dim> center (face->center());
+ const Point<spacedim> center (face->center());
if (std::abs(center(0)-p1[0]) < epsilon)
- face->set_boundary_indicator(0);
+ face->set_boundary_id(0);
else if (std::abs(center(0) - p2[0]) < epsilon)
- face->set_boundary_indicator(1);
+ face->set_boundary_id(1);
else if (dim > 1 && std::abs(center(1) - p1[1]) < epsilon)
- face->set_boundary_indicator(2);
+ face->set_boundary_id(2);
else if (dim > 1 && std::abs(center(1) - p2[1]) < epsilon)
- face->set_boundary_indicator(3);
+ face->set_boundary_id(3);
else if (dim > 2 && std::abs(center(2) - p1[2]) < epsilon)
- face->set_boundary_indicator(4);
+ face->set_boundary_id(4);
else if (dim > 2 && std::abs(center(2) - p2[2]) < epsilon)
- face->set_boundary_indicator(5);
+ face->set_boundary_id(5);
else
// triangulation says it
// is on the boundary,