ConstraintMatrix cm;
DoFTools::make_hanging_node_constraints(dh,cm);
- IndexSet support = DoFTools::extract_dofs_with_support_contained_within(dh, pred_d<dim>, cm);
+ IndexSet support = DoFTools::extract_dofs_with_support_contained_within(dh,
+ std::function<bool (const typename DoFHandler<dim>::active_cell_iterator &)>(&pred_d<dim>),
+ cm);
support.print(deallog);
// print grid and DoFs for visual inspection
DoFTools::make_hanging_node_constraints(dh,cm);
cm.close ();
- const IndexSet support = DoFTools::extract_dofs_with_support_contained_within(dh, pred_d<dim>, cm);
+ const IndexSet support = DoFTools::extract_dofs_with_support_contained_within(dh,
+ std::function<bool (const typename DoFHandler<dim>::active_cell_iterator &)>(&pred_d<dim>),
+ cm);
const IndexSet support_local = support & dh.locally_owned_dofs();
deallog << support.n_elements() << std::endl;
DoFTools::make_hanging_node_constraints(dh,cm);
IndexSet support = left ?
- DoFTools::extract_dofs_with_support_contained_within(dh, pred_left<dim>, cm) :
- DoFTools::extract_dofs_with_support_contained_within(dh, pred_right<dim>, cm);
+ DoFTools::extract_dofs_with_support_contained_within(dh,
+ std::function<bool (const typename DoFHandler<dim>::active_cell_iterator &)>(&pred_left<dim>),
+ cm) :
+ DoFTools::extract_dofs_with_support_contained_within(dh,
+ std::function<bool (const typename DoFHandler<dim>::active_cell_iterator &)>(&pred_right<dim>),
+ cm);
support.print(deallog);
// print grid and DoFs for visual inspection
std::vector<types::global_dof_index> local_dof_indices(fe.dofs_per_cell);
// get support on the predicate
- IndexSet support = DoFTools::extract_dofs_with_support_contained_within(dh, pred_d<dim>, cm);
+ IndexSet support = DoFTools::extract_dofs_with_support_contained_within(dh,
+ std::function<bool (const typename DoFHandler<dim>::active_cell_iterator &)>(&pred_d<dim>),
+ cm);
IndexSet local_support = support & locally_owned_set;
// rhs vectors: