4 #ifndef DUNE_LOCALFUNCTIONS_COMMON_LOCALTOGLOBALADAPTORS_HH
5 #define DUNE_LOCALFUNCTIONS_COMMON_LOCALTOGLOBALADAPTORS_HH
10 #include <dune/common/fmatrix.hh>
11 #include <dune/common/fvector.hh>
12 #include <dune/common/static_assert.hh>
13 #include <dune/common/typetraits.hh>
15 #include <dune/geometry/type.hh>
28 template<
class LocalBasisTraits, std::
size_t dimDomainGlobal_>
40 typedef FieldMatrix<RangeField, dimRange, dimDomainGlobal>
Jacobian;
65 template<
class LocalBasis,
class Geometry>
67 dune_static_assert(LocalBasis::Traits::dimRange == 1,
68 "ScalarLocalToGlobalBasisAdaptor can only wrap a "
69 "scalar local basis.");
70 dune_static_assert((is_same<
typename LocalBasis::Traits::DomainFieldType,
71 typename Geometry::ctype>::
value),
72 "ScalarLocalToGlobalBasisAdaptor: LocalBasis must use "
73 "the same ctype as Geometry");
75 ( static_cast<std::size_t>(LocalBasis::Traits::dimDomain) ==
76 static_cast<std::size_t>(Geometry::mydimension),
77 "ScalarLocalToGlobalBasisAdaptor: LocalBasis domain dimension must "
78 "match local dimension of Geometry");
80 const LocalBasis& localBasis;
81 const Geometry& geometry;
98 const Geometry& geometry_) :
99 localBasis(localBasis_), geometry(geometry_)
102 std::size_t
size()
const {
return localBasis.size(); }
112 if(geometry.affine())
114 return localBasis.order();
121 std::vector<typename Traits::Range>& out)
const
123 localBasis.evaluateFunction(in, out);
127 std::vector<typename Traits::Jacobian>& out)
const
129 std::vector<typename LocalBasis::Traits::JacobianType>
130 localJacobian(
size());
131 localBasis.evaluateJacobian(in, localJacobian);
133 const typename Geometry::Jacobian &geoJacobian =
134 geometry.jacobianInverseTransposed(in);
137 for(std::size_t i = 0; i <
size(); ++i)
138 geoJacobian.mv(localJacobian[i][0], out[i][0]);
149 template<
class LocalInterpolation,
class Traits_>
151 const LocalInterpolation& localInterpolation;
166 (
const LocalInterpolation& localInterpolation_) :
167 localInterpolation(localInterpolation_)
170 template<
class Function,
class Coeff>
171 void interpolate(
const Function&
function, std::vector<Coeff>& out)
const
172 { localInterpolation.interpolate(
function, out); }
187 template<
class LocalFiniteElement,
class Geometry>
194 Traits::LocalBasisType, Geometry>
Basis;
198 typedef typename LocalFiniteElement::Traits::LocalCoefficientsType
203 const LocalFiniteElement &localFE;
219 (
const LocalFiniteElement& localFE_,
const Geometry &geometry) :
221 basis_(localFE.localBasis(), geometry),
222 interpolation_(localFE.localInterpolation())
227 {
return interpolation_; }
229 {
return localFE.localCoefficients(); }
230 GeometryType
type()
const {
return localFE.type(); }
244 template<
class LocalFiniteElement,
class Geometry>
246 const LocalFiniteElement& localFE;
262 (
const LocalFiniteElement &localFE_) : localFE(localFE_) {}
282 #endif // DUNE_LOCALFUNCTIONS_COMMON_LOCALTOGLOBALADAPTORS_HH