From: Marc Fehling Date: Tue, 3 Nov 2020 01:53:06 +0000 (-0700) Subject: Functions of hp::Refinement now only available for DoFHandler in hp mode. X-Git-Tag: v9.3.0-rc1~916^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=4ddbd602e45d25b48432a483db0aa1414d2a2d74;p=dealii.git Functions of hp::Refinement now only available for DoFHandler in hp mode. --- diff --git a/source/hp/refinement.cc b/source/hp/refinement.cc index 7c9c02a494..dd976ef6ec 100644 --- a/source/hp/refinement.cc +++ b/source/hp/refinement.cc @@ -43,6 +43,11 @@ namespace hp void full_p_adaptivity(const dealii::DoFHandler &dof_handler) { + Assert( + dof_handler.has_hp_capabilities(), + (typename dealii::DoFHandler::ExcNotAvailableWithoutHP())); + std::vector p_flags( dof_handler.get_triangulation().n_active_cells(), true); @@ -57,6 +62,10 @@ namespace hp const dealii::DoFHandler &dof_handler, const std::vector & p_flags) { + Assert( + dof_handler.has_hp_capabilities(), + (typename dealii::DoFHandler::ExcNotAvailableWithoutHP())); AssertDimension(dof_handler.get_triangulation().n_active_cells(), p_flags.size()); @@ -99,6 +108,10 @@ namespace hp const ComparisonFunction::type> &compare_coarsen) { + Assert( + dof_handler.has_hp_capabilities(), + (typename dealii::DoFHandler::ExcNotAvailableWithoutHP())); AssertDimension(dof_handler.get_triangulation().n_active_cells(), criteria.size()); @@ -131,6 +144,10 @@ namespace hp const ComparisonFunction::type> &compare_coarsen) { + Assert( + dof_handler.has_hp_capabilities(), + (typename dealii::DoFHandler::ExcNotAvailableWithoutHP())); AssertDimension(dof_handler.get_triangulation().n_active_cells(), criteria.size()); Assert((p_refine_fraction >= 0) && (p_refine_fraction <= 1), @@ -221,6 +238,10 @@ namespace hp const ComparisonFunction::type> &compare_coarsen) { + Assert( + dof_handler.has_hp_capabilities(), + (typename dealii::DoFHandler::ExcNotAvailableWithoutHP())); AssertDimension(dof_handler.get_triangulation().n_active_cells(), criteria.size()); Assert((p_refine_fraction >= 0) && (p_refine_fraction <= 1), @@ -410,6 +431,10 @@ namespace hp const dealii::DoFHandler &dof_handler, const Vector & sobolev_indices) { + Assert( + dof_handler.has_hp_capabilities(), + (typename dealii::DoFHandler::ExcNotAvailableWithoutHP())); AssertDimension(dof_handler.get_triangulation().n_active_cells(), sobolev_indices.size()); @@ -465,6 +490,10 @@ namespace hp const ComparisonFunction::type> &compare_coarsen) { + Assert( + dof_handler.has_hp_capabilities(), + (typename dealii::DoFHandler::ExcNotAvailableWithoutHP())); AssertDimension(dof_handler.get_triangulation().n_active_cells(), criteria.size()); AssertDimension(dof_handler.get_triangulation().n_active_cells(), @@ -609,6 +638,11 @@ namespace hp void force_p_over_h(const dealii::DoFHandler &dof_handler) { + Assert( + dof_handler.has_hp_capabilities(), + (typename dealii::DoFHandler::ExcNotAvailableWithoutHP())); + for (const auto &cell : dof_handler.active_cell_iterators()) if (cell->is_locally_owned() && cell->future_fe_index_set()) { @@ -623,6 +657,11 @@ namespace hp void choose_p_over_h(const dealii::DoFHandler &dof_handler) { + Assert( + dof_handler.has_hp_capabilities(), + (typename dealii::DoFHandler::ExcNotAvailableWithoutHP())); + // Siblings of cells to be coarsened may not be owned by the same // processor. We will exchange coarsening flags on ghost cells and // temporarily store them.