4#ifndef DUNE_LOCALFUNCTIONS_META_POWER_INTERPOLATION_HH
5#define DUNE_LOCALFUNCTIONS_META_POWER_INTERPOLATION_HH
22 template<
class Backend,
class BasisTraits>
24 static_assert(Backend::Traits::dimRange == 1,
25 "PowerInterpolation works only with scalar backends");
27 const Backend *backend;
44 class ComponentEvaluator {
49 ComponentEvaluator(
const F &f_, std::size_t comp_) :
53 void evaluate(
const typename Backend::Traits::DomainLocal &x,
54 typename Backend::Traits::Range &y)
const
56 typename Traits::Range fy;
73 template<
typename F,
typename C>
77 for(std::size_t d = 0; d < Traits::dimRange; ++d) {
78 backend->interpolate(ComponentEvaluator<F>(f, d), cout);
80 out.resize(cout.size()*Traits::dimRange);
82 assert(out.size() == cout.size()*Traits::dimRange);
83 std::copy(cout.begin(), cout.end(), out.begin() + d*cout.size());
Meta-interpolation turning a scalar interpolation into vector-valued interpolation.
Definition: interpolation.hh:23
BasisTraits Traits
Export basis traits.
Definition: interpolation.hh:31
void interpolate(const F &f, std::vector< C > &out) const
Determine coefficients interpolating a given function.
Definition: interpolation.hh:74
PowerInterpolation(const Backend &backend_)
Construct a PowerInterpolation.
Definition: interpolation.hh:40
Dune namespace.
Definition: alignedallocator.hh:10