double mp_unexp_mp2d(int n, int radix, short int in[]) {
  int j;
  double d1_radix, dout;
  d1_radix = (double)1 / radix;
  dout = 0;
  for (j = n - 1; j >= 0; j--) {
    dout = d1_radix * dout + in[j];
  }
  return dout;
}


x86 assembly (clang, -O0):
<mp_unexp_mp2d>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
cvtsi2sdl -0x8(%rbp),%xmm1
movsd  0xee5(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x14(%rbp)
cmpl   $0x0,-0x14(%rbp)
jl     1174 <mp_unexp_mp2d+0x74>
movsd  -0x20(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm2
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movswl (%rax,%rcx,2),%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
jmp    1135 <mp_unexp_mp2d+0x35>
movsd  -0x28(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<mp_unexp_mp2d>:
test   %edi,%edi
jle    11c8 <mp_unexp_mp2d+0xc8>
cvtsi2sd %esi,%xmm0
movsd  0xeec(%rip),%xmm1
divsd  %xmm0,%xmm1
mov    %edi,%eax
lea    -0x1(%rax),%rcx
mov    %rax,%rsi
xorpd  %xmm0,%xmm0
and    $0x3,%rsi
je     114f <mp_unexp_mp2d+0x4f>
nopl   0x0(%rax,%rax,1)
add    $0xffffffffffffffff,%rax
mov    %eax,%edi
movswl (%rdx,%rdi,2),%edi
xorps  %xmm2,%xmm2
cvtsi2sd %edi,%xmm2
mulsd  %xmm1,%xmm0
addsd  %xmm2,%xmm0
add    $0xffffffffffffffff,%rsi
jne    1130 <mp_unexp_mp2d+0x30>
cmp    $0x3,%rcx
jb     11c7 <mp_unexp_mp2d+0xc7>
mov    $0xfffffffe,%ecx
nopw   0x0(%rax,%rax,1)
lea    (%rax,%rcx,1),%esi
add    $0x1,%esi
movswl (%rdx,%rsi,2),%esi
xorps  %xmm2,%xmm2
cvtsi2sd %esi,%xmm2
mulsd  %xmm1,%xmm0
addsd  %xmm2,%xmm0
lea    (%rax,%rcx,1),%esi
movswl (%rdx,%rsi,2),%esi
xorps  %xmm2,%xmm2
cvtsi2sd %esi,%xmm2
mulsd  %xmm1,%xmm0
addsd  %xmm2,%xmm0
lea    -0x3(%rax),%rsi
mov    %esi,%edi
movswl (%rdx,%rdi,2),%edi
xorps  %xmm2,%xmm2
cvtsi2sd %edi,%xmm2
mulsd  %xmm1,%xmm0
add    $0xfffffffffffffffc,%rax
mov    %eax,%edi
movswl (%rdx,%rdi,2),%edi
xorps  %xmm3,%xmm3
cvtsi2sd %edi,%xmm3
addsd  %xmm2,%xmm0
mulsd  %xmm1,%xmm0
addsd  %xmm3,%xmm0
cmp    $0x1,%rsi
ja     1160 <mp_unexp_mp2d+0x60>
ret
xorps  %xmm0,%xmm0
ret

LLVM IR (-O0):
define double @mp_unexp_mp2d(i32 %0, i32 %1, i16* %2) {
%4 = alloca i32
%5 = alloca i32
%6 = alloca i16*
%7 = alloca i32
%8 = alloca double
%9 = alloca double
store i32 %0, i32* %4
store i32 %1, i32* %5
store i16* %2, i16** %6
%10 = load i32, i32* %5
%11 = sitofp i32 %10 to double
%12 = fdiv double 1.000000e+00, %11
store double %12, double* %8
store double 0.000000e+00, double* %9
%13 = load i32, i32* %4
%14 = sub nsw i32 %13, 1
store i32 %14, i32* %7
br label %15
15: ; preds = %29, %3
%16 = load i32, i32* %7
%17 = icmp sge i32 %16, 0
br i1 %17, label %18, label %32
18: ; preds = %15
%19 = load double, double* %8
%20 = load double, double* %9
%21 = load i16*, i16** %6
%22 = load i32, i32* %7
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds i16, i16* %21, i64 %23
%25 = load i16, i16* %24
%26 = sext i16 %25 to i32
%27 = sitofp i32 %26 to double
%28 = call double @llvm.fmuladd.f64(double %19, double %20, double %27)
store double %28, double* %9
br label %29
29: ; preds = %18
%30 = load i32, i32* %7
%31 = add nsw i32 %30, -1
store i32 %31, i32* %7
br label %15
32: ; preds = %15
%33 = load double, double* %9
ret double %33
}

LLVM IR (-O2):
define double @mp_unexp_mp2d(i32 %0, i32 %1, i16* nocapture %2) {
%4 = sitofp i32 %1 to double
%5 = fdiv double 1.000000e+00, %4
%6 = icmp sgt i32 %0, 0
br i1 %6, label %7, label %57
7: ; preds = %3
%8 = zext i32 %0 to i64
%9 = add nsw i64 %8, -1
%10 = and i64 %8, 3
%11 = icmp eq i64 %10, 0
br i1 %11, label %24, label %12
12: ; preds = %7, %12
%13 = phi i64 [ %16, %12 ], [ %8, %7 ]
%14 = phi double [ %21, %12 ], [ 0.000000e+00, %7 ]
%15 = phi i64 [ %22, %12 ], [ 0, %7 ]
%16 = add nsw i64 %13, -1
%17 = and i64 %16, 4294967295
%18 = getelementptr inbounds i16, i16* %2, i64 %17
%19 = load i16, i16* %18
%20 = sitofp i16 %19 to double
%21 = tail call double @llvm.fmuladd.f64(double %5, double %14, double %20)
%22 = add i64 %15, 1
%23 = icmp eq i64 %22, %10
br i1 %23, label %24, label %12
24: ; preds = %12, %7
%25 = phi i64 [ %8, %7 ], [ %16, %12 ]
%26 = phi double [ 0.000000e+00, %7 ], [ %21, %12 ]
%27 = phi double [ undef, %7 ], [ %21, %12 ]
%28 = icmp ult i64 %9, 3
br i1 %28, label %57, label %29
29: ; preds = %24, %29
%30 = phi i64 [ %50, %29 ], [ %25, %24 ]
%31 = phi double [ %55, %29 ], [ %26, %24 ]
%32 = add i64 %30, 4294967295
%33 = and i64 %32, 4294967295
%34 = getelementptr inbounds i16, i16* %2, i64 %33
%35 = load i16, i16* %34
%36 = sitofp i16 %35 to double
%37 = tail call double @llvm.fmuladd.f64(double %5, double %31, double %36)
%38 = add i64 %30, 4294967294
%39 = and i64 %38, 4294967295
%40 = getelementptr inbounds i16, i16* %2, i64 %39
%41 = load i16, i16* %40
%42 = sitofp i16 %41 to double
%43 = tail call double @llvm.fmuladd.f64(double %5, double %37, double %42)
%44 = add nsw i64 %30, -3
%45 = and i64 %44, 4294967295
%46 = getelementptr inbounds i16, i16* %2, i64 %45
%47 = load i16, i16* %46
%48 = sitofp i16 %47 to double
%49 = tail call double @llvm.fmuladd.f64(double %5, double %43, double %48)
%50 = add nsw i64 %30, -4
%51 = and i64 %50, 4294967295
%52 = getelementptr inbounds i16, i16* %2, i64 %51
%53 = load i16, i16* %52
%54 = sitofp i16 %53 to double
%55 = tail call double @llvm.fmuladd.f64(double %5, double %49, double %54)
%56 = icmp ugt i64 %44, 1
br i1 %56, label %29, label %57
57: ; preds = %24, %29, %3
%58 = phi double [ 0.000000e+00, %3 ], [ %27, %24 ], [ %55, %29 ]
ret double %58
}