From: Marc Fehling Date: Thu, 28 Feb 2019 12:57:06 +0000 (+0100) Subject: hp::distribute_dofs: Move 'set_fe' call to the place where active_fe_indices are... X-Git-Tag: v9.1.0-rc1~310^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F7717%2Fhead;p=dealii.git hp::distribute_dofs: Move 'set_fe' call to the place where active_fe_indices are initialized. --- diff --git a/source/hp/dof_handler.cc b/source/hp/dof_handler.cc index 35e4a17666..f2cfe2686d 100644 --- a/source/hp/dof_handler.cc +++ b/source/hp/dof_handler.cc @@ -1431,9 +1431,6 @@ namespace hp ExcMessage("The Triangulation you are using is empty!")); Assert(ff.size() > 0, ExcMessage("The hp::FECollection given is empty!")); - // first, assign the fe_collection - set_fe(ff); - // If an underlying shared::Tria allows artificial cells, // then save the current set of subdomain ids, and set // subdomain ids to the "true" owner of each cell. we later @@ -1457,7 +1454,10 @@ namespace hp } } - // then allocate space for all tables + // assign the fe_collection and initialize all active_fe_indices + set_fe(ff); + + // then allocate space for all the other tables dealii::internal::hp::DoFHandlerImplementation::Implementation:: reserve_space(*this);