dune-pdelab  2.5-dev
leafgridviewordering.hh
Go to the documentation of this file.
1 // -*- tab-width: 8; indent-tabs-mode: nil; c-basic-offset: 2 -*-
2 // vi: set et ts=8 sw=2 sts=2:
3 
4 #ifndef DUNE_PDELAB_ORDERING_LEAFGRIDVIEWORDERING_HH
5 #define DUNE_PDELAB_ORDERING_LEAFGRIDVIEWORDERING_HH
6 
9 
10 namespace Dune {
11  namespace PDELab {
12 
15 
17  template<typename LocalOrdering>
19  : public LeafOrderingBase<LocalOrdering>
20  {
21  public:
22  typedef typename LocalOrdering::Traits Traits;
23 
24  private:
25 
26  using ES = typename Traits::EntitySet;
27 
29  typedef typename BaseT::NodeT NodeT;
30 
31  public:
32 
33  LeafGridViewOrdering(const typename NodeT::NodeStorage& local_ordering, bool container_blocked, typename BaseT::GFSData* gfs_data)
34  : BaseT(local_ordering, container_blocked, gfs_data)
35  , _es(this->template child<0>().entitySet())
36  {}
37 
38 #ifndef DOXYGEN
39 
40 // we need to override the default copy / move ctor to fix the delegate pointer, but that is
41 // hardly interesting to our users...
42 
44  : BaseT(r)
45  , _es(r._es)
46  {}
47 
49  : BaseT(std::move(r))
50  , _es(r._es)
51  {}
52 
53 #endif // DOXYGEN
54 
55  virtual void update() override
56  {
57  LocalOrdering& lo = this->localOrdering();
58  lo.update_a_priori_fixed_size();
59 
60  const std::size_t dim = ES::dimension;
61 
62  typename ES::CodimMask codims;
63  codims.set(0); // always need cells
64  lo.collect_used_codims(codims);
65 
66  for (typename ES::dim_type codim = 0; codim <= ES::dimension; ++codim)
67  if (codims.test(codim))
68  _es.addCodim(codim);
69 
70  _es.update();
71 
72  typedef typename Traits::SizeType size_type;
73  auto geom_types = _es.indexSet().types();
74 
75  if (lo._fixed_size)
76  {
77  lo.update_fixed_size(geom_types.begin(),geom_types.end());
78  }
79  else
80  {
81  lo.pre_collect_used_geometry_types_from_cell();
82 
83  for (const auto& element : elements(_es))
84  {
85  lo.collect_used_geometry_types_from_cell(element);
86  }
87 
88  lo.allocate_entity_offset_vector(geom_types.begin(),geom_types.end());
89 
90  for (const auto& element : elements(_es))
91  {
92  lo.extract_per_entity_sizes_from_cell(element);
93  }
94 
95  // FIXME: handling of blocked containers!
96  lo.finalize_non_fixed_size_update();
97  }
98 
99  // we need to re-test here, as the local ordering could have detected a fixed size ordering
100  // and switched its implementation
101  if (lo._fixed_size)
102  {
103  _gt_dof_offsets.assign(GlobalGeometryTypeIndex::size(dim) + 1,0);
104  _size = 0;
105 
106  for (auto gt : geom_types)
107  {
108  const size_type gt_index = GlobalGeometryTypeIndex::index(gt);
109  size_type gt_size = lo.size(gt_index,0);
110  size_type entity_count = _es.indexSet().size(gt);
111  _size += gt_size * entity_count;
112  if (_container_blocked)
113  gt_size = gt_size > 0;
114  _gt_dof_offsets[gt_index + 1] = gt_size * entity_count;
115  }
116  std::partial_sum(_gt_dof_offsets.begin(),_gt_dof_offsets.end(),_gt_dof_offsets.begin());
117  _block_count = _gt_dof_offsets.back();
118  _codim_fixed_size.set();
119  }
120  else
121  {
122  _block_count = _size = lo._entity_dof_offsets.back();
123  _codim_fixed_size.reset();
124  }
125 
126  _fixed_size = lo._fixed_size;
127  _max_local_size = lo.maxLocalSize();
128 
129  _codim_used = lo._codim_used;
130  _codim_fixed_size = lo._codim_fixed_size;
131 
132  }
133 
134  using BaseT::fixedSize;
135 
136  private:
137 
139  using BaseT::_size;
140  using BaseT::_block_count;
142  using BaseT::_fixed_size;
143  using BaseT::_codim_used;
146 
147  typename Traits::EntitySet _es;
148  };
149 
150 
151  template<typename GFS, typename Transformation>
153  {
154 
155  static const bool recursive = false;
156 
157  typedef DirectLeafLocalOrdering<typename GFS::Traits::OrderingTag,
158  typename GFS::Traits::FiniteElementMap,
159  typename GFS::Traits::EntitySet,
160  typename Transformation::DOFIndex,
161  typename Transformation::ContainerIndex
163 
165 
166  typedef GridViewOrdering transformed_type;
167  typedef std::shared_ptr<transformed_type> transformed_storage_type;
168 
169  static transformed_type transform(const GFS& gfs, const Transformation& t)
170  {
171  return transformed_type(make_tuple(std::make_shared<LocalOrdering>(gfs.finiteElementMapStorage(),gfs.entitySet())),gfs.backend().blocked(gfs),const_cast<GFS*>(&gfs));
172  }
173 
174  static transformed_storage_type transform_storage(std::shared_ptr<const GFS> gfs, const Transformation& t)
175  {
176  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()));
177  }
178 
179  };
180 
181 
182  template<typename GFS, typename Transformation, typename Params>
185 
187  } // namespace PDELab
188 } // namespace Dune
189 
190 #endif // DUNE_PDELAB_ORDERING_LEAFGRIDVIEWORDERING_HH
static const int dim
Definition: adaptivity.hh:84
virtual void update() override
Definition: leafgridviewordering.hh:55
static transformed_type transform(const GFS &gfs, const Transformation &t)
Definition: leafgridviewordering.hh:169
TypeTree::CompositeNode< LocalOrdering > NodeT
Definition: leaforderingbase.hh:37
bool fixedSize(typename Traits::SizeType codim) const
Definition: orderingbase.hh:209
DirectLeafLocalOrdering< typename GFS::Traits::OrderingTag, typename GFS::Traits::FiniteElementMap, typename GFS::Traits::EntitySet, typename Transformation::DOFIndex, typename Transformation::ContainerIndex > LocalOrdering
Definition: leafgridviewordering.hh:162
direct_leaf_gfs_to_gridview_ordering_descriptor< GFS, Transformation > register_leaf_gfs_to_ordering_descriptor(GFS *, Transformation *, LeafOrderingTag< Params > *)
Tag indicating a standard ordering for a leaf GridfunctionSpace.
Definition: gridfunctionspace/tags.hh:183
LocalOrdering::Traits Traits
Definition: leafgridviewordering.hh:22
std::vector< typename Traits::SizeType > _gt_dof_offsets
Definition: leaforderingbase.hh:279
For backward compatibility – Do not use this!
Definition: adaptivity.hh:28
GridViewOrdering transformed_type
Definition: leafgridviewordering.hh:166
LeafGridViewOrdering(const typename NodeT::NodeStorage &local_ordering, bool container_blocked, typename BaseT::GFSData *gfs_data)
Definition: leafgridviewordering.hh:33
Generic infrastructure for orderings for leaf spaces.
Definition: leaforderingbase.hh:21
Gridview ordering for leaf spaces.
Definition: leafgridviewordering.hh:18
Definition: directleaflocalordering.hh:25
LocalOrdering & localOrdering()
Definition: leaforderingbase.hh:44
LeafGridViewOrdering< LocalOrdering > GridViewOrdering
Definition: leafgridviewordering.hh:164
Dune::PDELab::impl::GridFunctionSpaceOrderingData< typename Traits::SizeType > GFSData
Definition: orderingbase.hh:32
std::shared_ptr< transformed_type > transformed_storage_type
Definition: leafgridviewordering.hh:167
static transformed_storage_type transform_storage(std::shared_ptr< const GFS > gfs, const Transformation &t)
Definition: leafgridviewordering.hh:174
std::size_t index
Definition: interpolate.hh:118