double frexp(d, exp)
double d;
int *exp;
{
  int e = 0;
  if (d > 0) {
    while (d < 0.5)
      d *= 2.0, e--;
    while (d >= 1.0)
      d /= 2.0, e++;
  }
  *exp = e;
  return d;
}


x86 assembly (clang, -O0):
<frexp>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
mov    %rdi,-0x10(%rbp)
movl   $0x0,-0x14(%rbp)
movsd  -0x8(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    11a3 <frexp+0xa3>
jmp    112b <frexp+0x2b>
movsd  0xecd(%rip),%xmm0
ucomisd -0x8(%rbp),%xmm0
jbe    115e <frexp+0x5e>
movsd  0xeca(%rip),%xmm0
mulsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
jmp    112b <frexp+0x2b>
jmp    1163 <frexp+0x63>
movsd  -0x8(%rbp),%xmm0
movsd  0xe98(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     119e <frexp+0x9e>
movsd  -0x8(%rbp),%xmm0
movsd  0xe89(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1163 <frexp+0x63>
jmp    11a3 <frexp+0xa3>
mov    -0x14(%rbp),%ecx
mov    -0x10(%rbp),%rax
mov    %ecx,(%rax)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<frexp>:
xor    %eax,%eax
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    114d <frexp+0x4d>
xor    %eax,%eax
movsd  0xeea(%rip),%xmm1
ucomisd %xmm0,%xmm1
jbe    112d <frexp+0x2d>
xor    %eax,%eax
xchg   %ax,%ax
addsd  %xmm0,%xmm0
add    $0xffffffff,%eax
ucomisd %xmm0,%xmm1
ja     1120 <frexp+0x20>
ucomisd 0xed3(%rip),%xmm0
jb     114d <frexp+0x4d>
movsd  0xec9(%rip),%xmm2
nop
mulsd  %xmm1,%xmm0
add    $0x1,%eax
ucomisd %xmm2,%xmm0
jae    1140 <frexp+0x40>
mov    %eax,(%rdi)
ret

LLVM IR (-O0):
define double @frexp(double %0, i32* %1) {
%3 = alloca double
%4 = alloca i32*
%5 = alloca i32
store double %0, double* %3
store i32* %1, i32** %4
store i32 0, i32* %5
%6 = load double, double* %3
%7 = fcmp ogt double %6, 0.000000e+00
br i1 %7, label %8, label %27
8: ; preds = %2
br label %9
9: ; preds = %12, %8
%10 = load double, double* %3
%11 = fcmp olt double %10, 5.000000e-01
br i1 %11, label %12, label %17
12: ; preds = %9
%13 = load double, double* %3
%14 = fmul double %13, 2.000000e+00
store double %14, double* %3
%15 = load i32, i32* %5
%16 = add nsw i32 %15, -1
store i32 %16, i32* %5
br label %9
17: ; preds = %9
br label %18
18: ; preds = %21, %17
%19 = load double, double* %3
%20 = fcmp oge double %19, 1.000000e+00
br i1 %20, label %21, label %26
21: ; preds = %18
%22 = load double, double* %3
%23 = fdiv double %22, 2.000000e+00
store double %23, double* %3
%24 = load i32, i32* %5
%25 = add nsw i32 %24, 1
store i32 %25, i32* %5
br label %18
26: ; preds = %18
br label %27
27: ; preds = %26, %2
%28 = load i32, i32* %5
%29 = load i32*, i32** %4
store i32 %28, i32* %29
%30 = load double, double* %3
ret double %30
}

LLVM IR (-O2):
define double @frexp(double %0, i32* nocapture writeonly %1) {
%3 = fcmp ogt double %0, 0.000000e+00
br i1 %3, label %4, label %22
4: ; preds = %2
%5 = fcmp olt double %0, 5.000000e-01
br i1 %5, label %10, label %6
6: ; preds = %10, %4
%7 = phi double [ %0, %4 ], [ %13, %10 ]
%8 = phi i32 [ 0, %4 ], [ %14, %10 ]
%9 = fcmp ult double %7, 1.000000e+00
br i1 %9, label %22, label %16
10: ; preds = %4, %10
%11 = phi i32 [ %14, %10 ], [ 0, %4 ]
%12 = phi double [ %13, %10 ], [ %0, %4 ]
%13 = fmul double %12, 2.000000e+00
%14 = add nsw i32 %11, -1
%15 = fcmp olt double %13, 5.000000e-01
br i1 %15, label %10, label %6
16: ; preds = %6, %16
%17 = phi i32 [ %20, %16 ], [ %8, %6 ]
%18 = phi double [ %19, %16 ], [ %7, %6 ]
%19 = fmul double %18, 5.000000e-01
%20 = add nsw i32 %17, 1
%21 = fcmp ult double %19, 1.000000e+00
br i1 %21, label %22, label %16
22: ; preds = %16, %6, %2
%23 = phi double [ %0, %2 ], [ %7, %6 ], [ %19, %16 ]
%24 = phi i32 [ 0, %2 ], [ %8, %6 ], [ %20, %16 ]
store i32 %24, i32* %1
ret double %23
}