4 #ifndef DUNE_LOCALFUNCTIONS_META_POWER_BASIS_HH 5 #define DUNE_LOCALFUNCTIONS_META_POWER_BASIS_HH 10 #include <dune/common/fmatrix.hh> 11 #include <dune/common/fvector.hh> 22 template<
class Backend, std::
size_t dimR>
24 static_assert(Backend::Traits::dimRange == 1,
25 "PowerBasis works only with scalar backends");
28 const Backend *backend;
32 struct Traits :
public Backend::Traits
37 typedef FieldVector<typename Traits::RangeField, dimR>
Range;
44 typedef FieldMatrix<
typename Traits::RangeField, dimR,
49 Backend::Traits::diffOrder > 0 ? 1 : 0;
59 PowerBasis(
const Backend &backend_) : backend(&backend_) { }
62 std::size_t
size ()
const {
return backend->size()*dimR; }
64 std::size_t
order ()
const {
return backend->order(); }
68 std::vector<typename Traits::Range>& out)
const 70 std::vector<typename Backend::Traits::Range> backendValues;
71 backend->evaluateFunction(in, backendValues);
73 for(std::size_t d = 0; d < dimR; ++d)
74 for(std::size_t i = 0; i < backend->size(); ++i)
75 out[d*backend->size()+i][d] = backendValues[i][0];
80 std::vector<typename Traits::Jacobian>& out)
const 82 std::vector<typename Backend::Traits::Jacobian> backendValues;
83 backend->evaluateJacobian(in, backendValues);
85 for(std::size_t d = 0; d < dimR; ++d)
86 for(std::size_t i = 0; i < backend->size(); ++i)
87 out[d*backend->size()+i][d] = backendValues[i][0];
93 #endif // DUNE_LOCALFUNCTIONS_META_POWER_BASIS_HH 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:67
static const std::size_t diffOrder
maximum number of partial derivatives supported
Definition: meta/power/basis.hh:48
Definition: brezzidouglasmarini1cube2dlocalbasis.hh:14
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:79
std::size_t order() const
Polynomial order of the shape functions for quadrature.
Definition: meta/power/basis.hh:64
Meta-basis turning a scalar basis into vector-valued basis.
Definition: meta/power/basis.hh:23
std::size_t size() const
Number of shape functions.
Definition: meta/power/basis.hh:62
FieldVector< typename Traits::RangeField, dimR > Range
Type used for range values.
Definition: meta/power/basis.hh:37
types of domain and range
Definition: meta/power/basis.hh:32
static const std::size_t dimRange
Dimension of the range values.
Definition: meta/power/basis.hh:35
PowerBasis(const Backend &backend_)
Construct a PowerBasis.
Definition: meta/power/basis.hh:59
FieldMatrix< typename Traits::RangeField, dimR, Traits::dimDomainGlobal > Jacobian
Jacobian properties.
Definition: meta/power/basis.hh:45