5#ifndef DUNE_LAGRANGEBASIS_INTERPOLATION_HH
6#define DUNE_LAGRANGEBASIS_INTERPOLATION_HH
12#include <dune/common/typeutilities.hh>
20 template<
template <
class,
unsigned int>
class LP,
21 unsigned int dim,
class F >
22 struct LagrangeInterpolationFactory;
27 template<
template <
class,
unsigned int>
class LP,
unsigned int dim,
class F >
34 typedef typename LagrangePointSet::Field
Field;
36 static const unsigned int dimension = LagrangePointSet::dimension;
48 template<
class Fn,
class Vector >
49 auto interpolate (
const Fn &fn, Vector &coefficients, PriorityTag< 1 > )
const
50 -> std::enable_if_t< std::is_invocable_v<
const Fn &,
decltype( this->lagrangePoints_.begin()->point() ) > >
52 unsigned int index = 0;
53 for(
const auto &lp : lagrangePoints_ )
54 field_cast( fn( lp.point() ), coefficients[ index++ ] );
56 template<
class Fn,
class Vector >
57 auto interpolate (
const Fn &fn, Vector &coefficients, PriorityTag< 0 > )
const
58 -> std::enable_if_t< models<Impl::FunctionWithEvaluate< typename Fn::DomainType, typename Fn::RangeType >, Fn>(),
void>
60 unsigned int index = 0;
61 for(
const auto &lp : lagrangePoints_ )
63 typename Fn::RangeType val;
64 fn.evaluate( field_cast< typename Fn::DomainType::field_type >( lp.point() ), val );
70 template<
class Fn,
class Vector >
72 -> std::enable_if_t< std::is_same< decltype(std::declval<Vector>().resize(1) ),
void >
::value,
void>
74 coefficients.resize( lagrangePoints_.size() );
75 interpolate( fn, coefficients, PriorityTag< 42 >() );
78 template<
class Basis,
class Matrix >
79 auto interpolate (
const Basis &basis, Matrix &coefficients )
const
80 -> std::enable_if_t< std::is_same<
81 decltype(std::declval<Matrix>().rowPtr(0)),
typename Matrix::Field* >
::value,
void>
83 coefficients.resize( lagrangePoints_.size(), basis.size( ) );
85 unsigned int index = 0;
86 for(
const auto &lp : lagrangePoints_ )
87 basis.template evaluate< 0 >( lp.point(), coefficients.rowPtr( index++ ) );
97 template<
template <
class,
unsigned int>
class LP,
98 unsigned int dim,
class F >
107 template< GeometryType::Id geometryId >
111 = LagrangePointSetFactory::template create< geometryId >( key );
112 if ( lagrangeCoeff == 0 )
115 return new Object( *lagrangeCoeff );
117 template< GeometryType::Id geometryId >
Definition: bdfmcube.hh:18
void field_cast(const F1 &f1, F2 &f2)
a helper class to cast from one field to another
Definition: field.hh:159
@ value
Definition: tensor.hh:168
Definition: lagrange/interpolation.hh:100
LagrangeCoefficientsFactory< LP, dim, F > LagrangePointSetFactory
Definition: lagrange/interpolation.hh:101
LagrangePointSetFactory::Object LagrangePointSet
Definition: lagrange/interpolation.hh:102
static bool supports(const Key &key)
Definition: lagrange/interpolation.hh:118
LagrangePointSetFactory::Key Key
Definition: lagrange/interpolation.hh:104
static Object * create(const Key &key)
Definition: lagrange/interpolation.hh:108
static void release(Object *object)
Definition: lagrange/interpolation.hh:122
const LocalLagrangeInterpolation< LP, dim, F > Object
Definition: lagrange/interpolation.hh:105
Definition: lagrange/interpolation.hh:29
LagrangePointSet::Field Field
Definition: lagrange/interpolation.hh:34
const LagrangePointSet & lagrangePoints() const
Definition: lagrange/interpolation.hh:90
static const unsigned int dimension
Definition: lagrange/interpolation.hh:36
LP< F, dim > LagrangePointSet
Definition: lagrange/interpolation.hh:33
auto interpolate(const Fn &fn, Vector &coefficients) const -> std::enable_if_t< std::is_same< decltype(std::declval< Vector >().resize(1)), void >::value, void >
Definition: lagrange/interpolation.hh:71
auto interpolate(const Basis &basis, Matrix &coefficients) const -> std::enable_if_t< std::is_same< decltype(std::declval< Matrix >().rowPtr(0)), typename Matrix::Field * >::value, void >
Definition: lagrange/interpolation.hh:79
Definition: lagrangecoefficients.hh:23
static void release(Object *object)
Definition: lagrangecoefficients.hh:42
const typedef LP< F, dim > Object
Definition: lagrangecoefficients.hh:25
std::size_t Key
Definition: lagrangecoefficients.hh:26