\n
\n
\n
\n
\n
\n-
7#ifndef DUNE_FUNCTIONS_COMMON_DEFAULT_DERIVATIVE_TRAITS_HH \n-
8#define DUNE_FUNCTIONS_COMMON_DEFAULT_DERIVATIVE_TRAITS_HH \n+
7#ifndef DUNE_FUNCTIONS_FUNCTIONSPACEBASES_REFINEDLAGRANGEBASIS_HH \n+
8#define DUNE_FUNCTIONS_FUNCTIONSPACEBASES_REFINEDLAGRANGEBASIS_HH \n
\n
\n-
\n-
\n-
13#include <dune/common/fvector.hh> \n-
14#include <dune/common/fmatrix.hh> \n-
\n-
\n-
\n-
\n-
\n+
\n+
12#include <dune/common/exceptions.hh> \n+
13#include <dune/common/math.hh> \n+
\n+
15#include <dune/localfunctions/refined.hh> \n+
\n+
\n+
\n+
\n
\n-
\n+
21#include <dune/geometry/type.hh> \n+
\n+
23#include <dune/grid/common/mcmgmapper.hh> \n+
\n+
\n+
\n+
\n
\n-
\n-
39template<
class Signature>
\n-
\n-
\n-
\n-
\n-
\n-
\n-
70template<
typename K,
int n>
\n-
\n+
29template<
typename GV,
int k,
typename R>
\n+
30class RefinedLagrangeNode;
\n+
\n+
43template <
typename GV,
int k,
typename R =
double>
\n+
\n+
\n+
\n+
\n+
\n+
\n+
49 static const int dim = GV::dimension;
\n+
\n+
\n+
52 static_assert(k == 0 || k == 1);
\n+
\n+
\n+
55 static MCMGLayout dofLayout()
\n+
\n+
\n+
\n+
59 return [](GeometryType gt,
int) ->
size_t {
\n+
60 return (gt.dim() == dim) ? (1 << dim) : 0;
\n+
\n+
62 else if constexpr(k == 1)
\n+
\n+
64 return [](GeometryType gt,
int) ->
size_t {
\n+
65 return Dune::binomial(
int(k),
int(gt.dim()));
\n+
\n+
\n+
68 DUNE_THROW(Dune::NotImplemented,
\n+
69 "Refined basis not implemented for higher-order Lagrange (k>=2) elements.");
\n+
\n+
\n+
\n+
\n+
\n
\n-
86template<
typename K,
int n,
int m>
\n-
\n-
\n-
\n-
90 typedef FieldMatrix<K,m,n>
Range;
\n-
\n-
\n-
\n-
102template<
typename K,
int n,
int m>
\n-
\n-
\n-
\n-
\n-
\n-
\n-
\n+
\n+
\n+
\n+
\n+
87 :
Base(gv, dofLayout())
\n+
\n+
89 for (
auto gt : gv.indexSet().types(0)) {
\n+
\n+
91 DUNE_THROW(Dune::NotImplemented,
\n+
92 "Refined Lagrange basis only implemented for simplex grids.");
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
109 static constexpr unsigned int order()
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
117template <
typename GV,
int k,
typename R>
\n+
\n+
\n+
\n+
\n+
121 static constexpr int dim = GV::dimension;
\n+
\n+
\n+
124 static_assert(k == 0 || k == 1);
\n+
\n+
\n+
128 using Element =
typename GV::template Codim<0>::Entity;
\n+
\n+
\n+
132 Dune::RefinedP0LocalFiniteElement<typename GV::ctype,R,dim>,
\n+
133 Dune::RefinedP1LocalFiniteElement<typename GV::ctype,R,dim>>;
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
180 static constexpr unsigned int order()
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
192namespace BasisFactory {
\n+
\n+
201template <
int k,
typename R=
double>
\n+
\n+
\n+
\n+
204 return [](
const auto& gridView) {
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
225template <
typename GV,
int k,
typename R=
double>
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
auto refinedLagrange()
Create a pre-basis factory that can create a RefinedLagrange pre-basis.
Definition refinedlagrangebasis.hh:202
\n
Definition polynomial.hh:17
\n-
Dummy range class to be used if no proper type is available.
Definition defaultderivativetraits.hh:27
\n-
Default implementation for derivative traits.
Definition defaultderivativetraits.hh:41
\n-
InvalidRange Range
Range of derivative for function with given signature.
Definition defaultderivativetraits.hh:43
\n-
double Range
Range of derivative for function with given signature.
Definition defaultderivativetraits.hh:58
\n-
FieldVector< K, n > Range
Range of derivative for function with given signature.
Definition defaultderivativetraits.hh:74
\n-
FieldMatrix< K, m, n > Range
Range of derivative for function with given signature.
Definition defaultderivativetraits.hh:90
\n-
FieldMatrix< K, m, n > Range
Range of derivative for function with given signature.
Definition defaultderivativetraits.hh:106
\n+
Global basis for given pre-basis.
Definition defaultglobalbasis.hh:50
\n+
A generic MixIn class for PreBasis with flat indices computed from a mapper.
Definition leafprebasismappermixin.hh:62
\n+
void setSize(const size_type size)
Definition nodes.hh:169
\n+
\n+
Definition refinedlagrangebasis.hh:120
\n+
const FiniteElement & finiteElement() const
Return the LocalFiniteElement for the element we are bound to.
Definition refinedlagrangebasis.hh:164
\n+
const Element * element_
Definition refinedlagrangebasis.hh:187
\n+
void bind(const Element &e)
Bind the node to the element e.
Definition refinedlagrangebasis.hh:170
\n+
const FiniteElement finiteElement_
Definition refinedlagrangebasis.hh:186
\n+
static constexpr unsigned int order()
Polynomial order used in the local Lagrange finite-elements in subdomains of the element.
Definition refinedlagrangebasis.hh:180
\n+
typename GV::template Codim< 0 >::Entity Element
Type of the element in the GridView.
Definition refinedlagrangebasis.hh:128
\n+
std::conditional_t<(k==0), Dune::RefinedP0LocalFiniteElement< typename GV::ctype, R, dim >, Dune::RefinedP1LocalFiniteElement< typename GV::ctype, R, dim > > FiniteElement
Type of the local finite-element.
Definition refinedlagrangebasis.hh:133
\n+
RefinedLagrangeNode()
The default constructor initializes all members to their default.
Definition refinedlagrangebasis.hh:144
\n+
const Element & element() const
Return current element. The behavior is undefined if the node is not bound to any element.
Definition refinedlagrangebasis.hh:153
\n+
A pre-basis for a refined Lagrange bases.
Definition refinedlagrangebasis.hh:46
\n+
static constexpr unsigned int order()
Polynomial order used in the local Lagrange finite-elements.
Definition refinedlagrangebasis.hh:109
\n+
Node makeNode() const
Create tree node.
Definition refinedlagrangebasis.hh:97
\n+
GV GridView
The grid view that the FE basis is defined on.
Definition refinedlagrangebasis.hh:75
\n+
RefinedLagrangePreBasis(const GridView &gv)
Constructor for a given grid view object.
Definition refinedlagrangebasis.hh:86
\n
![\"doxygen\"/](\"doxygen.svg\")
1.9.8\n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -1,30 +1,42 @@\n dune-functions\u00a02.10\n Loading...\n Searching...\n No Matches\n * _\bd_\bu_\bn_\be\n * _\bf_\bu_\bn_\bc_\bt_\bi_\bo_\bn_\bs\n- * _\ba_\bn_\ba_\bl_\by_\bt_\bi_\bc_\bf_\bu_\bn_\bc_\bt_\bi_\bo_\bn_\bs\n+ * _\bf_\bu_\bn_\bc_\bt_\bi_\bo_\bn_\bs_\bp_\ba_\bc_\be_\bb_\ba_\bs_\be_\bs\n _\bC_\bl_\ba_\bs_\bs_\be_\bs | _\bN_\ba_\bm_\be_\bs_\bp_\ba_\bc_\be_\bs | _\bF_\bu_\bn_\bc_\bt_\bi_\bo_\bn_\bs\n-trigonometricfunction.hh File Reference\n-#include