From: Pavan Nikhilesh <pbhagavat...@marvell.com> Fix off by one error in 64bit reciprocal division when divisor is 32bit.
Fixes: 6d45659eacb8 ("eal: add u64-bit variant for reciprocal divide") Cc: sta...@dpdk.org Signed-off-by: Pavan Nikhilesh <pbhagavat...@marvell.com> --- Example: Division failed, 17358247066007716387/244 = expected 71140356827900476 result 71140356827900477 Division failed, 17541123788887206374/41475 = expected 422932460250444 result 422932460250445 lib/librte_eal/common/rte_reciprocal.c | 7 +++++-- 1 file changed, 5 insertions(+), 2 deletions(-) diff --git a/lib/librte_eal/common/rte_reciprocal.c b/lib/librte_eal/common/rte_reciprocal.c index f017d0c28..1c6d10e73 100644 --- a/lib/librte_eal/common/rte_reciprocal.c +++ b/lib/librte_eal/common/rte_reciprocal.c @@ -133,12 +133,15 @@ rte_reciprocal_value_u64(uint64_t d) { struct rte_reciprocal_u64 R; uint64_t m; + uint64_t r; int l; l = 63 - __builtin_clzll(d); - m = divide_128_div_64_to_64((1ULL << l), 0, d, NULL) << 1; - m = (1ULL << l) - d ? m + 2 : 1; + m = divide_128_div_64_to_64((1ULL << l), 0, d, &r) << 1; + if (r << 1 < r || r << 1 >= d) + m++; + m = (1ULL << l) - d ? m + 1 : 1; R.m = m; R.sh1 = l > 1 ? 1 : l; -- 2.23.0