вычисление квадратного корня для реализации функции с фиксированной точкой

Я пытаюсь найти квадратный корень из фиксированной точки, и я использовал следующий расчет, чтобы найти аппроксимацию квадратного корня с помощью целочисленного алгоритма. Алгоритм описан в Википедии:
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

Большое спасибо