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
101
102
103
|
// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
// SPDX-FileCopyrightText: Copyright © DUNE Project contributors, see file AUTHORS.md
// SPDX-License-Identifier: LicenseRef-GPL-2.0-only-with-DUNE-exception OR LGPL-3.0-or-later
#ifndef DUNE_FUNCTIONS_COMMON_DIFFERENTIABLE_FUNCTION_IMP_HH
#define DUNE_FUNCTIONS_COMMON_DIFFERENTIABLE_FUNCTION_IMP_HH
#include <dune/common/exceptions.hh>
#include <dune/common/concept.hh>
#include <dune/functions/common/type_traits.hh>
namespace Dune {
namespace Functions {
namespace Imp {
/**
* A concept describing types that have a derivative() method found by ADL
*/
struct HasFreeDerivative
{
template<class F>
auto require(F&& f) -> decltype(
derivative(f)
);
};
template<class Dummy, class F,
std::enable_if_t<
models< HasFreeDerivative, F>() , int> = 0>
auto derivativeIfImplemented(const F& f) -> decltype(derivative(f))
{
return derivative(f);
}
template<class Dummy, class F,
std::enable_if_t<
not(models< HasFreeDerivative, F>()) , int> = 0>
Dummy derivativeIfImplemented(const F& f)
{
DUNE_THROW(Dune::NotImplemented, "Derivative not implemented");
}
template<class Signature, class DerivativeInterface>
class DifferentiableFunctionWrapperInterface
{};
// Interface of type erasure wrapper
//
// Notice that the basic interface of polymorphic classes (destructor, clone, ...)
// will be added by the type erasure foundation classes.
template<class Range, class Domain, class DerivativeInterface>
class DifferentiableFunctionWrapperInterface<Range(Domain), DerivativeInterface>
{
public:
virtual Range operator() (const Domain& x) const = 0;
virtual DerivativeInterface derivative() const = 0;
};
template<class Signature, class DerivativeInterface, class B>
class DifferentiableFunctionWrapperImplementation
{};
// Implementation of type erasure wrapper
template<class Range, class Domain, class DerivativeInterface, class B>
class DifferentiableFunctionWrapperImplementation< Range(Domain), DerivativeInterface, B> :
public B
{
public:
using B::B;
using Wrapped = typename B::Wrapped;
virtual Range operator() (const Domain& x) const
{
return this->get()(x);
}
virtual DerivativeInterface derivative() const
{
return derivativeIfImplemented<DerivativeInterface, Wrapped>(this->get());
}
};
}}} // namespace Dune::Functions::Imp
#endif // DUNE_FUNCTIONS_COMMON_DIFFERENTIABLE_FUNCTION_IMP_HH
|
