From da66aa81aa17700ee83c722d1af3eb58150d3d74 Mon Sep 17 00:00:00 2001 From: Timo Heister Date: Sun, 4 Oct 2015 17:03:59 -0400 Subject: [PATCH] step-50: avoid get_refinement_edge_boundary_indices --- examples/step-50/step-50.cc | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/examples/step-50/step-50.cc b/examples/step-50/step-50.cc index 205ef37210..05964f5c92 100644 --- a/examples/step-50/step-50.cc +++ b/examples/step-50/step-50.cc @@ -625,8 +625,12 @@ namespace Step50 boundary_constraints[level].close (); + boundary_interface_constraints[level].reinit(dofset); + // compute indices that are refinement edge and boundary (aka refinement_edge_boundary_indices): + IndexSet idxset = mg_constrained_dofs.get_refinement_edge_indices(level) & mg_constrained_dofs.get_boundary_indices(level); + boundary_interface_constraints[level] - .add_lines (mg_constrained_dofs.get_refinement_edge_boundary_indices()[ (level)]); + .add_lines (idxset); boundary_interface_constraints[level].close (); } -- 2.39.5