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
|
/* Copyright (C) 2008 Tobi Vollebregt */
/* Conditionally include streflop depending on STREFLOP_* #defines:
If one of those is present, #include "streflop.h", otherwise #include <math.h>
When faced with ambiguous call errors with e.g. fabs, use math::function.
Add it to math namespace if it doesn't exist there yet. */
#ifndef STREFLOP_COND_H
#define STREFLOP_COND_H
#if (defined(STREFLOP_X87) || defined(STREFLOP_SSE) || defined(STREFLOP_SOFT)) && (!defined(NOT_USING_STREFLOP))
#include "streflop.h"
namespace math {
using namespace streflop;
}
#else
#include <cmath>
#ifdef __APPLE__
// macosx's cmath doesn't include c++11's std::hypot yet (tested 2013)
namespace std {
template<typename T> T hypot(T x, T y);
}
#endif
namespace math {
using std::fabs;
// We are using fastmath::sqrt_sse instead!
// using std::sqrt;
using std::sin;
using std::cos;
using std::sinh;
using std::cosh;
using std::tan;
using std::tanh;
using std::asin;
using std::acos;
using std::atan;
using std::atan2;
using std::ceil;
using std::floor;
using std::fmod;
using std::hypot;
using std::pow;
using std::log;
using std::log10;
using std::exp;
using std::frexp;
using std::ldexp;
// the following are C99 functions -> not supported by VS C
#if !defined(_MSC_VER) || _MSC_VER < 1500
using std::isnan;
using std::isinf;
using std::isfinite;
#elif __cplusplus
template<typename T> inline bool isnan(T value) {
return value != value;
}
// requires include <limits>
template<typename T> inline bool isinf(T value) {
return std::numeric_limits<T>::has_infinity && value == std::numeric_limits<T>::infinity();
}
// requires include <limits>
template<typename T> inline bool isfinite(T value) {
return !isinf<T>(value);
}
#endif
}
#ifdef __APPLE__
// see above
#include <algorithm>
namespace std {
template<typename T>
T hypot(T x, T y) {
x = std::abs(x);
y = std::abs(y);
auto t = std::min(x,y);
x = std::max(x,y);
t = t / x;
return x * sqrtf(1.f + t*t);
}
}
#endif
#endif
#endif // STREFLOP_COND_H
|
