double ldexp(double v, int e) {
  if (e < 0) {
    for (;;) {
      v /= 2;
      if (++e == 0) {
        return (v);
      }
    }
  }
  while (e > 0) {
    v *= 2;
    e -= 1;
  }
  return (v);
}


x86 assembly (clang, -O0):
<ldexp>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
mov    %edi,-0x14(%rbp)
cmpl   $0x0,-0x14(%rbp)
jge    1157 <ldexp+0x57>
jmp    111b <ldexp+0x1b>
movsd  -0x10(%rbp),%xmm0
movsd  0xed8(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x10(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
cmp    $0x0,%eax
jne    1152 <ldexp+0x52>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1190 <ldexp+0x90>
jmp    111b <ldexp+0x1b>
jmp    115c <ldexp+0x5c>
cmpl   $0x0,-0x14(%rbp)
jle    1186 <ldexp+0x86>
movsd  0xe92(%rip),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x10(%rbp)
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    115c <ldexp+0x5c>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<ldexp>:
test   %edi,%edi
js     116a <ldexp+0x6a>
je     11d5 <ldexp+0xd5>
lea    -0x1(%rdi),%eax
mov    %edi,%ecx
and    $0x7,%ecx
je     112d <ldexp+0x2d>
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
addsd  %xmm0,%xmm0
add    $0x1,%edx
cmp    %edx,%ecx
jne    1120 <ldexp+0x20>
sub    %edx,%edi
cmp    $0x7,%eax
jb     11d5 <ldexp+0xd5>
add    $0xffffffff,%edi
nopl   0x0(%rax)
addsd  %xmm0,%xmm0
addsd  %xmm0,%xmm0
addsd  %xmm0,%xmm0
addsd  %xmm0,%xmm0
addsd  %xmm0,%xmm0
addsd  %xmm0,%xmm0
addsd  %xmm0,%xmm0
addsd  %xmm0,%xmm0
add    $0xfffffff8,%edi
cmp    $0xfffffffe,%edi
jb     1140 <ldexp+0x40>
jmp    11d5 <ldexp+0xd5>
mov    %edi,%eax
neg    %eax
and    $0x7,%eax
je     1196 <ldexp+0x96>
neg    %eax
xor    %ecx,%ecx
movsd  0xe81(%rip),%xmm1
nop
mulsd  %xmm1,%xmm0
add    $0xffffffff,%ecx
cmp    %ecx,%eax
jne    1180 <ldexp+0x80>
mov    %edi,%eax
sub    %ecx,%eax
cmp    $0xfffffff8,%edi
jbe    119d <ldexp+0x9d>
jmp    11d5 <ldexp+0xd5>
mov    %edi,%eax
cmp    $0xfffffff8,%edi
ja     11d5 <ldexp+0xd5>
movsd  0xe5b(%rip),%xmm1
cs nopw 0x0(%rax,%rax,1)
nop
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
add    $0x8,%eax
jne    11b0 <ldexp+0xb0>
ret

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

LLVM IR (-O2):
define double @ldexp(double %0, i32 %1) {
%3 = icmp slt i32 %1, 0
br i1 %3, label %4, label %21
4: ; preds = %2
%5 = sub i32 0, %1
%6 = and i32 %5, 7
%7 = icmp eq i32 %6, 0
br i1 %7, label %16, label %8
8: ; preds = %4, %8
%9 = phi double [ %12, %8 ], [ %0, %4 ]
%10 = phi i32 [ %13, %8 ], [ %1, %4 ]
%11 = phi i32 [ %14, %8 ], [ 0, %4 ]
%12 = fmul double %9, 5.000000e-01
%13 = add nsw i32 %10, 1
%14 = add i32 %11, 1
%15 = icmp eq i32 %14, %6
br i1 %15, label %16, label %8
16: ; preds = %8, %4
%17 = phi double [ undef, %4 ], [ %12, %8 ]
%18 = phi double [ %0, %4 ], [ %12, %8 ]
%19 = phi i32 [ %1, %4 ], [ %13, %8 ]
%20 = icmp ugt i32 %1, -8
br i1 %20, label %67, label %40
21: ; preds = %2
%22 = icmp eq i32 %1, 0
br i1 %22, label %67, label %23
23: ; preds = %21
%24 = add i32 %1, -1
%25 = and i32 %1, 7
%26 = icmp eq i32 %25, 0
br i1 %26, label %35, label %27
27: ; preds = %23, %27
%28 = phi i32 [ %32, %27 ], [ %1, %23 ]
%29 = phi double [ %31, %27 ], [ %0, %23 ]
%30 = phi i32 [ %33, %27 ], [ 0, %23 ]
%31 = fmul double %29, 2.000000e+00
%32 = add nsw i32 %28, -1
%33 = add i32 %30, 1
%34 = icmp eq i32 %33, %25
br i1 %34, label %35, label %27
35: ; preds = %27, %23
%36 = phi i32 [ %1, %23 ], [ %32, %27 ]
%37 = phi double [ %0, %23 ], [ %31, %27 ]
%38 = phi double [ undef, %23 ], [ %31, %27 ]
%39 = icmp ult i32 %24, 7
br i1 %39, label %67, label %53
40: ; preds = %16, %40
%41 = phi double [ %50, %40 ], [ %18, %16 ]
%42 = phi i32 [ %51, %40 ], [ %19, %16 ]
%43 = fmul double %41, 5.000000e-01
%44 = fmul double %43, 5.000000e-01
%45 = fmul double %44, 5.000000e-01
%46 = fmul double %45, 5.000000e-01
%47 = fmul double %46, 5.000000e-01
%48 = fmul double %47, 5.000000e-01
%49 = fmul double %48, 5.000000e-01
%50 = fmul double %49, 5.000000e-01
%51 = add nsw i32 %42, 8
%52 = icmp eq i32 %51, 0
br i1 %52, label %67, label %40
53: ; preds = %35, %53
%54 = phi i32 [ %64, %53 ], [ %36, %35 ]
%55 = phi double [ %63, %53 ], [ %37, %35 ]
%56 = fmul double %55, 2.000000e+00
%57 = fmul double %56, 2.000000e+00
%58 = fmul double %57, 2.000000e+00
%59 = fmul double %58, 2.000000e+00
%60 = fmul double %59, 2.000000e+00
%61 = fmul double %60, 2.000000e+00
%62 = fmul double %61, 2.000000e+00
%63 = fmul double %62, 2.000000e+00
%64 = add nsw i32 %54, -8
%65 = add i32 %54, -9
%66 = icmp ult i32 %65, -2
br i1 %66, label %53, label %67
67: ; preds = %35, %53, %16, %40, %21
%68 = phi double [ %0, %21 ], [ %17, %16 ], [ %50, %40 ], [ %38, %35 ], [ %63, %53 ]
ret double %68
}