57 template <
int dim,
int spacedim>
60 const ::DoFHandler<dim, spacedim> &dof_handler,
62 std::vector<std::vector<
63 std::pair<types::global_dof_index, types::global_dof_index>>>
65 std::vector<std::vector<
66 std::pair<types::global_dof_index, types::global_dof_index>>>
67 ©_indices_global_mine,
68 std::vector<std::vector<
69 std::pair<types::global_dof_index, types::global_dof_index>>>
70 & copy_indices_level_mine,
71 const bool skip_interface_dofs)
84 std::vector<DoFPair> send_data_temp;
86 const unsigned int n_levels =
87 dof_handler.get_triangulation().n_global_levels();
88 copy_indices.resize(n_levels);
89 copy_indices_global_mine.resize(n_levels);
90 copy_indices_level_mine.resize(n_levels);
94 const unsigned int dofs_per_cell = dof_handler.get_fe().n_dofs_per_cell();
95 std::vector<types::global_dof_index> global_dof_indices(dofs_per_cell);
96 std::vector<types::global_dof_index> level_dof_indices(dofs_per_cell);
100 std::vector<bool> dof_touched(globally_relevant.
n_elements(),
false);
106 std::vector<types::global_dof_index> unrolled_copy_indices;
108 copy_indices_level_mine[
level].clear();
109 copy_indices_global_mine[
level].clear();
111 for (
const auto &level_cell :
112 dof_handler.active_cell_iterators_on_level(
level))
114 if (dof_handler.get_triangulation().locally_owned_subdomain() !=
116 (level_cell->level_subdomain_id() ==
118 level_cell->subdomain_id() ==
122 unrolled_copy_indices.resize(
123 dof_handler.locally_owned_dofs().n_elements(),
128 level_cell->get_dof_indices(global_dof_indices);
129 level_cell->get_mg_dof_indices(level_dof_indices);
131 for (
unsigned int i = 0; i < dofs_per_cell; ++i)
135 if (skip_interface_dofs && mg_constrained_dofs !=
nullptr &&
137 level, level_dof_indices[i]))
145 dof_handler.locally_owned_dofs().is_element(
146 global_dof_indices[i]);
148 dof_handler.locally_owned_mg_dofs(
level).is_element(
149 level_dof_indices[i]);
151 if (global_mine && level_mine)
156 unrolled_copy_indices[dof_handler.locally_owned_dofs()
158 global_dof_indices[i])] =
159 level_dof_indices[i];
168 global_dof_indices[i]);
172 if (dof_touched[relevant_idx] ==
false)
176 copy_indices_global_mine[
level].emplace_back(
177 global_dof_indices[i], level_dof_indices[i]);
180 send_data_temp.emplace_back(
level,
181 global_dof_indices[i],
182 level_dof_indices[i]);
189 dof_touched[relevant_idx] =
true;
197 if (!unrolled_copy_indices.empty())
199 copy_indices[
level].clear();
203 if (copy_indices_global_mine[
level].empty())
204 copy_indices[
level].reserve(unrolled_copy_indices.size());
209 for (
unsigned int i = 0; i < unrolled_copy_indices.size(); ++i)
211 copy_indices[
level].emplace_back(
212 dof_handler.locally_owned_dofs().nth_index_in_set(i),
213 unrolled_copy_indices[i]);
217 const ::parallel::TriangulationBase<dim, spacedim> *
tria =
218 (
dynamic_cast<const ::parallel::TriangulationBase<dim, spacedim>
219 *
>(&dof_handler.get_triangulation()));
221 send_data_temp.size() == 0 ||
tria !=
nullptr,
223 "We should only be sending information with a parallel Triangulation!"));
225#ifdef DEAL_II_WITH_MPI
229 const std::set<types::subdomain_id> &neighbors =
230 tria->level_ghost_owners();
231 std::map<int, std::vector<DoFPair>> send_data;
233 std::sort(send_data_temp.begin(),
234 send_data_temp.end(),
236 if (lhs.level < rhs.level)
238 if (lhs.level > rhs.level)
241 if (lhs.level_dof_index < rhs.level_dof_index)
243 if (lhs.level_dof_index > rhs.level_dof_index)
246 if (lhs.global_dof_index < rhs.global_dof_index)
251 send_data_temp.erase(
252 std::unique(send_data_temp.begin(),
253 send_data_temp.end(),
255 return (lhs.level == rhs.level) &&
256 (lhs.level_dof_index == rhs.level_dof_index) &&
257 (lhs.global_dof_index == rhs.global_dof_index);
259 send_data_temp.end());
264 dof_handler.locally_owned_mg_dofs(
level);
266 std::vector<types::global_dof_index> level_dof_indices;
267 std::vector<types::global_dof_index> global_dof_indices;
268 for (
const auto &dofpair : send_data_temp)
269 if (dofpair.level ==
level)
271 level_dof_indices.push_back(dofpair.level_dof_index);
272 global_dof_indices.push_back(dofpair.global_dof_index);
277 level_dof_indices.end());
282 const auto index_owner =
287 AssertThrow(level_dof_indices.size() == index_owner.size(),
290 for (
unsigned int i = 0; i < index_owner.size(); ++i)
291 send_data[index_owner[i]].emplace_back(
level,
292 global_dof_indices[i],
293 level_dof_indices[i]);
306 std::vector<MPI_Request> requests;
308 for (
const auto dest : neighbors)
310 requests.push_back(MPI_Request());
311 std::vector<DoFPair> &data = send_data[dest];
317 const int ierr = MPI_Isend(data.data(),
318 data.size() *
sizeof(data[0]),
323 &*requests.rbegin());
328 const int ierr = MPI_Isend(
nullptr,
334 &*requests.rbegin());
343 std::vector<DoFPair> receive_buffer;
344 for (
unsigned int counter = 0; counter < neighbors.size();
348 int ierr = MPI_Probe(MPI_ANY_SOURCE,
354 ierr = MPI_Get_count(&status, MPI_BYTE, &len);
359 ierr = MPI_Recv(
nullptr,
370 int count = len /
sizeof(
DoFPair);
373 receive_buffer.resize(count);
375 void *ptr = receive_buffer.data();
385 for (
const auto &dof_pair : receive_buffer)
387 copy_indices_level_mine[dof_pair.level].emplace_back(
388 dof_pair.global_dof_index, dof_pair.level_dof_index);
394 if (requests.size() > 0)
396 const int ierr = MPI_Waitall(requests.size(),
398 MPI_STATUSES_IGNORE);
416 std::less<std::pair<types::global_dof_index, types::global_dof_index>>
418 for (
auto &level_indices : copy_indices_level_mine)
419 std::sort(level_indices.begin(), level_indices.end(), compare);
420 for (
auto &level_indices : copy_indices_global_mine)
421 std::sort(level_indices.begin(), level_indices.end(), compare);
431 std::vector<types::global_dof_index> &ghosted_level_dofs,
432 const std::shared_ptr<const Utilities::MPI::Partitioner>
433 & external_partitioner,
435 std::shared_ptr<const Utilities::MPI::Partitioner> &target_partitioner,
438 std::sort(ghosted_level_dofs.begin(), ghosted_level_dofs.end());
440 ghosted_dofs.
add_indices(ghosted_level_dofs.begin(),
441 std::unique(ghosted_level_dofs.begin(),
442 ghosted_level_dofs.end()));
446 if (target_partitioner.get() !=
nullptr &&
447 target_partitioner->size() == locally_owned.
size())
449 ghosted_dofs.
add_indices(target_partitioner->ghost_indices());
454 const int ghosts_locally_contained =
455 (external_partitioner.get() !=
nullptr &&
456 (external_partitioner->ghost_indices() & ghosted_dofs) ==
460 if (external_partitioner.get() !=
nullptr &&
466 if (target_partitioner.get() !=
nullptr &&
467 target_partitioner->size() == locally_owned.
size())
468 for (
unsigned int i = 0; i < copy_indices_global_mine.n_cols(); ++i)
469 copy_indices_global_mine(1, i) =
470 external_partitioner->global_to_local(
471 target_partitioner->local_to_global(
472 copy_indices_global_mine(1, i)));
473 target_partitioner = external_partitioner;
477 if (target_partitioner.get() !=
nullptr &&
478 target_partitioner->size() == locally_owned.
size())
479 for (
unsigned int i = 0; i < copy_indices_global_mine.n_cols(); ++i)
480 copy_indices_global_mine(1, i) =
483 target_partitioner->local_to_global(
484 copy_indices_global_mine(1, i)));
486 locally_owned, ghosted_dofs, communicator));
496 const std::vector<types::global_dof_index> &mine,
497 const std::vector<types::global_dof_index> &remote,
498 std::vector<unsigned int> & localized_indices)
500 localized_indices.resize(mine.size() + remote.size(),
502 for (
unsigned int i = 0; i < mine.size(); ++i)
506 for (
unsigned int i = 0; i < remote.size(); ++i)
518 const unsigned int fe_shift_1d,
519 const unsigned int fe_degree)
523 unsigned int c_tensor_index[dim];
524 unsigned int tmp = child;
525 for (
unsigned int d = 0; d < dim; ++d)
527 c_tensor_index[d] = tmp % 2;
530 const unsigned int n_child_dofs_1d = fe_degree + 1 + fe_shift_1d;
531 unsigned int factor = 1;
532 unsigned int shift = fe_shift_1d * c_tensor_index[0];
533 for (
unsigned int d = 1; d < dim; ++d)
535 factor *= n_child_dofs_1d;
536 shift = shift + factor * fe_shift_1d * c_tensor_index[d];
548 const unsigned int child,
549 const unsigned int fe_shift_1d,
550 const unsigned int fe_degree,
551 const std::vector<unsigned int> & lexicographic_numbering,
552 const std::vector<types::global_dof_index> &local_dof_indices,
555 const unsigned int n_child_dofs_1d = fe_degree + 1 + fe_shift_1d;
556 const unsigned int shift =
557 compute_shift_within_children<dim>(child, fe_shift_1d, fe_degree);
558 const unsigned int n_components =
559 local_dof_indices.size() / Utilities::fixed_power<dim>(fe_degree + 1);
561 const unsigned int n_scalar_cell_dofs =
562 Utilities::fixed_power<dim>(n_child_dofs_1d);
563 for (
unsigned int c = 0, m = 0; c < n_components; ++c)
564 for (
unsigned int k = 0; k < (dim > 2 ? (fe_degree + 1) : 1); ++k)
565 for (
unsigned int j = 0; j < (dim > 1 ? (fe_degree + 1) : 1); ++j)
566 for (
unsigned int i = 0; i < (fe_degree + 1); ++i, ++m)
568 const unsigned int index =
569 c * n_scalar_cell_dofs +
570 k * n_child_dofs_1d * n_child_dofs_1d + j * n_child_dofs_1d +
574 local_dof_indices[lexicographic_numbering[m]],
576 indices[
index] = local_dof_indices[lexicographic_numbering[m]];
582 template <
int dim,
typename Number>
586 const ::DoFHandler<dim> &dof_handler)
589 elem_info.
n_components = dof_handler.get_fe().element_multiplicity(0);
592 dof_handler.get_fe().n_dofs_per_cell());
618 const unsigned int n_child_dofs_1d =
623 elem_info.
n_components * Utilities::fixed_power<dim>(n_child_dofs_1d);
627 shape_info.
reinit(dummy_quadrature, dof_handler.get_fe(), 0);
634 for (
unsigned int c = 0; c < GeometryInfo<1>::max_children_per_cell; ++c)
638 .prolongation_matrix_1d[i * n_child_dofs_1d + j + c * shift] =
646 template <
int dim,
typename Number>
649 const ::DoFHandler<dim> &dof_handler,
651 const std::vector<std::shared_ptr<const Utilities::MPI::Partitioner>>
652 & external_partitioners,
654 std::vector<std::vector<unsigned int>> &level_dof_indices,
655 std::vector<std::vector<std::pair<unsigned int, unsigned int>>>
656 & parent_child_connect,
657 std::vector<unsigned int> &n_owned_level_cells,
658 std::vector<std::vector<std::vector<unsigned short>>> &dirichlet_indices,
659 std::vector<std::vector<Number>> & weights_on_refined,
661 MGLevelObject<std::shared_ptr<const Utilities::MPI::Partitioner>>
662 &target_partitioners)
664 level_dof_indices.clear();
665 parent_child_connect.clear();
666 n_owned_level_cells.clear();
667 dirichlet_indices.clear();
668 weights_on_refined.clear();
680 std::string fe_name = dof_handler.get_fe().base_element(0).get_name();
682 const std::size_t template_starts = fe_name.find_first_of(
'<');
683 Assert(fe_name[template_starts + 1] ==
684 (dim == 1 ?
'1' : (dim == 2 ?
'2' :
'3')),
686 fe_name[template_starts + 1] =
'1';
688 const std::unique_ptr<FiniteElement<1>> fe(
689 FETools::get_fe_by_name<1, 1>(fe_name));
697 level_dof_indices.resize(n_levels);
698 parent_child_connect.resize(n_levels - 1);
699 n_owned_level_cells.resize(n_levels - 1);
700 std::vector<std::vector<unsigned int>> coarse_level_indices(n_levels - 1);
701 for (
unsigned int level = 0;
706 std::vector<types::global_dof_index> local_dof_indices(
707 dof_handler.get_fe().n_dofs_per_cell());
708 dirichlet_indices.resize(n_levels - 1);
711 Assert(external_partitioners.empty() ||
712 external_partitioners.size() == n_levels,
717 unsigned int counter = 0;
718 std::vector<types::global_dof_index> global_level_dof_indices;
719 std::vector<types::global_dof_index> global_level_dof_indices_remote;
720 std::vector<types::global_dof_index> ghosted_level_dofs;
721 std::vector<types::global_dof_index> global_level_dof_indices_l0;
722 std::vector<types::global_dof_index> ghosted_level_dofs_l0;
725 typename ::DoFHandler<dim>::cell_iterator cell,
726 endc = dof_handler.end(
level - 1);
727 for (cell = dof_handler.begin(
level - 1); cell != endc; ++cell)
731 if (!cell->has_children())
742 const bool cell_is_remote = !consider_cell;
743 for (
unsigned int c = 0;
744 c < GeometryInfo<dim>::max_children_per_cell;
746 if (cell->child(c)->level_subdomain_id() ==
749 consider_cell =
true;
764 std::vector<types::global_dof_index> &next_indices =
765 cell_is_remote ? global_level_dof_indices_remote :
766 global_level_dof_indices;
767 const std::size_t start_index = next_indices.size();
770 for (
unsigned int c = 0;
771 c < GeometryInfo<dim>::max_children_per_cell;
774 if (cell_is_remote && cell->child(c)->level_subdomain_id() !=
777 cell->child(c)->get_mg_dof_indices(local_dof_indices);
784 dof_handler.locally_owned_mg_dofs(
level);
785 for (
const auto local_dof_index : local_dof_indices)
786 if (!owned_level_dofs.
is_element(local_dof_index))
787 ghosted_level_dofs.push_back(local_dof_index);
789 add_child_indices<dim>(c,
790 fe->n_dofs_per_cell() -
791 fe->n_dofs_per_vertex(),
795 &next_indices[start_index]);
798 if (cell->child(c)->has_children() &&
801 cell->child(c)->level_subdomain_id() ==
804 const unsigned int child_index =
805 coarse_level_indices[
level][cell->child(c)->index()];
807 parent_child_connect[
level].size());
808 unsigned int parent_index = counter;
819 parent_child_connect[
level][child_index] =
820 std::make_pair(parent_index, c);
822 static_cast<unsigned short>(-1));
826 if (mg_constrained_dofs !=
nullptr)
827 for (
unsigned int i = 0;
828 i < dof_handler.get_fe().n_dofs_per_cell();
834 dirichlet_indices[
level][child_index].push_back(i);
840 coarse_level_indices[
level - 1].size());
841 coarse_level_indices[
level - 1][cell->index()] = counter++;
848 if (
level == 1 && !cell_is_remote)
850 cell->get_mg_dof_indices(local_dof_indices);
856 const IndexSet &owned_level_dofs_l0 =
857 dof_handler.locally_owned_mg_dofs(0);
858 for (
const auto local_dof_index : local_dof_indices)
859 if (!owned_level_dofs_l0.
is_element(local_dof_index))
860 ghosted_level_dofs_l0.push_back(local_dof_index);
862 const std::size_t start_index =
863 global_level_dof_indices_l0.
size();
864 global_level_dof_indices_l0.resize(
867 add_child_indices<dim>(
869 fe->n_dofs_per_cell() - fe->n_dofs_per_vertex(),
873 &global_level_dof_indices_l0[start_index]);
875 dirichlet_indices[0].emplace_back();
876 if (mg_constrained_dofs !=
nullptr)
877 for (
unsigned int i = 0;
878 i < dof_handler.get_fe().n_dofs_per_cell();
882 local_dof_indices[elem_info
884 dirichlet_indices[0].back().push_back(i);
892 global_level_dof_indices.size());
893 n_owned_level_cells[
level - 1] = counter;
894 dirichlet_indices[
level - 1].resize(counter);
895 parent_child_connect[
level - 1].resize(
903 if (
level < n_levels - 1)
904 for (std::vector<std::pair<unsigned int, unsigned int>>::iterator
905 i = parent_child_connect[
level].begin();
906 i != parent_child_connect[
level].end();
922 external_partitioners.empty() ?
924 external_partitioners[
level],
926 target_partitioners[
level],
927 copy_indices_global_mine[
level]);
930 global_level_dof_indices,
931 global_level_dof_indices_remote,
932 level_dof_indices[
level]);
936 for (
unsigned int i = 0; i < parent_child_connect[0].size(); ++i)
937 parent_child_connect[0][i] = std::make_pair(i, 0U);
940 ghosted_level_dofs_l0,
941 external_partitioners.empty() ?
943 external_partitioners[0],
945 target_partitioners[0],
946 copy_indices_global_mine[0]);
949 *target_partitioners[0],
950 global_level_dof_indices_l0,
951 std::vector<types::global_dof_index>(),
952 level_dof_indices[0]);
958 const unsigned int n_child_dofs_1d =
959 fe->degree + 1 + fe->n_dofs_per_cell() - fe->n_dofs_per_vertex();
963 weights_on_refined.resize(n_levels - 1);
967 target_partitioners[
level]);
968 for (
unsigned int c = 0; c < n_owned_level_cells[
level - 1]; ++c)
976 std::vector<unsigned int> degree_to_3(n_child_dofs_1d);
978 for (
unsigned int i = 1; i < n_child_dofs_1d - 1; ++i)
980 degree_to_3.back() = 2;
984 weights_on_refined[
level - 1].resize(n_owned_level_cells[
level - 1] *
985 Utilities::fixed_power<dim>(3));
986 for (
unsigned int c = 0; c < n_owned_level_cells[
level - 1]; ++c)
987 for (
unsigned int k = 0, m = 0; k < (dim > 2 ? n_child_dofs_1d : 1);
989 for (
unsigned int j = 0; j < (dim > 1 ? n_child_dofs_1d : 1); ++j)
991 unsigned int shift = 9 * degree_to_3[k] + 3 * degree_to_3[j];
992 for (
unsigned int i = 0; i < n_child_dofs_1d; ++i, ++m)
993 weights_on_refined[
level -
994 1][c * Utilities::fixed_power<dim>(3) +
995 shift + degree_to_3[i]] =
998 level_dof_indices[
level]
1009 const unsigned int level,
1010 std::vector<types::global_dof_index> &dof_indices)
1012 if (mg_constrained_dofs !=
nullptr &&
1014 for (
auto &ind : dof_indices)
1034#include "mg_transfer_internal.inst"
const std::vector< std::pair< size_type, number > > * get_constraint_entries(const size_type line_n) const
bool is_identity_constrained(const size_type line_n) const
size_type n_constraints() const
unsigned int n_dofs_per_vertex() const
const unsigned int degree
unsigned int n_dofs_per_cell() const
unsigned int n_dofs_per_line() const
virtual const FullMatrix< double > & get_prolongation_matrix(const unsigned int child, const RefinementCase< dim > &refinement_case=RefinementCase< dim >::isotropic_refinement) const
size_type index_within_set(const size_type global_index) const
size_type n_elements() const
bool is_element(const size_type index) const
void add_indices(const ForwardIterator &begin, const ForwardIterator &end)
bool at_refinement_edge(const unsigned int level, const types::global_dof_index index) const
bool is_boundary_index(const unsigned int level, const types::global_dof_index index) const
const AffineConstraints< double > & get_level_constraints(const unsigned int level) const
virtual MPI_Comm get_communicator() const
virtual types::subdomain_id locally_owned_subdomain() const
unsigned int n_levels() const
unsigned int n_raw_cells(const unsigned int level) const
virtual unsigned int n_global_levels() const
unsigned int n_cells() const
Triangulation< dim, spacedim > & get_triangulation()
unsigned int global_to_local(const types::global_dof_index global_index) const
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
static ::ExceptionBase & ExcNotImplemented()
#define Assert(cond, exc)
#define AssertDimension(dim1, dim2)
#define AssertThrowMPI(error_code)
#define AssertIndexRange(index, range)
static ::ExceptionBase & ExcInternalError()
static ::ExceptionBase & ExcDimensionMismatch(std::size_t arg1, std::size_t arg2)
static ::ExceptionBase & ExcMessage(std::string arg1)
#define AssertThrow(cond, exc)
Number local_element(const size_type local_index) const
void update_ghost_values() const
virtual void compress(::VectorOperation::values operation) override
std::vector< unsigned int > compute_index_owner(const IndexSet &owned_indices, const IndexSet &indices_to_look_up, const MPI_Comm &comm)
T min(const T &t, const MPI_Comm &mpi_communicator)
T sum(const T &t, const MPI_Comm &mpi_communicator)
void reinit_level_partitioner(const IndexSet &locally_owned, std::vector< types::global_dof_index > &ghosted_level_dofs, const std::shared_ptr< const Utilities::MPI::Partitioner > &external_partitioner, const MPI_Comm &communicator, std::shared_ptr< const Utilities::MPI::Partitioner > &target_partitioner, Table< 2, unsigned int > ©_indices_global_mine)
void resolve_identity_constraints(const MGConstrainedDoFs *mg_constrained_dofs, const unsigned int level, std::vector< types::global_dof_index > &dof_indices)
unsigned int compute_shift_within_children(const unsigned int child, const unsigned int fe_shift_1d, const unsigned int fe_degree)
void setup_element_info(ElementInfo< Number > &elem_info, const FiniteElement< 1 > &fe, const ::DoFHandler< dim > &dof_handler)
void add_child_indices(const unsigned int child, const unsigned int fe_shift_1d, const unsigned int fe_degree, const std::vector< unsigned int > &lexicographic_numbering, const std::vector< types::global_dof_index > &local_dof_indices, types::global_dof_index *target_indices)
void setup_transfer(const DoFHandler< dim > &dof_handler, const MGConstrainedDoFs *mg_constrained_dofs, const std::vector< std::shared_ptr< const Utilities::MPI::Partitioner > > &external_partitioners, ElementInfo< Number > &elem_info, std::vector< std::vector< unsigned int > > &level_dof_indices, std::vector< std::vector< std::pair< unsigned int, unsigned int > > > &parent_child_connect, std::vector< unsigned int > &n_owned_level_cells, std::vector< std::vector< std::vector< unsigned short > > > &dirichlet_indices, std::vector< std::vector< Number > > &weights_on_refined, std::vector< Table< 2, unsigned int > > ©_indices_global_mine, MGLevelObject< std::shared_ptr< const Utilities::MPI::Partitioner > > &vector_partitioners)
void fill_copy_indices(const DoFHandler< dim, spacedim > &dof_handler, const MGConstrainedDoFs *mg_constrained_dofs, std::vector< std::vector< std::pair< types::global_dof_index, types::global_dof_index > > > ©_indices, std::vector< std::vector< std::pair< types::global_dof_index, types::global_dof_index > > > ©_indices_global_mine, std::vector< std::vector< std::pair< types::global_dof_index, types::global_dof_index > > > ©_indices_level_mine, const bool skip_interface_dofs=true)
void copy_indices_to_mpi_local_numbers(const Utilities::MPI::Partitioner &part, const std::vector< types::global_dof_index > &mine, const std::vector< types::global_dof_index > &remote, std::vector< unsigned int > &localized_indices)
const types::subdomain_id artificial_subdomain_id
const types::subdomain_id invalid_subdomain_id
static const unsigned int invalid_unsigned_int
const types::global_dof_index invalid_dof_index
::VectorizedArray< Number, width > min(const ::VectorizedArray< Number, width > &, const ::VectorizedArray< Number, width > &)
types::global_dof_index global_dof_index
types::global_dof_index level_dof_index
DoFPair(const unsigned int level, const types::global_dof_index global_dof_index, const types::global_dof_index level_dof_index)
std::vector< Number > prolongation_matrix_1d
unsigned int n_child_cell_dofs
std::vector< unsigned int > lexicographic_numbering
bool element_is_continuous
unsigned int n_components
void reinit(const Quadrature< dim_q > &quad, const FiniteElement< dim, spacedim > &fe_dim, const unsigned int base_element=0)
std::vector< unsigned int > lexicographic_numbering
const ::Triangulation< dim, spacedim > & tria