void generate(int k, int n, int *r, double *ranvec, double **b) {
  double unif(), cum;
  int i, l, nleft;
  for (i = 1; i <= k - 1; i++)
    ranvec[i] = unif();
  nleft = n;
  for (l = 1; l < k; l++) {
    cum = 0.0;
    for (i = 1; i <= nleft; i++) {
      cum += b[k - l][nleft - i] / (i * b[k - l + 1][nleft]);
      if (cum >= ranvec[l])
        break;
    }
    r[l] = i;
    nleft -= i;
  }
  r[k] = nleft;
}


x86 assembly (clang, -O0):
<generate>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     1164 <generate+0x54>
mov    $0x0,%al
call   1030 <unif@plt>
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1131 <generate+0x21>
mov    -0x8(%rbp),%eax
mov    %eax,-0x34(%rbp)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    123e <generate+0x12e>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jg     1217 <generate+0x107>
mov    -0x20(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    -0x30(%rbp),%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x34(%rbp),%ecx
sub    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x2c(%rbp),%xmm1
mov    -0x20(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    -0x30(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
divsd  %xmm1,%xmm0
addsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
jb     1204 <generate+0xf4>
jmp    1217 <generate+0x107>
jmp    1209 <generate+0xf9>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    118c <generate+0x7c>
mov    -0x2c(%rbp),%edx
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x2c(%rbp),%ecx
mov    -0x34(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x34(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1171 <generate+0x61>
mov    -0x34(%rbp),%edx
mov    -0x10(%rbp),%rax
movslq -0x4(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<generate>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    %rdx,%r14
mov    %esi,%ebx
cmp    $0x2,%edi
jl     11e0 <generate+0xd0>
mov    %r8,%r15
mov    %rcx,%r12
mov    %edi,0x4(%rsp)
mov    %edi,%ebp
mov    $0x1,%r13d
nopl   0x0(%rax,%rax,1)
xor    %eax,%eax
call   1030 <unif@plt>
movsd  %xmm0,(%r12,%r13,8)
add    $0x1,%r13
cmp    %r13,%rbp
jne    1140 <generate+0x30>
mov    0x4(%rsp),%edi
cmp    $0x1,%edi
jle    11e0 <generate+0xd0>
mov    $0x1,%eax
jmp    1184 <generate+0x74>
nopw   0x0(%rax,%rax,1)
mov    $0x1,%ecx
mov    %ecx,(%r14,%rax,4)
sub    %ecx,%ebx
add    $0x1,%rax
cmp    %rbp,%rax
je     11e3 <generate+0xd3>
test   %ebx,%ebx
jle    1170 <generate+0x60>
mov    %rbp,%rcx
sub    %rax,%rcx
mov    (%r15,%rcx,8),%rdx
mov    0x8(%r15,%rcx,8),%rcx
mov    %ebx,%esi
movsd  (%rcx,%rsi,8),%xmm0
movsd  (%r12,%rax,8),%xmm1
lea    0x1(%rbx),%edi
xorpd  %xmm2,%xmm2
mov    $0x1,%ecx
movapd %xmm2,%xmm3
movsd  -0x8(%rdx,%rsi,8),%xmm2
xorps  %xmm4,%xmm4
cvtsi2sd %ecx,%xmm4
mulsd  %xmm0,%xmm4
divsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
ucomisd %xmm1,%xmm2
jae    1175 <generate+0x65>
add    $0x1,%ecx
add    $0xffffffffffffffff,%rsi
jne    11b0 <generate+0xa0>
mov    %edi,%ecx
jmp    1175 <generate+0x65>
movslq %edi,%rbp
mov    %ebx,(%r14,%rbp,4)
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @generate(i32 %0, i32 %1, i32* %2, double* %3, double** %4) {
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32*
%9 = alloca double*
%10 = alloca double**
%11 = alloca double
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
store i32 %0, i32* %6
store i32 %1, i32* %7
store i32* %2, i32** %8
store double* %3, double** %9
store double** %4, double*** %10
store i32 1, i32* %12
br label %15
15: ; preds = %26, %5
%16 = load i32, i32* %12
%17 = load i32, i32* %6
%18 = sub nsw i32 %17, 1
%19 = icmp sle i32 %16, %18
br i1 %19, label %20, label %29
20: ; preds = %15
%21 = call double (...) @unif()
%22 = load double*, double** %9
%23 = load i32, i32* %12
%24 = sext i32 %23 to i64
%25 = getelementptr inbounds double, double* %22, i64 %24
store double %21, double* %25
br label %26
26: ; preds = %20
%27 = load i32, i32* %12
%28 = add nsw i32 %27, 1
store i32 %28, i32* %12
br label %15
29: ; preds = %15
%30 = load i32, i32* %7
store i32 %30, i32* %14
store i32 1, i32* %13
br label %31
31: ; preds = %93, %29
%32 = load i32, i32* %13
%33 = load i32, i32* %6
%34 = icmp slt i32 %32, %33
br i1 %34, label %35, label %96
35: ; preds = %31
store double 0.000000e+00, double* %11
store i32 1, i32* %12
br label %36
36: ; preds = %81, %35
%37 = load i32, i32* %12
%38 = load i32, i32* %14
%39 = icmp sle i32 %37, %38
br i1 %39, label %40, label %84
40: ; preds = %36
%41 = load double**, double*** %10
%42 = load i32, i32* %6
%43 = load i32, i32* %13
%44 = sub nsw i32 %42, %43
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds double*, double** %41, i64 %45
%47 = load double*, double** %46
%48 = load i32, i32* %14
%49 = load i32, i32* %12
%50 = sub nsw i32 %48, %49
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds double, double* %47, i64 %51
%53 = load double, double* %52
%54 = load i32, i32* %12
%55 = sitofp i32 %54 to double
%56 = load double**, double*** %10
%57 = load i32, i32* %6
%58 = load i32, i32* %13
%59 = sub nsw i32 %57, %58
%60 = add nsw i32 %59, 1
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds double*, double** %56, i64 %61
%63 = load double*, double** %62
%64 = load i32, i32* %14
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds double, double* %63, i64 %65
%67 = load double, double* %66
%68 = fmul double %55, %67
%69 = fdiv double %53, %68
%70 = load double, double* %11
%71 = fadd double %70, %69
store double %71, double* %11
%72 = load double, double* %11
%73 = load double*, double** %9
%74 = load i32, i32* %13
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds double, double* %73, i64 %75
%77 = load double, double* %76
%78 = fcmp oge double %72, %77
br i1 %78, label %79, label %80
79: ; preds = %40
br label %84
80: ; preds = %40
br label %81
81: ; preds = %80
%82 = load i32, i32* %12
%83 = add nsw i32 %82, 1
store i32 %83, i32* %12
br label %36
84: ; preds = %79, %36
%85 = load i32, i32* %12
%86 = load i32*, i32** %8
%87 = load i32, i32* %13
%88 = sext i32 %87 to i64
%89 = getelementptr inbounds i32, i32* %86, i64 %88
store i32 %85, i32* %89
%90 = load i32, i32* %12
%91 = load i32, i32* %14
%92 = sub nsw i32 %91, %90
store i32 %92, i32* %14
br label %93
93: ; preds = %84
%94 = load i32, i32* %13
%95 = add nsw i32 %94, 1
store i32 %95, i32* %13
br label %31
96: ; preds = %31
%97 = load i32, i32* %14
%98 = load i32*, i32** %8
%99 = load i32, i32* %6
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds i32, i32* %98, i64 %100
store i32 %97, i32* %101
ret void
}

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