void mp_unsgn_imul(int n, double dradix, int in1[], double din2, int out[]) {
  int j, carry, shift;
  double x, d1_radix;
  short int *in1dgt, *outdgt;
  in1dgt = (short int *)&in1[1];
  outdgt = (short int *)&out[1];
  d1_radix = (double)1 / dradix;
  carry = 0;
  for (j = n - 1; j >= 0; j--) {
    x = din2 * in1dgt[j] + carry + ((double)0.5);
    carry = (int)(d1_radix * x);
    outdgt[j] = (short int)(x - dradix * carry);
  }
  shift = 0;
  x = carry + ((double)0.5);
  while (x > 1) {
    x *= d1_radix;
    shift++;
  }
  out[0] = in1[0] + shift;
  if (shift > 0) {
    while (shift > n) {
      carry = (int)(d1_radix * carry + ((double)0.5));
      shift--;
    }
    for (j = n - 1; j >= shift; j--) {
      outdgt[j] = outdgt[j - shift];
    }
    for (j = shift - 1; j >= 0; j--) {
      x = carry + ((double)0.5);
      carry = (int)(d1_radix * x);
      outdgt[j] = (short int)(x - dradix * carry);
    }
  }
}


x86 assembly (clang, -O0):
<mp_unsgn_imul>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
mov    %rdx,-0x28(%rbp)
mov    -0x18(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x50(%rbp)
mov    -0x28(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x58(%rbp)
movsd  0xec7(%rip),%xmm0
divsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x0,-0x2c(%rbp)
jl     11dd <mp_unsgn_imul+0xdd>
movsd  -0x20(%rbp),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movswl (%rax,%rcx,2),%eax
cvtsi2sd %eax,%xmm2
cvtsi2sdl -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  0xe81(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x30(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  -0x10(%rbp),%xmm1
cvtsi2sdl -0x30(%rbp),%xmm2
mulsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %ax,%dx
mov    -0x58(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
mov    -0x2c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x2c(%rbp)
jmp    1153 <mp_unsgn_imul+0x53>
movl   $0x0,-0x34(%rbp)
cvtsi2sdl -0x30(%rbp),%xmm0
movsd  0xe17(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  0xdf9(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    122e <mp_unsgn_imul+0x12e>
movsd  -0x48(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11fa <mp_unsgn_imul+0xfa>
mov    -0x18(%rbp),%rax
mov    (%rax),%ecx
add    -0x34(%rbp),%ecx
mov    -0x28(%rbp),%rax
mov    %ecx,(%rax)
cmpl   $0x0,-0x34(%rbp)
jle    1342 <mp_unsgn_imul+0x242>
jmp    124c <mp_unsgn_imul+0x14c>
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jle    1287 <mp_unsgn_imul+0x187>
movsd  -0x48(%rbp),%xmm0
cvtsi2sdl -0x30(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xd9a(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x30(%rbp)
mov    -0x34(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x34(%rbp)
jmp    124c <mp_unsgn_imul+0x14c>
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jl     12c7 <mp_unsgn_imul+0x1c7>
mov    -0x58(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x34(%rbp),%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,2),%dx
mov    -0x58(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
mov    -0x2c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x2c(%rbp)
jmp    1290 <mp_unsgn_imul+0x190>
mov    -0x34(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x0,-0x2c(%rbp)
jl     133d <mp_unsgn_imul+0x23d>
cvtsi2sdl -0x30(%rbp),%xmm0
movsd  0xd21(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x30(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  -0x10(%rbp),%xmm1
cvtsi2sdl -0x30(%rbp),%xmm2
mulsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %ax,%dx
mov    -0x58(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
mov    -0x2c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x2c(%rbp)
jmp    12d0 <mp_unsgn_imul+0x1d0>
jmp    1342 <mp_unsgn_imul+0x242>
pop    %rbp
ret

x86 assembly (clang, -O2):
<mp_unsgn_imul>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
lea    0x4(%rdx),%r8
movsd  0xeec(%rip),%xmm3
movapd %xmm3,%xmm2
divsd  %xmm0,%xmm2
test   %edi,%edi
jle    118e <mp_unsgn_imul+0x8e>
lea    -0x1(%rdi),%eax
movapd 0xee5(%rip),%xmm4
xorpd  %xmm0,%xmm4
xor    %r10d,%r10d
movsd  0xece(%rip),%xmm5
nopw   0x0(%rax,%rax,1)
movswl 0x4(%rsi,%rax,2),%ecx
xorps  %xmm6,%xmm6
cvtsi2sd %ecx,%xmm6
xorps  %xmm7,%xmm7
cvtsi2sd %r10d,%xmm7
mulsd  %xmm1,%xmm6
addsd  %xmm7,%xmm6
addsd  %xmm5,%xmm6
movapd %xmm2,%xmm7
mulsd  %xmm6,%xmm7
cvttsd2si %xmm7,%r10d
cvttpd2dq %xmm7,%xmm7
cvtdq2pd %xmm7,%xmm7
mulsd  %xmm4,%xmm7
addsd  %xmm6,%xmm7
cvttsd2si %xmm7,%ecx
mov    %cx,(%r8,%rax,2)
add    $0xffffffffffffffff,%rax
jb     1140 <mp_unsgn_imul+0x40>
jmp    1191 <mp_unsgn_imul+0x91>
xor    %r10d,%r10d
xorps  %xmm1,%xmm1
cvtsi2sd %r10d,%xmm1
addsd  0xe67(%rip),%xmm1
ucomisd 0xe57(%rip),%xmm1
jbe    123b <mp_unsgn_imul+0x13b>
movslq %edi,%r9
lea    (%rdx,%r9,2),%rax
add    $0x2,%rax
xor    %r11d,%r11d
mov    %r8,%rcx
mulsd  %xmm2,%xmm1
add    $0x1,%r11
add    $0xfffffffffffffffe,%rax
add    $0xfffffffffffffffe,%rcx
ucomisd %xmm3,%xmm1
ja     11c0 <mp_unsgn_imul+0xc0>
mov    (%rsi),%esi
add    %r11d,%esi
mov    %esi,(%rdx)
lea    -0x1(%r11),%esi
cmp    %edi,%esi
jge    1244 <mp_unsgn_imul+0x144>
cmp    %edi,%r11d
jge    1236 <mp_unsgn_imul+0x136>
lea    -0x1(%r9),%rdi
cmp    %r11,%rdi
mov    %r11,%rsi
cmovl  %rdi,%rsi
mov    %r9,%r12
sub    %rsi,%r12
cmp    $0x10,%r12
jae    12d8 <mp_unsgn_imul+0x1d8>
mov    %r9,%rax
sub    %r11,%rax
lea    (%rdx,%rax,2),%rax
add    $0x2,%rax
cs nopw 0x0(%rax,%rax,1)
movzwl (%rax),%ecx
mov    %cx,0x2(%rdx,%r9,2)
add    $0xffffffffffffffff,%r9
add    $0xfffffffffffffffe,%rax
cmp    %r11,%r9
jg     1220 <mp_unsgn_imul+0x120>
mov    %r11d,%edi
jmp    1272 <mp_unsgn_imul+0x172>
mov    (%rsi),%eax
mov    %eax,(%rdx)
jmp    12cf <mp_unsgn_imul+0x1cf>
movsd  0xdbc(%rip),%xmm1
nopl   0x0(%rax)
xorps  %xmm3,%xmm3
cvtsi2sd %r10d,%xmm3
mulsd  %xmm2,%xmm3
addsd  %xmm1,%xmm3
cvttsd2si %xmm3,%r10d
add    $0xffffffff,%r11d
cmp    %edi,%r11d
jg     1250 <mp_unsgn_imul+0x150>
test   %edi,%edi
jle    12cf <mp_unsgn_imul+0x1cf>
xorpd  0xd96(%rip),%xmm0
mov    %edi,%eax
add    $0x1,%rax
movsd  0xd80(%rip),%xmm1
nopl   0x0(%rax,%rax,1)
xorps  %xmm3,%xmm3
cvtsi2sd %r10d,%xmm3
addsd  %xmm1,%xmm3
movapd %xmm2,%xmm4
mulsd  %xmm3,%xmm4
cvttsd2si %xmm4,%r10d
cvttpd2dq %xmm4,%xmm4
cvtdq2pd %xmm4,%xmm4
mulsd  %xmm0,%xmm4
addsd  %xmm3,%xmm4
cvttsd2si %xmm4,%ecx
lea    -0x2(%rax),%edx
mov    %cx,(%r8,%rdx,2)
add    $0xffffffffffffffff,%rax
cmp    $0x1,%rax
ja     1290 <mp_unsgn_imul+0x190>
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret
cmp    %r11,%rdi
mov    %r11,%rsi
cmovl  %rdi,%rsi
not    %rsi
add    %r9,%rsi
lea    (%rdx,%r9,2),%rbp
add    $0x2,%rbp
add    %rsi,%rsi
setb   %r14b
mov    %rbp,%rbx
sub    %rsi,%rbx
cmp    %rdi,%r11
mov    %rdi,%rsi
cmovl  %r11,%rsi
cmp    %rbp,%rbx
ja     1208 <mp_unsgn_imul+0x108>
test   %r14b,%r14b
jne    1208 <mp_unsgn_imul+0x108>
lea    (%rcx,%rsi,2),%rsi
cmp    %rax,%rsi
ja     1208 <mp_unsgn_imul+0x108>
test   %r14b,%r14b
jne    1208 <mp_unsgn_imul+0x108>
cmp    %r11,%rdi
mov    %r11,%rsi
cmovl  %rdi,%rsi
lea    0x4(%rdx,%rsi,2),%rsi
mov    %rdi,%rbx
cmovg  %r11,%rbx
add    $0x2,%rax
cmp    %rax,%rsi
jae    1363 <mp_unsgn_imul+0x263>
lea    (%rdx,%r9,2),%rax
add    $0x4,%rax
lea    (%rcx,%rbx,2),%rcx
cmp    %rax,%rcx
jb     1208 <mp_unsgn_imul+0x108>
mov    %r12,%r14
and    $0xfffffffffffffff0,%r14
lea    -0x10(%r14),%rax
mov    %rax,%r15
shr    $0x4,%r15
add    $0x1,%r15
test   %rax,%rax
je     13fe <mp_unsgn_imul+0x2fe>
mov    $0xfffffffffffffffa,%rbx
sub    %r11,%rbx
lea    (%r9,%r9,1),%rax
add    %rdx,%rax
cmp    %rdi,%r11
cmovl  %r11,%rdi
mov    %r9,%rcx
sub    %rdi,%rcx
and    $0xfffffffffffffff0,%rcx
add    $0xfffffffffffffff0,%rcx
shr    $0x4,%rcx
add    $0x1,%rcx
and    $0xfffffffffffffffe,%rcx
mov    $0xfffffffffffffffa,%rdi
nopl   0x0(%rax,%rax,1)
movups -0x10(%rax,%rbx,2),%xmm1
movups (%rax,%rbx,2),%xmm3
movups %xmm3,(%rax,%rdi,2)
movups %xmm1,-0x10(%rax,%rdi,2)
movupd -0x30(%rax,%rbx,2),%xmm1
movupd -0x20(%rax,%rbx,2),%xmm3
movupd %xmm3,-0x20(%rax,%rdi,2)
movupd %xmm1,-0x30(%rax,%rdi,2)
add    $0xffffffffffffffe0,%rdi
add    $0xffffffffffffffe0,%rbx
add    $0xfffffffffffffffe,%rcx
jne    13c0 <mp_unsgn_imul+0x2c0>
add    $0x5,%rdi
jmp    1405 <mp_unsgn_imul+0x305>
mov    $0xffffffffffffffff,%rdi
test   $0x1,%r15b
je     1430 <mp_unsgn_imul+0x330>
add    %r9,%rdi
mov    %rdi,%rax
sub    %r11,%rax
movupd -0x1e(%r8,%rax,2),%xmm1
movupd -0xe(%r8,%rax,2),%xmm3
movupd %xmm3,-0xe(%r8,%rdi,2)
movupd %xmm1,-0x1e(%r8,%rdi,2)
cmp    %r14,%r12
je     1236 <mp_unsgn_imul+0x136>
sub    %r14,%r9
jmp    1208 <mp_unsgn_imul+0x108>

LLVM IR (-O0):
define void @mp_unsgn_imul(i32 %0, double %1, i32* %2, double %3, i32* %4) {
%6 = alloca i32
%7 = alloca double
%8 = alloca i32*
%9 = alloca double
%10 = alloca i32*
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca double
%15 = alloca double
%16 = alloca i16*
%17 = alloca i16*
store i32 %0, i32* %6
store double %1, double* %7
store i32* %2, i32** %8
store double %3, double* %9
store i32* %4, i32** %10
%18 = load i32*, i32** %8
%19 = getelementptr inbounds i32, i32* %18, i64 1
%20 = bitcast i32* %19 to i16*
store i16* %20, i16** %16
%21 = load i32*, i32** %10
%22 = getelementptr inbounds i32, i32* %21, i64 1
%23 = bitcast i32* %22 to i16*
store i16* %23, i16** %17
%24 = load double, double* %7
%25 = fdiv double 1.000000e+00, %24
store double %25, double* %15
store i32 0, i32* %12
%26 = load i32, i32* %6
%27 = sub nsw i32 %26, 1
store i32 %27, i32* %11
br label %28
28: ; preds = %59, %5
%29 = load i32, i32* %11
%30 = icmp sge i32 %29, 0
br i1 %30, label %31, label %62
31: ; preds = %28
%32 = load double, double* %9
%33 = load i16*, i16** %16
%34 = load i32, i32* %11
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds i16, i16* %33, i64 %35
%37 = load i16, i16* %36
%38 = sext i16 %37 to i32
%39 = sitofp i32 %38 to double
%40 = load i32, i32* %12
%41 = sitofp i32 %40 to double
%42 = call double @llvm.fmuladd.f64(double %32, double %39, double %41)
%43 = fadd double %42, 5.000000e-01
store double %43, double* %14
%44 = load double, double* %15
%45 = load double, double* %14
%46 = fmul double %44, %45
%47 = fptosi double %46 to i32
store i32 %47, i32* %12
%48 = load double, double* %14
%49 = load double, double* %7
%50 = load i32, i32* %12
%51 = sitofp i32 %50 to double
%52 = fneg double %49
%53 = call double @llvm.fmuladd.f64(double %52, double %51, double %48)
%54 = fptosi double %53 to i16
%55 = load i16*, i16** %17
%56 = load i32, i32* %11
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds i16, i16* %55, i64 %57
store i16 %54, i16* %58
br label %59
59: ; preds = %31
%60 = load i32, i32* %11
%61 = add nsw i32 %60, -1
store i32 %61, i32* %11
br label %28
62: ; preds = %28
store i32 0, i32* %13
%63 = load i32, i32* %12
%64 = sitofp i32 %63 to double
%65 = fadd double %64, 5.000000e-01
store double %65, double* %14
br label %66
66: ; preds = %69, %62
%67 = load double, double* %14
%68 = fcmp ogt double %67, 1.000000e+00
br i1 %68, label %69, label %75
69: ; preds = %66
%70 = load double, double* %15
%71 = load double, double* %14
%72 = fmul double %71, %70
store double %72, double* %14
%73 = load i32, i32* %13
%74 = add nsw i32 %73, 1
store i32 %74, i32* %13
br label %66
75: ; preds = %66
%76 = load i32*, i32** %8
%77 = getelementptr inbounds i32, i32* %76, i64 0
%78 = load i32, i32* %77
%79 = load i32, i32* %13
%80 = add nsw i32 %78, %79
%81 = load i32*, i32** %10
%82 = getelementptr inbounds i32, i32* %81, i64 0
store i32 %80, i32* %82
%83 = load i32, i32* %13
%84 = icmp sgt i32 %83, 0
br i1 %84, label %85, label %149
85: ; preds = %75
br label %86
86: ; preds = %90, %85
%87 = load i32, i32* %13
%88 = load i32, i32* %6
%89 = icmp sgt i32 %87, %88
br i1 %89, label %90, label %98
90: ; preds = %86
%91 = load double, double* %15
%92 = load i32, i32* %12
%93 = sitofp i32 %92 to double
%94 = call double @llvm.fmuladd.f64(double %91, double %93, double 5.000000e-01)
%95 = fptosi double %94 to i32
store i32 %95, i32* %12
%96 = load i32, i32* %13
%97 = add nsw i32 %96, -1
store i32 %97, i32* %13
br label %86
98: ; preds = %86
%99 = load i32, i32* %6
%100 = sub nsw i32 %99, 1
store i32 %100, i32* %11
br label %101
101: ; preds = %117, %98
%102 = load i32, i32* %11
%103 = load i32, i32* %13
%104 = icmp sge i32 %102, %103
br i1 %104, label %105, label %120
105: ; preds = %101
%106 = load i16*, i16** %17
%107 = load i32, i32* %11
%108 = load i32, i32* %13
%109 = sub nsw i32 %107, %108
%110 = sext i32 %109 to i64
%111 = getelementptr inbounds i16, i16* %106, i64 %110
%112 = load i16, i16* %111
%113 = load i16*, i16** %17
%114 = load i32, i32* %11
%115 = sext i32 %114 to i64
%116 = getelementptr inbounds i16, i16* %113, i64 %115
store i16 %112, i16* %116
br label %117
117: ; preds = %105
%118 = load i32, i32* %11
%119 = add nsw i32 %118, -1
store i32 %119, i32* %11
br label %101
120: ; preds = %101
%121 = load i32, i32* %13
%122 = sub nsw i32 %121, 1
store i32 %122, i32* %11
br label %123
123: ; preds = %145, %120
%124 = load i32, i32* %11
%125 = icmp sge i32 %124, 0
br i1 %125, label %126, label %148
126: ; preds = %123
%127 = load i32, i32* %12
%128 = sitofp i32 %127 to double
%129 = fadd double %128, 5.000000e-01
store double %129, double* %14
%130 = load double, double* %15
%131 = load double, double* %14
%132 = fmul double %130, %131
%133 = fptosi double %132 to i32
store i32 %133, i32* %12
%134 = load double, double* %14
%135 = load double, double* %7
%136 = load i32, i32* %12
%137 = sitofp i32 %136 to double
%138 = fneg double %135
%139 = call double @llvm.fmuladd.f64(double %138, double %137, double %134)
%140 = fptosi double %139 to i16
%141 = load i16*, i16** %17
%142 = load i32, i32* %11
%143 = sext i32 %142 to i64
%144 = getelementptr inbounds i16, i16* %141, i64 %143
store i16 %140, i16* %144
br label %145
145: ; preds = %126
%146 = load i32, i32* %11
%147 = add nsw i32 %146, -1
store i32 %147, i32* %11
br label %123
148: ; preds = %123
br label %149
149: ; preds = %148, %75
ret void
}

LLVM IR (-O2):
define void @mp_unsgn_imul(i32 %0, double %1, i32* nocapture %2, double %3, i32* nocapture %4) {
%6 = bitcast i32* %4 to i8*
%7 = getelementptr inbounds i32, i32* %2, i64 1
%8 = bitcast i32* %7 to i16*
%9 = getelementptr inbounds i32, i32* %4, i64 1
%10 = bitcast i32* %9 to i16*
%11 = fdiv double 1.000000e+00, %1
%12 = icmp sgt i32 %0, 0
br i1 %12, label %13, label %34
13: ; preds = %5
%14 = add nsw i32 %0, -1
%15 = fneg double %1
%16 = zext i32 %14 to i64
br label %17
17: ; preds = %13, %17
%18 = phi i64 [ %16, %13 ], [ %32, %17 ]
%19 = phi i32 [ 0, %13 ], [ %27, %17 ]
%20 = getelementptr inbounds i16, i16* %8, i64 %18
%21 = load i16, i16* %20
%22 = sitofp i16 %21 to double
%23 = sitofp i32 %19 to double
%24 = tail call double @llvm.fmuladd.f64(double %3, double %22, double %23)
%25 = fadd double %24, 5.000000e-01
%26 = fmul double %11, %25
%27 = fptosi double %26 to i32
%28 = sitofp i32 %27 to double
%29 = tail call double @llvm.fmuladd.f64(double %15, double %28, double %25)
%30 = fptosi double %29 to i16
%31 = getelementptr inbounds i16, i16* %10, i64 %18
store i16 %30, i16* %31
%32 = add nsw i64 %18, -1
%33 = icmp eq i64 %18, 0
br i1 %33, label %34, label %17
34: ; preds = %17, %5
%35 = phi i32 [ 0, %5 ], [ %27, %17 ]
%36 = sitofp i32 %35 to double
%37 = fadd double %36, 5.000000e-01
%38 = fcmp ogt double %37, 1.000000e+00
br i1 %38, label %41, label %39
39: ; preds = %34
%40 = load i32, i32* %2
store i32 %40, i32* %4
br label %223
41: ; preds = %34, %41
%42 = phi i32 [ %45, %41 ], [ 0, %34 ]
%43 = phi double [ %44, %41 ], [ %37, %34 ]
%44 = fmul double %11, %43
%45 = add nuw nsw i32 %42, 1
%46 = fcmp ogt double %44, 1.000000e+00
br i1 %46, label %41, label %47
47: ; preds = %41
%48 = load i32, i32* %2
%49 = add nsw i32 %48, %45
store i32 %49, i32* %4
%50 = icmp slt i32 %42, %0
br i1 %50, label %51, label %186
51: ; preds = %47
%52 = icmp slt i32 %45, %0
br i1 %52, label %53, label %196
53: ; preds = %51
%54 = sext i32 %0 to i64
%55 = zext i32 %45 to i64
%56 = add nsw i64 %54, -1
%57 = call i64 @llvm.smin.i64(i64 %56, i64 %55)
%58 = sub i64 %54, %57
%59 = icmp ult i64 %58, 16
br i1 %59, label %184, label %60
60: ; preds = %53
%61 = add nsw i64 %54, -1
%62 = call i64 @llvm.smin.i64(i64 %61, i64 %55)
%63 = xor i64 %62, -1
%64 = add i64 %63, %54
%65 = shl nsw i64 %54, 1
%66 = add nsw i64 %65, 2
%67 = getelementptr i8, i8* %6, i64 %66
%68 = call { i64, i1 } @llvm.umul.with.overflow.i64(i64 %64, i64 2)
%69 = extractvalue { i64, i1 } %68, 0
%70 = extractvalue { i64, i1 } %68, 1
%71 = sub i64 0, %69
%72 = getelementptr i8, i8* %67, i64 %71
%73 = icmp ugt i8* %72, %67
%74 = or i1 %73, %70
%75 = add nsw i64 %65, 2
%76 = mul nsw i64 %55, -2
%77 = add nsw i64 %76, %75
%78 = getelementptr i8, i8* %6, i64 %77
%79 = call { i64, i1 } @llvm.umul.with.overflow.i64(i64 %64, i64 2)
%80 = extractvalue { i64, i1 } %79, 0
%81 = extractvalue { i64, i1 } %79, 1
%82 = sub i64 0, %80
%83 = getelementptr i8, i8* %78, i64 %82
%84 = icmp ugt i8* %83, %78
%85 = or i1 %84, %81
%86 = or i1 %74, %85
br i1 %86, label %184, label %87
87: ; preds = %60
%88 = getelementptr i32, i32* %4, i64 1
%89 = bitcast i32* %88 to i8*
%90 = add nsw i64 %54, -1
%91 = call i64 @llvm.smin.i64(i64 %90, i64 %55)
%92 = shl nsw i64 %91, 1
%93 = getelementptr i8, i8* %89, i64 %92
%94 = getelementptr i32, i32* %4, i64 1
%95 = bitcast i32* %94 to i8*
%96 = shl nsw i64 %54, 1
%97 = getelementptr i8, i8* %95, i64 %96
%98 = getelementptr i32, i32* %4, i64 1
%99 = bitcast i32* %98 to i8*
%100 = shl nuw nsw i64 %55, 1
%101 = sub i64 %92, %100
%102 = getelementptr i8, i8* %99, i64 %101
%103 = getelementptr i32, i32* %4, i64 1
%104 = bitcast i32* %103 to i8*
%105 = sub nsw i64 %96, %100
%106 = getelementptr i8, i8* %104, i64 %105
%107 = icmp ult i8* %93, %106
%108 = icmp ult i8* %102, %97
%109 = and i1 %107, %108
br i1 %109, label %184, label %110
110: ; preds = %87
%111 = and i64 %58, -16
%112 = sub i64 %54, %111
%113 = add i64 %111, -16
%114 = lshr exact i64 %113, 4
%115 = add nuw nsw i64 %114, 1
%116 = and i64 %115, 1
%117 = icmp eq i64 %113, 0
br i1 %117, label %162, label %118
118: ; preds = %110
%119 = and i64 %115, 2305843009213693950
br label %120
120: ; preds = %120, %118
%121 = phi i64 [ 0, %118 ], [ %157, %120 ]
%122 = phi i64 [ 0, %118 ], [ %158, %120 ]
%123 = xor i64 %121, -1
%124 = add i64 %123, %54
%125 = sub nsw i64 %124, %55
%126 = getelementptr inbounds i16, i16* %10, i64 %125
%127 = getelementptr inbounds i16, i16* %126, i64 -7
%128 = bitcast i16* %127 to <8 x i16>*
%129 = load <8 x i16>, <8 x i16>* %128
%130 = getelementptr inbounds i16, i16* %126, i64 -8
%131 = getelementptr inbounds i16, i16* %130, i64 -7
%132 = bitcast i16* %131 to <8 x i16>*
%133 = load <8 x i16>, <8 x i16>* %132
%134 = getelementptr inbounds i16, i16* %10, i64 %124
%135 = getelementptr inbounds i16, i16* %134, i64 -7
%136 = bitcast i16* %135 to <8 x i16>*
store <8 x i16> %129, <8 x i16>* %136
%137 = getelementptr inbounds i16, i16* %134, i64 -8
%138 = getelementptr inbounds i16, i16* %137, i64 -7
%139 = bitcast i16* %138 to <8 x i16>*
store <8 x i16> %133, <8 x i16>* %139
%140 = sub nuw nsw i64 -17, %121
%141 = add i64 %140, %54
%142 = sub nsw i64 %141, %55
%143 = getelementptr inbounds i16, i16* %10, i64 %142
%144 = getelementptr inbounds i16, i16* %143, i64 -7
%145 = bitcast i16* %144 to <8 x i16>*
%146 = load <8 x i16>, <8 x i16>* %145
%147 = getelementptr inbounds i16, i16* %143, i64 -8
%148 = getelementptr inbounds i16, i16* %147, i64 -7
%149 = bitcast i16* %148 to <8 x i16>*
%150 = load <8 x i16>, <8 x i16>* %149
%151 = getelementptr inbounds i16, i16* %10, i64 %141
%152 = getelementptr inbounds i16, i16* %151, i64 -7
%153 = bitcast i16* %152 to <8 x i16>*
store <8 x i16> %146, <8 x i16>* %153
%154 = getelementptr inbounds i16, i16* %151, i64 -8
%155 = getelementptr inbounds i16, i16* %154, i64 -7
%156 = bitcast i16* %155 to <8 x i16>*
store <8 x i16> %150, <8 x i16>* %156
%157 = add nuw i64 %121, 32
%158 = add i64 %122, 2
%159 = icmp eq i64 %158, %119
br i1 %159, label %160, label %120
160: ; preds = %120
%161 = sub i64 -33, %121
br label %162
162: ; preds = %160, %110
%163 = phi i64 [ -1, %110 ], [ %161, %160 ]
%164 = icmp eq i64 %116, 0
br i1 %164, label %182, label %165
165: ; preds = %162
%166 = add i64 %163, %54
%167 = sub nsw i64 %166, %55
%168 = getelementptr inbounds i16, i16* %10, i64 %167
%169 = getelementptr inbounds i16, i16* %168, i64 -7
%170 = bitcast i16* %169 to <8 x i16>*
%171 = load <8 x i16>, <8 x i16>* %170
%172 = getelementptr inbounds i16, i16* %168, i64 -8
%173 = getelementptr inbounds i16, i16* %172, i64 -7
%174 = bitcast i16* %173 to <8 x i16>*
%175 = load <8 x i16>, <8 x i16>* %174
%176 = getelementptr inbounds i16, i16* %10, i64 %166
%177 = getelementptr inbounds i16, i16* %176, i64 -7
%178 = bitcast i16* %177 to <8 x i16>*
store <8 x i16> %171, <8 x i16>* %178
%179 = getelementptr inbounds i16, i16* %176, i64 -8
%180 = getelementptr inbounds i16, i16* %179, i64 -7
%181 = bitcast i16* %180 to <8 x i16>*
store <8 x i16> %175, <8 x i16>* %181
br label %182
182: ; preds = %162, %165
%183 = icmp eq i64 %58, %111
br i1 %183, label %196, label %184
184: ; preds = %87, %60, %53, %182
%185 = phi i64 [ %54, %87 ], [ %54, %60 ], [ %54, %53 ], [ %112, %182 ]
br label %201
186: ; preds = %47, %186
%187 = phi i32 [ %191, %186 ], [ %35, %47 ]
%188 = phi i32 [ %192, %186 ], [ %45, %47 ]
%189 = sitofp i32 %187 to double
%190 = tail call double @llvm.fmuladd.f64(double %11, double %189, double 5.000000e-01)
%191 = fptosi double %190 to i32
%192 = add nsw i32 %188, -1
%193 = icmp sgt i32 %192, %0
br i1 %193, label %186, label %194
194: ; preds = %186
%195 = icmp sgt i32 %0, 0
br i1 %195, label %196, label %223
196: ; preds = %201, %182, %51, %194
%197 = phi i32 [ %0, %194 ], [ %45, %51 ], [ %45, %182 ], [ %45, %201 ]
%198 = phi i32 [ %191, %194 ], [ %35, %51 ], [ %35, %182 ], [ %35, %201 ]
%199 = fneg double %1
%200 = zext i32 %197 to i64
br label %209
201: ; preds = %184, %201
%202 = phi i64 [ %203, %201 ], [ %185, %184 ]
%203 = add nsw i64 %202, -1
%204 = sub nsw i64 %203, %55
%205 = getelementptr inbounds i16, i16* %10, i64 %204
%206 = load i16, i16* %205
%207 = getelementptr inbounds i16, i16* %10, i64 %203
store i16 %206, i16* %207
%208 = icmp sgt i64 %203, %55
br i1 %208, label %201, label %196
209: ; preds = %196, %209
%210 = phi i64 [ %200, %196 ], [ %212, %209 ]
%211 = phi i32 [ %198, %196 ], [ %216, %209 ]
%212 = add nsw i64 %210, -1
%213 = sitofp i32 %211 to double
%214 = fadd double %213, 5.000000e-01
%215 = fmul double %11, %214
%216 = fptosi double %215 to i32
%217 = sitofp i32 %216 to double
%218 = tail call double @llvm.fmuladd.f64(double %199, double %217, double %214)
%219 = fptosi double %218 to i16
%220 = and i64 %212, 4294967295
%221 = getelementptr inbounds i16, i16* %10, i64 %220
store i16 %219, i16* %221
%222 = icmp ugt i64 %210, 1
br i1 %222, label %209, label %223
223: ; preds = %209, %39, %194
ret void
}