dune-localfunctions  2.9.0
meta/power/basis.hh
Go to the documentation of this file.
1 // -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
2 // vi: set et ts=4 sw=2 sts=2:
3 // SPDX-FileCopyrightInfo: Copyright (C) DUNE Project contributors, see file LICENSE.md in module root
4 // SPDX-License-Identifier: LicenseRef-GPL-2.0-only-with-DUNE-exception
5 
6 #ifndef DUNE_LOCALFUNCTIONS_META_POWER_BASIS_HH
7 #define DUNE_LOCALFUNCTIONS_META_POWER_BASIS_HH
8 
9 #include <numeric>
10 #include <cstddef>
11 #include <vector>
12 
13 #include <dune/common/fmatrix.hh>
14 #include <dune/common/fvector.hh>
15 
16 namespace Dune {
17 
19 
25  template<class Backend, std::size_t dimR>
26  class PowerBasis {
27  static_assert(Backend::Traits::dimRange == 1,
28  "PowerBasis works only with scalar backends");
29 
30  // don't use a reference here so this class stays copyable
31  const Backend *backend;
32 
33  public:
35  struct Traits : public Backend::Traits
36  {
38  static const std::size_t dimRange = dimR;
40  typedef FieldVector<typename Traits::RangeField, dimR> Range;
41 
43 
47  typedef FieldMatrix<typename Traits::RangeField, dimR,
48  Traits::dimDomainGlobal> Jacobian;
49  };
50 
52 
58  PowerBasis(const Backend &backend_) : backend(&backend_) { }
59 
61  std::size_t size () const { return backend->size()*dimR; }
63  std::size_t order () const { return backend->order(); }
64 
66  void evaluateFunction(const typename Traits::DomainLocal& in,
67  std::vector<typename Traits::Range>& out) const
68  {
69  std::vector<typename Backend::Traits::Range> backendValues;
70  backend->evaluateFunction(in, backendValues);
71  out.assign(size(), typename Traits::Range(0));
72  for(std::size_t d = 0; d < dimR; ++d)
73  for(std::size_t i = 0; i < backend->size(); ++i)
74  out[d*backend->size()+i][d] = backendValues[i][0];
75  }
76 
78  void evaluateJacobian(const typename Traits::DomainLocal& in,
79  std::vector<typename Traits::Jacobian>& out) const
80  {
81  std::vector<typename Backend::Traits::Jacobian> backendValues;
82  backend->evaluateJacobian(in, backendValues);
83  out.assign(size(), typename Traits::Jacobian(0));
84  for(std::size_t d = 0; d < dimR; ++d)
85  for(std::size_t i = 0; i < backend->size(); ++i)
86  out[d*backend->size()+i][d] = backendValues[i][0];
87  }
88 
90  void partial (const std::array<unsigned int, Backend::Traits::dimDomainGlobal>& order,
91  const typename Traits::DomainLocal& in, // position
92  std::vector<typename Traits::Range>& out) const // return value
93  {
94  auto totalOrder = std::accumulate(order.begin(), order.end(), 0);
95  if (totalOrder == 0) {
96  evaluateFunction(in, out);
97  } else {
98  DUNE_THROW(NotImplemented, "Desired derivative order is not implemented");
99  }
100  }
101  };
102 
103 } // namespace Dune
104 
105 #endif // DUNE_LOCALFUNCTIONS_META_POWER_BASIS_HH
Definition: bdfmcube.hh:18
Meta-basis turning a scalar basis into vector-valued basis.
Definition: meta/power/basis.hh:26
PowerBasis(const Backend &backend_)
Construct a PowerBasis.
Definition: meta/power/basis.hh:58
std::size_t size() const
Number of shape functions.
Definition: meta/power/basis.hh:61
void partial(const std::array< unsigned int, Backend::Traits::dimDomainGlobal > &order, const typename Traits::DomainLocal &in, std::vector< typename Traits::Range > &out) const
Evaluate partial derivatives of all shape functions.
Definition: meta/power/basis.hh:90
void evaluateJacobian(const typename Traits::DomainLocal &in, std::vector< typename Traits::Jacobian > &out) const
Evaluate Jacobian of all shape functions at given position.
Definition: meta/power/basis.hh:78
void evaluateFunction(const typename Traits::DomainLocal &in, std::vector< typename Traits::Range > &out) const
Evaluate all shape functions at given position.
Definition: meta/power/basis.hh:66
std::size_t order() const
Polynomial order of the shape functions for quadrature.
Definition: meta/power/basis.hh:63
types of domain and range
Definition: meta/power/basis.hh:36
FieldMatrix< typename Traits::RangeField, dimR, Traits::dimDomainGlobal > Jacobian
Jacobian properties.
Definition: meta/power/basis.hh:48
FieldVector< typename Traits::RangeField, dimR > Range
Type used for range values.
Definition: meta/power/basis.hh:40
static const std::size_t dimRange
Dimension of the range values.
Definition: meta/power/basis.hh:38