double 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[0];
  i = *q >> 4;
  i &= 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):
<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
mov    %rax,-0x30(%rbp)
mov    -0x30(%rbp),%rax
movswl (%rax),%eax
sar    $0x4,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
and    $0x7ff,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
je     1146 <frexp+0x46>
jmp    11ba <frexp+0xba>
movsd  -0x20(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1175 <frexp+0x75>
jp     1175 <frexp+0x75>
mov    -0x18(%rbp),%rax
movl   $0x0,(%rax)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    11f8 <frexp+0xf8>
jmp    117a <frexp+0x7a>
movsd  0xe7e(%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     117a <frexp+0x7a>
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):
<frexp>:
movq   %xmm0,%rax
mov    %eax,%edx
shr    $0x4,%edx
and    $0x7ff,%edx
je     112c <frexp+0x2c>
add    $0xfffffc02,%edx
mov    %edx,(%rdi)
and    $0xffffffffffff800f,%rax
or     $0x3fe0,%rax
movq   %rax,%xmm0
ret
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1143 <frexp+0x43>
jp     1143 <frexp+0x43>
movl   $0x0,(%rdi)
movapd %xmm1,%xmm0
ret
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nop
addsd  %xmm0,%xmm0
add    $0xffffffff,%ecx
movq   %xmm0,%rax
mov    %eax,%edx
shr    $0x4,%edx
test   $0x7ff,%edx
je     1150 <frexp+0x50>
mov    %eax,%edx
shr    $0x4,%edx
and    $0x7ff,%edx
add    %ecx,%edx
jmp    1112 <frexp+0x12>

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

LLVM IR (-O2):
define double @frexp(double %0, i32* nocapture writeonly %1) {
%3 = bitcast double %0 to i64
%4 = trunc i64 %3 to i32
%5 = lshr i32 %4, 4
%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 = trunc i64 %16 to i32
%18 = lshr i32 %17, 4
%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, -32753
%28 = or i64 %27, 16352
%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
}