double floor(double x) {
  int i0, i1, j0;
  unsigned i, j;
  i0 = *(1 + (int *)&x);
  i1 = *(int *)&x;
  j0 = ((i0 >> 20) & 0x7ff) - 0x3ff;
  if (j0 < 20) {
    if (j0 < 0) {
      if (1.0e300 + x > 0.0) {
        if (i0 >= 0) {
          i0 = i1 = 0;
        } else if (((i0 & 0x7fffffff) | i1) != 0) {
          i0 = 0xbff00000;
          i1 = 0;
        }
      }
    } else {
      i = (0x000fffff) >> j0;
      if (((i0 & i) | i1) == 0) {
        return x;
      }
      if (1.0e300 + x > 0.0) {
        if (i0 < 0) {
          i0 += (0x00100000) >> j0;
        }
        i0 &= (~i);
        i1 = 0;
      }
    }
  } else if (j0 > 51) {
    if (j0 == 0x400) {
      return x + x;
    } else {
      return x;
    }
  } else {
    i = ((unsigned)(0xffffffff)) >> (j0 - 20);
    if ((i1 & i) == 0) {
      return x;
    }
    if (1.0e300 + x > 0.0) {
      if (i0 < 0) {
        if (j0 == 20) {
          i0 += 1;
        } else {
          j = i1 + (1 << (52 - j0));
          if (j < i1) {
            i0 += 1;
          }
          i1 = j;
        }
      }
      i1 &= (~i);
    }
  }
  *(1 + (int *)&x) = i0;
  *(int *)&x = i1;
  return x;
}


x86 assembly (clang, -O0):
<floor>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
mov    -0xc(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x10(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x14(%rbp),%eax
sar    $0x14,%eax
and    $0x7ff,%eax
sub    $0x3ff,%eax
mov    %eax,-0x1c(%rbp)
cmpl   $0x14,-0x1c(%rbp)
jge    1223 <floor+0x123>
cmpl   $0x0,-0x1c(%rbp)
jge    11a4 <floor+0xa4>
movsd  0xebc(%rip),%xmm0
addsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    119f <floor+0x9f>
cmpl   $0x0,-0x14(%rbp)
jl     1173 <floor+0x73>
movl   $0x0,-0x18(%rbp)
movl   $0x0,-0x14(%rbp)
jmp    119a <floor+0x9a>
mov    -0x14(%rbp),%eax
and    $0x7fffffff,%eax
or     -0x18(%rbp),%eax
cmp    $0x0,%eax
je     1195 <floor+0x95>
movl   $0xbff00000,-0x14(%rbp)
movl   $0x0,-0x18(%rbp)
jmp    119a <floor+0x9a>
jmp    119f <floor+0x9f>
jmp    121e <floor+0x11e>
mov    -0x1c(%rbp),%ecx
mov    $0xfffff,%eax
sar    %cl,%eax
mov    %eax,-0x20(%rbp)
mov    -0x14(%rbp),%eax
and    -0x20(%rbp),%eax
or     -0x18(%rbp),%eax
cmp    $0x0,%eax
jne    11d2 <floor+0xd2>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    132c <floor+0x22c>
movsd  0xe26(%rip),%xmm0
addsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1219 <floor+0x119>
cmpl   $0x0,-0x14(%rbp)
jge    1206 <floor+0x106>
mov    -0x1c(%rbp),%ecx
mov    $0x100000,%eax
sar    %cl,%eax
add    -0x14(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x20(%rbp),%eax
xor    $0xffffffff,%eax
and    -0x14(%rbp),%eax
mov    %eax,-0x14(%rbp)
movl   $0x0,-0x18(%rbp)
jmp    121e <floor+0x11e>
jmp    1316 <floor+0x216>
cmpl   $0x33,-0x1c(%rbp)
jle    125d <floor+0x15d>
cmpl   $0x400,-0x1c(%rbp)
jne    124e <floor+0x14e>
movsd  -0x10(%rbp),%xmm0
addsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    132c <floor+0x22c>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    132c <floor+0x22c>
mov    -0x1c(%rbp),%ecx
sub    $0x14,%ecx
mov    $0xffffffff,%eax
shr    %cl,%eax
mov    %eax,-0x20(%rbp)
mov    -0x18(%rbp),%eax
and    -0x20(%rbp),%eax
cmp    $0x0,%eax
jne    128b <floor+0x18b>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    132c <floor+0x22c>
movsd  0xd6d(%rip),%xmm0
addsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    130c <floor+0x20c>
cmpl   $0x0,-0x14(%rbp)
jge    1300 <floor+0x200>
cmpl   $0x14,-0x1c(%rbp)
jne    12c7 <floor+0x1c7>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    12fb <floor+0x1fb>
mov    -0x18(%rbp),%eax
mov    $0x34,%ecx
sub    -0x1c(%rbp),%ecx
mov    $0x1,%edx
shl    %cl,%edx
mov    %edx,%ecx
add    %ecx,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x18(%rbp),%eax
jae    12f5 <floor+0x1f5>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
mov    -0x24(%rbp),%eax
mov    %eax,-0x18(%rbp)
jmp    1300 <floor+0x200>
mov    -0x20(%rbp),%eax
xor    $0xffffffff,%eax
and    -0x18(%rbp),%eax
mov    %eax,-0x18(%rbp)
jmp    1311 <floor+0x211>
jmp    1316 <floor+0x216>
mov    -0x14(%rbp),%eax
mov    %eax,-0xc(%rbp)
mov    -0x18(%rbp),%eax
mov    %eax,-0x10(%rbp)
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<floor>:
movq   %xmm0,%rsi
mov    %rsi,%rax
shr    $0x20,%rax
mov    %rsi,%rdx
shr    $0x34,%rdx
and    $0x7ff,%edx
lea    -0x3ff(%rdx),%r8d
cmp    $0x412,%edx
ja     1172 <floor+0x72>
cmp    $0x3fe,%edx
ja     1194 <floor+0x94>
addsd  0xec8(%rip),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1209 <floor+0x109>
xor    %edx,%edx
mov    $0x0,%edi
test   %rsi,%rsi
jns    1227 <floor+0x127>
mov    %eax,%ecx
and    $0x7fffffff,%ecx
xor    %edi,%edi
or     %esi,%ecx
mov    $0xbff00000,%edx
cmove  %eax,%edx
cmove  %esi,%edi
jmp    1227 <floor+0x127>
cmp    $0x433,%edx
jb     11d9 <floor+0xd9>
movapd %xmm0,%xmm1
cmp    $0x400,%r8d
je     118f <floor+0x8f>
movsd  0xe79(%rip),%xmm1
addsd  %xmm1,%xmm0
ret
mov    $0xfffff,%edx
mov    %r8d,%ecx
shr    %cl,%edx
and    %eax,%edx
or     %esi,%edx
je     1235 <floor+0x135>
addsd  0xe50(%rip),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1209 <floor+0x109>
mov    $0x100000,%edi
mov    %r8d,%ecx
shr    %cl,%edi
sar    $0x3f,%rsi
mov    $0xfff00000,%edx
sar    %cl,%edx
and    %edi,%esi
add    %eax,%esi
and    %esi,%edx
xor    %edi,%edi
jmp    1227 <floor+0x127>
lea    -0x13(%rdx),%ecx
mov    $0xffffffff,%edi
shr    %cl,%edi
test   %esi,%edi
je     1235 <floor+0x135>
addsd  0xe11(%rip),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1209 <floor+0x109>
test   %rsi,%rsi
jns    1221 <floor+0x121>
cmp    $0x14,%r8d
jne    120f <floor+0x10f>
add    $0x1,%eax
jmp    1221 <floor+0x121>
mov    %eax,%edx
mov    %esi,%edi
jmp    1227 <floor+0x127>
mov    $0x33,%cl
sub    %dl,%cl
mov    $0x1,%edx
shl    %cl,%edx
add    %esi,%edx
adc    $0x0,%eax
mov    %edx,%esi
not    %edi
and    %esi,%edi
mov    %eax,%edx
shl    $0x20,%rdx
mov    %edi,%eax
or     %rdx,%rax
movq   %rax,%xmm0
ret

LLVM IR (-O0):
define double @floor(double %0) {
%2 = alloca double
%3 = alloca double
%4 = alloca i32
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
store double %0, double* %3
%9 = bitcast double* %3 to i32*
%10 = getelementptr inbounds i32, i32* %9, i64 1
%11 = load i32, i32* %10
store i32 %11, i32* %4
%12 = bitcast double* %3 to i32*
%13 = load i32, i32* %12
store i32 %13, i32* %5
%14 = load i32, i32* %4
%15 = ashr i32 %14, 20
%16 = and i32 %15, 2047
%17 = sub nsw i32 %16, 1023
store i32 %17, i32* %6
%18 = load i32, i32* %6
%19 = icmp slt i32 %18, 20
br i1 %19, label %20, label %71
20: ; preds = %1
%21 = load i32, i32* %6
%22 = icmp slt i32 %21, 0
br i1 %22, label %23, label %41
23: ; preds = %20
%24 = load double, double* %3
%25 = fadd double 1.000000e+300, %24
%26 = fcmp ogt double %25, 0.000000e+00
br i1 %26, label %27, label %40
27: ; preds = %23
%28 = load i32, i32* %4
%29 = icmp sge i32 %28, 0
br i1 %29, label %30, label %31
30: ; preds = %27
store i32 0, i32* %5
store i32 0, i32* %4
br label %39
31: ; preds = %27
%32 = load i32, i32* %4
%33 = and i32 %32, 2147483647
%34 = load i32, i32* %5
%35 = or i32 %33, %34
%36 = icmp ne i32 %35, 0
br i1 %36, label %37, label %38
37: ; preds = %31
store i32 -1074790400, i32* %4
store i32 0, i32* %5
br label %38
38: ; preds = %37, %31
br label %39
39: ; preds = %38, %30
br label %40
40: ; preds = %39, %23
br label %70
41: ; preds = %20
%42 = load i32, i32* %6
%43 = ashr i32 1048575, %42
store i32 %43, i32* %7
%44 = load i32, i32* %4
%45 = load i32, i32* %7
%46 = and i32 %44, %45
%47 = load i32, i32* %5
%48 = or i32 %46, %47
%49 = icmp eq i32 %48, 0
br i1 %49, label %50, label %52
50: ; preds = %41
%51 = load double, double* %3
store double %51, double* %2
br label %135
52: ; preds = %41
%53 = load double, double* %3
%54 = fadd double 1.000000e+300, %53
%55 = fcmp ogt double %54, 0.000000e+00
br i1 %55, label %56, label %69
56: ; preds = %52
%57 = load i32, i32* %4
%58 = icmp slt i32 %57, 0
br i1 %58, label %59, label %64
59: ; preds = %56
%60 = load i32, i32* %6
%61 = ashr i32 1048576, %60
%62 = load i32, i32* %4
%63 = add nsw i32 %62, %61
store i32 %63, i32* %4
br label %64
64: ; preds = %59, %56
%65 = load i32, i32* %7
%66 = xor i32 %65, -1
%67 = load i32, i32* %4
%68 = and i32 %67, %66
store i32 %68, i32* %4
store i32 0, i32* %5
br label %69
69: ; preds = %64, %52
br label %70
70: ; preds = %69, %40
br label %128
71: ; preds = %1
%72 = load i32, i32* %6
%73 = icmp sgt i32 %72, 51
br i1 %73, label %74, label %83
74: ; preds = %71
%75 = load i32, i32* %6
%76 = icmp eq i32 %75, 1024
br i1 %76, label %77, label %81
77: ; preds = %74
%78 = load double, double* %3
%79 = load double, double* %3
%80 = fadd double %78, %79
store double %80, double* %2
br label %135
81: ; preds = %74
%82 = load double, double* %3
store double %82, double* %2
br label %135
83: ; preds = %71
%84 = load i32, i32* %6
%85 = sub nsw i32 %84, 20
%86 = lshr i32 -1, %85
store i32 %86, i32* %7
%87 = load i32, i32* %5
%88 = load i32, i32* %7
%89 = and i32 %87, %88
%90 = icmp eq i32 %89, 0
br i1 %90, label %91, label %93
91: ; preds = %83
%92 = load double, double* %3
store double %92, double* %2
br label %135
93: ; preds = %83
%94 = load double, double* %3
%95 = fadd double 1.000000e+300, %94
%96 = fcmp ogt double %95, 0.000000e+00
br i1 %96, label %97, label %126
97: ; preds = %93
%98 = load i32, i32* %4
%99 = icmp slt i32 %98, 0
br i1 %99, label %100, label %121
100: ; preds = %97
%101 = load i32, i32* %6
%102 = icmp eq i32 %101, 20
br i1 %102, label %103, label %106
103: ; preds = %100
%104 = load i32, i32* %4
%105 = add nsw i32 %104, 1
store i32 %105, i32* %4
br label %120
106: ; preds = %100
%107 = load i32, i32* %5
%108 = load i32, i32* %6
%109 = sub nsw i32 52, %108
%110 = shl i32 1, %109
%111 = add nsw i32 %107, %110
store i32 %111, i32* %8
%112 = load i32, i32* %8
%113 = load i32, i32* %5
%114 = icmp ult i32 %112, %113
br i1 %114, label %115, label %118
115: ; preds = %106
%116 = load i32, i32* %4
%117 = add nsw i32 %116, 1
store i32 %117, i32* %4
br label %118
118: ; preds = %115, %106
%119 = load i32, i32* %8
store i32 %119, i32* %5
br label %120
120: ; preds = %118, %103
br label %121
121: ; preds = %120, %97
%122 = load i32, i32* %7
%123 = xor i32 %122, -1
%124 = load i32, i32* %5
%125 = and i32 %124, %123
store i32 %125, i32* %5
br label %126
126: ; preds = %121, %93
br label %127
127: ; preds = %126
br label %128
128: ; preds = %127, %70
%129 = load i32, i32* %4
%130 = bitcast double* %3 to i32*
%131 = getelementptr inbounds i32, i32* %130, i64 1
store i32 %129, i32* %131
%132 = load i32, i32* %5
%133 = bitcast double* %3 to i32*
store i32 %132, i32* %133
%134 = load double, double* %3
store double %134, double* %2
br label %135
135: ; preds = %128, %91, %81, %77, %50
%136 = load double, double* %2
ret double %136
}

LLVM IR (-O2):
define double @floor(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 = lshr i32 %4, 20
%7 = and i32 %6, 2047
%8 = add nsw i32 %7, -1023
%9 = icmp ult i32 %7, 1043
br i1 %9, label %10, label %38
10: ; preds = %1
%11 = icmp ult i32 %7, 1023
br i1 %11, label %12, label %23
12: ; preds = %10
%13 = fadd double %0, 1.000000e+300
%14 = fcmp ogt double %13, 0.000000e+00
br i1 %14, label %15, label %70
15: ; preds = %12
%16 = icmp sgt i64 %2, -1
br i1 %16, label %70, label %17
17: ; preds = %15
%18 = and i32 %4, 2147483647
%19 = or i32 %18, %5
%20 = icmp eq i32 %19, 0
%21 = select i1 %20, i32 %4, i32 -1074790400
%22 = select i1 %20, i32 %5, i32 0
br label %70
23: ; preds = %10
%24 = lshr i32 1048575, %8
%25 = and i32 %24, %4
%26 = or i32 %25, %5
%27 = icmp eq i32 %26, 0
br i1 %27, label %78, label %28
28: ; preds = %23
%29 = fadd double %0, 1.000000e+300
%30 = fcmp ogt double %29, 0.000000e+00
br i1 %30, label %31, label %70
31: ; preds = %28
%32 = icmp slt i64 %2, 0
%33 = lshr i32 1048576, %8
%34 = select i1 %32, i32 %33, i32 0
%35 = add nsw i32 %34, %4
%36 = ashr i32 -1048576, %8
%37 = and i32 %35, %36
br label %70
38: ; preds = %1
%39 = icmp ugt i32 %7, 1074
br i1 %39, label %40, label %44
40: ; preds = %38
%41 = icmp eq i32 %8, 1024
%42 = select i1 %41, double %0, double -0.000000e+00
%43 = fadd double %42, %0
br label %78
44: ; preds = %38
%45 = add nsw i32 %7, -1043
%46 = lshr i32 -1, %45
%47 = and i32 %46, %5
%48 = icmp eq i32 %47, 0
br i1 %48, label %78, label %49
49: ; preds = %44
%50 = fadd double %0, 1.000000e+300
%51 = fcmp ogt double %50, 0.000000e+00
br i1 %51, label %52, label %70
52: ; preds = %49
%53 = icmp slt i64 %2, 0
br i1 %53, label %54, label %65
54: ; preds = %52
%55 = icmp eq i32 %8, 20
br i1 %55, label %56, label %58
56: ; preds = %54
%57 = add nsw i32 %4, 1
br label %65
58: ; preds = %54
%59 = sub nuw nsw i32 1075, %7
%60 = shl nuw i32 1, %59
%61 = add nsw i32 %60, %5
%62 = icmp ult i32 %61, %5
%63 = zext i1 %62 to i32
%64 = add nsw i32 %63, %4
br label %65
65: ; preds = %56, %58, %52
%66 = phi i32 [ %57, %56 ], [ %64, %58 ], [ %4, %52 ]
%67 = phi i32 [ %5, %56 ], [ %61, %58 ], [ %5, %52 ]
%68 = xor i32 %46, -1
%69 = and i32 %67, %68
br label %70
70: ; preds = %17, %15, %65, %49, %12, %31, %28
%71 = phi i32 [ %4, %12 ], [ %37, %31 ], [ %4, %28 ], [ %66, %65 ], [ %4, %49 ], [ 0, %15 ], [ %21, %17 ]
%72 = phi i32 [ %5, %12 ], [ 0, %31 ], [ %5, %28 ], [ %69, %65 ], [ %5, %49 ], [ 0, %15 ], [ %22, %17 ]
%73 = zext i32 %71 to i64
%74 = shl nuw i64 %73, 32
%75 = zext i32 %72 to i64
%76 = or i64 %74, %75
%77 = bitcast i64 %76 to double
br label %78
78: ; preds = %40, %44, %23, %70
%79 = phi double [ %77, %70 ], [ %0, %23 ], [ %0, %44 ], [ %43, %40 ]
ret double %79
}