double cephes_frexp(x, pw2)
double x;
int *pw2;
{
  union {
    double y;
    unsigned short sh[4];
  } u;
  int i;
  int k;
  short *q;
  u.y = x;
  q = (short *)&u.sh[3];
  i = (*q >> 4) & 0x7ff;
  if (i != 0)
    goto ieeedon;
  if (u.y == 0.0) {
    *pw2 = 0;
    return (0.0);
  }
  do {
    u.y *= 2.0;
    i -= 1;
    k = (*q >> 4) & 0x7ff;
  } while (k == 0);
  i = i + k;
ieeedon:
  i -= 0x3fe;
  *pw2 = i;
  *q &= 0x800f;
  *q |= 0x3fe0;
  return (u.y);
}


x86 assembly (clang, -O0):
<cephes_frexp>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
mov    %rdi,-0x18(%rbp)
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
lea    -0x20(%rbp),%rax
add    $0x6,%rax
mov    %rax,-0x30(%rbp)
mov    -0x30(%rbp),%rax
movswl (%rax),%eax
sar    $0x4,%eax
and    $0x7ff,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
je     1144 <cephes_frexp+0x44>
jmp    11b8 <cephes_frexp+0xb8>
movsd  -0x20(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1173 <cephes_frexp+0x73>
jp     1173 <cephes_frexp+0x73>
mov    -0x18(%rbp),%rax
movl   $0x0,(%rax)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    11f6 <cephes_frexp+0xf6>
jmp    1178 <cephes_frexp+0x78>
movsd  0xe80(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x30(%rbp),%rax
movswl (%rax),%eax
sar    $0x4,%eax
and    $0x7ff,%eax
mov    %eax,-0x28(%rbp)
cmpl   $0x0,-0x28(%rbp)
je     1178 <cephes_frexp+0x78>
mov    -0x24(%rbp),%eax
add    -0x28(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
sub    $0x3fe,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x30(%rbp),%rax
movswl (%rax),%ecx
and    $0x800f,%ecx
mov    %cx,(%rax)
mov    -0x30(%rbp),%rax
movswl (%rax),%ecx
or     $0x3fe0,%ecx
mov    %cx,(%rax)
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<cephes_frexp>:
movq   %xmm0,%rax
mov    %rax,%rcx
shr    $0x34,%rcx
and    $0x7ff,%ecx
je     113c <cephes_frexp+0x3c>
add    $0xfffffc02,%ecx
mov    %ecx,(%rdi)
movabs $0x800fffffffffffff,%rcx
ff 0f 80
and    %rax,%rcx
movabs $0x3fe0000000000000,%rax
or     %rcx,%rax
movq   %rax,%xmm0
ret
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1153 <cephes_frexp+0x53>
jp     1153 <cephes_frexp+0x53>
movl   $0x0,(%rdi)
movapd %xmm1,%xmm0
ret
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nop
addsd  %xmm0,%xmm0
add    $0xffffffff,%edx
movq   %xmm0,%rax
mov    %rax,%rcx
shr    $0x34,%rcx
test   $0x7ff,%ecx
je     1160 <cephes_frexp+0x60>
and    $0x7ff,%ecx
add    %edx,%ecx
jmp    1114 <cephes_frexp+0x14>

LLVM IR (-O0):
%union.anon = type { double }
define double @cephes_frexp(double %0, i32* %1) {
%3 = alloca double
%4 = alloca double
%5 = alloca i32*
%6 = alloca %union.anon
%7 = alloca i32
%8 = alloca i32
%9 = alloca i16*
store double %0, double* %4
store i32* %1, i32** %5
%10 = load double, double* %4
%11 = bitcast %union.anon* %6 to double*
store double %10, double* %11
%12 = bitcast %union.anon* %6 to [4 x i16]*
%13 = getelementptr inbounds [4 x i16], [4 x i16]* %12, i64 0, i64 3
store i16* %13, i16** %9
%14 = load i16*, i16** %9
%15 = load i16, i16* %14
%16 = sext i16 %15 to i32
%17 = ashr i32 %16, 4
%18 = and i32 %17, 2047
store i32 %18, i32* %7
%19 = load i32, i32* %7
%20 = icmp ne i32 %19, 0
br i1 %20, label %21, label %22
21: ; preds = %2
br label %47
22: ; preds = %2
%23 = bitcast %union.anon* %6 to double*
%24 = load double, double* %23
%25 = fcmp oeq double %24, 0.000000e+00
br i1 %25, label %26, label %28
26: ; preds = %22
%27 = load i32*, i32** %5
store i32 0, i32* %27
store double 0.000000e+00, double* %3
br label %64
28: ; preds = %22
br label %29
29: ; preds = %40, %28
%30 = bitcast %union.anon* %6 to double*
%31 = load double, double* %30
%32 = fmul double %31, 2.000000e+00
store double %32, double* %30
%33 = load i32, i32* %7
%34 = sub nsw i32 %33, 1
store i32 %34, i32* %7
%35 = load i16*, i16** %9
%36 = load i16, i16* %35
%37 = sext i16 %36 to i32
%38 = ashr i32 %37, 4
%39 = and i32 %38, 2047
store i32 %39, i32* %8
br label %40
40: ; preds = %29
%41 = load i32, i32* %8
%42 = icmp eq i32 %41, 0
br i1 %42, label %29, label %43
43: ; preds = %40
%44 = load i32, i32* %7
%45 = load i32, i32* %8
%46 = add nsw i32 %44, %45
store i32 %46, i32* %7
br label %47
47: ; preds = %43, %21
%48 = load i32, i32* %7
%49 = sub nsw i32 %48, 1022
store i32 %49, i32* %7
%50 = load i32, i32* %7
%51 = load i32*, i32** %5
store i32 %50, i32* %51
%52 = load i16*, i16** %9
%53 = load i16, i16* %52
%54 = sext i16 %53 to i32
%55 = and i32 %54, 32783
%56 = trunc i32 %55 to i16
store i16 %56, i16* %52
%57 = load i16*, i16** %9
%58 = load i16, i16* %57
%59 = sext i16 %58 to i32
%60 = or i32 %59, 16352
%61 = trunc i32 %60 to i16
store i16 %61, i16* %57
%62 = bitcast %union.anon* %6 to double*
%63 = load double, double* %62
store double %63, double* %3
br label %64
64: ; preds = %47, %26
%65 = load double, double* %3
ret double %65
}

LLVM IR (-O2):
define double @cephes_frexp(double %0, i32* nocapture writeonly %1) {
%3 = bitcast double %0 to i64
%4 = lshr i64 %3, 52
%5 = trunc i64 %4 to i32
%6 = and i32 %5, 2047
%7 = icmp eq i32 %6, 0
br i1 %7, label %8, label %23
8: ; preds = %2
%9 = fcmp oeq double %0, 0.000000e+00
br i1 %9, label %10, label %11
10: ; preds = %8
store i32 0, i32* %1
br label %30
11: ; preds = %8, %11
%12 = phi double [ %14, %11 ], [ %0, %8 ]
%13 = phi i32 [ %15, %11 ], [ 0, %8 ]
%14 = fmul double %12, 2.000000e+00
%15 = add nsw i32 %13, -1
%16 = bitcast double %14 to i64
%17 = lshr i64 %16, 52
%18 = trunc i64 %17 to i32
%19 = and i32 %18, 2047
%20 = icmp eq i32 %19, 0
br i1 %20, label %11, label %21
21: ; preds = %11
%22 = add nsw i32 %19, %15
br label %23
23: ; preds = %2, %21
%24 = phi i64 [ %3, %2 ], [ %16, %21 ]
%25 = phi i32 [ %6, %2 ], [ %22, %21 ]
%26 = add nsw i32 %25, -1022
store i32 %26, i32* %1
%27 = and i64 %24, -9218868437227405313
%28 = or i64 %27, 4602678819172646912
%29 = bitcast i64 %28 to double
br label %30
30: ; preds = %23, %10
%31 = phi double [ %29, %23 ], [ 0.000000e+00, %10 ]
ret double %31
}