int mySqrt(int x) {
  unsigned short root = 0, i = 0;
  unsigned int rem = 0, divisor = 0, uX = x;
  double fX = uX;
  int lg2_n = (*((unsigned long long *)&fX) >> 52) - 1023;
  uX <<= (15 - (lg2_n >> 1)) << 1;
  i = (lg2_n >> 1) + 2;
  while (--i) {
    root <<= 1;
    rem = ((rem << 2) | uX >> 30);
    uX <<= 2;
    divisor = root << 1 | 1;
    if (divisor <= rem) {
      rem -= divisor;
      ++root;
    }
  }
  return (int)root;
}


x86 assembly (clang, -O0):
<mySqrt>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
movw   $0x0,-0x6(%rbp)
movw   $0x0,-0x8(%rbp)
movl   $0x0,-0xc(%rbp)
movl   $0x0,-0x10(%rbp)
mov    -0x4(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cvtsi2sd %rax,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x20(%rbp),%rax
shr    $0x34,%rax
sub    $0x3ff,%rax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
sar    $0x1,%eax
mov    $0xf,%ecx
sub    %eax,%ecx
shl    $0x1,%ecx
mov    -0x14(%rbp),%eax
shl    %cl,%eax
mov    %eax,-0x14(%rbp)
mov    -0x24(%rbp),%eax
sar    $0x1,%eax
add    $0x2,%eax
mov    %ax,-0x8(%rbp)
mov    -0x8(%rbp),%ax
add    $0xffff,%ax
mov    %ax,-0x8(%rbp)
cmp    $0x0,%ax
je     11da <mySqrt+0xda>
movzwl -0x6(%rbp),%eax
shl    $0x1,%eax
mov    %ax,-0x6(%rbp)
mov    -0xc(%rbp),%eax
shl    $0x2,%eax
mov    -0x14(%rbp),%ecx
shr    $0x1e,%ecx
or     %ecx,%eax
mov    %eax,-0xc(%rbp)
mov    -0x14(%rbp),%eax
shl    $0x2,%eax
mov    %eax,-0x14(%rbp)
movzwl -0x6(%rbp),%eax
shl    $0x1,%eax
or     $0x1,%eax
mov    %eax,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
ja     11d5 <mySqrt+0xd5>
mov    -0x10(%rbp),%ecx
mov    -0xc(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0xc(%rbp)
mov    -0x6(%rbp),%ax
add    $0x1,%ax
mov    %ax,-0x6(%rbp)
jmp    116a <mySqrt+0x6a>
movzwl -0x6(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<mySqrt>:
mov    %edi,%eax
cvtsi2sd %rax,%xmm0
movq   %xmm0,%rsi
shr    $0x34,%rsi
add    $0xfffffc01,%esi
mov    %esi,%r9d
shr    %r9d
xor    %eax,%eax
inc    %r9w
je     1170 <mySqrt+0x70>
and    $0xfe,%sil
mov    $0x1e,%cl
sub    %sil,%cl
shl    %cl,%edi
xor    %r8d,%r8d
xor    %eax,%eax
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
shld   $0x2,%edi,%esi
shl    $0x2,%edi
mov    %eax,%ecx
and    $0x7fff,%ecx
lea    0x1(,%rcx,4),%ecx
xor    %edx,%edx
cmp    %ecx,%esi
cmovb  %r8d,%ecx
setae  %dl
sub    %ecx,%esi
lea    (%rdx,%rax,2),%eax
add    $0xffff,%r9w
jne    1140 <mySqrt+0x40>
movzwl %ax,%eax
ret

LLVM IR (-O0):
define i32 @mySqrt(i32 %0) {
%2 = alloca i32
%3 = alloca i16
%4 = alloca i16
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca double
%9 = alloca i32
store i32 %0, i32* %2
store i16 0, i16* %3
store i16 0, i16* %4
store i32 0, i32* %5
store i32 0, i32* %6
%10 = load i32, i32* %2
store i32 %10, i32* %7
%11 = load i32, i32* %7
%12 = uitofp i32 %11 to double
store double %12, double* %8
%13 = bitcast double* %8 to i64*
%14 = load i64, i64* %13
%15 = lshr i64 %14, 52
%16 = sub i64 %15, 1023
%17 = trunc i64 %16 to i32
store i32 %17, i32* %9
%18 = load i32, i32* %9
%19 = ashr i32 %18, 1
%20 = sub nsw i32 15, %19
%21 = shl i32 %20, 1
%22 = load i32, i32* %7
%23 = shl i32 %22, %21
store i32 %23, i32* %7
%24 = load i32, i32* %9
%25 = ashr i32 %24, 1
%26 = add nsw i32 %25, 2
%27 = trunc i32 %26 to i16
store i16 %27, i16* %4
br label %28
28: ; preds = %57, %1
%29 = load i16, i16* %4
%30 = add i16 %29, -1
store i16 %30, i16* %4
%31 = icmp ne i16 %30, 0
br i1 %31, label %32, label %58
32: ; preds = %28
%33 = load i16, i16* %3
%34 = zext i16 %33 to i32
%35 = shl i32 %34, 1
%36 = trunc i32 %35 to i16
store i16 %36, i16* %3
%37 = load i32, i32* %5
%38 = shl i32 %37, 2
%39 = load i32, i32* %7
%40 = lshr i32 %39, 30
%41 = or i32 %38, %40
store i32 %41, i32* %5
%42 = load i32, i32* %7
%43 = shl i32 %42, 2
store i32 %43, i32* %7
%44 = load i16, i16* %3
%45 = zext i16 %44 to i32
%46 = shl i32 %45, 1
%47 = or i32 %46, 1
store i32 %47, i32* %6
%48 = load i32, i32* %6
%49 = load i32, i32* %5
%50 = icmp ule i32 %48, %49
br i1 %50, label %51, label %57
51: ; preds = %32
%52 = load i32, i32* %6
%53 = load i32, i32* %5
%54 = sub i32 %53, %52
store i32 %54, i32* %5
%55 = load i16, i16* %3
%56 = add i16 %55, 1
store i16 %56, i16* %3
br label %57
57: ; preds = %51, %32
br label %28
58: ; preds = %28
%59 = load i16, i16* %3
%60 = zext i16 %59 to i32
ret i32 %60
}

LLVM IR (-O2):
define i32 @mySqrt(i32 %0) {
%2 = uitofp i32 %0 to double
%3 = bitcast double %2 to i64
%4 = lshr i64 %3, 52
%5 = trunc i64 %4 to i32
%6 = add nsw i32 %5, -1023
%7 = ashr i32 %6, 1
%8 = trunc i32 %7 to i16
%9 = add nsw i16 %8, 1
%10 = icmp eq i16 %9, 0
br i1 %10, label %36, label %11
11: ; preds = %1
%12 = and i32 %6, -2
%13 = sub nsw i32 30, %12
%14 = shl i32 %0, %13
br label %15
15: ; preds = %11, %15
%16 = phi i16 [ %32, %15 ], [ %9, %11 ]
%17 = phi i16 [ %31, %15 ], [ 0, %11 ]
%18 = phi i32 [ %22, %15 ], [ %14, %11 ]
%19 = phi i32 [ %28, %15 ], [ 0, %11 ]
%20 = shl i16 %17, 1
%21 = tail call i32 @llvm.fshl.i32(i32 %19, i32 %18, i32 2)
%22 = shl i32 %18, 2
%23 = zext i16 %20 to i32
%24 = shl nuw nsw i32 %23, 1
%25 = or i32 %24, 1
%26 = icmp ult i32 %21, %25
%27 = select i1 %26, i32 0, i32 %25
%28 = sub i32 %21, %27
%29 = xor i1 %26, true
%30 = zext i1 %29 to i16
%31 = or i16 %20, %30
%32 = add i16 %16, -1
%33 = icmp eq i16 %32, 0
br i1 %33, label %34, label %15
34: ; preds = %15
%35 = zext i16 %31 to i32
br label %36
36: ; preds = %34, %1
%37 = phi i32 [ 0, %1 ], [ %35, %34 ]
ret i32 %37
}