dune-localfunctions  2.6-git
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 
4 #ifndef DUNE_LOCALFUNCTIONS_META_POWER_BASIS_HH
5 #define DUNE_LOCALFUNCTIONS_META_POWER_BASIS_HH
6 
7 #include <numeric>
8 #include <cstddef>
9 #include <vector>
10 
11 #include <dune/common/fmatrix.hh>
12 #include <dune/common/fvector.hh>
13 
14 namespace Dune {
15 
17 
23  template<class Backend, std::size_t dimR>
24  class PowerBasis {
25  static_assert(Backend::Traits::dimRange == 1,
26  "PowerBasis works only with scalar backends");
27 
28  // don't use a reference here so this class stays copyable
29  const Backend *backend;
30 
31  public:
33  struct Traits : public Backend::Traits
34  {
36  static const std::size_t dimRange = dimR;
38  typedef FieldVector<typename Traits::RangeField, dimR> Range;
39 
41 
45  typedef FieldMatrix<typename Traits::RangeField, dimR,
46  Traits::dimDomainGlobal> Jacobian;
47  };
48 
50 
56  PowerBasis(const Backend &backend_) : backend(&backend_) { }
57 
59  std::size_t size () const { return backend->size()*dimR; }
61  std::size_t order () const { return backend->order(); }
62 
64  void evaluateFunction(const typename Traits::DomainLocal& in,
65  std::vector<typename Traits::Range>& out) const
66  {
67  std::vector<typename Backend::Traits::Range> backendValues;
68  backend->evaluateFunction(in, backendValues);
69  out.assign(size(), typename Traits::Range(0));
70  for(std::size_t d = 0; d < dimR; ++d)
71  for(std::size_t i = 0; i < backend->size(); ++i)
72  out[d*backend->size()+i][d] = backendValues[i][0];
73  }
74 
76  void evaluateJacobian(const typename Traits::DomainLocal& in,
77  std::vector<typename Traits::Jacobian>& out) const
78  {
79  std::vector<typename Backend::Traits::Jacobian> backendValues;
80  backend->evaluateJacobian(in, backendValues);
81  out.assign(size(), typename Traits::Jacobian(0));
82  for(std::size_t d = 0; d < dimR; ++d)
83  for(std::size_t i = 0; i < backend->size(); ++i)
84  out[d*backend->size()+i][d] = backendValues[i][0];
85  }
86 
88  void partial (const std::array<unsigned int, Backend::Traits::dimDomainGlobal>& order,
89  const typename Traits::DomainLocal& in, // position
90  std::vector<typename Traits::Range>& out) const // return value
91  {
92  auto totalOrder = std::accumulate(order.begin(), order.end(), 0);
93  if (totalOrder == 0) {
94  evaluateFunction(in, out);
95  } else {
96  DUNE_THROW(NotImplemented, "Desired derivative order is not implemented");
97  }
98  }
99  };
100 
101 } // namespace Dune
102 
103 #endif // DUNE_LOCALFUNCTIONS_META_POWER_BASIS_HH
Dune::PowerBasis::Traits::Range
FieldVector< typename Traits::RangeField, dimR > Range
Type used for range values.
Definition: meta/power/basis.hh:38
Dune::PowerBasis::evaluateFunction
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:64
Dune::PowerBasis::Traits
types of domain and range
Definition: meta/power/basis.hh:33
Dune::PowerBasis::PowerBasis
PowerBasis(const Backend &backend_)
Construct a PowerBasis.
Definition: meta/power/basis.hh:56
Dune::PowerBasis::order
std::size_t order() const
Polynomial order of the shape functions for quadrature.
Definition: meta/power/basis.hh:61
Dune::PowerBasis
Meta-basis turning a scalar basis into vector-valued basis.
Definition: meta/power/basis.hh:24
Dune::PowerBasis::Traits::Jacobian
FieldMatrix< typename Traits::RangeField, dimR, Traits::dimDomainGlobal > Jacobian
Jacobian properties.
Definition: meta/power/basis.hh:46
Dune::PowerBasis::partial
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:88
Dune::PowerBasis::evaluateJacobian
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:76
Dune::PowerBasis::Traits::dimRange
static const std::size_t dimRange
Dimension of the range values.
Definition: meta/power/basis.hh:36
Dune
Definition: brezzidouglasmarini1cube2dlocalbasis.hh:15
Dune::PowerBasis::size
std::size_t size() const
Number of shape functions.
Definition: meta/power/basis.hh:59