#define iter1(N)      \
    try = root + (1 << (N));   \
    if (n >= try << (N))   \
    {            \
   n -= try << (N);   \
        root |= 2 << (N);   \
    }

u32 sqrt_u32 (u32 n)
{
    u32 root = 0, try;

    iter1 (15);
    iter1 (14);
    iter1 (13);
    iter1 (12);
    iter1 (11);
    iter1 (10);
    iter1 ( 9);
    iter1 ( 8);
    iter1 ( 7);
    iter1 ( 6);
    iter1 ( 5);
    iter1 ( 4);
    iter1 ( 3);
    iter1 ( 2);
    iter1 ( 1);
    iter1 ( 0);

    return root >> 1;
}