计算平方根以实现定点函数

Question

我试图找到一个固定点的平方根，我使用下面的计算来找到使用整数算法的平方根的近似值。该算法在维基百科中描述：  http://en.wikipedia.org/wiki/Methods_of_computing_square_roots

uint32_t SquareRoot(uint32_t a_nInput)
{
    uint32_t op  = a_nInput;
    uint32_t res = 0;
    uint32_t one = 1uL << 30; // The second-to-top bit is set: use 1u << 14 for uint16_t type; use 1uL<<30 for uint32_t type


    // "one" starts at the highest power of four <= than the argument.
    while (one > op)
    {
        one >>= 2;
    }

    while (one != 0)
    {
        if (op >= res + one)
        {
            op = op - (res + one);
            res = res +  2 * one;
        }
        res >>= 1;
        one >>= 2;
    }
    return res;
}

但是我无法跟踪代码中发生的内容// "one" starts at the highest power of four <= than the argument.的确切意味着什么。有人可以提示我在代码中发生什么来计算参数的平方根a_nInput

非常感谢

Answer 1

请注意one的初始化方式。

uint32_t one = 1uL << 30;

那2 ³⁰ 或1073741824。这也是4 ¹⁵ 。

这一行：

    one >>= 2;

相当于

    one = one / 4;

所以发生了什么的伪代码是：

• one = 4 ¹⁵

• 如果one超过a_nInput

  • one = 4 ¹⁴

• 如果one仍然超过a_nInput

  • one = 4 ¹³

• （依此类推......）

最终，one 不将超过a_nInput。

// "one" starts at the highest power of four less than or equal to a_nInput