double ceil(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 = 0x80000000;
          i1 = 0;
        } else if ((i0 | i1) != 0) {
          i0 = 0x3ff00000;
          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):
<ceil>:
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    121e <ceil+0x11e>
cmpl   $0x0,-0x1c(%rbp)
jge    119f <ceil+0x9f>
movsd  0xebc(%rip),%xmm0
addsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    119a <ceil+0x9a>
cmpl   $0x0,-0x14(%rbp)
jge    1173 <ceil+0x73>
movl   $0x80000000,-0x14(%rbp)
movl   $0x0,-0x18(%rbp)
jmp    1195 <ceil+0x95>
mov    -0x14(%rbp),%eax
or     -0x18(%rbp),%eax
cmp    $0x0,%eax
je     1190 <ceil+0x90>
movl   $0x3ff00000,-0x14(%rbp)
movl   $0x0,-0x18(%rbp)
jmp    1195 <ceil+0x95>
jmp    119a <ceil+0x9a>
jmp    1219 <ceil+0x119>
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    11cd <ceil+0xcd>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1327 <ceil+0x227>
movsd  0xe2b(%rip),%xmm0
addsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1214 <ceil+0x114>
cmpl   $0x0,-0x14(%rbp)
jle    1201 <ceil+0x101>
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    1219 <ceil+0x119>
jmp    1311 <ceil+0x211>
cmpl   $0x33,-0x1c(%rbp)
jle    1258 <ceil+0x158>
cmpl   $0x400,-0x1c(%rbp)
jne    1249 <ceil+0x149>
movsd  -0x10(%rbp),%xmm0
addsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1327 <ceil+0x227>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1327 <ceil+0x227>
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    1286 <ceil+0x186>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1327 <ceil+0x227>
movsd  0xd72(%rip),%xmm0
addsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1307 <ceil+0x207>
cmpl   $0x0,-0x14(%rbp)
jle    12fb <ceil+0x1fb>
cmpl   $0x14,-0x1c(%rbp)
jne    12c2 <ceil+0x1c2>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    12f6 <ceil+0x1f6>
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    12f0 <ceil+0x1f0>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
mov    -0x24(%rbp),%eax
mov    %eax,-0x18(%rbp)
jmp    12fb <ceil+0x1fb>
mov    -0x20(%rbp),%eax
xor    $0xffffffff,%eax
and    -0x18(%rbp),%eax
mov    %eax,-0x18(%rbp)
jmp    130c <ceil+0x20c>
jmp    1311 <ceil+0x211>
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):
<ceil>:
movq   %xmm0,%rdi
mov    %rdi,%rdx
shr    $0x20,%rdx
mov    %rdi,%rsi
shr    $0x34,%rsi
and    $0x7ff,%esi
lea    -0x3ff(%rsi),%eax
cmp    $0x412,%esi
ja     1168 <ceil+0x68>
cmp    $0x3fe,%esi
ja     1188 <ceil+0x88>
addsd  0xec9(%rip),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1220 <ceil+0x120>
test   %rdi,%rdi
js     11fb <ceil+0xfb>
xor    %eax,%eax
mov    %edx,%ecx
or     %edi,%ecx
mov    $0x3ff00000,%ecx
cmove  %edx,%ecx
cmove  %edi,%eax
mov    %ecx,%edx
mov    %eax,%edi
jmp    1220 <ceil+0x120>
cmp    $0x433,%esi
jb     11ca <ceil+0xca>
movapd %xmm0,%xmm1
cmp    $0x400,%eax
je     1183 <ceil+0x83>
movsd  0xe85(%rip),%xmm1
addsd  %xmm1,%xmm0
ret
mov    $0xfffff,%esi
mov    %eax,%ecx
shr    %cl,%esi
and    %edx,%esi
or     %edi,%esi
je     122e <ceil+0x12e>
addsd  0xe5d(%rip),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1220 <ceil+0x120>
mov    $0x100000,%esi
mov    %eax,%ecx
shr    %cl,%esi
xor    %edi,%edi
test   %edx,%edx
cmovle %edi,%esi
add    %edx,%esi
mov    $0xfff00000,%edx
sar    %cl,%edx
and    %esi,%edx
jmp    1220 <ceil+0x120>
lea    -0x13(%rsi),%ecx
mov    $0xffffffff,%r8d
shr    %cl,%r8d
test   %edi,%r8d
je     122e <ceil+0x12e>
addsd  0xe1d(%rip),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1220 <ceil+0x120>
test   %edx,%edx
jle    1217 <ceil+0x117>
cmp    $0x14,%eax
jne    1204 <ceil+0x104>
add    $0x1,%edx
jmp    1217 <ceil+0x117>
xor    %edi,%edi
mov    $0x80000000,%edx
jmp    1220 <ceil+0x120>
mov    $0x33,%cl
sub    %sil,%cl
mov    $0x1,%eax
shl    %cl,%eax
add    %edi,%eax
adc    $0x0,%edx
mov    %eax,%edi
not    %r8d
and    %edi,%r8d
mov    %r8d,%edi
shl    $0x20,%rdx
mov    %edi,%eax
or     %rdx,%rax
movq   %rax,%xmm0
ret

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

LLVM IR (-O2):
define double @ceil(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 slt i64 %2, 0
br i1 %16, label %70, label %17
17: ; preds = %15
%18 = or i64 %3, %2
%19 = trunc i64 %18 to i32
%20 = icmp eq i32 %19, 0
%21 = select i1 %20, i32 %4, i32 1072693248
%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 sgt i32 %4, 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 sgt i32 %4, 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 nuw 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 nuw 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 ], [ -2147483648, %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
}