4#ifndef DUNE_PDELAB_ORDERING_LEAFGRIDVIEWORDERING_HH
5#define DUNE_PDELAB_ORDERING_LEAFGRIDVIEWORDERING_HH
7#include <dune/pdelab/ordering/directleaflocalordering.hh>
8#include <dune/pdelab/ordering/leaforderingbase.hh>
17 template<
typename LocalOrdering>
22 typedef typename LocalOrdering::Traits Traits;
26 using ES =
typename Traits::EntitySet;
34 :
BaseT(local_ordering, container_blocked, gfs_data)
35 , _es(this->
template child<0>().entitySet())
59 LocalOrdering& lo = this->localOrdering();
60 lo.update_a_priori_fixed_size();
62 const std::size_t dim = ES::dimension;
64 typename ES::CodimMask codims;
66 lo.collect_used_codims(codims);
68 for (
typename ES::dim_type codim = 0; codim <= ES::dimension; ++codim)
69 if (codims.test(codim))
74 typedef typename Traits::SizeType size_type;
75 auto geom_types = _es.indexSet().types();
79 lo.update_fixed_size(geom_types.begin(),geom_types.end());
83 lo.pre_collect_used_geometry_types_from_cell();
85 for (
const auto& element : elements(_es))
87 lo.collect_used_geometry_types_from_cell(element);
90 lo.allocate_entity_offset_vector(geom_types.begin(),geom_types.end());
92 for (
const auto& element : elements(_es))
94 lo.extract_per_entity_sizes_from_cell(element);
98 lo.finalize_non_fixed_size_update();
108 for (
auto gt : geom_types)
111 size_type gt_size = lo.size(gt_index,0);
112 size_type entity_count = _es.indexSet().size(
gt);
113 _size += gt_size * entity_count;
114 if (_container_blocked)
115 gt_size = gt_size > 0;
116 _gt_dof_offsets[gt_index + 1] = gt_size * entity_count;
118 std::partial_sum(_gt_dof_offsets.begin(),_gt_dof_offsets.end(),_gt_dof_offsets.begin());
119 _block_count = _gt_dof_offsets.back();
120 _codim_fixed_size.set();
124 _block_count = _size = lo._entity_dof_offsets.back();
125 _codim_fixed_size.reset();
128 _fixed_size = lo._fixed_size;
129 _max_local_size = lo.maxLocalSize();
131 _codim_used = lo._codim_used;
132 _codim_fixed_size = lo._codim_fixed_size;
136 using BaseT::fixedSize;
140 using BaseT::_max_local_size;
142 using BaseT::_block_count;
143 using BaseT::_container_blocked;
144 using BaseT::_fixed_size;
145 using BaseT::_codim_used;
146 using BaseT::_codim_fixed_size;
147 using BaseT::_gt_dof_offsets;
149 typename Traits::EntitySet _es;
153 template<
typename GFS,
typename Transformation>
154 struct direct_leaf_gfs_to_gridview_ordering_descriptor
157 static const bool recursive =
false;
159 typedef DirectLeafLocalOrdering<
typename GFS::Traits::OrderingTag,
160 typename GFS::Traits::FiniteElementMap,
161 typename GFS::Traits::EntitySet,
162 typename Transformation::DOFIndex,
163 typename Transformation::ContainerIndex
166 typedef LeafGridViewOrdering<LocalOrdering> GridViewOrdering;
168 typedef GridViewOrdering transformed_type;
169 typedef std::shared_ptr<transformed_type> transformed_storage_type;
171 static transformed_type transform(
const GFS& gfs,
const Transformation& t)
173 return transformed_type(make_tuple(std::make_shared<LocalOrdering>(gfs.finiteElementMapStorage(),gfs.entitySet())),gfs.backend().blocked(gfs),
const_cast<GFS*
>(&gfs));
176 static transformed_storage_type transform_storage(std::shared_ptr<const GFS> gfs,
const Transformation& t)
178 return std::make_shared<transformed_type>(make_tuple(std::make_shared<LocalOrdering>(gfs->finiteElementMapStorage(),gfs->entitySet())),gfs->backend().blocked(*gfs),
const_cast<GFS*
>(gfs.get()));
184 template<
typename GFS,
typename Transformation,
typename Params>
185 direct_leaf_gfs_to_gridview_ordering_descriptor<GFS,Transformation>
186 register_leaf_gfs_to_ordering_descriptor(GFS*,Transformation*,LeafOrderingTag<Params>*);
static constexpr std::size_t index(const GeometryType >)
Compute the index for the given geometry type over all dimensions.
Definition: typeindex.hh:133
static constexpr std::size_t size(std::size_t maxdim)
Compute total number of geometry types up to and including the given dimension.
Definition: typeindex.hh:120
Gridview ordering for leaf spaces.
Definition: leafgridviewordering.hh:20
virtual void update() override
Definition: leafgridviewordering.hh:57
Generic infrastructure for orderings for leaf spaces.
Definition: leaforderingbase.hh:27
std::tuple< std::shared_ptr< Children >... > NodeStorage
The type used for storing the children.
Definition: compositenode.hh:33
bool gt(const T &first, const T &second, typename EpsilonType< T >::Type epsilon)
test if first greater than second
Definition: float_cmp.cc:156
Dune namespace.
Definition: alignedallocator.hh:14