1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
|
// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_FUNCTIONS_FUNCTIONSPACEBASES_SIZEINFO_HH
#define DUNE_FUNCTIONS_FUNCTIONSPACEBASES_SIZEINFO_HH
#include <array>
namespace Dune {
namespace Functions {
/**
* \brief A class encapsulating size information
*
* This class encapsulates size information of a basis.
*
* We may want to return this or a similar class when
* calling basis.size() without arguments. Until we
* decided on this, we can use sizeInfo(basis) and
* use the result as size provider for the HierarchicVectorWrapper.
*/
template<class B>
class SizeInfo
{
public:
using Basis = B;
using size_type = typename Basis::size_type;
using SizePrefix = typename Basis::SizePrefix;
/**
* \brief Construct from basis
*/
SizeInfo(const Basis& basis) :
basis_(&basis)
{}
/**
* \brief Return number possible values for next position in multi index
*/
size_type operator()(const SizePrefix& prefix) const
{
return basis_->size(prefix);
}
/**
* \brief Return number possible values for next position in multi index
*
* This shall vanish. It's just here such that this can be used
* as size provider n place of the basis.
*/
size_type size(const SizePrefix& prefix) const
{
return basis_->size(prefix);
}
operator size_type () const
{
return basis_->dimension();
}
protected:
const Basis* basis_;
};
template<class Basis>
SizeInfo<Basis> sizeInfo(const Basis& basis)
{
return SizeInfo<Basis>(basis);
}
} // end namespace Functions
} // end namespace Dune
#endif // DUNE_FUNCTIONS_FUNCTIONSPACEBASES_SIZEINFO_HH
|
