long double modfd(long double x, long double *i) {
  long double y;
  long double fract;
  long double abs_x;
  long int exp;
  long int bitsToClear;
  long int *yPtr = (long int *)&y;
  y = x;
  if (x < 0.0L) {
    y = -y;
  }
  abs_x = y;
  exp = (yPtr[1] & 0x7ff00000) >> 20;
  exp = exp - 1023;
  if (exp < 0) {
    *i = 0.0L;
    fract = x;
    return fract;
  }
  if (exp > 52) {
    *i = x;
    return 0.0L;
  }
  bitsToClear = 52 - exp;
  if (bitsToClear > 31) {
    bitsToClear = bitsToClear - 32;
    yPtr[0] = 0;
    yPtr[1] = yPtr[1] & (0xffffffff << bitsToClear);
  } else {
    yPtr[0] = yPtr[0] & (0xffffffff << bitsToClear);
  }
  fract = ((long double)(abs_x) - (long double)(y));
  if (x < 0.0L) {
    y = -y;
    fract = -fract;
  }
  *i = y;
  return fract;
}


x86 assembly (clang, -O0):
<modfd>:
push   %rbp
mov    %rsp,%rbp
sub    $0x10,%rsp
fldt   0x10(%rbp)
fstpt  -0x20(%rbp)
mov    %rdi,-0x28(%rbp)
lea    -0x40(%rbp),%rax
mov    %rax,-0x78(%rbp)
fldt   -0x20(%rbp)
fstpt  -0x40(%rbp)
fldt   -0x20(%rbp)
fldz
fucomip %st(1),%st
fstp   %st(0)
jbe    113c <modfd+0x3c>
jmp    1134 <modfd+0x34>
fldt   -0x40(%rbp)
fchs
fstpt  -0x40(%rbp)
fldt   -0x40(%rbp)
fstpt  -0x60(%rbp)
mov    -0x78(%rbp),%rax
mov    0x8(%rax),%rax
and    $0x7ff00000,%rax
sar    $0x14,%rax
mov    %rax,-0x68(%rbp)
mov    -0x68(%rbp),%rax
sub    $0x3ff,%rax
mov    %rax,-0x68(%rbp)
cmpq   $0x0,-0x68(%rbp)
jge    118a <modfd+0x8a>
mov    -0x28(%rbp),%rax
fldz
fstpt  (%rax)
fldt   -0x20(%rbp)
fstpt  -0x50(%rbp)
fldt   -0x50(%rbp)
fstpt  -0x10(%rbp)
jmp    1272 <modfd+0x172>
cmpq   $0x34,-0x68(%rbp)
jle    11a8 <modfd+0xa8>
fldt   -0x20(%rbp)
mov    -0x28(%rbp),%rax
fstpt  (%rax)
fldz
fstpt  -0x10(%rbp)
jmp    1272 <modfd+0x172>
mov    $0x34,%eax
sub    -0x68(%rbp),%rax
mov    %rax,-0x70(%rbp)
cmpq   $0x1f,-0x70(%rbp)
jle    1206 <modfd+0x106>
mov    -0x70(%rbp),%rax
sub    $0x20,%rax
mov    %rax,-0x70(%rbp)
mov    -0x78(%rbp),%rax
movq   $0x0,(%rax)
mov    -0x78(%rbp),%rax
mov    0x8(%rax),%rax
mov    %rax,-0x80(%rbp)
mov    -0x70(%rbp),%rax
mov    %eax,%ecx
mov    $0xffffffff,%eax
shl    %cl,%eax
mov    -0x80(%rbp),%rcx
mov    %eax,%eax
and    %rax,%rcx
mov    -0x78(%rbp),%rax
mov    %rcx,0x8(%rax)
jmp    1234 <modfd+0x134>
mov    -0x78(%rbp),%rax
mov    (%rax),%rax
mov    %rax,-0x88(%rbp)
mov    -0x70(%rbp),%rax
mov    %eax,%ecx
mov    $0xffffffff,%eax
shl    %cl,%eax
mov    -0x88(%rbp),%rcx
mov    %eax,%eax
and    %rax,%rcx
mov    -0x78(%rbp),%rax
mov    %rcx,(%rax)
fldt   -0x60(%rbp)
fldt   -0x40(%rbp)
fsubrp %st,%st(1)
fstpt  -0x50(%rbp)
fldt   -0x20(%rbp)
fldz
fucomip %st(1),%st
fstp   %st(0)
jbe    1263 <modfd+0x163>
jmp    1253 <modfd+0x153>
fldt   -0x40(%rbp)
fchs
fstpt  -0x40(%rbp)
fldt   -0x50(%rbp)
fchs
fstpt  -0x50(%rbp)
fldt   -0x40(%rbp)
mov    -0x28(%rbp),%rax
fstpt  (%rax)
fldt   -0x50(%rbp)
fstpt  -0x10(%rbp)
fldt   -0x10(%rbp)
add    $0x10,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<modfd>:
fldt   0x8(%rsp)
fldz
fucomi %st(1),%st
fld    %st(1)
fchs
fld    %st(2)
fcmovnbe %st(1),%st
fstp   %st(1)
fld    %st(0)
fstpt  -0x18(%rsp)
mov    -0x18(%rsp),%rax
fld    %st(0)
fstpt  -0x28(%rsp)
mov    -0x20(%rsp),%rdx
mov    %edx,%esi
shr    $0x14,%esi
and    $0x7ff,%esi
cmp    $0x3ff,%esi
jae    1140 <modfd+0x40>
fstp   %st(0)
fstpt  (%rdi)
ret
fstp   %st(1)
fldz
cmp    $0x433,%esi
jbe    1153 <modfd+0x53>
fstp   %st(1)
fxch   %st(1)
fstpt  (%rdi)
ret
fstp   %st(0)
cmp    $0x413,%esi
ja     1176 <modfd+0x76>
movq   $0x0,-0x28(%rsp)
mov    $0x13,%cl
sub    %sil,%cl
shr    %cl,%edx
shl    %cl,%edx
mov    %rdx,-0x20(%rsp)
jmp    1184 <modfd+0x84>
mov    $0x33,%cl
sub    %sil,%cl
shr    %cl,%eax
shl    %cl,%eax
mov    %rax,-0x28(%rsp)
fldz
fucomip %st(2),%st
fstp   %st(1)
fldt   -0x28(%rsp)
fsubr  %st,%st(1)
fld    %st(0)
fchs
fxch   %st(1)
fcmovnbe %st(1),%st
fstp   %st(1)
fld    %st(1)
fchs
fxch   %st(2)
fcmovnbe %st(2),%st
fstp   %st(2)
fstpt  (%rdi)
ret

LLVM IR (-O0):
define x86_fp80 @modfd(x86_fp80 %0, x86_fp80* %1) {
%3 = alloca x86_fp80
%4 = alloca x86_fp80
%5 = alloca x86_fp80*
%6 = alloca x86_fp80
%7 = alloca x86_fp80
%8 = alloca x86_fp80
%9 = alloca i64
%10 = alloca i64
%11 = alloca i64*
store x86_fp80 %0, x86_fp80* %4
store x86_fp80* %1, x86_fp80** %5
%12 = bitcast x86_fp80* %6 to i64*
store i64* %12, i64** %11
%13 = load x86_fp80, x86_fp80* %4
store x86_fp80 %13, x86_fp80* %6
%14 = load x86_fp80, x86_fp80* %4
%15 = fcmp olt x86_fp80 %14, 0xK00000000000000000000
br i1 %15, label %16, label %19
16: ; preds = %2
%17 = load x86_fp80, x86_fp80* %6
%18 = fneg x86_fp80 %17
store x86_fp80 %18, x86_fp80* %6
br label %19
19: ; preds = %16, %2
%20 = load x86_fp80, x86_fp80* %6
store x86_fp80 %20, x86_fp80* %8
%21 = load i64*, i64** %11
%22 = getelementptr inbounds i64, i64* %21, i64 1
%23 = load i64, i64* %22
%24 = and i64 %23, 2146435072
%25 = ashr i64 %24, 20
store i64 %25, i64* %9
%26 = load i64, i64* %9
%27 = sub nsw i64 %26, 1023
store i64 %27, i64* %9
%28 = load i64, i64* %9
%29 = icmp slt i64 %28, 0
br i1 %29, label %30, label %34
30: ; preds = %19
%31 = load x86_fp80*, x86_fp80** %5
store x86_fp80 0xK00000000000000000000, x86_fp80* %31
%32 = load x86_fp80, x86_fp80* %4
store x86_fp80 %32, x86_fp80* %7
%33 = load x86_fp80, x86_fp80* %7
store x86_fp80 %33, x86_fp80* %3
br label %86
34: ; preds = %19
%35 = load i64, i64* %9
%36 = icmp sgt i64 %35, 52
br i1 %36, label %37, label %40
37: ; preds = %34
%38 = load x86_fp80, x86_fp80* %4
%39 = load x86_fp80*, x86_fp80** %5
store x86_fp80 %38, x86_fp80* %39
store x86_fp80 0xK00000000000000000000, x86_fp80* %3
br label %86
40: ; preds = %34
%41 = load i64, i64* %9
%42 = sub nsw i64 52, %41
store i64 %42, i64* %10
%43 = load i64, i64* %10
%44 = icmp sgt i64 %43, 31
br i1 %44, label %45, label %60
45: ; preds = %40
%46 = load i64, i64* %10
%47 = sub nsw i64 %46, 32
store i64 %47, i64* %10
%48 = load i64*, i64** %11
%49 = getelementptr inbounds i64, i64* %48, i64 0
store i64 0, i64* %49
%50 = load i64*, i64** %11
%51 = getelementptr inbounds i64, i64* %50, i64 1
%52 = load i64, i64* %51
%53 = load i64, i64* %10
%54 = trunc i64 %53 to i32
%55 = shl i32 -1, %54
%56 = zext i32 %55 to i64
%57 = and i64 %52, %56
%58 = load i64*, i64** %11
%59 = getelementptr inbounds i64, i64* %58, i64 1
store i64 %57, i64* %59
br label %71
60: ; preds = %40
%61 = load i64*, i64** %11
%62 = getelementptr inbounds i64, i64* %61, i64 0
%63 = load i64, i64* %62
%64 = load i64, i64* %10
%65 = trunc i64 %64 to i32
%66 = shl i32 -1, %65
%67 = zext i32 %66 to i64
%68 = and i64 %63, %67
%69 = load i64*, i64** %11
%70 = getelementptr inbounds i64, i64* %69, i64 0
store i64 %68, i64* %70
br label %71
71: ; preds = %60, %45
%72 = load x86_fp80, x86_fp80* %8
%73 = load x86_fp80, x86_fp80* %6
%74 = fsub x86_fp80 %72, %73
store x86_fp80 %74, x86_fp80* %7
%75 = load x86_fp80, x86_fp80* %4
%76 = fcmp olt x86_fp80 %75, 0xK00000000000000000000
br i1 %76, label %77, label %82
77: ; preds = %71
%78 = load x86_fp80, x86_fp80* %6
%79 = fneg x86_fp80 %78
store x86_fp80 %79, x86_fp80* %6
%80 = load x86_fp80, x86_fp80* %7
%81 = fneg x86_fp80 %80
store x86_fp80 %81, x86_fp80* %7
br label %82
82: ; preds = %77, %71
%83 = load x86_fp80, x86_fp80* %6
%84 = load x86_fp80*, x86_fp80** %5
store x86_fp80 %83, x86_fp80* %84
%85 = load x86_fp80, x86_fp80* %7
store x86_fp80 %85, x86_fp80* %3
br label %86
86: ; preds = %82, %37, %30
%87 = load x86_fp80, x86_fp80* %3
ret x86_fp80 %87
}

LLVM IR (-O2):
define x86_fp80 @modfd(x86_fp80 %0, x86_fp80* nocapture writeonly %1) {
%3 = alloca x86_fp80
%4 = bitcast x86_fp80* %3 to i8*
%5 = fcmp olt x86_fp80 %0, 0xK00000000000000000000
%6 = fneg x86_fp80 %0
%7 = select i1 %5, x86_fp80 %6, x86_fp80 %0
store x86_fp80 %7, x86_fp80* %3
%8 = getelementptr inbounds i8, i8* %4, i64 8
%9 = bitcast i8* %8 to i64*
%10 = load i64, i64* %9
%11 = lshr i64 %10, 20
%12 = and i64 %11, 2047
%13 = icmp ult i64 %12, 1023
%14 = bitcast x86_fp80 %7 to i80
%15 = trunc i80 %14 to i64
br i1 %13, label %40, label %16
16: ; preds = %2
%17 = icmp ugt i64 %12, 1075
br i1 %17, label %40, label %18
18: ; preds = %16
%19 = icmp ult i64 %12, 1044
%20 = bitcast x86_fp80* %3 to i64*
br i1 %19, label %21, label %27
21: ; preds = %18
store i64 0, i64* %20
%22 = trunc i64 %12 to i32
%23 = sub nuw nsw i32 1043, %22
%24 = shl nsw i32 -1, %23
%25 = zext i32 %24 to i64
%26 = and i64 %10, %25
store i64 %26, i64* %9
br label %33
27: ; preds = %18
%28 = trunc i64 %12 to i32
%29 = sub nuw nsw i32 1075, %28
%30 = shl nsw i32 -1, %29
%31 = zext i32 %30 to i64
%32 = and i64 %15, %31
store i64 %32, i64* %20
br label %33
33: ; preds = %27, %21
%34 = load x86_fp80, x86_fp80* %3
%35 = fsub x86_fp80 %7, %34
%36 = fneg x86_fp80 %34
%37 = fneg x86_fp80 %35
%38 = select i1 %5, x86_fp80 %36, x86_fp80 %34
%39 = select i1 %5, x86_fp80 %37, x86_fp80 %35
br label %40
40: ; preds = %16, %2, %33
%41 = phi x86_fp80 [ %38, %33 ], [ 0xK00000000000000000000, %2 ], [ %0, %16 ]
%42 = phi x86_fp80 [ %39, %33 ], [ %0, %2 ], [ 0xK00000000000000000000, %16 ]
store x86_fp80 %41, x86_fp80* %1
ret x86_fp80 %42
}