double sqrt(double x) {
  double z;
  int sign = (int)0x80000000;
  unsigned r, t1, s1, ix1, q1;
  int ix0, s0, q, m, t, i;
  ix0 = *(1 + (int *)&x);
  ix1 = *(int *)&x;
  if ((ix0 & 0x7ff00000) == 0x7ff00000) {
    return x * x + x;
  }
  if (ix0 <= 0) {
    if (((ix0 & (~sign)) | ix1) == 0) {
      return x;
    } else if (ix0 < 0) {
      return (x - x) / (x - x);
    }
  }
  m = (ix0 >> 20);
  if (m == 0) {
    while (ix0 == 0) {
      m -= 21;
      ix0 |= (ix1 >> 11);
      ix1 <<= 21;
    }
    for (i = 0; (ix0 & 0x00100000) == 0; i++) {
      ix0 <<= 1;
    }
    m -= i - 1;
    ix0 |= (ix1 >> (32 - i));
    ix1 <<= i;
  }
  m -= 1023;
  ix0 = (ix0 & 0x000fffff) | 0x00100000;
  if (m & 1) {
    ix0 += ix0 + ((ix1 & sign) >> 31);
    ix1 += ix1;
  }
  m >>= 1;
  ix0 += ix0 + ((ix1 & sign) >> 31);
  ix1 += ix1;
  q = q1 = s0 = s1 = 0;
  r = 0x00200000;
  while (r != 0) {
    t = s0 + r;
    if (t <= ix0) {
      s0 = t + r;
      ix0 -= t;
      q += r;
    }
    ix0 += ix0 + ((ix1 & sign) >> 31);
    ix1 += ix1;
    r >>= 1;
  }
  r = sign;
  while (r != 0) {
    t1 = s1 + r;
    t = s0;
    if ((t < ix0) || ((t == ix0) && (t1 <= ix1))) {
      s1 = t1 + r;
      if (((t1 & sign) == sign) && (s1 & sign) == 0) {
        s0 += 1;
      }
      ix0 -= t;
      if (ix1 < t1) {
        ix0 -= 1;
      }
      ix1 -= t1;
      q1 += r;
    }
    ix0 += ix0 + ((ix1 & sign) >> 31);
    ix1 += ix1;
    r >>= 1;
  }
  if ((ix0 | ix1) != 0) {
    z = 1.0 - 1.0e-300;
    if (z >= 1.0) {
      z = 1.0 + 1.0e-300;
      if (q1 == (unsigned)0xffffffff) {
        q1 = 0;
        q += 1;
      } else if (z > 1.0) {
        if (q1 == (unsigned)0xfffffffe) {
          q += 1;
        }
        q1 += 2;
      } else {
        q1 += (q1 & 1);
      }
    }
  }
  ix0 = (q >> 1) + 0x3fe00000;
  ix1 = q1 >> 1;
  if ((q & 1) == 1) {
    ix1 |= sign;
  }
  ix0 += (m << 20);
  *(1 + (int *)&z) = ix0;
  *(int *)&z = ix1;
  return z;
}


x86 assembly (clang, -O0):
<sqrt>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
movl   $0x80000000,-0x1c(%rbp)
mov    -0xc(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x10(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x34(%rbp),%eax
and    $0x7ff00000,%eax
cmp    $0x7ff00000,%eax
jne    1150 <sqrt+0x50>
movsd  -0x10(%rbp),%xmm0
movsd  -0x10(%rbp),%xmm2
movsd  -0x10(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1522 <sqrt+0x422>
cmpl   $0x0,-0x34(%rbp)
jg     11b6 <sqrt+0xb6>
mov    -0x34(%rbp),%eax
mov    -0x1c(%rbp),%ecx
xor    $0xffffffff,%ecx
and    %ecx,%eax
or     -0x2c(%rbp),%eax
cmp    $0x0,%eax
jne    1180 <sqrt+0x80>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1522 <sqrt+0x422>
cmpl   $0x0,-0x34(%rbp)
jge    11ac <sqrt+0xac>
movsd  -0x10(%rbp),%xmm0
subsd  -0x10(%rbp),%xmm0
movsd  -0x10(%rbp),%xmm1
subsd  -0x10(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1522 <sqrt+0x422>
jmp    11b1 <sqrt+0xb1>
jmp    11b6 <sqrt+0xb6>
mov    -0x34(%rbp),%eax
sar    $0x14,%eax
mov    %eax,-0x40(%rbp)
cmpl   $0x0,-0x40(%rbp)
jne    1256 <sqrt+0x156>
jmp    11ce <sqrt+0xce>
cmpl   $0x0,-0x34(%rbp)
jne    11fb <sqrt+0xfb>
mov    -0x40(%rbp),%eax
sub    $0x15,%eax
mov    %eax,-0x40(%rbp)
mov    -0x2c(%rbp),%eax
shr    $0xb,%eax
or     -0x34(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x2c(%rbp),%eax
shl    $0x15,%eax
mov    %eax,-0x2c(%rbp)
jmp    11ce <sqrt+0xce>
movl   $0x0,-0x48(%rbp)
mov    -0x34(%rbp),%eax
and    $0x100000,%eax
cmp    $0x0,%eax
jne    122a <sqrt+0x12a>
mov    -0x34(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1202 <sqrt+0x102>
mov    -0x48(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x40(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x40(%rbp)
mov    -0x2c(%rbp),%eax
mov    $0x20,%ecx
sub    -0x48(%rbp),%ecx
shr    %cl,%eax
or     -0x34(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x48(%rbp),%ecx
mov    -0x2c(%rbp),%eax
shl    %cl,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x40(%rbp),%eax
sub    $0x3ff,%eax
mov    %eax,-0x40(%rbp)
mov    -0x34(%rbp),%eax
and    $0xfffff,%eax
or     $0x100000,%eax
mov    %eax,-0x34(%rbp)
mov    -0x40(%rbp),%eax
and    $0x1,%eax
cmp    $0x0,%eax
je     129d <sqrt+0x19d>
mov    -0x34(%rbp),%eax
mov    -0x2c(%rbp),%ecx
and    -0x1c(%rbp),%ecx
shr    $0x1f,%ecx
add    %ecx,%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x2c(%rbp),%eax
add    -0x2c(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x40(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x40(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x2c(%rbp),%ecx
and    -0x1c(%rbp),%ecx
shr    $0x1f,%ecx
add    %ecx,%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x2c(%rbp),%eax
add    -0x2c(%rbp),%eax
mov    %eax,-0x2c(%rbp)
movl   $0x0,-0x28(%rbp)
movl   $0x0,-0x38(%rbp)
movl   $0x0,-0x30(%rbp)
movl   $0x0,-0x3c(%rbp)
movl   $0x200000,-0x20(%rbp)
cmpl   $0x0,-0x20(%rbp)
je     134d <sqrt+0x24d>
mov    -0x38(%rbp),%eax
add    -0x20(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x34(%rbp),%eax
jg     1322 <sqrt+0x222>
mov    -0x44(%rbp),%eax
add    -0x20(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x44(%rbp),%ecx
mov    -0x34(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x34(%rbp)
mov    -0x20(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x2c(%rbp),%ecx
and    -0x1c(%rbp),%ecx
shr    $0x1f,%ecx
add    %ecx,%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x2c(%rbp),%eax
add    -0x2c(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x20(%rbp),%eax
shr    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    12e6 <sqrt+0x1e6>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x20(%rbp)
cmpl   $0x0,-0x20(%rbp)
je     141f <sqrt+0x31f>
mov    -0x28(%rbp),%eax
add    -0x20(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x38(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x34(%rbp),%eax
jl     1390 <sqrt+0x290>
mov    -0x44(%rbp),%eax
cmp    -0x34(%rbp),%eax
jne    13f4 <sqrt+0x2f4>
mov    -0x24(%rbp),%eax
cmp    -0x2c(%rbp),%eax
ja     13f4 <sqrt+0x2f4>
mov    -0x24(%rbp),%eax
add    -0x20(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x24(%rbp),%eax
and    -0x1c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jne    13c0 <sqrt+0x2c0>
mov    -0x28(%rbp),%eax
and    -0x1c(%rbp),%eax
cmp    $0x0,%eax
jne    13c0 <sqrt+0x2c0>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
mov    -0x44(%rbp),%ecx
mov    -0x34(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x34(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jae    13e0 <sqrt+0x2e0>
mov    -0x34(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x24(%rbp),%ecx
mov    -0x2c(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x20(%rbp),%eax
add    -0x30(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x2c(%rbp),%ecx
and    -0x1c(%rbp),%ecx
shr    $0x1f,%ecx
add    %ecx,%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x2c(%rbp),%eax
add    -0x2c(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x20(%rbp),%eax
shr    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1353 <sqrt+0x253>
mov    -0x34(%rbp),%eax
or     -0x2c(%rbp),%eax
cmp    $0x0,%eax
je     14d1 <sqrt+0x3d1>
movsd  0xbca(%rip),%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  0xbb8(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     14cc <sqrt+0x3cc>
movsd  0xba6(%rip),%xmm0
movsd  %xmm0,-0x18(%rbp)
cmpl   $0xffffffff,-0x30(%rbp)
jne    147e <sqrt+0x37e>
movl   $0x0,-0x30(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    14c7 <sqrt+0x3c7>
movsd  -0x18(%rbp),%xmm0
movsd  0xb75(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    14b6 <sqrt+0x3b6>
cmpl   $0xfffffffe,-0x30(%rbp)
jne    14a8 <sqrt+0x3a8>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x30(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x30(%rbp)
jmp    14c2 <sqrt+0x3c2>
mov    -0x30(%rbp),%eax
and    $0x1,%eax
add    -0x30(%rbp),%eax
mov    %eax,-0x30(%rbp)
jmp    14c7 <sqrt+0x3c7>
jmp    14cc <sqrt+0x3cc>
jmp    14d1 <sqrt+0x3d1>
mov    -0x3c(%rbp),%eax
sar    $0x1,%eax
add    $0x3fe00000,%eax
mov    %eax,-0x34(%rbp)
mov    -0x30(%rbp),%eax
shr    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x3c(%rbp),%eax
and    $0x1,%eax
cmp    $0x1,%eax
jne    1500 <sqrt+0x400>
mov    -0x1c(%rbp),%eax
or     -0x2c(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x40(%rbp),%eax
shl    $0x14,%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<sqrt>:
push   %r14
push   %rbx
movq   %xmm0,%r14
mov    %r14,%rdi
shr    $0x20,%rdi
mov    %edi,%eax
not    %eax
test   $0x7ff00000,%eax
jne    112f <sqrt+0x2f>
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
addsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
jmp    12e3 <sqrt+0x1e3>
test   %edi,%edi
jle    119b <sqrt+0x9b>
cmp    $0xfffff,%edi
ja     1192 <sqrt+0x92>
mov    $0x1,%r8d
test   %edi,%edi
jne    1166 <sqrt+0x66>
mov    $0x1,%r8d
nopl   0x0(%rax,%rax,1)
mov    %r14d,%edi
shl    $0x15,%r14d
add    $0xffffffeb,%r8d
cmp    $0x800,%edi
jb     1150 <sqrt+0x50>
shr    $0xb,%edi
xor    %eax,%eax
test   $0x100000,%edi
jne    11bd <sqrt+0xbd>
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
lea    (%rdi,%rdi,1),%edx
add    $0x1,%eax
test   $0x80000,%edi
mov    %edx,%edi
je     1180 <sqrt+0x80>
jmp    11bf <sqrt+0xbf>
mov    %edi,%r8d
sar    $0x14,%r8d
jmp    11d2 <sqrt+0xd2>
mov    %edi,%eax
and    $0x7fffffff,%eax
or     %r14d,%eax
je     12e3 <sqrt+0x1e3>
test   %r14,%r14
jns    1133 <sqrt+0x33>
subsd  %xmm0,%xmm0
divsd  %xmm0,%xmm0
jmp    12e3 <sqrt+0x1e3>
mov    %edi,%edx
sub    %eax,%r8d
mov    %eax,%ecx
neg    %cl
mov    %r14d,%edi
shr    %cl,%edi
or     %edx,%edi
mov    %eax,%ecx
shl    %cl,%r14d
add    $0xfffffc01,%r8d
and    $0xfffff,%edi
or     $0x100000,%edi
test   $0x1,%r8b
je     11f3 <sqrt+0xf3>
shld   $0x1,%r14d,%edi
add    %r14d,%r14d
shld   $0x1,%r14d,%edi
add    %r14d,%r14d
mov    $0x200000,%ecx
xor    %r9d,%r9d
xor    %r11d,%r11d
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
mov    %edx,%eax
lea    (%rax,%rcx,1),%esi
lea    (%rsi,%rcx,1),%edx
cmp    %esi,%edi
cmovl  %r9d,%esi
cmovl  %eax,%edx
mov    %ecx,%eax
cmovl  %r9d,%eax
sub    %esi,%edi
shld   $0x1,%r14d,%edi
add    %eax,%r11d
mov    %ecx,%eax
shr    %eax
add    %r14d,%r14d
cmp    $0x1,%ecx
mov    %eax,%ecx
ja     1210 <sqrt+0x110>
mov    $0x80000000,%ecx
xor    %r9d,%r9d
xor    %r10d,%r10d
jmp    1287 <sqrt+0x187>
nopl   0x0(%rax)
test   %eax,%eax
sets   %sil
mov    %eax,%r10d
add    %ecx,%r10d
setns  %bl
and    %sil,%bl
movzbl %bl,%esi
add    %edx,%esi
sub    %eax,%r14d
sbb    $0x0,%edi
sub    %edx,%edi
add    %ecx,%r9d
mov    %esi,%edx
shld   $0x1,%r14d,%edi
add    %r14d,%r14d
mov    %ecx,%eax
shr    %eax
cmp    $0x1,%ecx
mov    %eax,%ecx
jbe    1298 <sqrt+0x198>
lea    (%r10,%rcx,1),%eax
cmp    %edx,%edi
jg     1250 <sqrt+0x150>
jne    1274 <sqrt+0x174>
cmp    %eax,%r14d
jae    1250 <sqrt+0x150>
jmp    1274 <sqrt+0x174>
or     %r14d,%edi
je     12b8 <sqrt+0x1b8>
cmp    $0xffffffff,%r9d
je     12b1 <sqrt+0x1b1>
mov    %r9d,%eax
and    $0x1,%eax
add    %r9d,%eax
mov    %eax,%r9d
jmp    12b8 <sqrt+0x1b8>
add    $0x1,%r11d
xor    %r9d,%r9d
mov    %r11d,%eax
shrd   $0x1,%r11d,%r9d
sar    %eax
shl    $0x13,%r8d
add    $0x3fe00000,%r8d
and    $0xfff00000,%r8d
add    %eax,%r8d
shl    $0x20,%r8
or     %r8,%r9
movq   %r9,%xmm0
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define double @sqrt(double %0) {
%2 = alloca double
%3 = alloca double
%4 = alloca double
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
store double %0, double* %3
store i32 -2147483648, i32* %5
%17 = bitcast double* %3 to i32*
%18 = getelementptr inbounds i32, i32* %17, i64 1
%19 = load i32, i32* %18
store i32 %19, i32* %11
%20 = bitcast double* %3 to i32*
%21 = load i32, i32* %20
store i32 %21, i32* %9
%22 = load i32, i32* %11
%23 = and i32 %22, 2146435072
%24 = icmp eq i32 %23, 2146435072
br i1 %24, label %25, label %30
25: ; preds = %1
%26 = load double, double* %3
%27 = load double, double* %3
%28 = load double, double* %3
%29 = call double @llvm.fmuladd.f64(double %26, double %27, double %28)
store double %29, double* %2
br label %295
30: ; preds = %1
%31 = load i32, i32* %11
%32 = icmp sle i32 %31, 0
br i1 %32, label %33, label %56
33: ; preds = %30
%34 = load i32, i32* %11
%35 = load i32, i32* %5
%36 = xor i32 %35, -1
%37 = and i32 %34, %36
%38 = load i32, i32* %9
%39 = or i32 %37, %38
%40 = icmp eq i32 %39, 0
br i1 %40, label %41, label %43
41: ; preds = %33
%42 = load double, double* %3
store double %42, double* %2
br label %295
43: ; preds = %33
%44 = load i32, i32* %11
%45 = icmp slt i32 %44, 0
br i1 %45, label %46, label %54
46: ; preds = %43
%47 = load double, double* %3
%48 = load double, double* %3
%49 = fsub double %47, %48
%50 = load double, double* %3
%51 = load double, double* %3
%52 = fsub double %50, %51
%53 = fdiv double %49, %52
store double %53, double* %2
br label %295
54: ; preds = %43
br label %55
55: ; preds = %54
br label %56
56: ; preds = %55, %30
%57 = load i32, i32* %11
%58 = ashr i32 %57, 20
store i32 %58, i32* %14
%59 = load i32, i32* %14
%60 = icmp eq i32 %59, 0
br i1 %60, label %61, label %99
61: ; preds = %56
br label %62
62: ; preds = %65, %61
%63 = load i32, i32* %11
%64 = icmp eq i32 %63, 0
br i1 %64, label %65, label %74
65: ; preds = %62
%66 = load i32, i32* %14
%67 = sub nsw i32 %66, 21
store i32 %67, i32* %14
%68 = load i32, i32* %9
%69 = lshr i32 %68, 11
%70 = load i32, i32* %11
%71 = or i32 %70, %69
store i32 %71, i32* %11
%72 = load i32, i32* %9
%73 = shl i32 %72, 21
store i32 %73, i32* %9
br label %62
74: ; preds = %62
store i32 0, i32* %16
br label %75
75: ; preds = %82, %74
%76 = load i32, i32* %11
%77 = and i32 %76, 1048576
%78 = icmp eq i32 %77, 0
br i1 %78, label %79, label %85
79: ; preds = %75
%80 = load i32, i32* %11
%81 = shl i32 %80, 1
store i32 %81, i32* %11
br label %82
82: ; preds = %79
%83 = load i32, i32* %16
%84 = add nsw i32 %83, 1
store i32 %84, i32* %16
br label %75
85: ; preds = %75
%86 = load i32, i32* %16
%87 = sub nsw i32 %86, 1
%88 = load i32, i32* %14
%89 = sub nsw i32 %88, %87
store i32 %89, i32* %14
%90 = load i32, i32* %9
%91 = load i32, i32* %16
%92 = sub nsw i32 32, %91
%93 = lshr i32 %90, %92
%94 = load i32, i32* %11
%95 = or i32 %94, %93
store i32 %95, i32* %11
%96 = load i32, i32* %16
%97 = load i32, i32* %9
%98 = shl i32 %97, %96
store i32 %98, i32* %9
br label %99
99: ; preds = %85, %56
%100 = load i32, i32* %14
%101 = sub nsw i32 %100, 1023
store i32 %101, i32* %14
%102 = load i32, i32* %11
%103 = and i32 %102, 1048575
%104 = or i32 %103, 1048576
store i32 %104, i32* %11
%105 = load i32, i32* %14
%106 = and i32 %105, 1
%107 = icmp ne i32 %106, 0
br i1 %107, label %108, label %120
108: ; preds = %99
%109 = load i32, i32* %11
%110 = load i32, i32* %9
%111 = load i32, i32* %5
%112 = and i32 %110, %111
%113 = lshr i32 %112, 31
%114 = add i32 %109, %113
%115 = load i32, i32* %11
%116 = add i32 %115, %114
store i32 %116, i32* %11
%117 = load i32, i32* %9
%118 = load i32, i32* %9
%119 = add i32 %118, %117
store i32 %119, i32* %9
br label %120
120: ; preds = %108, %99
%121 = load i32, i32* %14
%122 = ashr i32 %121, 1
store i32 %122, i32* %14
%123 = load i32, i32* %11
%124 = load i32, i32* %9
%125 = load i32, i32* %5
%126 = and i32 %124, %125
%127 = lshr i32 %126, 31
%128 = add i32 %123, %127
%129 = load i32, i32* %11
%130 = add i32 %129, %128
store i32 %130, i32* %11
%131 = load i32, i32* %9
%132 = load i32, i32* %9
%133 = add i32 %132, %131
store i32 %133, i32* %9
store i32 0, i32* %8
store i32 0, i32* %12
store i32 0, i32* %10
store i32 0, i32* %13
store i32 2097152, i32* %6
br label %134
134: ; preds = %154, %120
%135 = load i32, i32* %6
%136 = icmp ne i32 %135, 0
br i1 %136, label %137, label %168
137: ; preds = %134
%138 = load i32, i32* %12
%139 = load i32, i32* %6
%140 = add i32 %138, %139
store i32 %140, i32* %15
%141 = load i32, i32* %15
%142 = load i32, i32* %11
%143 = icmp sle i32 %141, %142
br i1 %143, label %144, label %154
144: ; preds = %137
%145 = load i32, i32* %15
%146 = load i32, i32* %6
%147 = add i32 %145, %146
store i32 %147, i32* %12
%148 = load i32, i32* %15
%149 = load i32, i32* %11
%150 = sub nsw i32 %149, %148
store i32 %150, i32* %11
%151 = load i32, i32* %6
%152 = load i32, i32* %13
%153 = add i32 %152, %151
store i32 %153, i32* %13
br label %154
154: ; preds = %144, %137
%155 = load i32, i32* %11
%156 = load i32, i32* %9
%157 = load i32, i32* %5
%158 = and i32 %156, %157
%159 = lshr i32 %158, 31
%160 = add i32 %155, %159
%161 = load i32, i32* %11
%162 = add i32 %161, %160
store i32 %162, i32* %11
%163 = load i32, i32* %9
%164 = load i32, i32* %9
%165 = add i32 %164, %163
store i32 %165, i32* %9
%166 = load i32, i32* %6
%167 = lshr i32 %166, 1
store i32 %167, i32* %6
br label %134
168: ; preds = %134
%169 = load i32, i32* %5
store i32 %169, i32* %6
br label %170
170: ; preds = %223, %168
%171 = load i32, i32* %6
%172 = icmp ne i32 %171, 0
br i1 %172, label %173, label %237
173: ; preds = %170
%174 = load i32, i32* %8
%175 = load i32, i32* %6
%176 = add i32 %174, %175
store i32 %176, i32* %7
%177 = load i32, i32* %12
store i32 %177, i32* %15
%178 = load i32, i32* %15
%179 = load i32, i32* %11
%180 = icmp slt i32 %178, %179
br i1 %180, label %189, label %181
181: ; preds = %173
%182 = load i32, i32* %15
%183 = load i32, i32* %11
%184 = icmp eq i32 %182, %183
br i1 %184, label %185, label %223
185: ; preds = %181
%186 = load i32, i32* %7
%187 = load i32, i32* %9
%188 = icmp ule i32 %186, %187
br i1 %188, label %189, label %223
189: ; preds = %185, %173
%190 = load i32, i32* %7
%191 = load i32, i32* %6
%192 = add i32 %190, %191
store i32 %192, i32* %8
%193 = load i32, i32* %7
%194 = load i32, i32* %5
%195 = and i32 %193, %194
%196 = load i32, i32* %5
%197 = icmp eq i32 %195, %196
br i1 %197, label %198, label %206
198: ; preds = %189
%199 = load i32, i32* %8
%200 = load i32, i32* %5
%201 = and i32 %199, %200
%202 = icmp eq i32 %201, 0
br i1 %202, label %203, label %206
203: ; preds = %198
%204 = load i32, i32* %12
%205 = add nsw i32 %204, 1
store i32 %205, i32* %12
br label %206
206: ; preds = %203, %198, %189
%207 = load i32, i32* %15
%208 = load i32, i32* %11
%209 = sub nsw i32 %208, %207
store i32 %209, i32* %11
%210 = load i32, i32* %9
%211 = load i32, i32* %7
%212 = icmp ult i32 %210, %211
br i1 %212, label %213, label %216
213: ; preds = %206
%214 = load i32, i32* %11
%215 = sub nsw i32 %214, 1
store i32 %215, i32* %11
br label %216
216: ; preds = %213, %206
%217 = load i32, i32* %7
%218 = load i32, i32* %9
%219 = sub i32 %218, %217
store i32 %219, i32* %9
%220 = load i32, i32* %6
%221 = load i32, i32* %10
%222 = add i32 %221, %220
store i32 %222, i32* %10
br label %223
223: ; preds = %216, %185, %181
%224 = load i32, i32* %11
%225 = load i32, i32* %9
%226 = load i32, i32* %5
%227 = and i32 %225, %226
%228 = lshr i32 %227, 31
%229 = add i32 %224, %228
%230 = load i32, i32* %11
%231 = add i32 %230, %229
store i32 %231, i32* %11
%232 = load i32, i32* %9
%233 = load i32, i32* %9
%234 = add i32 %233, %232
store i32 %234, i32* %9
%235 = load i32, i32* %6
%236 = lshr i32 %235, 1
store i32 %236, i32* %6
br label %170
237: ; preds = %170
%238 = load i32, i32* %11
%239 = load i32, i32* %9
%240 = or i32 %238, %239
%241 = icmp ne i32 %240, 0
br i1 %241, label %242, label %271
242: ; preds = %237
store double 1.000000e+00, double* %4
%243 = load double, double* %4
%244 = fcmp oge double %243, 1.000000e+00
br i1 %244, label %245, label %270
245: ; preds = %242
store double 1.000000e+00, double* %4
%246 = load i32, i32* %10
%247 = icmp eq i32 %246, -1
br i1 %247, label %248, label %251
248: ; preds = %245
store i32 0, i32* %10
%249 = load i32, i32* %13
%250 = add nsw i32 %249, 1
store i32 %250, i32* %13
br label %269
251: ; preds = %245
%252 = load double, double* %4
%253 = fcmp ogt double %252, 1.000000e+00
br i1 %253, label %254, label %263
254: ; preds = %251
%255 = load i32, i32* %10
%256 = icmp eq i32 %255, -2
br i1 %256, label %257, label %260
257: ; preds = %254
%258 = load i32, i32* %13
%259 = add nsw i32 %258, 1
store i32 %259, i32* %13
br label %260
260: ; preds = %257, %254
%261 = load i32, i32* %10
%262 = add i32 %261, 2
store i32 %262, i32* %10
br label %268
263: ; preds = %251
%264 = load i32, i32* %10
%265 = and i32 %264, 1
%266 = load i32, i32* %10
%267 = add i32 %266, %265
store i32 %267, i32* %10
br label %268
268: ; preds = %263, %260
br label %269
269: ; preds = %268, %248
br label %270
270: ; preds = %269, %242
br label %271
271: ; preds = %270, %237
%272 = load i32, i32* %13
%273 = ashr i32 %272, 1
%274 = add nsw i32 %273, 1071644672
store i32 %274, i32* %11
%275 = load i32, i32* %10
%276 = lshr i32 %275, 1
store i32 %276, i32* %9
%277 = load i32, i32* %13
%278 = and i32 %277, 1
%279 = icmp eq i32 %278, 1
br i1 %279, label %280, label %284
280: ; preds = %271
%281 = load i32, i32* %5
%282 = load i32, i32* %9
%283 = or i32 %282, %281
store i32 %283, i32* %9
br label %284
284: ; preds = %280, %271
%285 = load i32, i32* %14
%286 = shl i32 %285, 20
%287 = load i32, i32* %11
%288 = add nsw i32 %287, %286
store i32 %288, i32* %11
%289 = load i32, i32* %11
%290 = bitcast double* %4 to i32*
%291 = getelementptr inbounds i32, i32* %290, i64 1
store i32 %289, i32* %291
%292 = load i32, i32* %9
%293 = bitcast double* %4 to i32*
store i32 %292, i32* %293
%294 = load double, double* %4
store double %294, double* %2
br label %295
295: ; preds = %284, %46, %41, %25
%296 = load double, double* %2
ret double %296
}

LLVM IR (-O2):
define double @sqrt(double %0) {
%2 = bitcast double %0 to i64
%3 = lshr i64 %2, 32
%4 = trunc i64 %3 to i32
%5 = trunc i64 %2 to i32
%6 = and i32 %4, 2146435072
%7 = icmp eq i32 %6, 2146435072
br i1 %7, label %8, label %10
8: ; preds = %1
%9 = tail call double @llvm.fmuladd.f64(double %0, double %0, double %0)
br label %151
10: ; preds = %1
%11 = icmp slt i32 %4, 1
br i1 %11, label %12, label %21
12: ; preds = %10
%13 = and i32 %4, 2147483647
%14 = or i32 %13, %5
%15 = icmp eq i32 %14, 0
br i1 %15, label %151, label %16
16: ; preds = %12
%17 = icmp slt i64 %2, 0
br i1 %17, label %18, label %21
18: ; preds = %16
%19 = fsub double %0, %0
%20 = fdiv double %19, %19
br label %151
21: ; preds = %16, %10
%22 = ashr i32 %4, 20
%23 = icmp ult i32 %4, 1048576
br i1 %23, label %24, label %56
24: ; preds = %21
%25 = icmp eq i32 %4, 0
br i1 %25, label %35, label %29
26: ; preds = %35
%27 = lshr i32 %37, 11
%28 = add i32 %36, -20
br label %29
29: ; preds = %26, %24
%30 = phi i32 [ %5, %24 ], [ %39, %26 ]
%31 = phi i32 [ %4, %24 ], [ %27, %26 ]
%32 = phi i32 [ 1, %24 ], [ %28, %26 ]
%33 = and i32 %31, 1048576
%34 = icmp eq i32 %33, 0
br i1 %34, label %41, label %48
35: ; preds = %24, %35
%36 = phi i32 [ %38, %35 ], [ 0, %24 ]
%37 = phi i32 [ %39, %35 ], [ %5, %24 ]
%38 = add nsw i32 %36, -21
%39 = shl i32 %37, 21
%40 = icmp ult i32 %37, 2048
br i1 %40, label %35, label %26
41: ; preds = %29, %41
%42 = phi i32 [ %45, %41 ], [ 0, %29 ]
%43 = phi i32 [ %44, %41 ], [ %31, %29 ]
%44 = shl i32 %43, 1
%45 = add nuw nsw i32 %42, 1
%46 = and i32 %43, 524288
%47 = icmp eq i32 %46, 0
br i1 %47, label %41, label %48
48: ; preds = %41, %29
%49 = phi i32 [ %31, %29 ], [ %44, %41 ]
%50 = phi i32 [ 0, %29 ], [ %45, %41 ]
%51 = sub i32 %32, %50
%52 = sub nsw i32 32, %50
%53 = lshr i32 %30, %52
%54 = or i32 %53, %49
%55 = shl i32 %30, %50
br label %56
56: ; preds = %48, %21
%57 = phi i32 [ %55, %48 ], [ %5, %21 ]
%58 = phi i32 [ %54, %48 ], [ %4, %21 ]
%59 = phi i32 [ %51, %48 ], [ %22, %21 ]
%60 = add nsw i32 %59, -1023
%61 = and i32 %58, 1048575
%62 = or i32 %61, 1048576
%63 = and i32 %60, 1
%64 = icmp eq i32 %63, 0
br i1 %64, label %68, label %65
65: ; preds = %56
%66 = tail call i32 @llvm.fshl.i32(i32 %62, i32 %57, i32 1)
%67 = shl i32 %57, 1
br label %68
68: ; preds = %65, %56
%69 = phi i32 [ %67, %65 ], [ %57, %56 ]
%70 = phi i32 [ %66, %65 ], [ %62, %56 ]
%71 = tail call i32 @llvm.fshl.i32(i32 %70, i32 %69, i32 1)
%72 = shl i32 %69, 1
br label %73
73: ; preds = %68, %73
%74 = phi i32 [ %72, %68 ], [ %89, %73 ]
%75 = phi i32 [ 0, %68 ], [ %86, %73 ]
%76 = phi i32 [ 0, %68 ], [ %84, %73 ]
%77 = phi i32 [ %71, %68 ], [ %87, %73 ]
%78 = phi i32 [ 2097152, %68 ], [ %88, %73 ]
%79 = add i32 %76, %78
%80 = icmp slt i32 %77, %79
%81 = add i32 %79, %78
%82 = select i1 %80, i32 0, i32 %79
%83 = sub nsw i32 %77, %82
%84 = select i1 %80, i32 %76, i32 %81
%85 = select i1 %80, i32 0, i32 %78
%86 = add i32 %85, %75
%87 = tail call i32 @llvm.fshl.i32(i32 %83, i32 %74, i32 1)
%88 = lshr i32 %78, 1
%89 = shl i32 %74, 1
%90 = icmp ult i32 %78, 2
br i1 %90, label %91, label %73
91: ; preds = %73, %117
%92 = phi i32 [ %122, %117 ], [ %84, %73 ]
%93 = phi i32 [ %123, %117 ], [ %87, %73 ]
%94 = phi i32 [ %120, %117 ], [ 0, %73 ]
%95 = phi i32 [ %124, %117 ], [ %89, %73 ]
%96 = phi i32 [ %118, %117 ], [ 0, %73 ]
%97 = phi i32 [ %125, %117 ], [ -2147483648, %73 ]
%98 = add i32 %96, %97
%99 = icmp sgt i32 %93, %92
br i1 %99, label %104, label %100
100: ; preds = %91
%101 = icmp ne i32 %93, %92
%102 = icmp ult i32 %95, %98
%103 = select i1 %101, i1 true, i1 %102
br i1 %103, label %117, label %104
104: ; preds = %100, %91
%105 = add i32 %98, %97
%106 = icmp slt i32 %98, 0
%107 = icmp sgt i32 %105, -1
%108 = and i1 %106, %107
%109 = zext i1 %108 to i32
%110 = add nsw i32 %92, %109
%111 = icmp ult i32 %95, %98
%112 = sext i1 %111 to i32
%113 = add i32 %93, %112
%114 = sub i32 %113, %92
%115 = sub i32 %95, %98
%116 = add i32 %94, %97
br label %117
117: ; preds = %104, %100
%118 = phi i32 [ %105, %104 ], [ %96, %100 ]
%119 = phi i32 [ %115, %104 ], [ %95, %100 ]
%120 = phi i32 [ %116, %104 ], [ %94, %100 ]
%121 = phi i32 [ %114, %104 ], [ %93, %100 ]
%122 = phi i32 [ %110, %104 ], [ %92, %100 ]
%123 = tail call i32 @llvm.fshl.i32(i32 %121, i32 %119, i32 1)
%124 = shl i32 %119, 1
%125 = lshr i32 %97, 1
%126 = icmp ult i32 %97, 2
br i1 %126, label %127, label %91
127: ; preds = %117
%128 = or i32 %123, %124
%129 = icmp eq i32 %128, 0
br i1 %129, label %137, label %130
130: ; preds = %127
%131 = icmp eq i32 %120, -1
br i1 %131, label %132, label %134
132: ; preds = %130
%133 = add nsw i32 %86, 1
br label %137
134: ; preds = %130
%135 = and i32 %120, 1
%136 = add nuw i32 %135, %120
br label %137
137: ; preds = %134, %132, %127
%138 = phi i32 [ 0, %132 ], [ %136, %134 ], [ %120, %127 ]
%139 = phi i32 [ %133, %132 ], [ %86, %134 ], [ %86, %127 ]
%140 = ashr i32 %139, 1
%141 = tail call i32 @llvm.fshl.i32(i32 %139, i32 %138, i32 31)
%142 = shl i32 %60, 19
%143 = add i32 %142, 1071644672
%144 = and i32 %143, -1048576
%145 = add i32 %144, %140
%146 = zext i32 %145 to i64
%147 = shl nuw i64 %146, 32
%148 = zext i32 %141 to i64
%149 = or i64 %147, %148
%150 = bitcast i64 %149 to double
br label %151
151: ; preds = %12, %137, %18, %8
%152 = phi double [ %9, %8 ], [ %20, %18 ], [ %150, %137 ], [ %0, %12 ]
ret double %152
}