void mp_mul_i2d(int n, int radix, int nfft, int shift, int in[],
                double dout[]) {
  int j, x, carry, ndata, radix_2, topdgt;
  ndata = 0;
  topdgt = 0;
  if (n > shift) {
    topdgt = in[shift + 2];
    ndata = (nfft >> 1) + 1;
    if (ndata > n - shift) {
      ndata = n - shift;
    }
  }
  dout[nfft + 1] = in[0] * topdgt;
  for (j = nfft; j > ndata; j--) {
    dout[j] = 0;
  }
  if (ndata > 1) {
    radix_2 = radix / 2;
    carry = 0;
    for (j = ndata + 1; j > 3; j--) {
      x = in[j + shift] - carry;
      carry = x >= radix_2 ? -1 : 0;
      dout[j - 1] = x - (radix & carry);
    }
    dout[2] = in[shift + 3] - carry;
  }
  dout[1] = topdgt;
  dout[0] = in[1] - shift;
}


x86 assembly (clang, -O0):
<mp_mul_i2d>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %ecx,-0x10(%rbp)
mov    %r8,-0x18(%rbp)
mov    %r9,-0x20(%rbp)
movl   $0x0,-0x30(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x4(%rbp),%eax
cmp    -0x10(%rbp),%eax
jle    1170 <mp_mul_i2d+0x70>
mov    -0x18(%rbp),%rax
mov    -0x10(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x38(%rbp)
mov    -0xc(%rbp),%eax
sar    $0x1,%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    -0x10(%rbp),%ecx
cmp    %ecx,%eax
jle    116b <mp_mul_i2d+0x6b>
mov    -0x4(%rbp),%eax
sub    -0x10(%rbp),%eax
mov    %eax,-0x30(%rbp)
jmp    1170 <mp_mul_i2d+0x70>
mov    -0x18(%rbp),%rax
mov    (%rax),%eax
imul   -0x38(%rbp),%eax
cvtsi2sd %eax,%xmm0
mov    -0x20(%rbp),%rax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xc(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x30(%rbp),%eax
jle    11c0 <mp_mul_i2d+0xc0>
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    1196 <mp_mul_i2d+0x96>
cmpl   $0x1,-0x30(%rbp)
jle    126c <mp_mul_i2d+0x16c>
mov    -0x8(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x34(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x3,-0x24(%rbp)
jle    124c <mp_mul_i2d+0x14c>
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x10(%rbp),%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
sub    -0x2c(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%edx
mov    -0x34(%rbp),%esi
xor    %eax,%eax
mov    $0xffffffff,%ecx
cmp    %esi,%edx
cmovge %ecx,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x8(%rbp),%ecx
and    -0x2c(%rbp),%ecx
sub    %ecx,%eax
cvtsi2sd %eax,%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    11e8 <mp_mul_i2d+0xe8>
mov    -0x18(%rbp),%rax
mov    -0x10(%rbp),%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
sub    -0x2c(%rbp),%eax
cvtsi2sd %eax,%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,0x10(%rax)
cvtsi2sdl -0x38(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,0x8(%rax)
mov    -0x18(%rbp),%rax
mov    0x4(%rax),%eax
sub    -0x10(%rbp),%eax
cvtsi2sd %eax,%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<mp_mul_i2d>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    %r9,%r15
mov    %esi,%r12d
xor    %r13d,%r13d
movslq %ecx,%r14
mov    $0x0,%r10d
sub    %ecx,%edi
jle    1144 <mp_mul_i2d+0x34>
mov    0x8(%r8,%r14,4),%r10d
mov    %edx,%eax
sar    %eax
lea    0x1(%rax),%r13d
cmp    %edi,%eax
cmovge %edi,%r13d
mov    (%r8),%eax
imul   %r10d,%eax
cvtsi2sd %eax,%xmm0
movslq %edx,%rax
movsd  %xmm0,0x8(%r15,%rax,8)
cmp    %edx,%r13d
jge    118f <mp_mul_i2d+0x7f>
mov    %r13d,%esi
not    %esi
add    %edx,%esi
sub    %rsi,%rax
lea    (%r15,%rax,8),%rdi
lea    0x8(,%rsi,8),%rdx
xor    %esi,%esi
mov    %r10d,%ebp
mov    %ecx,0x4(%rsp)
mov    %r8,%rbx
call   1030 <memset@plt>
mov    %rbx,%r8
mov    0x4(%rsp),%ecx
mov    %ebp,%r10d
cmp    $0x2,%r13d
jl     1202 <mp_mul_i2d+0xf2>
xor    %ebx,%ebx
cmp    $0x2,%r13d
je     11ee <mp_mul_i2d+0xde>
mov    %r12d,%eax
shr    $0x1f,%eax
add    %r12d,%eax
sar    %eax
add    $0x1,%r13d
movslq %r13d,%rdx
add    $0x1,%rdx
lea    (%r8,%r14,4),%rsi
xor    %ebx,%ebx
nopl   0x0(%rax)
mov    -0x4(%rsi,%rdx,4),%edi
sub    %ebx,%edi
xor    %ebx,%ebx
cmp    %eax,%edi
setge  %bl
neg    %ebx
mov    %ebx,%ebp
and    %r12d,%ebp
sub    %ebp,%edi
xorps  %xmm0,%xmm0
cvtsi2sd %edi,%xmm0
movsd  %xmm0,-0x10(%r15,%rdx,8)
add    $0xffffffffffffffff,%rdx
cmp    $0x4,%rdx
jg     11c0 <mp_mul_i2d+0xb0>
mov    0xc(%r8,%r14,4),%eax
sub    %ebx,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,0x10(%r15)
xorps  %xmm0,%xmm0
cvtsi2sd %r10d,%xmm0
movsd  %xmm0,0x8(%r15)
mov    0x4(%r8),%eax
sub    %ecx,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%r15)
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @mp_mul_i2d(i32 %0, i32 %1, i32 %2, i32 %3, i32* %4, double* %5) {
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32*
%12 = alloca double*
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
store i32 %0, i32* %7
store i32 %1, i32* %8
store i32 %2, i32* %9
store i32 %3, i32* %10
store i32* %4, i32** %11
store double* %5, double** %12
store i32 0, i32* %16
store i32 0, i32* %18
%19 = load i32, i32* %7
%20 = load i32, i32* %10
%21 = icmp sgt i32 %19, %20
br i1 %21, label %22, label %42
22: ; preds = %6
%23 = load i32*, i32** %11
%24 = load i32, i32* %10
%25 = add nsw i32 %24, 2
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds i32, i32* %23, i64 %26
%28 = load i32, i32* %27
store i32 %28, i32* %18
%29 = load i32, i32* %9
%30 = ashr i32 %29, 1
%31 = add nsw i32 %30, 1
store i32 %31, i32* %16
%32 = load i32, i32* %16
%33 = load i32, i32* %7
%34 = load i32, i32* %10
%35 = sub nsw i32 %33, %34
%36 = icmp sgt i32 %32, %35
br i1 %36, label %37, label %41
37: ; preds = %22
%38 = load i32, i32* %7
%39 = load i32, i32* %10
%40 = sub nsw i32 %38, %39
store i32 %40, i32* %16
br label %41
41: ; preds = %37, %22
br label %42
42: ; preds = %41, %6
%43 = load i32*, i32** %11
%44 = getelementptr inbounds i32, i32* %43, i64 0
%45 = load i32, i32* %44
%46 = load i32, i32* %18
%47 = mul nsw i32 %45, %46
%48 = sitofp i32 %47 to double
%49 = load double*, double** %12
%50 = load i32, i32* %9
%51 = add nsw i32 %50, 1
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds double, double* %49, i64 %52
store double %48, double* %53
%54 = load i32, i32* %9
store i32 %54, i32* %13
br label %55
55: ; preds = %64, %42
%56 = load i32, i32* %13
%57 = load i32, i32* %16
%58 = icmp sgt i32 %56, %57
br i1 %58, label %59, label %67
59: ; preds = %55
%60 = load double*, double** %12
%61 = load i32, i32* %13
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds double, double* %60, i64 %62
store double 0.000000e+00, double* %63
br label %64
64: ; preds = %59
%65 = load i32, i32* %13
%66 = add nsw i32 %65, -1
store i32 %66, i32* %13
br label %55
67: ; preds = %55
%68 = load i32, i32* %16
%69 = icmp sgt i32 %68, 1
br i1 %69, label %70, label %119
70: ; preds = %67
%71 = load i32, i32* %8
%72 = sdiv i32 %71, 2
store i32 %72, i32* %17
store i32 0, i32* %15
%73 = load i32, i32* %16
%74 = add nsw i32 %73, 1
store i32 %74, i32* %13
br label %75
75: ; preds = %104, %70
%76 = load i32, i32* %13
%77 = icmp sgt i32 %76, 3
br i1 %77, label %78, label %107
78: ; preds = %75
%79 = load i32*, i32** %11
%80 = load i32, i32* %13
%81 = load i32, i32* %10
%82 = add nsw i32 %80, %81
%83 = sext i32 %82 to i64
%84 = getelementptr inbounds i32, i32* %79, i64 %83
%85 = load i32, i32* %84
%86 = load i32, i32* %15
%87 = sub nsw i32 %85, %86
store i32 %87, i32* %14
%88 = load i32, i32* %14
%89 = load i32, i32* %17
%90 = icmp sge i32 %88, %89
%91 = zext i1 %90 to i64
%92 = select i1 %90, i32 -1, i32 0
store i32 %92, i32* %15
%93 = load i32, i32* %14
%94 = load i32, i32* %8
%95 = load i32, i32* %15
%96 = and i32 %94, %95
%97 = sub nsw i32 %93, %96
%98 = sitofp i32 %97 to double
%99 = load double*, double** %12
%100 = load i32, i32* %13
%101 = sub nsw i32 %100, 1
%102 = sext i32 %101 to i64
%103 = getelementptr inbounds double, double* %99, i64 %102
store double %98, double* %103
br label %104
104: ; preds = %78
%105 = load i32, i32* %13
%106 = add nsw i32 %105, -1
store i32 %106, i32* %13
br label %75
107: ; preds = %75
%108 = load i32*, i32** %11
%109 = load i32, i32* %10
%110 = add nsw i32 %109, 3
%111 = sext i32 %110 to i64
%112 = getelementptr inbounds i32, i32* %108, i64 %111
%113 = load i32, i32* %112
%114 = load i32, i32* %15
%115 = sub nsw i32 %113, %114
%116 = sitofp i32 %115 to double
%117 = load double*, double** %12
%118 = getelementptr inbounds double, double* %117, i64 2
store double %116, double* %118
br label %119
119: ; preds = %107, %67
%120 = load i32, i32* %18
%121 = sitofp i32 %120 to double
%122 = load double*, double** %12
%123 = getelementptr inbounds double, double* %122, i64 1
store double %121, double* %123
%124 = load i32*, i32** %11
%125 = getelementptr inbounds i32, i32* %124, i64 1
%126 = load i32, i32* %125
%127 = load i32, i32* %10
%128 = sub nsw i32 %126, %127
%129 = sitofp i32 %128 to double
%130 = load double*, double** %12
%131 = getelementptr inbounds double, double* %130, i64 0
store double %129, double* %131
ret void
}

LLVM IR (-O2):
define void @mp_mul_i2d(i32 %0, i32 %1, i32 %2, i32 %3, i32* nocapture %4, double* nocapture writeonly %5) {
%7 = icmp sgt i32 %0, %3
br i1 %7, label %8, label %18
8: ; preds = %6
%9 = add nsw i32 %3, 2
%10 = sext i32 %9 to i64
%11 = getelementptr inbounds i32, i32* %4, i64 %10
%12 = load i32, i32* %11
%13 = ashr i32 %2, 1
%14 = add nsw i32 %13, 1
%15 = sub nsw i32 %0, %3
%16 = icmp slt i32 %13, %15
%17 = select i1 %16, i32 %14, i32 %15
br label %18
18: ; preds = %8, %6
%19 = phi i32 [ 0, %6 ], [ %17, %8 ]
%20 = phi i32 [ 0, %6 ], [ %12, %8 ]
%21 = load i32, i32* %4
%22 = mul nsw i32 %21, %20
%23 = sitofp i32 %22 to double
%24 = add nsw i32 %2, 1
%25 = sext i32 %24 to i64
%26 = getelementptr inbounds double, double* %5, i64 %25
store double %23, double* %26
%27 = icmp slt i32 %19, %2
br i1 %27, label %28, label %38
28: ; preds = %18
%29 = sext i32 %2 to i64
%30 = xor i32 %19, -1
%31 = add i32 %30, %2
%32 = zext i32 %31 to i64
%33 = sub nsw i64 %29, %32
%34 = getelementptr double, double* %5, i64 %33
%35 = bitcast double* %34 to i8*
%36 = shl nuw nsw i64 %32, 3
%37 = add nuw nsw i64 %36, 8
br label %38
38: ; preds = %28, %18
%39 = icmp sgt i32 %19, 1
br i1 %39, label %40, label %71
40: ; preds = %38
%41 = sdiv i32 %1, 2
%42 = icmp eq i32 %19, 2
br i1 %42, label %62, label %43
43: ; preds = %40
%44 = add nuw i32 %19, 1
%45 = sext i32 %44 to i64
%46 = sext i32 %3 to i64
br label %47
47: ; preds = %43, %47
%48 = phi i64 [ %45, %43 ], [ %59, %47 ]
%49 = phi i32 [ 0, %43 ], [ %55, %47 ]
%50 = add nsw i64 %48, %46
%51 = getelementptr inbounds i32, i32* %4, i64 %50
%52 = load i32, i32* %51
%53 = sub nsw i32 %52, %49
%54 = icmp sge i32 %53, %41
%55 = sext i1 %54 to i32
%56 = and i32 %55, %1
%57 = sub nsw i32 %53, %56
%58 = sitofp i32 %57 to double
%59 = add nsw i64 %48, -1
%60 = getelementptr inbounds double, double* %5, i64 %59
store double %58, double* %60
%61 = icmp sgt i64 %48, 4
br i1 %61, label %47, label %62
62: ; preds = %47, %40
%63 = phi i32 [ 0, %40 ], [ %55, %47 ]
%64 = add nsw i32 %3, 3
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds i32, i32* %4, i64 %65
%67 = load i32, i32* %66
%68 = sub nsw i32 %67, %63
%69 = sitofp i32 %68 to double
%70 = getelementptr inbounds double, double* %5, i64 2
store double %69, double* %70
br label %71
71: ; preds = %62, %38
%72 = sitofp i32 %20 to double
%73 = getelementptr inbounds double, double* %5, i64 1
store double %72, double* %73
%74 = getelementptr inbounds i32, i32* %4, i64 1
%75 = load i32, i32* %74
%76 = sub nsw i32 %75, %3
%77 = sitofp i32 %76 to double
store double %77, double* %5
ret void
}