-
Notifications
You must be signed in to change notification settings - Fork 23
Expand file tree
/
Copy pathconfig.h
More file actions
92 lines (82 loc) · 2.26 KB
/
config.h
File metadata and controls
92 lines (82 loc) · 2.26 KB
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
/*
* Copyright (C) 2020-2026 MEmilio
*
* Authors: Daniel Abele
*
* Contact: Martin J. Kuehn <[email protected]>
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*/
/**
* Configuration of memilio library.
*/
#ifndef MIO_CONFIG_H
#define MIO_CONFIG_H
#include "memilio/config_internal.h"
using ScalarType = double;
namespace ad
{
namespace internal
{
// Forward declaration used for support of ad types in Limits.
// These templates are called AD_TAPE_REAL and DATA_HANDLER internally in AD.
template <class FP, class DataHandler>
struct active_type;
} // namespace internal
} // namespace ad
namespace mio
{
/**
* @brief Type specific limits for floating-points.
*
* Specializations are provided for `float` and `double` types. In order to use
* other floating point types, a new specialization for that type must be added.
*
* Member functions of the unspecialized class are disabled. Using them results
* in compile time errors.
* @{
*/
/**
* @tparam FP A floating point type.
*/
template <typename FP>
struct Limits {
static constexpr FP zero_tolerance() = delete;
};
template <>
struct Limits<float> {
/// @brief Returns the limit under which a float may be rounded down to zero.
static constexpr float zero_tolerance()
{
return 1e-6f;
}
};
template <>
struct Limits<double> {
/// @brief Returns the limit under which a double may be rounded down to zero.
static constexpr double zero_tolerance()
{
return 1e-12;
}
};
template <class FP, class DataHandler>
struct Limits<ad::internal::active_type<FP, DataHandler>> {
/// @brief Returns the limit under which an AD value may be rounded down to zero.
static constexpr FP zero_tolerance()
{
return Limits<FP>::zero_tolerance();
}
};
/** @} */
} // namespace mio
#endif