constraints.close ();
hanging_node_constraints.close ();
- check(constraints);
- check(hanging_node_constraints);
+ //check(constraints);
+ //check(hanging_node_constraints);
CompressedSimpleSparsityPattern csp(mg_dof_handler.n_dofs(), mg_dof_handler.n_dofs());
boundary_constraints[level].add_lines (mg_constrained_dofs.get_refinement_edge_indices(level));
boundary_constraints[level].add_lines (mg_constrained_dofs.get_boundary_indices()[level]);
- check(boundary_constraints[level]);
+
boundary_constraints[level].close ();
+ check(boundary_constraints[level]);
- // boundary_interface_constraints[level]
- //.add_lines (mg_constrained_dofs.get_refinement_edge_boundary_indices (level));
+ boundary_interface_constraints[level]
+ .add_lines (mg_constrained_dofs.get_refinement_edge_boundary_indices (level));
boundary_interface_constraints[level].close ();
+ check(boundary_interface_constraints[level]);
}
// Now that we're done with most of our preliminaries, let's start
<< "Aborting!" << std::endl
<< "----------------------------------------------------"
<< std::endl;
-
- return 1;
+ throw;
}
catch (...)
{
<< "Aborting!" << std::endl
<< "----------------------------------------------------"
<< std::endl;
- return 1;
+ throw;
}
return 0;