From a501894dc24a677998ba57fb6392547cb7d0fc71 Mon Sep 17 00:00:00 2001 From: Marc Fehling Date: Thu, 28 Feb 2019 13:57:06 +0100 Subject: [PATCH] hp::distribute_dofs: Move 'set_fe' call to the place where active_fe_indices are initialized. --- source/hp/dof_handler.cc | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) 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); -- 2.39.5