typename BaseClass::ExcVectorDoesNotMatch());
// this function really only makes
- // sense on non-active objects if
+ // sense if either a) there are
+ // degrees of freedom defined on
+ // the present object, or b) the
+ // object is non-active objects but
// all degrees of freedom are
// located on vertices, since
// otherwise there are degrees of
// freedom on sub-objects which are
// not allocated for this
// non-active thing
-// TODO: This assertion leads to problems with the do_make_hp_hanging_node_constraints
-// method!. Check how this assertion can be stated that it does not create a conflict
-// with the above mentioned method.
-/*
- Assert (!this->has_children() ||
+ Assert (this->fe_index_is_active (fe_index)
+ ||
(this->dof_handler->get_fe()[fe_index].dofs_per_cell ==
2*this->dof_handler->get_fe()[fe_index].dofs_per_vertex),
- typename BaseClass::ExcNotActive());
-*/
+ ExcInternalError());
const unsigned int dofs_per_vertex = this->dof_handler->get_fe()[fe_index].dofs_per_vertex,
dofs_per_line = this->dof_handler->get_fe()[fe_index].dofs_per_line;
typename BaseClass::ExcVectorDoesNotMatch());
// this function really only makes
- // sense on non-active objects if
+ // sense if either a) there are
+ // degrees of freedom defined on
+ // the present object, or b) the
+ // object is non-active objects but
// all degrees of freedom are
// located on vertices, since
// otherwise there are degrees of
// freedom on sub-objects which are
// not allocated for this
// non-active thing
-
-// TODO: This assertion leads to problems with the do_make_hp_hanging_node_constraints
-// method!. Check how this assertion can be stated that it does not create a conflict
-// with the above mentioned method.
-/*
- Assert (!this->has_children() ||
+ Assert (this->fe_index_is_active (fe_index)
+ ||
(this->dof_handler->get_fe()[fe_index].dofs_per_cell ==
4*this->dof_handler->get_fe()[fe_index].dofs_per_vertex),
- typename BaseClass::ExcNotActive());
-*/
+ ExcInternalError());
Assert (static_cast<unsigned int>(this->present_level) < this->dof_handler->levels.size(),
ExcMessage ("DoFHandler not initialized"));
ExcMessage ("DoFHandler not initialized"));
// this function really only makes
- // sense on non-active objects if
+ // sense if either a) there are
+ // degrees of freedom defined on
+ // the present object, or b) the
+ // object is non-active objects but
// all degrees of freedom are
// located on vertices, since
// otherwise there are degrees of
// freedom on sub-objects which are
// not allocated for this
// non-active thing
-// TODO: This assertion leads to problems with the do_make_hp_hanging_node_constraints
-// method!. Check how this assertion can be stated that it does not create a conflict
-// with the above mentioned method.
-/*
- Assert (!this->has_children() ||
+ Assert (this->fe_index_is_active (fe_index)
+ ||
(this->dof_handler->get_fe()[fe_index].dofs_per_cell ==
8*this->dof_handler->get_fe()[fe_index].dofs_per_vertex),
- typename BaseClass::ExcNotActive());
-*/
+ ExcInternalError());
+
Assert (static_cast<unsigned int>(this->present_level) < this->dof_handler->levels.size(),
ExcMessage ("DoFHandler not initialized"));