long double __rinfl(long double x) {
  union {
    long 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.0l : 1.0l;
  long mask = 0x000FFFFFFFFFFFFF >> exp;
  if (u.i & mask) {
    u.i &= ~mask;
    u.i += mask + 1;
  }
  return u.f;
}


x86 assembly (clang, -O0):
<__rinfl>:
push   %rbp
mov    %rsp,%rbp
fldt   0x10(%rbp)
fstpt  -0x20(%rbp)
fldt   -0x20(%rbp)
fstpt  -0x30(%rbp)
mov    -0x2c(%rbp),%eax
sar    $0x14,%eax
and    $0x7ff,%eax
sub    $0x3ff,%eax
mov    %eax,-0x34(%rbp)
cmpl   $0x35,-0x34(%rbp)
jge    113a <__rinfl+0x3a>
cmpl   $0xfffffc01,-0x34(%rbp)
jne    1145 <__rinfl+0x45>
fldt   -0x20(%rbp)
fstpt  -0x10(%rbp)
jmp    11b9 <__rinfl+0xb9>
cmpl   $0x0,-0x34(%rbp)
jge    116b <__rinfl+0x6b>
mov    -0x2c(%rbp),%eax
test   $0x80000000,%eax
fld1
fld    %st(0)
fchs
fxch   %st(1)
fcmovne %st(1),%st
fstp   %st(1)
fstpt  -0x10(%rbp)
jmp    11b9 <__rinfl+0xb9>
mov    -0x34(%rbp),%eax
mov    %eax,%ecx
movabs $0xfffffffffffff,%rax
sar    %cl,%rax
mov    %rax,-0x40(%rbp)
mov    -0x30(%rbp),%rax
and    -0x40(%rbp),%rax
cmp    $0x0,%rax
je     11b3 <__rinfl+0xb3>
mov    -0x40(%rbp),%rax
xor    $0xffffffffffffffff,%rax
and    -0x30(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
add    -0x30(%rbp),%rax
mov    %rax,-0x30(%rbp)
fldt   -0x30(%rbp)
fstpt  -0x10(%rbp)
fldt   -0x10(%rbp)
pop    %rbp
ret

x86 assembly (clang, -O2):
<__rinfl>:
fldt   0x8(%rsp)
fld    %st(0)
fstpt  -0x18(%rsp)
mov    -0x18(%rsp),%rax
mov    %rax,%rcx
shr    $0x34,%rcx
and    $0x7ff,%ecx
lea    -0x434(%rcx),%edx
cmp    $0xfffffbcd,%edx
jb     1190 <__rinfl+0x90>
cmp    $0x3fe,%ecx
ja     114a <__rinfl+0x4a>
fstp   %st(0)
shr    $0x20,%rax
test   %eax,%eax
setns  %al
test   %al,%al
fld1
fld    %st(0)
fchs
fcmovne %st(1),%st
fstp   %st(1)
ret
add    $0xfffffc01,%ecx
movabs $0xfffffffffffff,%rdx
shr    %cl,%rdx
test   %rax,%rdx
je     1190 <__rinfl+0x90>
fstp   %st(0)
mov    -0x10(%rsp),%esi
mov    %ecx,%ecx
movabs $0xfff0000000000000,%rdi
sar    %cl,%rdi
and    %rax,%rdi
lea    (%rdx,%rdi,1),%rax
add    $0x1,%rax
mov    %si,-0x20(%rsp)
mov    %rax,-0x28(%rsp)
fldt   -0x28(%rsp)
ret

LLVM IR (-O0):
%union.anon = type { x86_fp80 }
%struct.anon = type { i32, i32 }
define x86_fp80 @__rinfl(x86_fp80 %0) {
%2 = alloca x86_fp80
%3 = alloca x86_fp80
%4 = alloca %union.anon
%5 = alloca i32
%6 = alloca i64
store x86_fp80 %0, x86_fp80* %3
%7 = bitcast %union.anon* %4 to x86_fp80*
%8 = load x86_fp80, x86_fp80* %3
store x86_fp80 %8, x86_fp80* %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 x86_fp80, x86_fp80* %3
store x86_fp80 %21, x86_fp80* %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, x86_fp80 0xKBFFF8000000000000000, x86_fp80 0xK3FFF8000000000000000
store x86_fp80 %32, x86_fp80* %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 x86_fp80*
%55 = load x86_fp80, x86_fp80* %54
store x86_fp80 %55, x86_fp80* %2
br label %56
56: ; preds = %53, %25, %20
%57 = load x86_fp80, x86_fp80* %2
ret x86_fp80 %57
}

LLVM IR (-O2):
define x86_fp80 @__rinfl(x86_fp80 %0) {
%2 = bitcast x86_fp80 %0 to i80
%3 = lshr i80 %2, 32
%4 = trunc i80 %3 to i32
%5 = lshr i32 %4, 20
%6 = and i32 %5, 2047
%7 = add nsw i32 %6, -1023
%8 = add nsw i32 %6, -1076
%9 = icmp ult i32 %8, -1075
br i1 %9, label %30, label %10
10: ; preds = %1
%11 = icmp ult i32 %6, 1023
br i1 %11, label %12, label %15
12: ; preds = %10
%13 = icmp sgt i32 %4, -1
%14 = select i1 %13, x86_fp80 0xK3FFF8000000000000000, x86_fp80 0xKBFFF8000000000000000
br label %30
15: ; preds = %10
%16 = zext i32 %7 to i64
%17 = lshr i64 4503599627370495, %16
%18 = trunc i80 %2 to i64
%19 = and i64 %17, %18
%20 = icmp eq i64 %19, 0
br i1 %20, label %30, label %21
21: ; preds = %15
%22 = ashr i64 -4503599627370496, %16
%23 = and i64 %22, %18
%24 = and i80 %2, -18446744073709551616
%25 = add nuw nsw i64 %17, 1
%26 = add nsw i64 %25, %23
%27 = zext i64 %26 to i80
%28 = or i80 %24, %27
%29 = bitcast i80 %28 to x86_fp80
br label %30
30: ; preds = %15, %21, %1, %12
%31 = phi x86_fp80 [ %14, %12 ], [ %0, %1 ], [ %29, %21 ], [ %0, %15 ]
ret x86_fp80 %31
}