1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
|
// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_FUNCTIONS_FUNCTIONSPACEBASES_LAGRANGEBASIS_HH
#define DUNE_FUNCTIONS_FUNCTIONSPACEBASES_LAGRANGEBASIS_HH
#include <dune/functions/functionspacebases/flatmultiindex.hh>
#include <dune/functions/functionspacebases/defaultglobalbasis.hh>
#include <dune/functions/functionspacebases/pqknodalbasis.hh>
namespace Dune {
namespace Functions {
namespace BasisBuilder {
/**
* \brief Create a pre-basis factory that can create a PQ_k pre-basis
*
* \ingroup FunctionSpaceBasesImplementations
*
* \tparam k The polynomial order of ansatz functions
*/
template<std::size_t k>
auto lagrange()
{
return Imp::PQkPreBasisFactory<k>();
}
} // end namespace BasisBuilder
/** \brief Nodal basis of a scalar k-th-order Lagrangean finite element space
*
* \ingroup FunctionSpaceBasesImplementations
*
* \note This only works for certain grids. The following restrictions hold
* - If k is no larger than 2, then the grids can have any dimension
* - If k is larger than 3 then the grid must be two-dimensional
* - If k is 3, then the grid can be 3d *if* it is a simplex grid
*
* All arguments passed to the constructor will be forwarded to the constructor
* of PQkPreBasis.
*
* \tparam GV The GridView that the space is defined on
* \tparam k The order of the basis
*/
template<typename GV, int k>
using LagrangeBasis = DefaultGlobalBasis<PQkPreBasis<GV, k, FlatMultiIndex<std::size_t>> >;
} // end namespace Functions
} // end namespace Dune
#endif // DUNE_FUNCTIONS_FUNCTIONSPACEBASES_PQKNODALBASIS_HH
|
