int isqrt(i)
int i;
{
  register int a, b;
  if (i <= 0)
    return (0);
  for (a = i, b = 1; b <<= 1, a >>= 2;)
    ;
  while ((a = i / b) < b)
    b = (b + a) >> 1;
  return (b);
}


x86 assembly (clang, -O0):
<isqrt>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x8(%rbp)
cmpl   $0x0,-0x8(%rbp)
jg     111d <isqrt+0x1d>
movl   $0x0,-0x4(%rbp)
jmp    1179 <isqrt+0x79>
mov    -0x8(%rbp),%eax
mov    %eax,-0xc(%rbp)
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x10(%rbp)
mov    -0xc(%rbp),%eax
sar    $0x2,%eax
mov    %eax,-0xc(%rbp)
cmp    $0x0,%eax
je     114a <isqrt+0x4a>
jmp    112a <isqrt+0x2a>
jmp    114f <isqrt+0x4f>
mov    -0x8(%rbp),%eax
cltd
idivl  -0x10(%rbp)
mov    %eax,-0xc(%rbp)
cmp    -0x10(%rbp),%eax
jge    1173 <isqrt+0x73>
mov    -0x10(%rbp),%eax
add    -0xc(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    114f <isqrt+0x4f>
mov    -0x10(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<isqrt>:
test   %edi,%edi
jle    1130 <isqrt+0x30>
mov    $0x1,%ecx
mov    %edi,%eax
mov    %edi,%edx
nopl   (%rax)
add    %ecx,%ecx
sar    $0x2,%edx
cmp    $0x3,%eax
mov    %edx,%eax
ja     1110 <isqrt+0x10>
jmp    1124 <isqrt+0x24>
xchg   %ax,%ax
add    %eax,%ecx
sar    %ecx
mov    %edi,%eax
cltd
idiv   %ecx
cmp    %ecx,%eax
jl     1120 <isqrt+0x20>
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 sle i32 %6, 0
br i1 %7, label %8, label %9
8: ; preds = %1
store i32 0, i32* %2
br label %32
9: ; preds = %1
%10 = load i32, i32* %3
store i32 %10, i32* %4
store i32 1, i32* %5
br label %11
11: ; preds = %17, %9
%12 = load i32, i32* %5
%13 = shl i32 %12, 1
store i32 %13, i32* %5
%14 = load i32, i32* %4
%15 = ashr i32 %14, 2
store i32 %15, i32* %4
%16 = icmp ne i32 %15, 0
br i1 %16, label %17, label %18
17: ; preds = %11
br label %11
18: ; preds = %11
br label %19
19: ; preds = %25, %18
%20 = load i32, i32* %3
%21 = load i32, i32* %5
%22 = sdiv i32 %20, %21
store i32 %22, i32* %4
%23 = load i32, i32* %5
%24 = icmp slt i32 %22, %23
br i1 %24, label %25, label %30
25: ; preds = %19
%26 = load i32, i32* %5
%27 = load i32, i32* %4
%28 = add nsw i32 %26, %27
%29 = ashr i32 %28, 1
store i32 %29, i32* %5
br label %19
30: ; preds = %19
%31 = load i32, i32* %5
store i32 %31, i32* %2
br label %32
32: ; preds = %30, %8
%33 = load i32, i32* %2
ret i32 %33
}

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