dune-localfunctions
2.7.0
|
Go to the documentation of this file.
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/typetraits.hh>
14 #include <dune/geometry/type.hh>
27 template<
class LocalBasisTraits, std::
size_t dimDomainGlobal_>
39 typedef FieldMatrix<RangeField, dimRange, dimDomainGlobal>
Jacobian;
62 template<
class LocalBasis,
class Geometry>
64 static_assert(LocalBasis::Traits::dimRange == 1,
65 "ScalarLocalToGlobalBasisAdaptor can only wrap a "
66 "scalar local basis.");
67 static_assert((std::is_same<
typename LocalBasis::Traits::DomainFieldType,
68 typename Geometry::ctype>::
value),
69 "ScalarLocalToGlobalBasisAdaptor: LocalBasis must use "
70 "the same ctype as Geometry");
72 (
static_cast<std::size_t
>(LocalBasis::Traits::dimDomain) ==
73 static_cast<std::size_t
>(Geometry::mydimension),
74 "ScalarLocalToGlobalBasisAdaptor: LocalBasis domain dimension must "
75 "match local dimension of Geometry");
77 const LocalBasis& localBasis;
95 const Geometry& geometry_) :
96 localBasis(localBasis_), geometry(geometry_)
99 std::size_t
size()
const {
return localBasis.size(); }
109 if(geometry.affine())
111 return localBasis.order();
118 std::vector<typename Traits::Range>& out)
const
120 localBasis.evaluateFunction(in, out);
124 std::vector<typename Traits::Jacobian>& out)
const
126 std::vector<typename LocalBasis::Traits::JacobianType>
127 localJacobian(
size());
128 localBasis.evaluateJacobian(in, localJacobian);
130 const typename Geometry::JacobianInverseTransposed &geoJacobian =
131 geometry.jacobianInverseTransposed(in);
134 for(std::size_t i = 0; i <
size(); ++i)
135 geoJacobian.mv(localJacobian[i][0], out[i][0]);
146 template<
class LocalInterpolation,
class Traits_>
148 const LocalInterpolation& localInterpolation;
163 (
const LocalInterpolation& localInterpolation_) :
164 localInterpolation(localInterpolation_)
167 template<
class Function,
class Coeff>
168 void interpolate(
const Function&
function, std::vector<Coeff>& out)
const
169 { localInterpolation.interpolate(
function, out); }
184 template<
class LocalFiniteElement,
class Geometry>
191 Traits::LocalBasisType, Geometry>
Basis;
195 typedef typename LocalFiniteElement::Traits::LocalCoefficientsType
200 const LocalFiniteElement &localFE;
216 (
const LocalFiniteElement& localFE_,
const Geometry &geometry) :
218 basis_(localFE.localBasis(), geometry),
219 interpolation_(localFE.localInterpolation())
224 {
return interpolation_; }
226 {
return localFE.localCoefficients(); }
227 GeometryType
type()
const {
return localFE.type(); }
241 template<
class LocalFiniteElement,
class Geometry>
243 const LocalFiniteElement& localFE;
259 (
const LocalFiniteElement &localFE_) : localFE(localFE_) {}
279 #endif // DUNE_LOCALFUNCTIONS_COMMON_LOCALTOGLOBALADAPTORS_HH
std::size_t order() const
return maximum polynomial order of the base function
Definition: localtoglobaladaptors.hh:108
ScalarLocalToGlobalBasisAdaptor(const LocalBasis &localBasis_, const Geometry &geometry_)
construct a ScalarLocalToGlobalBasisAdaptor
Definition: localtoglobaladaptors.hh:94
ScalarLocalToGlobalFiniteElementAdaptor< LocalFiniteElement, Geometry > FiniteElement
Definition: localtoglobaladaptors.hh:247
const Traits::Basis & basis() const
Definition: localtoglobaladaptors.hh:222
GeometryType type() const
Definition: localtoglobaladaptors.hh:227
@ dimDomain
dimension of the domain
Definition: common/localbasis.hh:39
static const std::size_t dimRange
Definition: localtoglobaladaptors.hh:36
LocalBasisTraits::RangeFieldType RangeField
Definition: localtoglobaladaptors.hh:35
@ dimRange
dimension of the range
Definition: common/localbasis.hh:51
std::size_t size() const
Definition: localtoglobaladaptors.hh:99
LocalToGlobalInterpolationAdaptor(const LocalInterpolation &localInterpolation_)
construct a LocalToGlobalInterpolationAdaptor
Definition: localtoglobaladaptors.hh:163
void evaluateFunction(const typename Traits::DomainLocal &in, std::vector< typename Traits::Range > &out) const
Definition: localtoglobaladaptors.hh:117
Definition: bdfmcube.hh:15
const FiniteElement make(const Geometry &geometry)
construct ScalarLocalToGlobalFiniteElementAdaptor
Definition: localtoglobaladaptors.hh:272
ScalarLocalToGlobalFiniteElementAdaptorFactory(const LocalFiniteElement &localFE_)
construct a ScalarLocalToGlobalFiniteElementAdaptorFactory
Definition: localtoglobaladaptors.hh:259
ScalarLocalToGlobalBasisAdaptor< typename LocalFiniteElement::Traits::LocalBasisType, Geometry > Basis
Definition: localtoglobaladaptors.hh:191
@ value
Definition: tensor.hh:166
RF RangeFieldType
Export type for range field.
Definition: common/localbasis.hh:46
LocalToGlobalInterpolationAdaptor< typename LocalFiniteElement::Traits::LocalInterpolationType, typename Basis::Traits > Interpolation
Definition: localtoglobaladaptors.hh:194
LocalBasisTraits::DomainFieldType DomainField
Definition: localtoglobaladaptors.hh:29
LocalToGlobalBasisAdaptorTraits< typename LocalBasis::Traits, Geometry::coorddimension > Traits
Definition: localtoglobaladaptors.hh:82
D DomainType
domain type
Definition: common/localbasis.hh:43
Traits_ Traits
Definition: localtoglobaladaptors.hh:151
Convert a simple scalar local basis into a global basis.
Definition: localtoglobaladaptors.hh:63
ScalarLocalToGlobalFiniteElementAdaptor(const LocalFiniteElement &localFE_, const Geometry &geometry)
construct a ScalarLocalToGlobalFiniteElementAdaptor
Definition: localtoglobaladaptors.hh:216
void evaluateJacobian(const typename Traits::DomainLocal &in, std::vector< typename Traits::Jacobian > &out) const
Definition: localtoglobaladaptors.hh:123
LocalBasisTraits::RangeType Range
Definition: localtoglobaladaptors.hh:37
R RangeType
range type
Definition: common/localbasis.hh:55
ImplementationDefined FiniteElement
Type of the finite element.
Definition: interface.hh:115
FieldVector< DomainField, dimDomainGlobal > DomainGlobal
Definition: localtoglobaladaptors.hh:33
LocalFiniteElement::Traits::LocalCoefficientsType Coefficients
Definition: localtoglobaladaptors.hh:196
Factory for ScalarLocalToGlobalFiniteElementAdaptor objects.
Definition: localtoglobaladaptors.hh:242
void interpolate(const Function &function, std::vector< Coeff > &out) const
Definition: localtoglobaladaptors.hh:168
FieldMatrix< RangeField, dimRange, dimDomainGlobal > Jacobian
Definition: localtoglobaladaptors.hh:39
const Traits::Interpolation & interpolation() const
Definition: localtoglobaladaptors.hh:223
Convert a simple scalar local finite element into a global finite element.
Definition: localtoglobaladaptors.hh:185
Traits class for local-to-global basis adaptors.
Definition: localtoglobaladaptors.hh:28
const Traits::Coefficients & coefficients() const
Definition: localtoglobaladaptors.hh:225
Definition: localtoglobaladaptors.hh:189
DF DomainFieldType
Export type for domain field.
Definition: common/localbasis.hh:34
static const std::size_t dimDomainLocal
Definition: localtoglobaladaptors.hh:30
static const std::size_t dimDomainGlobal
Definition: localtoglobaladaptors.hh:31
Convert a local interpolation into a global interpolation.
Definition: localtoglobaladaptors.hh:147
LocalBasisTraits::DomainType DomainLocal
Definition: localtoglobaladaptors.hh:32