long double roundl(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)
    return x;
  if (exp < -1)
    return u.i &= 0x8000000000000000, u.f;
  if (exp == -1)
    return u.i1 & 0x80000000 ? -1.0l : 1.0l;
  if (u.i & (0x0008000000000000 >> exp))
    u.i += 0x0010000000000000 >> exp;
  return u.i &= ~(0x000FFFFFFFFFFFFF >> exp), u.f;
}


x86 assembly (clang, -O0):
<roundl>:
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)
jl     1138 <roundl+0x38>
fldt   -0x20(%rbp)
fstpt  -0x10(%rbp)
jmp    11ed <roundl+0xed>
cmpl   $0xffffffff,-0x34(%rbp)
jge    1162 <roundl+0x62>
mov    -0x30(%rbp),%rax
movabs $0x8000000000000000,%rcx
and    %rcx,%rax
mov    %rax,-0x30(%rbp)
fldt   -0x30(%rbp)
fstpt  -0x10(%rbp)
jmp    11ed <roundl+0xed>
cmpl   $0xffffffff,-0x34(%rbp)
jne    1188 <roundl+0x88>
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    11ed <roundl+0xed>
mov    -0x30(%rbp),%rax
mov    -0x34(%rbp),%ecx
movabs $0x8000000000000,%rdx
sar    %cl,%rdx
mov    %rdx,%rcx
and    %rcx,%rax
cmp    $0x0,%rax
je     11c6 <roundl+0xc6>
mov    -0x34(%rbp),%eax
mov    %eax,%ecx
movabs $0x10000000000000,%rax
sar    %cl,%rax
add    -0x30(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x34(%rbp),%cl
movabs $0xfffffffffffff,%rax
shr    %cl,%rax
mov    %rax,%rcx
not    %rcx
mov    -0x30(%rbp),%rax
and    %rcx,%rax
mov    %rax,-0x30(%rbp)
fldt   -0x30(%rbp)
fstpt  -0x10(%rbp)
fldt   -0x10(%rbp)
pop    %rbp
ret

x86 assembly (clang, -O2):
<roundl>:
fldt   0x8(%rsp)
fld    %st(0)
fstpt  -0x18(%rsp)
mov    -0x18(%rsp),%rax
mov    %rax,%rcx
shr    $0x34,%rcx
and    $0x7ff,%ecx
cmp    $0x433,%ecx
ja     119e <roundl+0x9e>
fstp   %st(0)
mov    -0x10(%rsp),%edx
cmp    $0x3fd,%ecx
ja     114e <roundl+0x4e>
movabs $0x8000000000000000,%rcx
and    %rcx,%rax
mov    %rax,-0x38(%rsp)
mov    %dx,-0x30(%rsp)
fldt   -0x38(%rsp)
ret
add    $0xfffffc01,%ecx
cmp    $0xffffffff,%ecx
je     119f <roundl+0x9f>
mov    %rax,%rsi
shl    %cl,%rsi
mov    %ecx,%ecx
movabs $0x8000000000000,%rdi
test   %rdi,%rsi
je     1180 <roundl+0x80>
movabs $0x10000000000000,%rsi
shr    %cl,%rsi
add    %rsi,%rax
movabs $0xfff0000000000000,%rsi
sar    %cl,%rsi
and    %rax,%rsi
mov    %dx,-0x20(%rsp)
mov    %rsi,-0x28(%rsp)
fldt   -0x28(%rsp)
ret
shr    $0x20,%rax
test   %eax,%eax
setns  %al
test   %al,%al
fld1
fld    %st(0)
fchs
fcmovne %st(1),%st
fstp   %st(1)
ret

LLVM IR (-O0):
%union.anon = type { x86_fp80 }
%struct.anon = type { i32, i32 }
define x86_fp80 @roundl(x86_fp80 %0) {
%2 = alloca x86_fp80
%3 = alloca x86_fp80
%4 = alloca %union.anon
%5 = alloca i32
store x86_fp80 %0, x86_fp80* %3
%6 = bitcast %union.anon* %4 to x86_fp80*
%7 = load x86_fp80, x86_fp80* %3
store x86_fp80 %7, x86_fp80* %6
%8 = bitcast %union.anon* %4 to %struct.anon*
%9 = getelementptr inbounds %struct.anon, %struct.anon* %8, i32 0, i32 1
%10 = load i32, i32* %9
%11 = ashr i32 %10, 20
%12 = and i32 %11, 2047
%13 = sub nsw i32 %12, 1023
store i32 %13, i32* %5
%14 = load i32, i32* %5
%15 = icmp sge i32 %14, 53
br i1 %15, label %16, label %18
16: ; preds = %1
%17 = load x86_fp80, x86_fp80* %3
store x86_fp80 %17, x86_fp80* %2
br label %63
18: ; preds = %1
%19 = load i32, i32* %5
%20 = icmp slt i32 %19, -1
br i1 %20, label %21, label %27
21: ; preds = %18
%22 = bitcast %union.anon* %4 to i64*
%23 = load i64, i64* %22
%24 = and i64 %23, -9223372036854775808
store i64 %24, i64* %22
%25 = bitcast %union.anon* %4 to x86_fp80*
%26 = load x86_fp80, x86_fp80* %25
store x86_fp80 %26, x86_fp80* %2
br label %63
27: ; preds = %18
%28 = load i32, i32* %5
%29 = icmp eq i32 %28, -1
br i1 %29, label %30, label %38
30: ; preds = %27
%31 = bitcast %union.anon* %4 to %struct.anon*
%32 = getelementptr inbounds %struct.anon, %struct.anon* %31, i32 0, i32 1
%33 = load i32, i32* %32
%34 = and i32 %33, -2147483648
%35 = icmp ne i32 %34, 0
%36 = zext i1 %35 to i64
%37 = select i1 %35, x86_fp80 0xKBFFF8000000000000000, x86_fp80 0xK3FFF8000000000000000
store x86_fp80 %37, x86_fp80* %2
br label %63
38: ; preds = %27
%39 = bitcast %union.anon* %4 to i64*
%40 = load i64, i64* %39
%41 = load i32, i32* %5
%42 = zext i32 %41 to i64
%43 = ashr i64 2251799813685248, %42
%44 = and i64 %40, %43
%45 = icmp ne i64 %44, 0
br i1 %45, label %46, label %53
46: ; preds = %38
%47 = load i32, i32* %5
%48 = zext i32 %47 to i64
%49 = ashr i64 4503599627370496, %48
%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 = %46, %38
%54 = load i32, i32* %5
%55 = zext i32 %54 to i64
%56 = ashr i64 4503599627370495, %55
%57 = xor i64 %56, -1
%58 = bitcast %union.anon* %4 to i64*
%59 = load i64, i64* %58
%60 = and i64 %59, %57
store i64 %60, i64* %58
%61 = bitcast %union.anon* %4 to x86_fp80*
%62 = load x86_fp80, x86_fp80* %61
store x86_fp80 %62, x86_fp80* %2
br label %63
63: ; preds = %53, %30, %21, %16
%64 = load x86_fp80, x86_fp80* %2
ret x86_fp80 %64
}

LLVM IR (-O2):
define x86_fp80 @roundl(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 = icmp ugt i32 %6, 1075
br i1 %8, label %40, label %9
9: ; preds = %1
%10 = icmp ult i32 %6, 1022
br i1 %10, label %11, label %14
11: ; preds = %9
%12 = and i80 %2, -9223372036854775808
%13 = bitcast i80 %12 to x86_fp80
br label %40
14: ; preds = %9
%15 = icmp eq i32 %7, -1
br i1 %15, label %16, label %19
16: ; preds = %14
%17 = icmp sgt i32 %4, -1
%18 = select i1 %17, x86_fp80 0xK3FFF8000000000000000, x86_fp80 0xKBFFF8000000000000000
br label %40
19: ; preds = %14
%20 = trunc i80 %2 to i64
%21 = zext i32 %7 to i64
%22 = lshr i64 2251799813685248, %21
%23 = and i64 %22, %20
%24 = icmp eq i64 %23, 0
br i1 %24, label %31, label %25
25: ; preds = %19
%26 = lshr i64 4503599627370496, %21
%27 = add nsw i64 %26, %20
%28 = zext i64 %27 to i80
%29 = and i80 %2, -18446744073709551616
%30 = or i80 %29, %28
br label %31
31: ; preds = %25, %19
%32 = phi i64 [ %27, %25 ], [ %20, %19 ]
%33 = phi i80 [ %30, %25 ], [ %2, %19 ]
%34 = ashr i64 -4503599627370496, %21
%35 = and i64 %34, %32
%36 = zext i64 %35 to i80
%37 = and i80 %33, -18446744073709551616
%38 = or i80 %37, %36
%39 = bitcast i80 %38 to x86_fp80
br label %40
40: ; preds = %1, %31, %16, %11
%41 = phi x86_fp80 [ %13, %11 ], [ %18, %16 ], [ %39, %31 ], [ %0, %1 ]
ret x86_fp80 %41
}