int isqrt(int x) {
  int s, t;
  if (x == 0)
    return 0;
  s = 1;
  t = x;
  while (s < t) {
    s <<= 1;
    t >>= 1;
  }
  do {
    t = s;
    s = (x / s + s) >> 1;
  } while (s < t);
  return (t);
}


x86 assembly (clang, -O0):
<isqrt>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x8(%rbp)
cmpl   $0x0,-0x8(%rbp)
jne    111d <isqrt+0x1d>
movl   $0x0,-0x4(%rbp)
jmp    117a <isqrt+0x7a>
movl   $0x1,-0xc(%rbp)
mov    -0x8(%rbp),%eax
mov    %eax,-0x10(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    114d <isqrt+0x4d>
mov    -0xc(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0xc(%rbp)
mov    -0x10(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    112a <isqrt+0x2a>
jmp    1152 <isqrt+0x52>
mov    -0xc(%rbp),%eax
mov    %eax,-0x10(%rbp)
mov    -0x8(%rbp),%eax
cltd
idivl  -0xc(%rbp)
add    -0xc(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0xc(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x10(%rbp),%eax
jl     1152 <isqrt+0x52>
mov    -0x10(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<isqrt>:
test   %edi,%edi
je     1142 <isqrt+0x42>
mov    $0x1,%eax
cmp    $0x2,%edi
jl     1130 <isqrt+0x30>
mov    $0x1,%eax
mov    %edi,%ecx
cs nopw 0x0(%rax,%rax,1)
nop
add    %eax,%eax
sar    %ecx
cmp    %ecx,%eax
jl     1120 <isqrt+0x20>
nopl   0x0(%rax,%rax,1)
mov    %eax,%ecx
mov    %edi,%eax
cltd
idiv   %ecx
add    %ecx,%eax
sar    %eax
cmp    %ecx,%eax
jl     1130 <isqrt+0x30>
mov    %ecx,%eax
ret
xor    %eax,%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 eq i32 %6, 0
br i1 %7, label %8, label %9
8: ; preds = %1
store i32 0, i32* %2
br label %35
9: ; preds = %1
store i32 1, i32* %4
%10 = load i32, i32* %3
store i32 %10, i32* %5
br label %11
11: ; preds = %15, %9
%12 = load i32, i32* %4
%13 = load i32, i32* %5
%14 = icmp slt i32 %12, %13
br i1 %14, label %15, label %20
15: ; preds = %11
%16 = load i32, i32* %4
%17 = shl i32 %16, 1
store i32 %17, i32* %4
%18 = load i32, i32* %5
%19 = ashr i32 %18, 1
store i32 %19, i32* %5
br label %11
20: ; preds = %11
br label %21
21: ; preds = %29, %20
%22 = load i32, i32* %4
store i32 %22, i32* %5
%23 = load i32, i32* %3
%24 = load i32, i32* %4
%25 = sdiv i32 %23, %24
%26 = load i32, i32* %4
%27 = add nsw i32 %25, %26
%28 = ashr i32 %27, 1
store i32 %28, i32* %4
br label %29
29: ; preds = %21
%30 = load i32, i32* %4
%31 = load i32, i32* %5
%32 = icmp slt i32 %30, %31
br i1 %32, label %21, label %33
33: ; preds = %29
%34 = load i32, i32* %5
store i32 %34, i32* %2
br label %35
35: ; preds = %33, %8
%36 = load i32, i32* %2
ret i32 %36
}

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