void mp_unsgn_idiv(int n, double dradix, int in1[], double din2, int out[]) {
  int j, ix, carry, shift;
  double x, d1_in2;
  d1_in2 = 1.0 / din2;
  shift = 0;
  x = 0;
  do {
    shift++;
    x *= dradix;
    if (shift <= n) {
      x += in1[shift];
    }
  } while (x < din2 - 0.5);
  x += 0.5;
  ix = (int)(d1_in2 * x);
  carry = (int)(x - din2 * ix);
  out[1] = ix;
  shift--;
  out[0] = in1[0] - shift;
  if (shift >= n) {
    shift = n - 1;
  }
  for (j = 2; j <= n - shift; j++) {
    x = in1[j + shift] + dradix * carry + 0.5;
    ix = (int)(d1_in2 * x);
    carry = (int)(x - din2 * ix);
    out[j] = ix;
  }
  for (j = n - shift + 1; j <= n; j++) {
    x = dradix * carry + 0.5;
    ix = (int)(d1_in2 * x);
    carry = (int)(x - din2 * ix);
    out[j] = ix;
  }
}


x86 assembly (clang, -O0):
<mp_unsgn_idiv>:
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)
movsd  0xedf(%rip),%xmm0
divsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
movsd  -0x10(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     1175 <mp_unsgn_idiv+0x75>
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    117a <mp_unsgn_idiv+0x7a>
movsd  -0x40(%rbp),%xmm1
movsd  -0x20(%rbp),%xmm0
movsd  0xe7c(%rip),%xmm2
subsd  %xmm2,%xmm0
ucomisd %xmm1,%xmm0
ja     113a <mp_unsgn_idiv+0x3a>
movsd  0xe66(%rip),%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x48(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x30(%rbp)
movsd  -0x40(%rbp),%xmm1
movsd  -0x20(%rbp),%xmm0
cvtsi2sdl -0x30(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x34(%rbp)
mov    -0x30(%rbp),%ecx
mov    -0x28(%rbp),%rax
mov    %ecx,0x4(%rax)
mov    -0x38(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x38(%rbp)
mov    -0x18(%rbp),%rax
mov    (%rax),%ecx
sub    -0x38(%rbp),%ecx
mov    -0x28(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x38(%rbp),%eax
cmp    -0x4(%rbp),%eax
jl     1229 <mp_unsgn_idiv+0x129>
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x38(%rbp)
movl   $0x2,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    -0x38(%rbp),%ecx
cmp    %ecx,%eax
jg     12d8 <mp_unsgn_idiv+0x1d8>
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm1
movsd  -0x10(%rbp),%xmm0
cvtsi2sdl -0x34(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  0xd9b(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x48(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x30(%rbp)
movsd  -0x40(%rbp),%xmm1
movsd  -0x20(%rbp),%xmm0
cvtsi2sdl -0x30(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x34(%rbp)
mov    -0x30(%rbp),%edx
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1230 <mp_unsgn_idiv+0x130>
mov    -0x4(%rbp),%eax
sub    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     1371 <mp_unsgn_idiv+0x271>
movsd  -0x10(%rbp),%xmm0
cvtsi2sdl -0x34(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xd02(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x48(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x30(%rbp)
movsd  -0x40(%rbp),%xmm1
movsd  -0x20(%rbp),%xmm0
cvtsi2sdl -0x30(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x34(%rbp)
mov    -0x30(%rbp),%edx
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    12e4 <mp_unsgn_idiv+0x1e4>
pop    %rbp
ret

x86 assembly (clang, -O2):
<mp_unsgn_idiv>:
movsd  0xef8(%rip),%xmm2
addsd  %xmm1,%xmm2
movslq %edi,%r11
xorpd  %xmm4,%xmm4
mov    $0xffffffff,%edi
xor    %ecx,%ecx
jmp    112d <mp_unsgn_idiv+0x2d>
nopl   0x0(%rax)
add    $0x1,%edi
lea    0x1(%rax),%rcx
ucomisd %xmm4,%xmm2
jbe    1148 <mp_unsgn_idiv+0x48>
mov    %rcx,%rax
mulsd  %xmm0,%xmm4
cmp    %r11,%rcx
jge    1120 <mp_unsgn_idiv+0x20>
xorps  %xmm3,%xmm3
cvtsi2sdl 0x4(%rsi,%rax,4),%xmm3
addsd  %xmm3,%xmm4
jmp    1120 <mp_unsgn_idiv+0x20>
movsd  0xeb8(%rip),%xmm2
addsd  0xeb8(%rip),%xmm4
divsd  %xmm1,%xmm2
movapd %xmm2,%xmm3
mulsd  %xmm4,%xmm3
cvttsd2si %xmm3,%ecx
mov    %ecx,0x4(%rdx)
mov    (%rsi),%ecx
sub    %edi,%ecx
lea    -0x1(%r11),%r8d
cvttpd2dq %xmm3,%xmm3
cvtdq2pd %xmm3,%xmm5
cmp    %r11,%rax
movapd 0xe9a(%rip),%xmm3
xorpd  %xmm1,%xmm3
mulsd  %xmm1,%xmm5
subsd  %xmm5,%xmm4
mov    %ecx,(%rdx)
cmovl  %edi,%r8d
mov    %r11d,%r9d
sub    %r8d,%r9d
cvttsd2si %xmm4,%ecx
cmp    $0x2,%r9d
jl     1207 <mp_unsgn_idiv+0x107>
movslq %r8d,%rax
lea    0x1(%r9),%r10d
lea    (%rsi,%rax,4),%rsi
mov    $0x2,%edi
movsd  0xe50(%rip),%xmm1
xorps  %xmm4,%xmm4
cvtsi2sdl (%rsi,%rdi,4),%xmm4
xorps  %xmm5,%xmm5
cvtsi2sd %ecx,%xmm5
mulsd  %xmm0,%xmm5
addsd  %xmm4,%xmm5
addsd  %xmm1,%xmm5
movapd %xmm2,%xmm4
mulsd  %xmm5,%xmm4
cvttsd2si %xmm4,%eax
cvttpd2dq %xmm4,%xmm4
cvtdq2pd %xmm4,%xmm4
mulsd  %xmm3,%xmm4
addsd  %xmm5,%xmm4
mov    %eax,(%rdx,%rdi,4)
add    $0x1,%rdi
cvttsd2si %xmm4,%ecx
cmp    %rdi,%r10
jne    11c0 <mp_unsgn_idiv+0xc0>
test   %r8d,%r8d
jle    125c <mp_unsgn_idiv+0x15c>
movslq %r9d,%rsi
movsd  0xdf9(%rip),%xmm1
nopw   0x0(%rax,%rax,1)
xorps  %xmm4,%xmm4
cvtsi2sd %ecx,%xmm4
mulsd  %xmm0,%xmm4
addsd  %xmm1,%xmm4
movapd %xmm2,%xmm5
mulsd  %xmm4,%xmm5
cvttsd2si %xmm5,%eax
cvttpd2dq %xmm5,%xmm5
cvtdq2pd %xmm5,%xmm5
mulsd  %xmm3,%xmm5
addsd  %xmm4,%xmm5
cvttsd2si %xmm5,%ecx
mov    %eax,0x4(%rdx,%rsi,4)
add    $0x1,%rsi
cmp    %r11,%rsi
jl     1220 <mp_unsgn_idiv+0x120>
ret

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

LLVM IR (-O2):
define void @mp_unsgn_idiv(i32 %0, double %1, i32* nocapture %2, double %3, i32* nocapture writeonly %4) {
%6 = fadd double %3, -5.000000e-01
%7 = sext i32 %0 to i64
br label %8
8: ; preds = %19, %5
%9 = phi i64 [ %11, %19 ], [ 0, %5 ]
%10 = phi double [ %20, %19 ], [ 0.000000e+00, %5 ]
%11 = add nuw i64 %9, 1
%12 = fmul double %10, %1
%13 = icmp slt i64 %9, %7
br i1 %13, label %14, label %19
14: ; preds = %8
%15 = getelementptr inbounds i32, i32* %2, i64 %11
%16 = load i32, i32* %15
%17 = sitofp i32 %16 to double
%18 = fadd double %12, %17
br label %19
19: ; preds = %8, %14
%20 = phi double [ %18, %14 ], [ %12, %8 ]
%21 = fcmp olt double %20, %6
br i1 %21, label %8, label %22
22: ; preds = %19
%23 = fdiv double 1.000000e+00, %3
%24 = trunc i64 %9 to i32
%25 = fadd double %20, 5.000000e-01
%26 = fmul double %23, %25
%27 = fptosi double %26 to i32
%28 = sitofp i32 %27 to double
%29 = fneg double %3
%30 = tail call double @llvm.fmuladd.f64(double %29, double %28, double %25)
%31 = getelementptr inbounds i32, i32* %4, i64 1
store i32 %27, i32* %31
%32 = load i32, i32* %2
%33 = sub nsw i32 %32, %24
store i32 %33, i32* %4
%34 = add nsw i32 %0, -1
%35 = select i1 %13, i32 %24, i32 %34
%36 = sub i32 %0, %35
%37 = fptosi double %30 to i32
%38 = icmp slt i32 %36, 2
br i1 %38, label %44, label %39
39: ; preds = %22
%40 = sext i32 %35 to i64
%41 = add i32 %0, 1
%42 = sub i32 %41, %35
%43 = zext i32 %42 to i64
br label %49
44: ; preds = %49, %22
%45 = phi i32 [ %37, %22 ], [ %65, %49 ]
%46 = icmp sgt i32 %35, 0
br i1 %46, label %47, label %80
47: ; preds = %44
%48 = sext i32 %36 to i64
br label %67
49: ; preds = %39, %49
%50 = phi i64 [ 2, %39 ], [ %64, %49 ]
%51 = phi i32 [ %37, %39 ], [ %65, %49 ]
%52 = add nsw i64 %50, %40
%53 = getelementptr inbounds i32, i32* %2, i64 %52
%54 = load i32, i32* %53
%55 = sitofp i32 %54 to double
%56 = sitofp i32 %51 to double
%57 = tail call double @llvm.fmuladd.f64(double %1, double %56, double %55)
%58 = fadd double %57, 5.000000e-01
%59 = fmul double %23, %58
%60 = fptosi double %59 to i32
%61 = sitofp i32 %60 to double
%62 = tail call double @llvm.fmuladd.f64(double %29, double %61, double %58)
%63 = getelementptr inbounds i32, i32* %4, i64 %50
store i32 %60, i32* %63
%64 = add nuw nsw i64 %50, 1
%65 = fptosi double %62 to i32
%66 = icmp eq i64 %64, %43
br i1 %66, label %44, label %49
67: ; preds = %47, %67
%68 = phi i64 [ %48, %47 ], [ %70, %67 ]
%69 = phi i32 [ %45, %47 ], [ %77, %67 ]
%70 = add nsw i64 %68, 1
%71 = sitofp i32 %69 to double
%72 = tail call double @llvm.fmuladd.f64(double %1, double %71, double 5.000000e-01)
%73 = fmul double %23, %72
%74 = fptosi double %73 to i32
%75 = sitofp i32 %74 to double
%76 = tail call double @llvm.fmuladd.f64(double %29, double %75, double %72)
%77 = fptosi double %76 to i32
%78 = getelementptr inbounds i32, i32* %4, i64 %70
store i32 %74, i32* %78
%79 = icmp slt i64 %70, %7
br i1 %79, label %67, label %80
80: ; preds = %67, %44
ret void
}