bi.print(deallog);
{
- // check that refinement_edge_boundary_indices (deprecated) is really
- // the intersection of the sets above:
- IndexSet intersect = rei & bi;
- std::vector<bool> ref = mg_constrained_dofs.get_refinement_edge_boundary_indices()[level];
-
- unsigned int idx = 0;
- for (std::vector<bool>::iterator it = ref.begin(); it != ref.end(); ++it, ++idx)
- {
- if ((*it == true) != intersect.is_element(idx))
- deallog << "mismatch idx=" << idx << " "
- << (*it == true) << " " << intersect.is_element(idx) << std::endl;
- }
+ // check that refinement_edge_boundary_indices (deprecated) is really
+ // the intersection of the sets above:
+ IndexSet intersect = rei & bi;
+ std::vector<bool> ref = mg_constrained_dofs.get_refinement_edge_boundary_indices()[level];
+
+ unsigned int idx = 0;
+ for (std::vector<bool>::iterator it = ref.begin(); it != ref.end(); ++it, ++idx)
+ {
+ if ((*it == true) != intersect.is_element(idx))
+ deallog << "mismatch idx=" << idx << " "
+ << (*it == true) << " " << intersect.is_element(idx) << std::endl;
+ }
}
-
-
+
+
IndexSet relevant;
DoFTools::extract_locally_relevant_mg_dofs (dofh,
relevant, level);
setup_tria(tr);
if (false)
- {
- DataOut<dim> data_out;
- Vector<float> subdomain (tr.n_active_cells());
- for (unsigned int i=0; i<subdomain.size(); ++i)
- subdomain(i) = tr.locally_owned_subdomain();
- data_out.attach_triangulation (tr);
- data_out.add_data_vector (subdomain, "subdomain");
- data_out.build_patches (0);
- const std::string filename = ("solution." +
- Utilities::int_to_string
- (tr.locally_owned_subdomain(), 4) +
- ".vtu");
- std::ofstream output (filename.c_str());
- data_out.write_vtu (output);
+ {
+ DataOut<dim> data_out;
+ Vector<float> subdomain (tr.n_active_cells());
+ for (unsigned int i=0; i<subdomain.size(); ++i)
+ subdomain(i) = tr.locally_owned_subdomain();
+ data_out.attach_triangulation (tr);
+ data_out.add_data_vector (subdomain, "subdomain");
+ data_out.build_patches (0);
+ const std::string filename = ("solution." +
+ Utilities::int_to_string
+ (tr.locally_owned_subdomain(), 4) +
+ ".vtu");
+ std::ofstream output (filename.c_str());
+ data_out.write_vtu (output);
- }
+ }
ZeroFunction<dim> zero;