mirror of
https://github.com/fastfloat/fast_float.git
synced 2026-10-28 18:07:16 +08:00
Merge pull request #426 from sahvx655-wq/power-negative-shift
Avoid right-shifting a negative value in detail::power()
This commit is contained in:
commit
638a8acd6d
@ -61,9 +61,15 @@ namespace detail {
|
||||
* -ceil(p) + q
|
||||
* where
|
||||
* p = log(5**-q)/log(2) = -q * log(5)/log(2)
|
||||
*
|
||||
* The shift must floor towards negative infinity, as an arithmetic right
|
||||
* shift does, but right-shifting a negative value is implementation-defined
|
||||
* before C++20. Biasing the product by 2048 * 65536 keeps the shifted value
|
||||
* non-negative without changing the result: the bias is a multiple of 2^16
|
||||
* and is removed again after the shift.
|
||||
*/
|
||||
constexpr fastfloat_really_inline int32_t power(int32_t q) noexcept {
|
||||
return (((152170 + 65536) * q) >> 16) + 63;
|
||||
return (((152170 + 65536) * q + (2048 * 65536)) >> 16) - 2048 + 63;
|
||||
}
|
||||
} // namespace detail
|
||||
|
||||
|
||||
Loading…
x
Reference in New Issue
Block a user