\n
\n
\n
\n
\n
\n-
7#ifndef DUNE_FUNCTIONS_COMMON_DEFAULT_DERIVATIVE_TRAITS_HH \n-
8#define DUNE_FUNCTIONS_COMMON_DEFAULT_DERIVATIVE_TRAITS_HH \n+
7#ifndef DUNE_FUNCTIONS_FUNCTIONSPACEBASES_ISTLVECTORBACKEND_HH \n+
8#define DUNE_FUNCTIONS_FUNCTIONSPACEBASES_ISTLVECTORBACKEND_HH \n
\n-
\n+
\n
\n-
\n-
13#include <dune/common/fvector.hh> \n-
14#include <dune/common/fmatrix.hh> \n-
\n-
\n-
\n+
\n+
\n+
14#include <dune/common/std/type_traits.hh> \n+
15#include <dune/common/indices.hh> \n+
16#include <dune/common/hybridutilities.hh> \n+
17#include <dune/common/concept.hh> \n
\n-
\n-
\n-
\n-
\n-
\n-
39template<
class Signature>
\n-
\n-
\n-
\n-
\n-
\n-
\n-
70template<
typename K,
int n>
\n-
\n-
\n-
86template<
typename K,
int n,
int m>
\n-
\n-
\n-
\n-
90 typedef FieldMatrix<K,m,n>
Range;
\n-
\n-
\n-
\n-
102template<
typename K,
int n,
int m>
\n-
\n-
\n-
\n-
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
29 std::enable_if_t<not Dune::models<Imp::Concept::HasStaticIndexAccess, V>() ,
int> = 0>
\n+
\n+
\n+
\n+
\n+
\n+
\n+
36 std::enable_if_t<Dune::models<Imp::Concept::HasStaticIndexAccess, V>(),
int> = 0>
\n+
\n+
\n+
39 if constexpr (Dune::models<Imp::Concept::HasDynamicIndexAccess<std::size_t>, V>())
\n+
\n+
\n+
\n+
43 auto indexRange =
typename decltype(range(Hybrid::size(v)))::integer_sequence();
\n+
44 return unpackIntegerSequence([&](
auto... i) {
\n+
45 return uniqueTypeList(std::tuple_cat(
fieldTypes(v[i])...));
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
69 return decltype(Impl::fieldTypes(std::declval<V>())){};
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
80 return std::tuple_size_v<std::decay_t<decltype(fieldTypes<V>())>> ==1;
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
105class ISTLVectorBackend
\n+
\n+
\n+
\n+
\n+
110 using dynamicIndexAccess_t =
decltype(std::declval<C>()[0]);
\n
\n-
\n-
\n+
\n+
113 using staticIndexAccess_t =
decltype(std::declval<C>()[Dune::Indices::_0]);
\n+
\n+
\n+
116 using resizeMethod_t =
decltype(std::declval<C>().resize(0));
\n+
\n+
\n+
\n+
\n+
\n+
122 using hasDynamicIndexAccess = Dune::Std::is_detected<dynamicIndexAccess_t, std::remove_reference_t<C>>;
\n+
\n+
\n+
125 using hasStaticIndexAccess = Dune::Std::is_detected<staticIndexAccess_t, std::remove_reference_t<C>>;
\n+
\n+
\n+
128 using hasResizeMethod = Dune::Std::is_detected<resizeMethod_t, std::remove_reference_t<C>>;
\n+
\n+
\n+
131 using isDynamicVector = Dune::Std::is_detected<dynamicIndexAccess_t, std::remove_reference_t<C>>;
\n+
\n+
\n+
134 using isStaticVector = Dune::Std::bool_constant<
\n+
135 Dune::Std::is_detected_v<staticIndexAccess_t, std::remove_reference_t<C>>
\n+
136 and not Dune::Std::is_detected_v<dynamicIndexAccess_t, std::remove_reference_t<C>>>;
\n+
\n+
\n+
139 using isScalar = Dune::Std::bool_constant<not Dune::Std::is_detected_v<staticIndexAccess_t, std::remove_reference_t<C>>>;
\n+
\n+
\n+
142 using isVector = Dune::Std::bool_constant<Dune::Std::is_detected_v<staticIndexAccess_t, std::remove_reference_t<C>>>;
\n+
\n+
\n+
\n+
146 template<
class... Args>
\n+
147 static void forwardToResize(Args&&... args)
\n+
\n+
149 resize(std::forward<Args>(args)...);
\n+
\n+
\n+
\n+
153 template<
class C,
class SizeProvider,
\n+
154 std::enable_if_t<hasResizeMethod<C>::value,
int> = 0>
\n+
155 static void resize(C&& c,
const SizeProvider& sizeProvider,
typename SizeProvider::SizePrefix prefix)
\n+
\n+
157 auto size = sizeProvider.size(prefix);
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
170 DUNE_THROW(RangeError,
"The vector entry v[" << prefix <<
"] should refer to a " \n+
171 <<
"scalar coefficient, but is a dynamically sized vector of size==0");
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
180 for(std::size_t i=0; i<size; ++i)
\n+
\n+
\n+
183 resize(c[i], sizeProvider, prefix);
\n+
\n+
\n+
\n+
187 template<
class C,
class SizeProvider,
\n+
188 std::enable_if_t<not hasResizeMethod<C>::value,
int> = 0,
\n+
189 std::enable_if_t<isVector<C>::value,
int> = 0>
\n+
190 static void resize(C&& c,
const SizeProvider& sizeProvider,
typename SizeProvider::SizePrefix prefix)
\n+
\n+
192 auto size = sizeProvider.size(prefix);
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
207 if (c.size() != size)
\n+
208 DUNE_THROW(RangeError,
"Can't resize non-resizable entry v[" << prefix <<
"] of size " << c.size() <<
" to size(" << prefix <<
")=" << size);
\n+
\n+
\n+
211 using namespace Dune::Hybrid;
\n+
\n+
213 forEach(integralRange(Hybrid::size(c)), [&](
auto&& i) {
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
221 ISTLVectorBackend<V>::forwardToResize(c[i], sizeProvider, prefix);
\n+
\n+
\n+
\n+
225 template<
class C,
class SizeProvider,
\n+
226 std::enable_if_t<not hasResizeMethod<C>::value,
int> = 0,
\n+
227 std::enable_if_t<isScalar<C>::value,
int> = 0>
\n+
228 static void resize(C&&,
const SizeProvider& sizeProvider,
typename SizeProvider::SizePrefix prefix)
\n+
\n+
230 auto size = sizeProvider.size(prefix);
\n+
\n+
232 DUNE_THROW(RangeError,
"Can't resize scalar vector entry v[" << prefix <<
"] to size(" << prefix <<
")=" << size);
\n+
\n+
\n+
235 template<
class C,
class T,
\n+
236 std::enable_if_t<std::is_assignable_v<C&,T>,
int> = 0>
\n+
237 void recursiveAssign(C& c,
const T& t)
\n+
\n+
\n+
\n+
\n+
242 template<
class C,
class T,
\n+
243 std::enable_if_t<not std::is_assignable_v<C&,T>,
int> = 0>
\n+
244 void recursiveAssign(C& c,
const T& t)
\n+
\n+
246 Dune::Hybrid::forEach(c, [&](
auto&& ci) {
\n+
247 recursiveAssign(ci, t);
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
255 ISTLVectorBackend(Vector& vector) :
\n+
\n+
\n+
\n+
259 template<
class SizeProv
ider>
\n+
260 void resize(
const SizeProvider& sizeProvider)
\n+
\n+
262 auto prefix =
typename SizeProvider::SizePrefix();
\n+
\n+
264 resize(*vector_, sizeProvider, prefix);
\n+
\n+
\n+
267 template<
class MultiIndex>
\n+
268 decltype(
auto)
operator[](
const MultiIndex& index)
const \n+
\n+
\n+
\n+
\n+
273 template<
class MultiIndex>
\n+
274 decltype(
auto)
operator[](
const MultiIndex& index)
\n+
\n+
\n+
\n+
\n+
\n+
288 void operator= (
const T& other)
\n+
\n+
290 recursiveAssign(vector(), other);
\n+
\n+
\n+
\n+
294 void operator= (
const ISTLVectorBackend<T>& other)
\n+
\n+
296 vector() = other.vector();
\n+
\n+
\n+
299 const Vector& vector()
const \n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
349template<
class Vector>
\n+
\n+
\n+
\n+
352 static_assert(hasUniqueFieldType<Vector&>(),
"Vector type passed to istlVectorBackend() does not have a unique field type.");
\n+
353 return Impl::ISTLVectorBackend<Vector>(v);
\n+
\n+
\n+
\n+
\n+
\n+
387template<
class Vector>
\n+
\n+
\n+
\n+
390 static_assert(hasUniqueFieldType<const Vector&>(),
"Vector type passed to istlVectorBackend() does not have a unique field type.");
\n+
391 return Impl::ISTLVectorBackend<const Vector>(v);
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
\n+
auto istlVectorBackend(Vector &v)
Return a vector backend wrapping non-const ISTL like containers.
Definition istlvectorbackend.hh:350
\n+
constexpr decltype(auto) resolveDynamicMultiIndex(C &&c, const MultiIndex &multiIndex, const IsFinal &isFinal)
Provide multi-index access by chaining operator[].
Definition indexaccess.hh:377
\n
Definition polynomial.hh:17
\n-
Dummy range class to be used if no proper type is available.
Definition defaultderivativetraits.hh:27
\n-
Default implementation for derivative traits.
Definition defaultderivativetraits.hh:41
\n-
InvalidRange Range
Range of derivative for function with given signature.
Definition defaultderivativetraits.hh:43
\n-
double Range
Range of derivative for function with given signature.
Definition defaultderivativetraits.hh:58
\n-
FieldVector< K, n > Range
Range of derivative for function with given signature.
Definition defaultderivativetraits.hh:74
\n-
FieldMatrix< K, m, n > Range
Range of derivative for function with given signature.
Definition defaultderivativetraits.hh:90
\n-
FieldMatrix< K, m, n > Range
Range of derivative for function with given signature.
Definition defaultderivativetraits.hh:106
\n+
constexpr auto fieldTypes()
Generate list of field types in container.
Definition istlvectorbackend.hh:67
\n+
constexpr bool hasUniqueFieldType()
Check if container has a unique field type.
Definition istlvectorbackend.hh:78
\n+
std::vector< Child > Vector
Descriptor for vectors with all children of the same type and dynamic size.
Definition containerdescriptors.hh:124
\n+
\n
![\"doxygen\"/](\"doxygen.svg\")
1.9.8\n \n \n", "details": [{"source1": "html2text {}", "source2": "html2text {}", "unified_diff": "@@ -1,44 +1,28 @@\n dune-functions\u00a02.10\n Loading...\n Searching...\n No Matches\n * _\bd_\bu_\bn_\be\n * _\bf_\bu_\bn_\bc_\bt_\bi_\bo_\bn_\bs\n- * _\bb_\ba_\bc_\bk_\be_\bn_\bd_\bs\n-_\bN_\ba_\bm_\be_\bs_\bp_\ba_\bc_\be_\bs | _\bF_\bu_\bn_\bc_\bt_\bi_\bo_\bn_\bs\n-istlvectorbackend.hh File Reference\n-#include