template <int dim>
{
public:
BoundaryIds()
{}
virtual void
evaluate_scalar_field(
{
// Get the cell and face we are currently dealing with:
inputs.template get_cell<dim>();
// Then fill the output fields with the boundary_id of the face
for (auto &output : computed_quantities)
{
output(0) = cell->face(face)->boundary_id();
}
}
};