double __rinf(double x) {
  union {
    double f;
    long i;
    struct {
      int i0, i1;
    };
  } u = {x};
  int exp = ((u.i1 >> 20) & 0x7FF) - 1023;
  if (exp >= 53 || exp == -1023)
    return x;
  if (exp < 0)
    return u.i1 & 0x80000000 ? -1.0 : 1.0;
  long mask = 0x000FFFFFFFFFFFFF >> exp;
  if (u.i & mask) {
    u.i &= ~mask;
    u.i += mask + 1;
  }
  return u.f;
}


x86 assembly (clang, -O0):
<__rinf>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x14(%rbp),%eax
sar    $0x14,%eax
and    $0x7ff,%eax
sub    $0x3ff,%eax
mov    %eax,-0x1c(%rbp)
cmpl   $0x35,-0x1c(%rbp)
jge    113d <__rinf+0x3d>
cmpl   $0xfffffc01,-0x1c(%rbp)
jne    114c <__rinf+0x4c>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    11ec <__rinf+0xec>
cmpl   $0x0,-0x1c(%rbp)
jge    119a <__rinf+0x9a>
mov    -0x14(%rbp),%eax
and    $0x80000000,%eax
movsd  0xe9a(%rip),%xmm0
movsd  0xe9a(%rip),%xmm1
movsd  %xmm1,-0x38(%rbp)
cmp    $0x0,%eax
movsd  %xmm0,-0x30(%rbp)
jne    118b <__rinf+0x8b>
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    11ec <__rinf+0xec>
mov    -0x1c(%rbp),%eax
mov    %eax,%ecx
movabs $0xfffffffffffff,%rax
sar    %cl,%rax
mov    %rax,-0x28(%rbp)
mov    -0x18(%rbp),%rax
and    -0x28(%rbp),%rax
cmp    $0x0,%rax
je     11e2 <__rinf+0xe2>
mov    -0x28(%rbp),%rax
xor    $0xffffffffffffffff,%rax
and    -0x18(%rbp),%rax
mov    %rax,-0x18(%rbp)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
add    -0x18(%rbp),%rax
mov    %rax,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<__rinf>:
movq   %xmm0,%rax
mov    %rax,%rcx
shr    $0x34,%rcx
and    $0x7ff,%ecx
lea    -0x434(%rcx),%edx
cmp    $0xfffffbcd,%edx
jb     1170 <__rinf+0x70>
cmp    $0x3fe,%ecx
ja     1139 <__rinf+0x39>
pand   0xed0(%rip),%xmm0
por    0xed8(%rip),%xmm0
ret
add    $0xfffffc01,%ecx
movabs $0xfffffffffffff,%rdx
shr    %cl,%rdx
test   %rax,%rdx
je     1170 <__rinf+0x70>
mov    %ecx,%ecx
movabs $0xfff0000000000000,%rsi
sar    %cl,%rsi
and    %rax,%rsi
lea    (%rdx,%rsi,1),%rax
add    $0x1,%rax
movq   %rax,%xmm0
ret

LLVM IR (-O0):
%union.anon = type { double }
%struct.anon = type { i32, i32 }
define double @__rinf(double %0) {
%2 = alloca double
%3 = alloca double
%4 = alloca %union.anon
%5 = alloca i32
%6 = alloca i64
store double %0, double* %3
%7 = bitcast %union.anon* %4 to double*
%8 = load double, double* %3
store double %8, double* %7
%9 = bitcast %union.anon* %4 to %struct.anon*
%10 = getelementptr inbounds %struct.anon, %struct.anon* %9, i32 0, i32 1
%11 = load i32, i32* %10
%12 = ashr i32 %11, 20
%13 = and i32 %12, 2047
%14 = sub nsw i32 %13, 1023
store i32 %14, i32* %5
%15 = load i32, i32* %5
%16 = icmp sge i32 %15, 53
br i1 %16, label %20, label %17
17: ; preds = %1
%18 = load i32, i32* %5
%19 = icmp eq i32 %18, -1023
br i1 %19, label %20, label %22
20: ; preds = %17, %1
%21 = load double, double* %3
store double %21, double* %2
br label %56
22: ; preds = %17
%23 = load i32, i32* %5
%24 = icmp slt i32 %23, 0
br i1 %24, label %25, label %33
25: ; preds = %22
%26 = bitcast %union.anon* %4 to %struct.anon*
%27 = getelementptr inbounds %struct.anon, %struct.anon* %26, i32 0, i32 1
%28 = load i32, i32* %27
%29 = and i32 %28, -2147483648
%30 = icmp ne i32 %29, 0
%31 = zext i1 %30 to i64
%32 = select i1 %30, double -1.000000e+00, double 1.000000e+00
store double %32, double* %2
br label %56
33: ; preds = %22
%34 = load i32, i32* %5
%35 = zext i32 %34 to i64
%36 = ashr i64 4503599627370495, %35
store i64 %36, i64* %6
%37 = bitcast %union.anon* %4 to i64*
%38 = load i64, i64* %37
%39 = load i64, i64* %6
%40 = and i64 %38, %39
%41 = icmp ne i64 %40, 0
br i1 %41, label %42, label %53
42: ; preds = %33
%43 = load i64, i64* %6
%44 = xor i64 %43, -1
%45 = bitcast %union.anon* %4 to i64*
%46 = load i64, i64* %45
%47 = and i64 %46, %44
store i64 %47, i64* %45
%48 = load i64, i64* %6
%49 = add nsw i64 %48, 1
%50 = bitcast %union.anon* %4 to i64*
%51 = load i64, i64* %50
%52 = add nsw i64 %51, %49
store i64 %52, i64* %50
br label %53
53: ; preds = %42, %33
%54 = bitcast %union.anon* %4 to double*
%55 = load double, double* %54
store double %55, double* %2
br label %56
56: ; preds = %53, %25, %20
%57 = load double, double* %2
ret double %57
}

LLVM IR (-O2):
define double @__rinf(double %0) {
%2 = bitcast double %0 to i64
%3 = lshr i64 %2, 52
%4 = trunc i64 %3 to i32
%5 = and i32 %4, 2047
%6 = add nsw i32 %5, -1023
%7 = add nsw i32 %5, -1076
%8 = icmp ult i32 %7, -1075
br i1 %8, label %24, label %9
9: ; preds = %1
%10 = icmp ult i32 %5, 1023
br i1 %10, label %11, label %13
11: ; preds = %9
%12 = tail call double @llvm.copysign.f64(double 1.000000e+00, double %0)
br label %24
13: ; preds = %9
%14 = zext i32 %6 to i64
%15 = lshr i64 4503599627370495, %14
%16 = and i64 %15, %2
%17 = icmp eq i64 %16, 0
br i1 %17, label %24, label %18
18: ; preds = %13
%19 = ashr i64 -4503599627370496, %14
%20 = and i64 %19, %2
%21 = add nuw nsw i64 %15, 1
%22 = add nsw i64 %21, %20
%23 = bitcast i64 %22 to double
br label %24
24: ; preds = %13, %18, %1, %11
%25 = phi double [ %12, %11 ], [ %0, %1 ], [ %23, %18 ], [ %0, %13 ]
ret double %25
}