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
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
|
// SPDX-FileCopyrightText: Copyright © DUNE Project contributors, see file LICENSE.md in module root
// SPDX-License-Identifier: LicenseRef-GPL-2.0-only-with-DUNE-exception
// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_GEOGRID_COORDFUNCTIONCALLER_HH
#define DUNE_GEOGRID_COORDFUNCTIONCALLER_HH
#include <dune/grid/geometrygrid/hostcorners.hh>
#include <dune/grid/geometrygrid/coordfunction.hh>
namespace Dune
{
namespace GeoGrid
{
// CoordFunctionCaller
// -------------------
template< class HostEntity, class CoordFunctionInterface >
class CoordFunctionCaller;
template< class HostEntity, class ct, unsigned int dimD, unsigned int dimR, class Impl >
class CoordFunctionCaller< HostEntity, AnalyticalCoordFunctionInterface< ct, dimD, dimR, Impl > >
{
typedef AnalyticalCoordFunctionInterface< ct, dimD, dimR, Impl > CoordFunctionInterface;
typedef CoordFunctionCaller< HostEntity, CoordFunctionInterface > This;
static const int codimension = HostEntity::codimension;
public:
typedef typename CoordFunctionInterface::RangeVector RangeVector;
CoordFunctionCaller ( const HostEntity &hostEntity,
const CoordFunctionInterface &coordFunction )
: hostCorners_( hostEntity ),
coordFunction_( coordFunction )
{}
void evaluate ( unsigned int i, RangeVector &y ) const
{
coordFunction_.evaluate( hostCorners_[ i ], y );
}
GeometryType type () const
{
return hostCorners_.type();
}
std::size_t size () const
{
return hostCorners_.size();
}
private:
const HostCorners< HostEntity > hostCorners_;
const CoordFunctionInterface &coordFunction_;
};
template< class HostEntity, class ct, unsigned int dimR, class Impl >
class CoordFunctionCaller< HostEntity, DiscreteCoordFunctionInterface< ct, dimR, Impl > >
{
typedef DiscreteCoordFunctionInterface< ct, dimR, Impl > CoordFunctionInterface;
typedef CoordFunctionCaller< HostEntity, CoordFunctionInterface > This;
typedef typename CoordFunctionInterface::RangeVector RangeVector;
public:
CoordFunctionCaller ( const HostEntity &hostEntity,
const CoordFunctionInterface &coordFunction )
: hostEntity_( hostEntity ),
coordFunction_( coordFunction )
{}
void evaluate ( unsigned int i, RangeVector &y ) const
{
coordFunction_.evaluate( hostEntity_, i, y );
}
GeometryType type () const
{
return hostEntity_.type();
}
std::size_t size () const
{
auto refElement = referenceElement< ct, HostEntity::mydimension >( type() );
return refElement.size( HostEntity::mydimension );
}
private:
const HostEntity &hostEntity_;
const CoordFunctionInterface &coordFunction_;
};
} // namespace GeoGrid
} // namespace Dune
#endif // #ifndef DUNE_GEOGRID_COORDFUNCTIONCALLER_HH
|
