double ieee_nextafter(double x, double y) {
  int hx, hy, ix, iy;
  unsigned lx, ly;
  hx = *(int *)&x;
  lx = *(1 + (int *)&x);
  hy = *(int *)&y;
  ly = *(1 + (int *)&y);
  ix = hx & 0x7fffffff;
  iy = hy & 0x7fffffff;
  if (((ix >= 0x7ff00000) && ((ix - 0x7ff00000) | lx) != 0) ||
      ((iy >= 0x7ff00000) && ((iy - 0x7ff00000) | ly) != 0))
    return x + y;
  if (x == y)
    return x;
  if ((ix | lx) == 0) {
    *(int *)&x = hy & 0x80000000;
    *(1 + (int *)&x) = 1;
    y = x * x;
    if (y == x)
      return y;
    else
      return x;
  }
  if (hx >= 0) {
    if (hx > hy || ((hx == hy) && (lx > ly))) {
      if (lx == 0)
        hx -= 1;
      lx -= 1;
    } else {
      lx += 1;
      if (lx == 0)
        hx += 1;
    }
  } else {
    if (hy >= 0 || hx > hy || ((hx == hy) && (lx > ly))) {
      if (lx == 0)
        hx -= 1;
      lx -= 1;
    } else {
      lx += 1;
      if (lx == 0)
        hx += 1;
    }
  }
  hy = hx & 0x7ff00000;
  if (hy >= 0x7ff00000)
    return x + x;
  if (hy < 0x00100000) {
    y = x * x;
    if (y != x) {
      *(int *)&y = hx;
      *(1 + (int *)&y) = lx;
      return y;
    }
  }
  *(int *)&x = hx;
  *(1 + (int *)&x) = lx;
  return x;
}


x86 assembly (clang, -O0):
<ieee_nextafter>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
mov    -0x10(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0xc(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x18(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x1c(%rbp),%eax
and    $0x7fffffff,%eax
mov    %eax,-0x24(%rbp)
mov    -0x20(%rbp),%eax
and    $0x7fffffff,%eax
mov    %eax,-0x28(%rbp)
cmpl   $0x7ff00000,-0x24(%rbp)
jl     115d <ieee_nextafter+0x5d>
mov    -0x24(%rbp),%eax
sub    $0x7ff00000,%eax
or     -0x2c(%rbp),%eax
cmp    $0x0,%eax
jne    117e <ieee_nextafter+0x7e>
cmpl   $0x7ff00000,-0x28(%rbp)
jl     1192 <ieee_nextafter+0x92>
mov    -0x28(%rbp),%eax
sub    $0x7ff00000,%eax
or     -0x30(%rbp),%eax
cmp    $0x0,%eax
je     1192 <ieee_nextafter+0x92>
movsd  -0x10(%rbp),%xmm0
addsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    139e <ieee_nextafter+0x29e>
movsd  -0x10(%rbp),%xmm0
ucomisd -0x18(%rbp),%xmm0
jne    11b7 <ieee_nextafter+0xb7>
jp     11b7 <ieee_nextafter+0xb7>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    139e <ieee_nextafter+0x29e>
mov    -0x24(%rbp),%eax
or     -0x2c(%rbp),%eax
cmp    $0x0,%eax
jne    121b <ieee_nextafter+0x11b>
mov    -0x20(%rbp),%eax
and    $0x80000000,%eax
mov    %eax,-0x10(%rbp)
movl   $0x1,-0xc(%rbp)
movsd  -0x10(%rbp),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm0
ucomisd -0x10(%rbp),%xmm0
jne    120c <ieee_nextafter+0x10c>
jp     120c <ieee_nextafter+0x10c>
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    139e <ieee_nextafter+0x29e>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    139e <ieee_nextafter+0x29e>
cmpl   $0x0,-0x1c(%rbp)
jl     1290 <ieee_nextafter+0x190>
mov    -0x1c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jg     1249 <ieee_nextafter+0x149>
mov    -0x1c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jne    126a <ieee_nextafter+0x16a>
mov    -0x2c(%rbp),%eax
cmp    -0x30(%rbp),%eax
jbe    126a <ieee_nextafter+0x16a>
cmpl   $0x0,-0x2c(%rbp)
jne    125c <ieee_nextafter+0x15c>
mov    -0x1c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x2c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    128b <ieee_nextafter+0x18b>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x0,-0x2c(%rbp)
jne    1286 <ieee_nextafter+0x186>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    128b <ieee_nextafter+0x18b>
jmp    1305 <ieee_nextafter+0x205>
cmpl   $0x0,-0x20(%rbp)
jge    12be <ieee_nextafter+0x1be>
mov    -0x1c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jg     12be <ieee_nextafter+0x1be>
mov    -0x1c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jne    12df <ieee_nextafter+0x1df>
mov    -0x2c(%rbp),%eax
cmp    -0x30(%rbp),%eax
jbe    12df <ieee_nextafter+0x1df>
cmpl   $0x0,-0x2c(%rbp)
jne    12d1 <ieee_nextafter+0x1d1>
mov    -0x1c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x2c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1300 <ieee_nextafter+0x200>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x0,-0x2c(%rbp)
jne    12fb <ieee_nextafter+0x1fb>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1300 <ieee_nextafter+0x200>
jmp    1305 <ieee_nextafter+0x205>
mov    -0x1c(%rbp),%eax
and    $0x7ff00000,%eax
mov    %eax,-0x20(%rbp)
cmpl   $0x7ff00000,-0x20(%rbp)
jl     1331 <ieee_nextafter+0x231>
movsd  -0x10(%rbp),%xmm0
addsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    139e <ieee_nextafter+0x29e>
cmpl   $0x100000,-0x20(%rbp)
jge    1388 <ieee_nextafter+0x288>
movsd  -0x10(%rbp),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm0
ucomisd -0x10(%rbp),%xmm0
jne    1368 <ieee_nextafter+0x268>
jp     1368 <ieee_nextafter+0x268>
jmp    1383 <ieee_nextafter+0x283>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x14(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    139e <ieee_nextafter+0x29e>
jmp    1388 <ieee_nextafter+0x288>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x10(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0xc(%rbp)
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<ieee_nextafter>:
movq   %xmm0,%r8
mov    %r8,%rax
shr    $0x20,%rax
mov    %r8d,%edi
and    $0x7fffffff,%edi
cmp    $0x7ff00000,%edi
jb     1127 <ieee_nextafter+0x27>
lea    -0x7ff00000(%rdi),%ecx
or     %eax,%ecx
jne    114d <ieee_nextafter+0x4d>
movq   %xmm1,%rsi
mov    %rsi,%rdx
shr    $0x20,%rdx
mov    %esi,%ecx
and    $0x7fffffff,%ecx
cmp    $0x7ff00000,%ecx
jb     1153 <ieee_nextafter+0x53>
add    $0x80100000,%ecx
or     %edx,%ecx
je     1153 <ieee_nextafter+0x53>
addsd  %xmm1,%xmm0
jmp    115b <ieee_nextafter+0x5b>
ucomisd %xmm1,%xmm0
jne    115c <ieee_nextafter+0x5c>
jp     115c <ieee_nextafter+0x5c>
ret
or     %eax,%edi
je     117b <ieee_nextafter+0x7b>
test   %r8d,%r8d
jns    1169 <ieee_nextafter+0x69>
test   %esi,%esi
jns    1174 <ieee_nextafter+0x74>
cmp    %esi,%r8d
jg     1174 <ieee_nextafter+0x74>
jne    11b1 <ieee_nextafter+0xb1>
cmp    %edx,%eax
jbe    11b1 <ieee_nextafter+0xb1>
sub    $0x1,%eax
sbb    %edx,%edx
jmp    11b8 <ieee_nextafter+0xb8>
and    $0x80000000,%esi
movabs $0x100000000,%rax
or     %rsi,%rax
movq   %rax,%xmm1
movdqa %xmm1,%xmm2
mulsd  %xmm1,%xmm2
movapd %xmm2,%xmm0
cmpeqsd %xmm1,%xmm0
andpd  %xmm0,%xmm2
andnpd %xmm1,%xmm0
orpd   %xmm2,%xmm0
ret
xor    %edx,%edx
inc    %eax
sete   %dl
add    %r8d,%edx
mov    %edx,%ecx
and    $0x7ff00000,%ecx
je     11d2 <ieee_nextafter+0xd2>
cmp    $0x7ff00000,%ecx
jne    11de <ieee_nextafter+0xde>
addsd  %xmm0,%xmm0
ret
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
ucomisd %xmm0,%xmm1
mov    %edx,%ecx
shl    $0x20,%rax
or     %rcx,%rax
movq   %rax,%xmm0
ret

LLVM IR (-O0):
define double @ieee_nextafter(double %0, double %1) {
%3 = alloca double
%4 = alloca double
%5 = alloca double
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
store double %0, double* %4
store double %1, double* %5
%12 = bitcast double* %4 to i32*
%13 = load i32, i32* %12
store i32 %13, i32* %6
%14 = bitcast double* %4 to i32*
%15 = getelementptr inbounds i32, i32* %14, i64 1
%16 = load i32, i32* %15
store i32 %16, i32* %10
%17 = bitcast double* %5 to i32*
%18 = load i32, i32* %17
store i32 %18, i32* %7
%19 = bitcast double* %5 to i32*
%20 = getelementptr inbounds i32, i32* %19, i64 1
%21 = load i32, i32* %20
store i32 %21, i32* %11
%22 = load i32, i32* %6
%23 = and i32 %22, 2147483647
store i32 %23, i32* %8
%24 = load i32, i32* %7
%25 = and i32 %24, 2147483647
store i32 %25, i32* %9
%26 = load i32, i32* %8
%27 = icmp sge i32 %26, 2146435072
br i1 %27, label %28, label %34
28: ; preds = %2
%29 = load i32, i32* %8
%30 = sub nsw i32 %29, 2146435072
%31 = load i32, i32* %10
%32 = or i32 %30, %31
%33 = icmp ne i32 %32, 0
br i1 %33, label %43, label %34
34: ; preds = %28, %2
%35 = load i32, i32* %9
%36 = icmp sge i32 %35, 2146435072
br i1 %36, label %37, label %47
37: ; preds = %34
%38 = load i32, i32* %9
%39 = sub nsw i32 %38, 2146435072
%40 = load i32, i32* %11
%41 = or i32 %39, %40
%42 = icmp ne i32 %41, 0
br i1 %42, label %43, label %47
43: ; preds = %37, %28
%44 = load double, double* %4
%45 = load double, double* %5
%46 = fadd double %44, %45
store double %46, double* %3
br label %176
47: ; preds = %37, %34
%48 = load double, double* %4
%49 = load double, double* %5
%50 = fcmp oeq double %48, %49
br i1 %50, label %51, label %53
51: ; preds = %47
%52 = load double, double* %4
store double %52, double* %3
br label %176
53: ; preds = %47
%54 = load i32, i32* %8
%55 = load i32, i32* %10
%56 = or i32 %54, %55
%57 = icmp eq i32 %56, 0
br i1 %57, label %58, label %74
58: ; preds = %53
%59 = load i32, i32* %7
%60 = and i32 %59, -2147483648
%61 = bitcast double* %4 to i32*
store i32 %60, i32* %61
%62 = bitcast double* %4 to i32*
%63 = getelementptr inbounds i32, i32* %62, i64 1
store i32 1, i32* %63
%64 = load double, double* %4
%65 = load double, double* %4
%66 = fmul double %64, %65
store double %66, double* %5
%67 = load double, double* %5
%68 = load double, double* %4
%69 = fcmp oeq double %67, %68
br i1 %69, label %70, label %72
70: ; preds = %58
%71 = load double, double* %5
store double %71, double* %3
br label %176
72: ; preds = %58
%73 = load double, double* %4
store double %73, double* %3
br label %176
74: ; preds = %53
%75 = load i32, i32* %6
%76 = icmp sge i32 %75, 0
br i1 %76, label %77, label %108
77: ; preds = %74
%78 = load i32, i32* %6
%79 = load i32, i32* %7
%80 = icmp sgt i32 %78, %79
br i1 %80, label %89, label %81
81: ; preds = %77
%82 = load i32, i32* %6
%83 = load i32, i32* %7
%84 = icmp eq i32 %82, %83
br i1 %84, label %85, label %98
85: ; preds = %81
%86 = load i32, i32* %10
%87 = load i32, i32* %11
%88 = icmp ugt i32 %86, %87
br i1 %88, label %89, label %98
89: ; preds = %85, %77
%90 = load i32, i32* %10
%91 = icmp eq i32 %90, 0
br i1 %91, label %92, label %95
92: ; preds = %89
%93 = load i32, i32* %6
%94 = sub nsw i32 %93, 1
store i32 %94, i32* %6
br label %95
95: ; preds = %92, %89
%96 = load i32, i32* %10
%97 = sub i32 %96, 1
store i32 %97, i32* %10
br label %107
98: ; preds = %85, %81
%99 = load i32, i32* %10
%100 = add i32 %99, 1
store i32 %100, i32* %10
%101 = load i32, i32* %10
%102 = icmp eq i32 %101, 0
br i1 %102, label %103, label %106
103: ; preds = %98
%104 = load i32, i32* %6
%105 = add nsw i32 %104, 1
store i32 %105, i32* %6
br label %106
106: ; preds = %103, %98
br label %107
107: ; preds = %106, %95
br label %142
108: ; preds = %74
%109 = load i32, i32* %7
%110 = icmp sge i32 %109, 0
br i1 %110, label %123, label %111
111: ; preds = %108
%112 = load i32, i32* %6
%113 = load i32, i32* %7
%114 = icmp sgt i32 %112, %113
br i1 %114, label %123, label %115
115: ; preds = %111
%116 = load i32, i32* %6
%117 = load i32, i32* %7
%118 = icmp eq i32 %116, %117
br i1 %118, label %119, label %132
119: ; preds = %115
%120 = load i32, i32* %10
%121 = load i32, i32* %11
%122 = icmp ugt i32 %120, %121
br i1 %122, label %123, label %132
123: ; preds = %119, %111, %108
%124 = load i32, i32* %10
%125 = icmp eq i32 %124, 0
br i1 %125, label %126, label %129
126: ; preds = %123
%127 = load i32, i32* %6
%128 = sub nsw i32 %127, 1
store i32 %128, i32* %6
br label %129
129: ; preds = %126, %123
%130 = load i32, i32* %10
%131 = sub i32 %130, 1
store i32 %131, i32* %10
br label %141
132: ; preds = %119, %115
%133 = load i32, i32* %10
%134 = add i32 %133, 1
store i32 %134, i32* %10
%135 = load i32, i32* %10
%136 = icmp eq i32 %135, 0
br i1 %136, label %137, label %140
137: ; preds = %132
%138 = load i32, i32* %6
%139 = add nsw i32 %138, 1
store i32 %139, i32* %6
br label %140
140: ; preds = %137, %132
br label %141
141: ; preds = %140, %129
br label %142
142: ; preds = %141, %107
%143 = load i32, i32* %6
%144 = and i32 %143, 2146435072
store i32 %144, i32* %7
%145 = load i32, i32* %7
%146 = icmp sge i32 %145, 2146435072
br i1 %146, label %147, label %151
147: ; preds = %142
%148 = load double, double* %4
%149 = load double, double* %4
%150 = fadd double %148, %149
store double %150, double* %3
br label %176
151: ; preds = %142
%152 = load i32, i32* %7
%153 = icmp slt i32 %152, 1048576
br i1 %153, label %154, label %169
154: ; preds = %151
%155 = load double, double* %4
%156 = load double, double* %4
%157 = fmul double %155, %156
store double %157, double* %5
%158 = load double, double* %5
%159 = load double, double* %4
%160 = fcmp une double %158, %159
br i1 %160, label %161, label %168
161: ; preds = %154
%162 = load i32, i32* %6
%163 = bitcast double* %5 to i32*
store i32 %162, i32* %163
%164 = load i32, i32* %10
%165 = bitcast double* %5 to i32*
%166 = getelementptr inbounds i32, i32* %165, i64 1
store i32 %164, i32* %166
%167 = load double, double* %5
store double %167, double* %3
br label %176
168: ; preds = %154
br label %169
169: ; preds = %168, %151
%170 = load i32, i32* %6
%171 = bitcast double* %4 to i32*
store i32 %170, i32* %171
%172 = load i32, i32* %10
%173 = bitcast double* %4 to i32*
%174 = getelementptr inbounds i32, i32* %173, i64 1
store i32 %172, i32* %174
%175 = load double, double* %4
store double %175, double* %3
br label %176
176: ; preds = %169, %161, %147, %72, %70, %51, %43
%177 = load double, double* %3
ret double %177
}

LLVM IR (-O2):
define double @ieee_nextafter(double %0, double %1) {
%3 = bitcast double %0 to i64
%4 = trunc i64 %3 to i32
%5 = lshr i64 %3, 32
%6 = trunc i64 %5 to i32
%7 = bitcast double %1 to i64
%8 = trunc i64 %7 to i32
%9 = lshr i64 %7, 32
%10 = trunc i64 %9 to i32
%11 = and i32 %4, 2147483647
%12 = and i32 %8, 2147483647
%13 = icmp ugt i32 %11, 2146435071
br i1 %13, label %14, label %18
14: ; preds = %2
%15 = add nsw i32 %11, -2146435072
%16 = or i32 %15, %6
%17 = icmp eq i32 %16, 0
br i1 %17, label %18, label %24
18: ; preds = %14, %2
%19 = icmp ugt i32 %12, 2146435071
br i1 %19, label %20, label %26
20: ; preds = %18
%21 = add nsw i32 %12, -2146435072
%22 = or i32 %21, %10
%23 = icmp eq i32 %22, 0
br i1 %23, label %26, label %24
24: ; preds = %20, %14
%25 = fadd double %0, %1
br label %92
26: ; preds = %20, %18
%27 = fcmp oeq double %0, %1
br i1 %27, label %92, label %28
28: ; preds = %26
%29 = or i32 %11, %6
%30 = icmp eq i32 %29, 0
br i1 %30, label %31, label %38
31: ; preds = %28
%32 = and i64 %7, 2147483648
%33 = or i64 %32, 4294967296
%34 = bitcast i64 %33 to double
%35 = fmul double %34, %34
%36 = fcmp oeq double %35, %34
%37 = select i1 %36, double %35, double %34
br label %92
38: ; preds = %28
%39 = icmp sgt i32 %4, -1
br i1 %39, label %40, label %54
40: ; preds = %38
%41 = icmp sgt i32 %4, %8
br i1 %41, label %46, label %42
42: ; preds = %40
%43 = icmp eq i32 %4, %8
%44 = icmp ugt i32 %6, %10
%45 = and i1 %43, %44
br i1 %45, label %46, label %50
46: ; preds = %42, %40
%47 = icmp eq i32 %6, 0
%48 = sext i1 %47 to i32
%49 = add i32 %6, -1
br label %70
50: ; preds = %42
%51 = add i32 %6, 1
%52 = icmp eq i32 %51, 0
%53 = zext i1 %52 to i32
br label %70
54: ; preds = %38
%55 = icmp sgt i32 %8, -1
%56 = icmp sgt i32 %4, %8
%57 = or i1 %55, %56
br i1 %57, label %62, label %58
58: ; preds = %54
%59 = icmp eq i32 %4, %8
%60 = icmp ugt i32 %6, %10
%61 = and i1 %59, %60
br i1 %61, label %62, label %66
62: ; preds = %58, %54
%63 = icmp eq i32 %6, 0
%64 = sext i1 %63 to i32
%65 = add i32 %6, -1
br label %70
66: ; preds = %58
%67 = add i32 %6, 1
%68 = icmp eq i32 %67, 0
%69 = zext i1 %68 to i32
br label %70
70: ; preds = %66, %50, %62, %46
%71 = phi i32 [ %48, %46 ], [ %64, %62 ], [ %53, %50 ], [ %69, %66 ]
%72 = phi i32 [ %49, %46 ], [ %65, %62 ], [ %51, %50 ], [ %67, %66 ]
%73 = add nsw i32 %71, %4
%74 = and i32 %73, 2146435072
switch i32 %74, label %86 [
i32 2146435072, label %75
i32 0, label %77
]
75: ; preds = %70
%76 = fadd double %0, %0
br label %92
77: ; preds = %70
%78 = fmul double %0, %0
%79 = fcmp une double %78, %0
br i1 %79, label %80, label %86
80: ; preds = %77
%81 = zext i32 %73 to i64
%82 = zext i32 %72 to i64
%83 = shl nuw i64 %82, 32
%84 = or i64 %83, %81
%85 = bitcast i64 %84 to double
br label %92
86: ; preds = %70, %77
%87 = zext i32 %73 to i64
%88 = zext i32 %72 to i64
%89 = shl nuw i64 %88, 32
%90 = or i64 %89, %87
%91 = bitcast i64 %90 to double
br label %92
92: ; preds = %31, %26, %86, %80, %75, %24
%93 = phi double [ %25, %24 ], [ %76, %75 ], [ %85, %80 ], [ %91, %86 ], [ %0, %26 ], [ %37, %31 ]
ret double %93
}