我怎样才能将分裂强度降低2 ^ n + 1?

时间:2010-10-25 17:32:19

标签: c division integer-division strength-reduction

我需要在代码的热路径中执行一些整数除法。我已经通过分析和循环计数确定了整数除法对我造成的损失。我希望我能做些什么来强化将分裂降低到更便宜的东西。

在这条路径中,我除以2 ^ n + 1,其中n是可变的。基本上我想优化这个函数来删除除法运算符:

unsigned long compute(unsigned long a, unsigned int n)
{
    return a / ((1 << n) + 1);
}

如果我除以2 ^ n,我只需用右移n替换div。如果我用常数除法,我会让编译器强度减少那个特定的除法,可能会把它变成乘法和一些变化。

是否存在适用于2 ^ n + 1的类似优化?

编辑:这里可以是任意64位整数。 n只需要10和25之间的几个值。我当然可以为每个n预先计算一些值,但不能用于。

2 个答案:

答案 0 :(得分:13)

由于你只能移动int这么多地方,你可以把所有这些情况用一个常数选择几个部分中的一个:

unsigned long compute(unsigned long a, unsigned int n)
{
    // assuming a 32-bit architecture (making this work for 64-bits 
    // is left as an exercise for the reader):
    switch (n) {
        case  0: return a / ((1 << 0) + 1);
        case  1: return a / ((1 << 1) + 1);
        case  2: return a / ((1 << 2) + 1);

            // cases 3 through 30...

        case 31: return a / ((1 << 31) + 1);
    }
}

所以现在每个除法都是常数,编译器通常会减少到一系列乘法/移位/加法指令(如上所述)。有关deatils,请参阅Does a c/c++ compiler optimize constant divisions by power-of-two value into shifts?

答案 1 :(得分:8)

你可以用一个常数替换整数除法,用一个幻数和一个移位乘以(模数词大小)。

可以为已知常数预先计算幻数。

因为n不能采用多个值,例如0..31为所有n预先计算这些幻数并“存储”在32个元素的表中是“容易的”。

Javascript Page for calculating the magic numbers

如果除数在编译时是常数,那么好的编译器可以计算幻数并通过乘法和移位替换整数除法。根据围绕性能关键代码构建其余代码的方式,您可以使用宏或内联技巧来展开n的所有可能值,并让编译器完成查找幻数的工作(类似于答案使用开关,但我会在常量区域放置更多代码,否则它可能是不值得的交易 - 分支也会花费你的性能

详细描述以及计算幻数的代码可以在Henry S. Warren,Jr。(强烈推荐必须有书!)pp.180ff的书“Hackers Delight”中获得资金。

链接到相关章节的Google图书:

Chapter 10-9 Unsigned Division by Divisors >= 1