│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
41 using LocalFiniteElement = LFE;
│ │ │ │ +
42 using LocalBasis =
typename LFE::Traits::LocalBasisType;
│ │ │ │ +
43 using LocalBasisTraits =
typename LocalBasis::Traits;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
47 using ctype =
typename LocalBasisTraits::DomainFieldType;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
75 inline static bool isPrism (
unsigned int topologyId,
int dim,
int codim = 0 ) noexcept
│ │ │ │ -
│ │ │ │ -
77 assert( (dim > 0) && (topologyId < numTopologies( dim )) );
│ │ │ │ -
78 assert( (0 <= codim) && (codim < dim) );
│ │ │ │ -
79 return (( (topologyId | 1) & (1u << (dim-codim-1))) != 0);
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
89 inline static unsigned int baseTopologyId (
unsigned int topologyId,
int dim,
int codim = 1 ) noexcept
│ │ │ │ -
│ │ │ │ -
91 assert( (dim >= 0) && (topologyId < numTopologies( dim )) );
│ │ │ │ -
92 assert( (0 <= codim) && (codim <= dim) );
│ │ │ │ -
93 return topologyId & ((1u << (dim-codim)) - 1);
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
65 using Jacobian = FieldMatrix<ctype, coorddimension, mydimension>;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
104 const LocalFiniteElement& localFE,
│ │ │ │ +
105 std::vector<GlobalCoordinate> vertices)
│ │ │ │ +
106 : refElement_(refElement)
│ │ │ │ +
│ │ │ │ +
108 , vertices_(
std::move(vertices))
│ │ │ │ +
│ │ │ │ +
110 assert(localFE_.size() == vertices_.size());
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
126 template <
class Param,
│ │ │ │ +
127 std::enable_if_t<std::is_invocable_r_v<GlobalCoordinate,Param,LocalCoordinate>,
int> = 0>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
129 const LocalFiniteElement& localFE,
│ │ │ │ +
130 Param&& parametrization)
│ │ │ │ +
131 : refElement_(refElement)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
134 localFE_.localInterpolation().interpolate(parametrization, vertices_);
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
138 unsigned int topologyId_;
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
146 enum class IdType : std::uint64_t
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
190 constexpr operator Id()
const
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
194 std::uint64_t
id = dim_ | (std::uint64_t(none_) << 8) | (std::uint64_t(topologyId_) << 32);
│ │ │ │ -
195 return static_cast<Id
>(
id);
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
212 return static_cast<Id
>(*this);
│ │ │ │ -
│ │ │ │ +
143 template <
class... Args>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
162 affine_.emplace(affineImpl());
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
169 return refElement_.type();
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
181 assert( (i >= 0) && (i <
corners()) );
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
188 return global(refElement_.position(0, 0));
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
204 thread_local std::vector<typename LocalBasisTraits::RangeType> shapeValues;
│ │ │ │ +
│ │ │ │ +
206 assert(shapeValues.size() == vertices_.size());
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
209 for (std::size_t i = 0; i < shapeValues.size(); ++i)
│ │ │ │ +
210 out.axpy(shapeValues[i], vertices_[i]);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
223 : dim_(static_cast<
std::uint64_t>(
id) & 0xFF)
│ │ │ │ -
224 , none_(static_cast<
std::uint64_t>(
id) & 0x100)
│ │ │ │ -
225 , topologyId_(static_cast<
std::uint64_t>(
id) >> 32)
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │
│ │ │ │ -
│ │ │ │ -
233 : dim_(0), none_(true), topologyId_(0)
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
243 : dim_(
dim), none_(
isNone), topologyId_(topologyId)
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
252 : dim_(
dim), none_(false), topologyId_(topologyId)
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
265 template<
class TopologyType,
│ │ │ │ -
266 class = std::void_t<
decltype(TopologyType::dimension),
decltype(TopologyType::id)>>
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
268 : dim_(TopologyType::dimension), none_(false), topologyId_(TopologyType::
id)
│ │ │ │ -
│ │ │ │ -
270 DUNE_UNUSED_PARAMETER(t);
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
235 Impl::GaussNewtonErrorCode err = Impl::gaussNewton(
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
241 if (err != Impl::GaussNewtonErrorCode::OK)
│ │ │ │ +
242 DUNE_THROW(Dune::Exception,
│ │ │ │ +
243 "Local coordinate can not be recovered from global coordinate, error code = " <<
int(err) <<
"\n"
│ │ │ │ +
244 <<
" (global(x) - y).two_norm() = " << (
global(x) - y).two_norm()
│ │ │ │ +
245 <<
" > tol = " << opts.absTol);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
283 for (
int p = 2; p < opts.maxIt; ++p) {
│ │ │ │ +
│ │ │ │ +
285 if (abs(vol1 - vol0) < opts.absTol)
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
290 return ! none_ && dim_==2 && (topologyId_ | 1) == 0b0001;
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │ -
│ │ │ │ -
295 return ! none_ && dim_==2 && (topologyId_ | 1) == 0b0011;
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
300 return ! none_ && dim_==3 && (topologyId_ | 1) == 0b0001;
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
305 return ! none_ && dim_==3 && (topologyId_ | 1) == 0b0011;
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
310 return ! none_ && dim_==3 && (topologyId_ | 1) == 0b0101;
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
297 for (
const auto& qp : quadRule)
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
309 thread_local std::vector<typename LocalBasisTraits::JacobianType> shapeJacobians;
│ │ │ │ +
│ │ │ │ +
311 assert(shapeJacobians.size() == vertices_.size());
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
315 return ! none_ && dim_==3 && (topologyId_ | 1) == 0b0111;
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
320 return ! none_ && (topologyId_ | 1) == 1;
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
325 return ! none_ && ((topologyId_ ^ ((1 << dim_)-1)) >> 1 == 0);
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
330 return ! none_ && (((topologyId_ & ~1) & (1u << (dim_-1))) == 0);
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
338 return ! none_ && (((topologyId_ & ~1) & (1u << step)) == 0);
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
343 return ! none_ && (( (topologyId_ | 1) & (1u << (dim_-1))) != 0);
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
314 for (std::size_t i = 0; i < shapeJacobians.size(); ++i) {
│ │ │ │ +
315 for (
int j = 0; j < Jacobian::rows; ++j) {
│ │ │ │ +
316 shapeJacobians[i].umtv(vertices_[i][j], out[j]);
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
351 return ! none_ && (( (topologyId_ | 1) & (1u << step)) != 0);
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
360 constexpr unsigned int dim()
const {
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
361 return geometry.refElement_;
│ │ │ │ +
│ │ │ │
│ │ │ │
│ │ │ │
│ │ │ │ -
365 constexpr unsigned int id()
const {
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
378 return ( ( none_ == other.none_ )
│ │ │ │ -
379 && ( ( none_ ==
true )
│ │ │ │ -
380 || ( ( dim_ == other.dim_ )
│ │ │ │ -
381 && ( (topologyId_ >> 1) == (other.topologyId_ >> 1) )
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
389 return ! ((*this)==other);
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
394 return ( ( none_ < other.none_ )
│ │ │ │ -
395 || ( !( other.none_ < none_ )
│ │ │ │ -
396 && ( ( dim_ < other.dim_ )
│ │ │ │ -
397 || ( (other.dim_ == dim_)
│ │ │ │ -
398 && ((topologyId_ >> 1) < (other.topologyId_ >> 1) )
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
414 s <<
"(simplex, " << a.
dim() <<
")";
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
419 s <<
"(cube, " << a.
dim() <<
")";
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
434 s <<
"(none, " << a.
dim() <<
")";
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
437 s <<
"(other [" << a.
id() <<
"], " << a.
dim() <<
")";
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
447 namespace GeometryTypes {
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
379 return localFE_.localBasis();
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
384 bool affineImpl ()
const
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
405 auto refSimplex = referenceElement<ctype,mydimension>(GeometryTypes::simplex(
mydimension));
│ │ │ │ +
406 auto simplexIndices = Dune::range(refSimplex.size(
mydimension));
│ │ │ │ +
407 auto simplexCorners = Dune::transformedRangeView(simplexIndices,
│ │ │ │ +
│ │ │ │ +
409 AffineGeometry<ctype,mydimension,coorddimension> affineGeo(refSimplex,simplexCorners);
│ │ │ │ +
│ │ │ │ +
411 const ctype tol = sqrt(std::numeric_limits<ctype>::epsilon());
│ │ │ │ +
412 for (
int i = 0; i <
corners(); ++i) {
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
427 LocalFiniteElement localFE_{};
│ │ │ │ +
│ │ │ │ +
430 std::vector<GlobalCoordinate> vertices_{};
│ │ │ │ +
│ │ │ │ +
432 mutable std::optional<bool> affine_ = std::nullopt;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
439using LocalCoordinate_t
│ │ │ │ +
440 = FieldVector<
typename LFE::Traits::LocalBasisType::Traits::DomainFieldType,
│ │ │ │ +
441 LFE::Traits::LocalBasisType::Traits::dimDomain>;
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
447template <
class I,
class LFE,
class GlobalCoordinate>
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
451template <
class I,
class LFE,
class F,
│ │ │ │ +
452 class Range = std::invoke_result_t<F,Impl::LocalCoordinate_t<LFE>>>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
456template <
class LFE,
class GlobalCoordinate>
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
464 return GeometryType(((dim>1) ? ((1 << dim) - 1) : 0),dim,
false);
│ │ │ │ -
│ │ │ │ +
460template <
class LFE,
class F,
│ │ │ │ +
461 class Range = std::invoke_result_t<F,Impl::LocalCoordinate_t<LFE>>>
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
550 template<
template< GeometryType::Id >
class Operation,
int dim,
GeometryType::Id geometryId = GeometryTypes::vertex >
│ │ │ │ -
551 struct IfGeometryType
│ │ │ │ -
│ │ │ │ -
553 static constexpr GeometryType geometry = geometryId;
│ │ │ │ -
554 template<
class... Args >
│ │ │ │ -
555 static auto apply ( GeometryType gt, Args &&... args )
│ │ │ │ -
│ │ │ │ -
557 GeometryType lowerGeometry(gt.id() >>1 , gt.dim()-1, gt.isNone());
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
560 return IfGeometryType< Operation, dim-1, GeometryTypes::prismaticExtension(geometry).toId() >::apply( lowerGeometry, std::forward< Args >( args )... );
│ │ │ │ -
│ │ │ │ -
562 return IfGeometryType< Operation, dim-1, GeometryTypes::conicalExtension(geometry).toId() >::apply( lowerGeometry, std::forward< Args >( args )... );
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
566 template<
template< GeometryType::Id >
class Operation,
GeometryType::Id geometryId >
│ │ │ │ -
567 struct IfGeometryType< Operation, 0, geometryId>
│ │ │ │ -
│ │ │ │ -
569 template<
class... Args >
│ │ │ │ -
570 static auto apply ([[maybe_unused]] GeometryType gt, Args &&... args )
│ │ │ │ -
│ │ │ │ -
572 return Operation< geometryId >::apply( std::forward< Args >( args )... );
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ -
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
An implementation of the Geometry interface for affine geometries.
│ │ │ │ +
│ │ │ │ +
│ │ │ │ +
A unique label for each type of element that can occur in a grid.
│ │ │ │
│ │ │ │
Definition affinegeometry.hh:21
│ │ │ │ -
std::ostream & operator<<(std::ostream &s, const GeometryType &a)
Prints the type to an output stream.
Definition type.hh:410
│ │ │ │ +
This class provides access to geometric and topological properties of a reference element.
Definition referenceelement.hh:52
│ │ │ │ +
Class providing access to the singletons of the reference elements.
Definition referenceelements.hh:128
│ │ │ │ +
typename Container::ReferenceElement ReferenceElement
The reference element type.
Definition referenceelements.hh:146
│ │ │ │ +
Geometry implementation based on local-basis function parametrization.
Definition localfiniteelementgeometry.hh:40
│ │ │ │ +
GeometryType type() const
Obtain the name of the reference element.
Definition localfiniteelementgeometry.hh:167
│ │ │ │ +
decltype(power(std::declval< ctype >(), mydimension)) Volume
type of volume
Definition localfiniteelementgeometry.hh:62
│ │ │ │ +
LocalCoordinate local(const GlobalCoordinate &y, Impl::GaussNewtonOptions< ctype > opts={}) const
Evaluate the inverse coordinate mapping.
Definition localfiniteelementgeometry.hh:232
│ │ │ │ +
LocalFiniteElementGeometry(const ReferenceElement &refElement, const LocalFiniteElement &localFE, std::vector< GlobalCoordinate > vertices)
Constructor from a vector of coefficients of the LocalBasis parametrizing the geometry.
Definition localfiniteelementgeometry.hh:103
│ │ │ │ +
JacobianTransposed jacobianTransposed(const LocalCoordinate &local) const
Obtain the transposed of the Jacobian.
Definition localfiniteelementgeometry.hh:327
│ │ │ │ +
const LocalBasis & localBasis() const
The local basis of the stored local finite-element.
Definition localfiniteelementgeometry.hh:377
│ │ │ │ +
typename LocalBasisTraits::DomainFieldType ctype
coordinate type
Definition localfiniteelementgeometry.hh:47
│ │ │ │ +
JacobianInverseTransposed jacobianInverseTransposed(const LocalCoordinate &local) const
Obtain the transposed of the Jacobian's inverse.
Definition localfiniteelementgeometry.hh:353
│ │ │ │ +
friend ReferenceElement referenceElement(const LocalFiniteElementGeometry &geometry)
Obtain the reference-element related to this geometry.
Definition localfiniteelementgeometry.hh:359
│ │ │ │ +
LocalFiniteElementGeometry()=default
Default constructed geometry results in an empty/invalid representation.
│ │ │ │ +
bool affine() const
Is this mapping affine? Geometries of order 1 might be affine, but it needs to be checked on non-simp...
Definition localfiniteelementgeometry.hh:159
│ │ │ │ +
Jacobian jacobian(const LocalCoordinate &local) const
Obtain the Jacobian.
Definition localfiniteelementgeometry.hh:307
│ │ │ │ +
Volume volume(const QuadratureRule< ctype, mydimension > &quadRule) const
Obtain the volume of the mapping's image by given quadrature rules.
Definition localfiniteelementgeometry.hh:294
│ │ │ │ +
const std::vector< GlobalCoordinate > & coefficients() const
Obtain the coefficients of the parametrization.
Definition localfiniteelementgeometry.hh:371
│ │ │ │ +
LocalFiniteElementGeometry(GeometryType gt, Args &&... args)
Constructor, forwarding to the other constructors that take a reference-element.
Definition localfiniteelementgeometry.hh:144
│ │ │ │ +
static const int coorddimension
coordinate dimension
Definition localfiniteelementgeometry.hh:53
│ │ │ │ +
int corners() const
Obtain number of corners of the corresponding reference element.
Definition localfiniteelementgeometry.hh:173
│ │ │ │ +
GlobalCoordinate center() const
Obtain the centroid of the mapping's image.
Definition localfiniteelementgeometry.hh:186
│ │ │ │ +
const LocalFiniteElement & finiteElement() const
Obtain the local finite-element.
Definition localfiniteelementgeometry.hh:365
│ │ │ │ +
GlobalCoordinate global(const LocalCoordinate &local) const
Evaluate the coordinate mapping.
Definition localfiniteelementgeometry.hh:202
│ │ │ │ +
Impl::FieldMatrixHelper< ctype > MatrixHelper
Definition localfiniteelementgeometry.hh:82
│ │ │ │ +
typename ReferenceElements::ReferenceElement ReferenceElement
Definition localfiniteelementgeometry.hh:79
│ │ │ │ +
GlobalCoordinate corner(int i) const
Obtain coordinates of the i-th corner.
Definition localfiniteelementgeometry.hh:179
│ │ │ │ +
Volume volume(Impl::ConvergenceOptions< ctype > opts={}) const
Obtain the volume of the mapping's image.
Definition localfiniteelementgeometry.hh:276
│ │ │ │ +
static const int mydimension
geometry dimension
Definition localfiniteelementgeometry.hh:50
│ │ │ │ +
JacobianInverse jacobianInverse(const LocalCoordinate &local) const
Obtain the Jacobian's inverse.
Definition localfiniteelementgeometry.hh:339
│ │ │ │ +
FieldVector< ctype, coorddimension > GlobalCoordinate
type of global coordinates
Definition localfiniteelementgeometry.hh:59
│ │ │ │ +
FieldMatrix< ctype, coorddimension, mydimension > Jacobian
type of jacobian
Definition localfiniteelementgeometry.hh:65
│ │ │ │ +
FieldMatrix< ctype, coorddimension, mydimension > JacobianInverseTransposed
type of jacobian inverse transposed
Definition localfiniteelementgeometry.hh:74
│ │ │ │ +
FieldMatrix< ctype, mydimension, coorddimension > JacobianInverse
type of jacobian inverse
Definition localfiniteelementgeometry.hh:71
│ │ │ │ +
ctype integrationElement(const LocalCoordinate &local) const
Obtain the integration element.
Definition localfiniteelementgeometry.hh:260
│ │ │ │ +
int order() const
Obtain the (highest) polynomial order of the parametrization.
Definition localfiniteelementgeometry.hh:149
│ │ │ │ +
FieldVector< ctype, mydimension > LocalCoordinate
type of local coordinates
Definition localfiniteelementgeometry.hh:56
│ │ │ │ +
FieldMatrix< ctype, mydimension, coorddimension > JacobianTransposed
type of jacobian transposed
Definition localfiniteelementgeometry.hh:68
│ │ │ │ +
LocalFiniteElementGeometry(const ReferenceElement &refElement, const LocalFiniteElement &localFE, Param &¶metrization)
Constructor from a local parametrization function, mapping local to (curved) global coordinates.
Definition localfiniteelementgeometry.hh:128
│ │ │ │ +
Abstract base class for quadrature rules.
Definition quadraturerules.hh:214
│ │ │ │ +
A container for all quadrature rules of dimension dim
Definition quadraturerules.hh:260
│ │ │ │
Unique label for each type of entities that can occur in DUNE grids.
Definition type.hh:114
│ │ │ │ -
constexpr GeometryType(unsigned int topologyId, unsigned int dim)
Constructor, using the topologyId (integer) and the dimension.
Definition type.hh:251
│ │ │ │ -
constexpr bool operator<(const GeometryType &other) const
less-than operation for use with maps
Definition type.hh:393
│ │ │ │ -
constexpr bool operator!=(const GeometryType &other) const
Check for inequality.
Definition type.hh:388
│ │ │ │ -
constexpr bool isPyramid() const
Return true if entity is a pyramid.
Definition type.hh:304
│ │ │ │ -
constexpr bool isTetrahedron() const
Return true if entity is a tetrahedron.
Definition type.hh:299
│ │ │ │ -
constexpr bool isPrism() const
Return true if entity is a prism.
Definition type.hh:309
│ │ │ │ -
constexpr bool isVertex() const
Return true if entity is a vertex.
Definition type.hh:279
│ │ │ │ -
constexpr bool operator==(const GeometryType &other) const
Check for equality. This method knows that in dimension 0 and 1 all BasicTypes are equal.
Definition type.hh:377
│ │ │ │ -
constexpr Id toId() const
Create an Id representation of this GeometryType.
Definition type.hh:210
│ │ │ │ -
constexpr bool isConical(const int &step) const
Return true if entity was constructed with a conical product in the chosen step.
Definition type.hh:337
│ │ │ │ -
constexpr unsigned int dim() const
Return dimension of the type.
Definition type.hh:360
│ │ │ │ -
constexpr bool isPrismatic(const int &step) const
Return true if entity was constructed with a prismatic product in the chosen step.
Definition type.hh:350
│ │ │ │ -
constexpr bool isTriangle() const
Return true if entity is a triangle.
Definition type.hh:289
│ │ │ │ -
GeometryType(TopologyType t)
Constructor from static TopologyType class.
Definition type.hh:267
│ │ │ │ -
constexpr GeometryType(unsigned int topologyId, unsigned int dim, bool isNone)
Constructor, using the topologyId (integer), the dimension and a flag for type none.
Definition type.hh:242
│ │ │ │ -
BasicType
Each entity can be tagged by one of these basic types plus its space dimension.
Definition type.hh:120
│ │ │ │ -
@ cube
Cube element in any nonnegative dimension.
Definition type.hh:122
│ │ │ │ -
@ simplex
Simplicial element in any nonnegative dimension.
Definition type.hh:121
│ │ │ │ -
@ pyramid
Four sided pyramid in three dimensions.
Definition type.hh:123
│ │ │ │ -
@ extended
Other, more general topology, representable as topologyId.
Definition type.hh:125
│ │ │ │ -
@ none
Even more general topology, cannot be specified by a topologyId. Two GeometryTypes with 'none' type a...
Definition type.hh:126
│ │ │ │ -
@ prism
Prism element in three dimensions.
Definition type.hh:124
│ │ │ │ -
constexpr GeometryType(Id id)
Reconstruct a Geometry type from a GeometryType::Id.
Definition type.hh:222
│ │ │ │ -
constexpr bool isCube() const
Return true if entity is a cube of any dimension.
Definition type.hh:324
│ │ │ │ -
constexpr GeometryType()
Default constructor, not initializing anything.
Definition type.hh:232
│ │ │ │ -
constexpr bool isConical() const
Return true if entity was constructed with a conical product in the last step.
Definition type.hh:329
│ │ │ │ -
constexpr bool isLine() const
Return true if entity is a line segment.
Definition type.hh:284
│ │ │ │ -
constexpr bool isQuadrilateral() const
Return true if entity is a quadrilateral.
Definition type.hh:294
│ │ │ │ -
constexpr bool isPrismatic() const
Return true if entity was constructed with a prismatic product in the last step.
Definition type.hh:342
│ │ │ │ -
constexpr unsigned int id() const
Return the topology id of the type.
Definition type.hh:365
│ │ │ │ -
constexpr bool isNone() const
Return true if entity is a singular of any dimension.
Definition type.hh:355
│ │ │ │ -
constexpr bool isHexahedron() const
Return true if entity is a hexahedron.
Definition type.hh:314
│ │ │ │
constexpr bool isSimplex() const
Return true if entity is a simplex of any dimension.
Definition type.hh:319
│ │ │ │ -
IdType Id
An integral id representing a GeometryType.
Definition type.hh:181
│ │ │ │
