unsigned int isqrt(unsigned int x) {
  unsigned int r;
  unsigned int v;
  if (x < 2)
    return (x);
  for (v = x, r = 1; v; v >>= 2, r <<= 1)
    ;
  while (1) {
    v = ((x / r) + r) / 2;
    if ((v == r) || (v == r + 1))
      return (r);
    r = v;
  }
}


x86 assembly (clang, -O0):
<isqrt>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x8(%rbp)
cmpl   $0x2,-0x8(%rbp)
jae    111c <isqrt+0x1c>
mov    -0x8(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    1198 <isqrt+0x98>
mov    -0x8(%rbp),%eax
mov    %eax,-0x10(%rbp)
movl   $0x1,-0xc(%rbp)
cmpl   $0x0,-0x10(%rbp)
je     114f <isqrt+0x4f>
jmp    1138 <isqrt+0x38>
mov    -0x10(%rbp),%eax
shr    $0x2,%eax
mov    %eax,-0x10(%rbp)
mov    -0xc(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    1129 <isqrt+0x29>
jmp    1154 <isqrt+0x54>
mov    -0x8(%rbp),%eax
xor    %edx,%edx
divl   -0xc(%rbp)
add    -0xc(%rbp),%eax
shr    $0x1,%eax
mov    %eax,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
je     1182 <isqrt+0x82>
mov    -0x10(%rbp),%eax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jne    118d <isqrt+0x8d>
mov    -0xc(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    1198 <isqrt+0x98>
mov    -0x10(%rbp),%eax
mov    %eax,-0xc(%rbp)
jmp    1154 <isqrt+0x54>
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<isqrt>:
cmp    $0x2,%edi
jae    110a <isqrt+0xa>
mov    %edi,%ecx
mov    %ecx,%eax
ret
mov    $0x1,%eax
mov    %edi,%ecx
mov    %edi,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
shr    $0x2,%edx
add    %eax,%eax
cmp    $0x3,%ecx
mov    %edx,%ecx
ja     1120 <isqrt+0x20>
jmp    1137 <isqrt+0x37>
xchg   %ax,%ax
lea    0x1(%rcx),%edx
cmp    %edx,%eax
je     1147 <isqrt+0x47>
mov    %eax,%ecx
mov    %edi,%eax
xor    %edx,%edx
div    %ecx
add    %ecx,%eax
shr    %eax
cmp    %ecx,%eax
jne    1130 <isqrt+0x30>
mov    %ecx,%eax
ret

LLVM IR (-O0):
define i32 @isqrt(i32 %0) {
%2 = alloca i32
%3 = alloca i32
%4 = alloca i32
%5 = alloca i32
store i32 %0, i32* %3
%6 = load i32, i32* %3
%7 = icmp ult i32 %6, 2
br i1 %7, label %8, label %10
8: ; preds = %1
%9 = load i32, i32* %3
store i32 %9, i32* %2
br label %41
10: ; preds = %1
%11 = load i32, i32* %3
store i32 %11, i32* %5
store i32 1, i32* %4
br label %12
12: ; preds = %16, %10
%13 = load i32, i32* %5
%14 = icmp ne i32 %13, 0
br i1 %14, label %15, label %21
15: ; preds = %12
br label %16
16: ; preds = %15
%17 = load i32, i32* %5
%18 = lshr i32 %17, 2
store i32 %18, i32* %5
%19 = load i32, i32* %4
%20 = shl i32 %19, 1
store i32 %20, i32* %4
br label %12
21: ; preds = %12
br label %22
22: ; preds = %21, %39
%23 = load i32, i32* %3
%24 = load i32, i32* %4
%25 = udiv i32 %23, %24
%26 = load i32, i32* %4
%27 = add i32 %25, %26
%28 = udiv i32 %27, 2
store i32 %28, i32* %5
%29 = load i32, i32* %5
%30 = load i32, i32* %4
%31 = icmp eq i32 %29, %30
br i1 %31, label %37, label %32
32: ; preds = %22
%33 = load i32, i32* %5
%34 = load i32, i32* %4
%35 = add i32 %34, 1
%36 = icmp eq i32 %33, %35
br i1 %36, label %37, label %39
37: ; preds = %32, %22
%38 = load i32, i32* %4
store i32 %38, i32* %2
br label %41
39: ; preds = %32
%40 = load i32, i32* %5
store i32 %40, i32* %4
br label %22
41: ; preds = %37, %8
%42 = load i32, i32* %2
ret i32 %42
}

LLVM IR (-O2):
define i32 @isqrt(i32 %0) {
%2 = icmp ult i32 %0, 2
br i1 %2, label %18, label %3
3: ; preds = %1, %3
%4 = phi i32 [ %6, %3 ], [ %0, %1 ]
%5 = phi i32 [ %7, %3 ], [ 1, %1 ]
%6 = lshr i32 %4, 2
%7 = shl i32 %5, 1
%8 = icmp ult i32 %4, 4
br i1 %8, label %9, label %3
9: ; preds = %3, %9
%10 = phi i32 [ %13, %9 ], [ %7, %3 ]
%11 = udiv i32 %0, %10
%12 = add i32 %11, %10
%13 = lshr i32 %12, 1
%14 = icmp eq i32 %13, %10
%15 = add nuw i32 %10, 1
%16 = icmp eq i32 %13, %15
%17 = or i1 %14, %16
br i1 %17, label %18, label %9
18: ; preds = %9, %1
%19 = phi i32 [ %0, %1 ], [ %10, %9 ]
ret i32 %19
}