4 #ifndef DUNE_PDELAB_ORDERING_LEAFGRIDVIEWORDERING_HH
5 #define DUNE_PDELAB_ORDERING_LEAFGRIDVIEWORDERING_HH
17 template<
typename LocalOrdering>
22 typedef typename LocalOrdering::Traits
Traits;
26 typedef typename Traits::GridView GV;
34 : BaseT(local_ordering, container_blocked, gfs_data)
35 , _gv(this->template child<0>().gridView())
58 lo.update_a_priori_fixed_size();
60 const std::size_t
dim = GV::dimension;
62 typedef typename Traits::SizeType size_type;
63 typedef std::vector<GeometryType> GTVector;
66 for (size_type cc = 0; cc <=
dim; ++cc)
68 auto per_codim_geom_types = _gv.indexSet().types(cc);
69 std::copy(per_codim_geom_types.begin(),per_codim_geom_types.end(),std::back_inserter(geom_types));
74 lo.update_fixed_size(geom_types.begin(),geom_types.end());
78 lo.pre_collect_used_geometry_types_from_cell();
80 typedef typename GV::template Codim<0>::Iterator CellIterator;
82 const CellIterator end_it = _gv.template end<0>();
83 for (CellIterator it = _gv.template begin<0>(); it != end_it; ++it)
85 lo.collect_used_geometry_types_from_cell(*it);
88 lo.allocate_entity_offset_vector(geom_types.begin(),geom_types.end());
90 for (CellIterator it = _gv.template begin<0>(); it != end_it; ++it)
92 lo.extract_per_entity_sizes_from_cell(*it);
96 lo.finalize_non_fixed_size_update();
106 const GTVector::const_iterator end_it = geom_types.end();
107 for (GTVector::const_iterator it = geom_types.begin(); it != end_it; ++it)
109 const size_type gt_index = GlobalGeometryTypeIndex::index(*it);
110 size_type gt_size = lo.size(gt_index,0);
111 size_type entity_count = _gv.indexSet().size(*it);
112 _size += gt_size * entity_count;
114 gt_size = gt_size > 0;
148 typename Traits::GridView _gv;
152 template<
typename GFS,
typename Transformation>
159 typename GFS::Traits::FiniteElementMap,
160 typename GFS::Traits::GridView,
161 typename Transformation::DOFIndex,
162 typename Transformation::ContainerIndex
170 static transformed_type
transform(
const GFS& gfs,
const Transformation& t)
172 return transformed_type(make_tuple(std::make_shared<LocalOrdering>(gfs.finiteElementMapStorage(),gfs.gridView())),gfs.backend().blocked(gfs),
const_cast<GFS*
>(&gfs));
175 static transformed_storage_type
transform_storage(std::shared_ptr<const GFS> gfs,
const Transformation& t)
177 return std::make_shared<transformed_type>(make_tuple(std::make_shared<LocalOrdering>(gfs->finiteElementMapStorage(),gfs->gridView())),gfs->backend().blocked(*gfs),
const_cast<GFS*
>(gfs.get()));
183 template<
typename GFS,
typename Transformation,
typename Params>
184 direct_leaf_gfs_to_gridview_ordering_descriptor<GFS,Transformation>
191 #endif // DUNE_PDELAB_ORDERING_LEAFORDERING_HH
Generic infrastructure for orderings for leaf spaces.
Definition: leaforderingbase.hh:21
TypeTree::CompositeNode< LocalOrdering > NodeT
Definition: leaforderingbase.hh:37
std::size_t _size
Definition: orderingbase.hh:289
bool fixedSize(typename Traits::SizeType codim) const
Definition: orderingbase.hh:215
LocalOrdering & localOrdering()
Definition: leaforderingbase.hh:44
std::size_t _max_local_size
Definition: orderingbase.hh:288
static const bool recursive
Definition: leafgridviewordering.hh:156
GridViewOrdering transformed_type
Definition: leafgridviewordering.hh:167
std::size_t _block_count
Definition: orderingbase.hh:290
direct_leaf_gfs_to_gridview_ordering_descriptor< GFS, Transformation > register_leaf_gfs_to_ordering_descriptor(GFS *, Transformation *, LeafOrderingTag< Params > *)
const bool _container_blocked
Definition: orderingbase.hh:276
Definition: directleaflocalordering.hh:27
std::vector< typename Traits::SizeType > _gt_dof_offsets
Definition: leaforderingbase.hh:282
Gridview ordering for leaf spaces.
Definition: leafgridviewordering.hh:18
bool _fixed_size
Definition: orderingbase.hh:275
LeafGridViewOrdering< LocalOrdering > GridViewOrdering
Definition: leafgridviewordering.hh:165
virtual void update()
Definition: leafgridviewordering.hh:55
Definition: leafgridviewordering.hh:153
Traits::CodimFlag _codim_used
Definition: orderingbase.hh:285
DirectLeafLocalOrdering< typename GFS::Traits::OrderingTag, typename GFS::Traits::FiniteElementMap, typename GFS::Traits::GridView, typename Transformation::DOFIndex, typename Transformation::ContainerIndex > LocalOrdering
Definition: leafgridviewordering.hh:163
Definition: adaptivity.hh:27
LocalOrdering::Traits Traits
Definition: leafgridviewordering.hh:22
Dune::PDELab::impl::GridFunctionSpaceOrderingData< typename Traits::SizeType > GFSData
Definition: orderingbase.hh:34
static transformed_storage_type transform_storage(std::shared_ptr< const GFS > gfs, const Transformation &t)
Definition: leafgridviewordering.hh:175
static const int dim
Definition: adaptivity.hh:83
static transformed_type transform(const GFS &gfs, const Transformation &t)
Definition: leafgridviewordering.hh:170
std::shared_ptr< transformed_type > transformed_storage_type
Definition: leafgridviewordering.hh:168
LeafGridViewOrdering(const typename NodeT::NodeStorage &local_ordering, bool container_blocked, typename BaseT::GFSData *gfs_data)
Definition: leafgridviewordering.hh:33
Traits::CodimFlag _codim_fixed_size
Definition: orderingbase.hh:286