GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-04-19 08:09:36 +08:00
Code Issues Packages Projects Releases Wiki Activity

Tweak bound for pow to account for floating-point types.

Browse Source
This commit is contained in:
Antonio Sánchez 2022-09-08 17:40:45 +00:00
parent f9dfda28ab
commit 3c37dd2a1d
1 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
test/array_cwise.cpp
View File

@ -136,8 +136,11 @@ Scalar calc_overflow_threshold(const ScalarExponent exponent) {
return NumTraits<Scalar>::highest(); return NumTraits<Scalar>::highest();
else { else {
// base^e <= highest ==> base <= 2^(log2(highest)/e) // base^e <= highest ==> base <= 2^(log2(highest)/e)
// For floating-point types, consider the bound for integer values that can be reproduced exactly = 2 ^ digits
double highest_bits = numext::mini(static_cast<double>(NumTraits<Scalar>::digits()),
log2(NumTraits<Scalar>::highest()));
return static_cast<Scalar>( return static_cast<Scalar>(
numext::floor(exp2(log2(NumTraits<Scalar>::highest()) / static_cast<double>(exponent)))); numext::floor(exp2(highest_bits / static_cast<double>(exponent))));
} }
} }