Finite Volume Solver  prototype
A framework to build finite volume solvers for the AG Klein at the Freie Universität Berlin.
algorithm.hpp
Go to the documentation of this file.
1 // Copyright (c) 2018 Maikel Nadolski
2 //
3 // Permission is hereby granted, free of charge, to any person obtaining a copy
4 // of this software and associated documentation files (the "Software"), to deal
5 // in the Software without restriction, including without limitation the rights
6 // to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
7 // copies of the Software, and to permit persons to whom the Software is
8 // furnished to do so, subject to the following conditions:
9 //
10 // The above copyright notice and this permission notice shall be included in
11 // all copies or substantial portions of the Software.
12 //
13 // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
14 // IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
15 // FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
16 // AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
17 // LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
18 // OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
19 // SOFTWARE.
20 
21 #ifndef FUB_CORE_ALGORITHM_HPP
22 #define FUB_CORE_ALGORITHM_HPP
23 
24 #include <cassert>
25 #include <cstdint>
26 
27 namespace fub {
28 
29 template <class T, class Compare>
30 constexpr const T& clamp(const T& v, const T& lo, const T& hi, Compare comp) {
31  return assert(!comp(hi, lo)), comp(v, lo) ? lo : comp(hi, v) ? hi : v;
32 }
33 
34 template <class T>
35 constexpr const T& clamp(const T& v, const T& lo, const T& hi) {
36  return clamp(v, lo, hi, std::less<>());
37 }
38 
39 constexpr std::ptrdiff_t ipow(int base, int exponent) {
40  std::ptrdiff_t prod{1};
41  while (exponent > 0) {
42  prod *= base;
43  exponent -= 1;
44  }
45  return prod;
46 }
47 
48 } // namespace fub
49 
50 #endif
The fub namespace.
Definition: AnyBoundaryCondition.hpp:31
constexpr std::ptrdiff_t ipow(int base, int exponent)
Definition: algorithm.hpp:39
constexpr const T & clamp(const T &v, const T &lo, const T &hi, Compare comp)
Definition: algorithm.hpp:30