3 #ifndef DUNE_FUNCTIONS_FUNCTIONSPACEBASES_HIERARCHICALLAGRANGEBASIS_HH
4 #define DUNE_FUNCTIONS_FUNCTIONSPACEBASES_HIERARCHICALLAGRANGEBASIS_HH
6 #include <dune/common/exceptions.hh>
7 #include <dune/localfunctions/hierarchical/hierarchicalp2.hh>
12 #include <dune/grid/common/mcmgmapper.hh>
33 template<
typename GV,
int k,
typename R=
double>
34 class HierarchicalLagrangeNode;
36 template<
typename GV,
int k,
class MI,
typename R=
double>
37 class HierarchicalLagrangePreBasis;
49 template<
typename GV,
int k,
class MI,
typename R>
52 static const int dim = GV::dimension;
66 using IndexSet = Impl::DefaultNodeIndexSet<HierarchicalLagrangePreBasis>;
126 assert(prefix.size() == 0 || prefix.size() == 1);
127 return (prefix.size() == 0) ?
size() : 0;
143 return Dune::binomial(std::size_t(
order() + (
unsigned int)GV::dimension),std::size_t(
order()));
146 template<
typename It>
151 Dune::LocalKey localKey = node.
finiteElement().localCoefficients().localKey(i);
152 const auto& element = node.
element();
174 static auto p2Layout()
176 return [](Dune::GeometryType type,
int gridDim)
182 if (type.isTriangle())
184 assert(type.isTetrahedron());
192 template<
typename GV,
int k,
typename R>
196 static const int dim = GV::dimension;
201 using Element =
typename GV::template Codim<0>::Entity;
202 using FiniteElement = HierarchicalP2LocalFiniteElement<typename GV::ctype,R,dim>;
230 DUNE_THROW(Dune::Exception,
231 "HierarchicalLagrange-elements do not exist for elements of type " << e.type());
249 namespace BasisFactory {
253 template<
int k,
typename R=
double>
254 class HierarchicalLagrangePreBasisFactory
258 static const std::size_t requiredMultiIndexSize = 1;
260 template<
class MultiIndex,
class Gr
idView>
261 auto makePreBasis(
const GridView& gridView)
const
278 template<std::
size_t k,
typename R=
double>
281 return Impl::HierarchicalLagrangePreBasisFactory<k,R>();
296 template<
typename GV,
int k,
typename R=
double>
auto hierarchicalLagrange()
Create a pre-basis factory that can create a HierarchicalLagrange pre-basis.
Definition: hierarchicallagrangebasis.hh:279
Definition: polynomial.hh:10
Global basis for given pre-basis.
Definition: defaultglobalbasis.hh:47
Definition: hierarchicallagrangebasis.hh:195
HierarchicalP2LocalFiniteElement< typename GV::ctype, R, dim > FiniteElement
Definition: hierarchicallagrangebasis.hh:202
const Element & element() const
Return current element, throw if unbound.
Definition: hierarchicallagrangebasis.hh:210
const FiniteElement & finiteElement() const
Return the LocalFiniteElement for the element we are bound to.
Definition: hierarchicallagrangebasis.hh:219
void bind(const Element &e)
Bind to element.
Definition: hierarchicallagrangebasis.hh:225
const Element * element_
Definition: hierarchicallagrangebasis.hh:244
const FiniteElement finiteElement_
Definition: hierarchicallagrangebasis.hh:243
unsigned int order() const
Definition: hierarchicallagrangebasis.hh:238
HierarchicalLagrangeNode()
Definition: hierarchicallagrangebasis.hh:204
std::size_t size_type
Definition: hierarchicallagrangebasis.hh:200
typename GV::template Codim< 0 >::Entity Element
Definition: hierarchicallagrangebasis.hh:201
A pre-basis for a hierarchical basis.
Definition: hierarchicallagrangebasis.hh:51
size_type size() const
Same as size(prefix) with empty prefix.
Definition: hierarchicallagrangebasis.hh:118
MI MultiIndex
Type used for global numbering of the basis vectors.
Definition: hierarchicallagrangebasis.hh:69
Impl::DefaultNodeIndexSet< HierarchicalLagrangePreBasis > IndexSet
Type of created tree node index set.
Definition: hierarchicallagrangebasis.hh:66
size_type maxNodeSize() const
Get the maximal number of DOFs associated to node for any element.
Definition: hierarchicallagrangebasis.hh:140
size_type size(const SizePrefix prefix) const
Return number of possible values for next position in multi index.
Definition: hierarchicallagrangebasis.hh:124
It indices(const Node &node, It it) const
Definition: hierarchicallagrangebasis.hh:147
Node makeNode() const
Create tree node.
Definition: hierarchicallagrangebasis.hh:101
unsigned int order() const
Definition: hierarchicallagrangebasis.hh:162
size_type dimension() const
Get the total dimension of the space spanned by this basis.
Definition: hierarchicallagrangebasis.hh:131
HierarchicalLagrangePreBasis(const GridView &gv)
Constructor for a given grid view object with layout for second order.
Definition: hierarchicallagrangebasis.hh:78
Dune::ReservedVector< size_type, 1 > SizePrefix
Type used for prefixes handed to the size() method.
Definition: hierarchicallagrangebasis.hh:72
GV GridView
The grid view that the FE basis is defined on.
Definition: hierarchicallagrangebasis.hh:57
IndexSet makeIndexSet() const
Create tree node index set.
Definition: hierarchicallagrangebasis.hh:112
void initializeIndices()
Initialize the global indices.
Definition: hierarchicallagrangebasis.hh:82
const GridView & gridView() const
Obtain the grid view that the basis is defined on.
Definition: hierarchicallagrangebasis.hh:86
std::size_t size_type
Type used for indices and size information.
Definition: hierarchicallagrangebasis.hh:60
void update(const GridView &gv)
Update the stored grid view & MultipleCodimMultipleGeomTypeMapper, to be called if the grid has chang...
Definition: hierarchicallagrangebasis.hh:92
GridView gridView_
Definition: hierarchicallagrangebasis.hh:160
MultipleCodimMultipleGeomTypeMapper< GridView > mcmgMapper_
Definition: hierarchicallagrangebasis.hh:167
void setSize(const size_type size)
Definition: nodes.hh:162