4#ifndef DUNE_LOCALFUNCTIONS_META_POWER_BASIS_HH
5#define DUNE_LOCALFUNCTIONS_META_POWER_BASIS_HH
11#include <dune/common/fmatrix.hh>
12#include <dune/common/fvector.hh>
23 template<
class Backend, std::
size_t dimR>
25 static_assert(Backend::Traits::dimRange == 1,
26 "PowerBasis works only with scalar backends");
29 const Backend *backend;
33 struct Traits :
public Backend::Traits
38 typedef FieldVector<typename Traits::RangeField, dimR>
Range;
45 typedef FieldMatrix<
typename Traits::RangeField, dimR,
50 Backend::Traits::diffOrder > 0 ? 1 : 0;
60 PowerBasis(
const Backend &backend_) : backend(&backend_) { }
63 std::size_t
size ()
const {
return backend->size()*dimR; }
65 std::size_t
order ()
const {
return backend->order(); }
69 std::vector<typename Traits::Range>& out)
const
71 std::vector<typename Backend::Traits::Range> backendValues;
72 backend->evaluateFunction(in, backendValues);
74 for(std::size_t d = 0; d < dimR; ++d)
75 for(std::size_t i = 0; i < backend->size(); ++i)
76 out[d*backend->size()+i][d] = backendValues[i][0];
81 std::vector<typename Traits::Jacobian>& out)
const
83 std::vector<typename Backend::Traits::Jacobian> backendValues;
84 backend->evaluateJacobian(in, backendValues);
86 for(std::size_t d = 0; d < dimR; ++d)
87 for(std::size_t i = 0; i < backend->size(); ++i)
88 out[d*backend->size()+i][d] = backendValues[i][0];
92 void partial (
const std::array<unsigned int, Backend::Traits::dimDomainGlobal>&
order,
93 const typename Traits::DomainLocal& in,
94 std::vector<typename Traits::Range>& out)
const
96 auto totalOrder = std::accumulate(
order.begin(),
order.end(), 0);
97 if (totalOrder == 0) {
100 DUNE_THROW(NotImplemented,
"Desired derivative order is not implemented");
Definition brezzidouglasmarini1cube2d.hh:14
Meta-basis turning a scalar basis into vector-valued basis.
Definition meta/power/basis.hh:24
PowerBasis(const Backend &backend_)
Construct a PowerBasis.
Definition meta/power/basis.hh:60
std::size_t size() const
Number of shape functions.
Definition meta/power/basis.hh:63
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:92
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:80
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:68
std::size_t order() const
Polynomial order of the shape functions for quadrature.
Definition meta/power/basis.hh:65
types of domain and range
Definition meta/power/basis.hh:34
FieldMatrix< typename Traits::RangeField, dimR, Traits::dimDomainGlobal > Jacobian
Jacobian properties.
Definition meta/power/basis.hh:46
static const std::size_t diffOrder
maximum number of partial derivatives supported
Definition meta/power/basis.hh:49
FieldVector< typename Traits::RangeField, dimR > Range
Type used for range values.
Definition meta/power/basis.hh:38
static const std::size_t dimRange
Dimension of the range values.
Definition meta/power/basis.hh:36