1 #ifndef DataFormatsMathAPPROX_LOG_H
2 #define DataFormatsMathAPPROX_LOG_H
39 return y * (float(0x1.0671c4p0) +
y * (float(-0x7.27744p-4)));
45 return y * (float(0x1.013354p0) +
y * (-float(0x8.33006p-4) +
y * float(0x4.0d16cp-4)));
52 (float(0xf.ff5bap-4) +
y * (-float(0x8.13e5ep-4) +
y * (float(0x5.826ep-4) +
y * (-float(0x2.e87fb8p-4)))));
58 return y * (float(0xf.ff652p-4) +
59 y * (-float(0x8.0048ap-4) +
60 y * (float(0x5.72782p-4) +
y * (-float(0x4.20904p-4) +
y * float(0x2.1d7fd8p-4)))));
66 return y * (float(0xf.fff14p-4) +
67 y * (-float(0x7.ff4bfp-4) +
68 y * (float(0x5.582f6p-4) +
69 y * (-float(0x4.1dcf2p-4) +
y * (float(0x3.3863f8p-4) +
y * (-float(0x1.9288d4p-4)))))));
75 return y * (float(0x1.000034p0) +
76 y * (-float(0x7.ffe57p-4) +
77 y * (float(0x5.5422ep-4) +
78 y * (-float(0x4.037a6p-4) +
79 y * (float(0x3.541c88p-4) +
y * (-float(0x2.af842p-4) +
y * float(0x1.48b3d8p-4)))))));
86 (float(0x1.00000cp0) +
87 y * (float(-0x8.0003p-4) +
88 y * (float(0x5.55087p-4) +
89 y * (float(-0x3.fedcep-4) +
90 y * (float(0x3.3a1dap-4) +
91 y * (float(-0x2.cb55fp-4) +
y * (float(0x2.38831p-4) +
y * (float(-0xf.e87cap-8)))))))));
96 using namespace approx_math;
102 int e = (((xx.i32) >> 23) & 0xFF) - 127;
103 m.i32 = (xx.i32 & 0x007FFFFF) | 0x3F800000;
105 int adjust = (xx.i32 >> 22) & 1;
106 m.i32 -= adjust << 23;
110 float y =
m.f - 1.0f;
114 float p = approx_logf_P<DEGREE>(y);
116 constexpr
float Log2 = 0xb.17218p-4;
117 return float(
e) * Log2 +
p;
120 #ifndef NO_APPROX_MATH
121 template <
int DEGREE>
123 return unsafe_logf_impl<DEGREE>(
x);
126 template <
int DEGREE>
128 using namespace approx_math;
130 constexpr
float MAXNUMF = 3.4028234663852885981170418348451692544e38f;
133 float res = unsafe_logf<DEGREE>(x);
135 return (x > 0) ? res : std::numeric_limits<float>::quiet_NaN();
139 template <
int DEGREE>
144 template <
int DEGREE>
149 #endif // NO_APPROX_MATH
static std::vector< std::string > checklist log
constexpr float approx_logf_P< 2 >(float y)
constexpr float approx_logf_P< 7 >(float y)
constexpr float unsafe_logf(float x)
constexpr float approx_logf_P< 4 >(float y)
constexpr float approx_logf_P< 5 >(float y)
constexpr float approx_logf_P(float p)
constexpr float unsafe_logf_impl(float x)
uint16_t const *__restrict__ x
constexpr float approx_logf_P< 8 >(float y)
constexpr float approx_logf_P< 6 >(float y)
constexpr float approx_logf_P< 3 >(float y)
constexpr float approx_logf(float x)