#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
float *malloc(int);
int mint_random_binomial(float, int);
void mint_random_multinomial(float *p, int *r, int len, int n) {
  int i;
  float *pr, sumpr, scale;
  pr = malloc(len * sizeof(int));
  if (p == 0) {
    for (i = 0; i < len; i++)
      pr[i] = 1. / len;
    sumpr = 1.;
  } else {
    sumpr = 0;
    for (i = 0; i < len; i++) {
      pr[i] = p[i];
      sumpr += p[i];
    }
  }
  scale = 1;
  for (i = 0; i < n; i++) {
    r[i] = mint_random_binomial(pr[i] / (sumpr * scale), n);
    n -= r[i];
    scale -= pr[i] / sumpr;
  }
}


x86 assembly (clang, -O0):
<mint_random_multinomial>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
movslq -0x14(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1040 <malloc@plt>
mov    %rax,-0x28(%rbp)
cmpq   $0x0,-0x8(%rbp)
jne    11a9 <mint_random_multinomial+0x89>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1197 <mint_random_multinomial+0x77>
cvtsi2sdl -0x14(%rbp),%xmm1
movsd  0xe94(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x28(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    115b <mint_random_multinomial+0x3b>
movss  0xe61(%rip),%xmm0
movss  %xmm0,-0x2c(%rbp)
jmp    1208 <mint_random_multinomial+0xe8>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x2c(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1203 <mint_random_multinomial+0xe3>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    11b8 <mint_random_multinomial+0x98>
jmp    1208 <mint_random_multinomial+0xe8>
movss  0xdf0(%rip),%xmm0
movss  %xmm0,-0x30(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    1299 <mint_random_multinomial+0x179>
mov    -0x28(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  -0x2c(%rbp),%xmm1
mulss  -0x30(%rbp),%xmm1
divss  %xmm1,%xmm0
mov    -0x18(%rbp),%edi
call   1030 <mint_random_binomial@plt>
mov    %eax,%edx
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    (%rax,%rcx,4),%ecx
mov    -0x18(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x18(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
divss  -0x2c(%rbp),%xmm1
movss  -0x30(%rbp),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    121c <mint_random_multinomial+0xfc>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<mint_random_multinomial>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %ecx,%ebp
mov    %edx,%r13d
mov    %rsi,%r14
mov    %rdi,%r12
lea    0x0(,%r13,4),%edi
call   1050 <malloc@plt>
mov    %rax,%r15
test   %r12,%r12
je     1193 <mint_random_multinomial+0x63>
test   %r13d,%r13d
jle    11cc <mint_random_multinomial+0x9c>
mov    %r13d,%r13d
lea    0x0(,%r13,4),%rdx
mov    %r15,%rdi
mov    %r12,%rsi
call   1040 <memcpy@plt>
lea    -0x1(%r13),%rcx
mov    %r13d,%eax
and    $0x7,%eax
cmp    $0x7,%rcx
jae    11d4 <mint_random_multinomial+0xa4>
xorps  %xmm2,%xmm2
xor    %ecx,%ecx
jmp    1220 <mint_random_multinomial+0xf0>
movss  0xe65(%rip),%xmm2
test   %r13d,%r13d
jle    134f <mint_random_multinomial+0x21f>
cvtsi2sd %r13d,%xmm0
movsd  0xe57(%rip),%xmm1
divsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
mov    %r13d,%eax
cmp    $0x8,%r13d
jae    1243 <mint_random_multinomial+0x113>
xor    %ecx,%ecx
jmp    1340 <mint_random_multinomial+0x210>
xorps  %xmm2,%xmm2
jmp    134f <mint_random_multinomial+0x21f>
and    $0xfffffff8,%r13d
xorps  %xmm2,%xmm2
xor    %ecx,%ecx
nopl   (%rax)
addss  (%r12,%rcx,4),%xmm2
addss  0x4(%r12,%rcx,4),%xmm2
addss  0x8(%r12,%rcx,4),%xmm2
addss  0xc(%r12,%rcx,4),%xmm2
addss  0x10(%r12,%rcx,4),%xmm2
addss  0x14(%r12,%rcx,4),%xmm2
addss  0x18(%r12,%rcx,4),%xmm2
addss  0x1c(%r12,%rcx,4),%xmm2
add    $0x8,%rcx
cmp    %rcx,%r13
jne    11e0 <mint_random_multinomial+0xb0>
test   %rax,%rax
je     134f <mint_random_multinomial+0x21f>
lea    (%r12,%rcx,4),%rcx
xor    %edx,%edx
nop
addss  (%rcx,%rdx,4),%xmm2
add    $0x1,%rdx
cmp    %rdx,%rax
jne    1230 <mint_random_multinomial+0x100>
jmp    134f <mint_random_multinomial+0x21f>
mov    %eax,%ecx
and    $0xfffffff8,%ecx
movaps %xmm0,%xmm1
shufps $0x0,%xmm0,%xmm1
lea    -0x8(%rcx),%rdi
mov    %rdi,%rsi
shr    $0x3,%rsi
add    $0x1,%rsi
mov    %esi,%edx
and    $0x7,%edx
cmp    $0x38,%rdi
jae    1270 <mint_random_multinomial+0x140>
xor    %edi,%edi
jmp    1305 <mint_random_multinomial+0x1d5>
and    $0xfffffffffffffff8,%rsi
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
movups %xmm1,(%r15,%rdi,4)
movups %xmm1,0x10(%r15,%rdi,4)
movups %xmm1,0x20(%r15,%rdi,4)
movups %xmm1,0x30(%r15,%rdi,4)
movups %xmm1,0x40(%r15,%rdi,4)
movups %xmm1,0x50(%r15,%rdi,4)
movups %xmm1,0x60(%r15,%rdi,4)
movups %xmm1,0x70(%r15,%rdi,4)
movups %xmm1,0x80(%r15,%rdi,4)
movups %xmm1,0x90(%r15,%rdi,4)
movups %xmm1,0xa0(%r15,%rdi,4)
movups %xmm1,0xb0(%r15,%rdi,4)
movups %xmm1,0xc0(%r15,%rdi,4)
movups %xmm1,0xd0(%r15,%rdi,4)
movups %xmm1,0xe0(%r15,%rdi,4)
movups %xmm1,0xf0(%r15,%rdi,4)
add    $0x40,%rdi
add    $0xfffffffffffffff8,%rsi
jne    1280 <mint_random_multinomial+0x150>
test   %rdx,%rdx
je     1332 <mint_random_multinomial+0x202>
lea    (%r15,%rdi,4),%rsi
add    $0x10,%rsi
shl    $0x5,%rdx
xor    %edi,%edi
nopl   0x0(%rax,%rax,1)
movups %xmm1,-0x10(%rsi,%rdi,1)
movups %xmm1,(%rsi,%rdi,1)
add    $0x20,%rdi
cmp    %rdi,%rdx
jne    1320 <mint_random_multinomial+0x1f0>
cmp    %rax,%rcx
je     134f <mint_random_multinomial+0x21f>
nopw   0x0(%rax,%rax,1)
movss  %xmm0,(%r15,%rcx,4)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1340 <mint_random_multinomial+0x210>
test   %ebp,%ebp
jle    13c3 <mint_random_multinomial+0x293>
movss  0xca5(%rip),%xmm0
xor    %ebx,%ebx
movss  %xmm2,0xc(%rsp)
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movss  %xmm0,0x14(%rsp)
movss  (%r15,%rbx,4),%xmm3
movss  %xmm3,0x10(%rsp)
movaps %xmm2,%xmm1
mulss  %xmm0,%xmm1
movaps %xmm3,%xmm0
divss  %xmm1,%xmm0
mov    %ebp,%edi
call   1030 <mint_random_binomial@plt>
movss  0x14(%rsp),%xmm0
movss  0xc(%rsp),%xmm2
mov    %eax,(%r14,%rbx,4)
sub    %eax,%ebp
movss  0x10(%rsp),%xmm1
divss  %xmm2,%xmm1
subss  %xmm1,%xmm0
add    $0x1,%rbx
movslq %ebp,%rax
cmp    %rax,%rbx
jl     1370 <mint_random_multinomial+0x240>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @mint_random_multinomial(float* %0, i32* nocapture writeonly %1, i32 %2, i32 %3) {
%5 = bitcast float* %0 to i8*
%6 = shl i32 %2, 2
%7 = tail call float* @malloc(i32 %6)
%8 = bitcast float* %7 to i8*
%9 = icmp eq float* %0, null
%10 = icmp sgt i32 %2, 0
br i1 %9, label %20, label %11
11: ; preds = %4
br i1 %10, label %12, label %160
12: ; preds = %11
%13 = zext i32 %2 to i64
%14 = shl nuw nsw i64 %13, 2
%15 = add nsw i64 %13, -1
%16 = and i64 %13, 7
%17 = icmp ult i64 %15, 7
br i1 %17, label %145, label %18
18: ; preds = %12
%19 = and i64 %13, 4294967288
br label %107
20: ; preds = %4
br i1 %10, label %21, label %160
21: ; preds = %20
%22 = sitofp i32 %2 to double
%23 = fdiv double 1.000000e+00, %22
%24 = fptrunc double %23 to float
%25 = zext i32 %2 to i64
%26 = icmp ult i32 %2, 8
br i1 %26, label %100, label %27
27: ; preds = %21
%28 = and i64 %25, 4294967288
%29 = insertelement <4 x float> poison, float %24, i64 0
%30 = shufflevector <4 x float> %29, <4 x float> poison, <4 x i32> zeroinitializer
%31 = insertelement <4 x float> poison, float %24, i64 0
%32 = shufflevector <4 x float> %31, <4 x float> poison, <4 x i32> zeroinitializer
%33 = add nsw i64 %28, -8
%34 = lshr exact i64 %33, 3
%35 = add nuw nsw i64 %34, 1
%36 = and i64 %35, 7
%37 = icmp ult i64 %33, 56
br i1 %37, label %85, label %38
38: ; preds = %27
%39 = and i64 %35, 4611686018427387896
br label %40
40: ; preds = %40, %38
%41 = phi i64 [ 0, %38 ], [ %82, %40 ]
%42 = phi i64 [ 0, %38 ], [ %83, %40 ]
%43 = getelementptr inbounds float, float* %7, i64 %41
%44 = bitcast float* %43 to <4 x float>*
store <4 x float> %30, <4 x float>* %44
%45 = getelementptr inbounds float, float* %43, i64 4
%46 = bitcast float* %45 to <4 x float>*
store <4 x float> %32, <4 x float>* %46
%47 = or i64 %41, 8
%48 = getelementptr inbounds float, float* %7, i64 %47
%49 = bitcast float* %48 to <4 x float>*
store <4 x float> %30, <4 x float>* %49
%50 = getelementptr inbounds float, float* %48, i64 4
%51 = bitcast float* %50 to <4 x float>*
store <4 x float> %32, <4 x float>* %51
%52 = or i64 %41, 16
%53 = getelementptr inbounds float, float* %7, i64 %52
%54 = bitcast float* %53 to <4 x float>*
store <4 x float> %30, <4 x float>* %54
%55 = getelementptr inbounds float, float* %53, i64 4
%56 = bitcast float* %55 to <4 x float>*
store <4 x float> %32, <4 x float>* %56
%57 = or i64 %41, 24
%58 = getelementptr inbounds float, float* %7, i64 %57
%59 = bitcast float* %58 to <4 x float>*
store <4 x float> %30, <4 x float>* %59
%60 = getelementptr inbounds float, float* %58, i64 4
%61 = bitcast float* %60 to <4 x float>*
store <4 x float> %32, <4 x float>* %61
%62 = or i64 %41, 32
%63 = getelementptr inbounds float, float* %7, i64 %62
%64 = bitcast float* %63 to <4 x float>*
store <4 x float> %30, <4 x float>* %64
%65 = getelementptr inbounds float, float* %63, i64 4
%66 = bitcast float* %65 to <4 x float>*
store <4 x float> %32, <4 x float>* %66
%67 = or i64 %41, 40
%68 = getelementptr inbounds float, float* %7, i64 %67
%69 = bitcast float* %68 to <4 x float>*
store <4 x float> %30, <4 x float>* %69
%70 = getelementptr inbounds float, float* %68, i64 4
%71 = bitcast float* %70 to <4 x float>*
store <4 x float> %32, <4 x float>* %71
%72 = or i64 %41, 48
%73 = getelementptr inbounds float, float* %7, i64 %72
%74 = bitcast float* %73 to <4 x float>*
store <4 x float> %30, <4 x float>* %74
%75 = getelementptr inbounds float, float* %73, i64 4
%76 = bitcast float* %75 to <4 x float>*
store <4 x float> %32, <4 x float>* %76
%77 = or i64 %41, 56
%78 = getelementptr inbounds float, float* %7, i64 %77
%79 = bitcast float* %78 to <4 x float>*
store <4 x float> %30, <4 x float>* %79
%80 = getelementptr inbounds float, float* %78, i64 4
%81 = bitcast float* %80 to <4 x float>*
store <4 x float> %32, <4 x float>* %81
%82 = add nuw i64 %41, 64
%83 = add nuw i64 %42, 8
%84 = icmp eq i64 %83, %39
br i1 %84, label %85, label %40
85: ; preds = %40, %27
%86 = phi i64 [ 0, %27 ], [ %82, %40 ]
%87 = icmp eq i64 %36, 0
br i1 %87, label %98, label %88
88: ; preds = %85, %88
%89 = phi i64 [ %95, %88 ], [ %86, %85 ]
%90 = phi i64 [ %96, %88 ], [ 0, %85 ]
%91 = getelementptr inbounds float, float* %7, i64 %89
%92 = bitcast float* %91 to <4 x float>*
store <4 x float> %30, <4 x float>* %92
%93 = getelementptr inbounds float, float* %91, i64 4
%94 = bitcast float* %93 to <4 x float>*
store <4 x float> %32, <4 x float>* %94
%95 = add nuw i64 %89, 8
%96 = add i64 %90, 1
%97 = icmp eq i64 %96, %36
br i1 %97, label %98, label %88
98: ; preds = %88, %85
%99 = icmp eq i64 %28, %25
br i1 %99, label %160, label %100
100: ; preds = %21, %98
%101 = phi i64 [ 0, %21 ], [ %28, %98 ]
br label %102
102: ; preds = %100, %102
%103 = phi i64 [ %105, %102 ], [ %101, %100 ]
%104 = getelementptr inbounds float, float* %7, i64 %103
store float %24, float* %104
%105 = add nuw nsw i64 %103, 1
%106 = icmp eq i64 %105, %25
br i1 %106, label %160, label %102
107: ; preds = %107, %18
%108 = phi i64 [ 0, %18 ], [ %142, %107 ]
%109 = phi float [ 0.000000e+00, %18 ], [ %141, %107 ]
%110 = phi i64 [ 0, %18 ], [ %143, %107 ]
%111 = getelementptr inbounds float, float* %0, i64 %108
%112 = load float, float* %111
%113 = fadd float %109, %112
%114 = or i64 %108, 1
%115 = getelementptr inbounds float, float* %0, i64 %114
%116 = load float, float* %115
%117 = fadd float %113, %116
%118 = or i64 %108, 2
%119 = getelementptr inbounds float, float* %0, i64 %118
%120 = load float, float* %119
%121 = fadd float %117, %120
%122 = or i64 %108, 3
%123 = getelementptr inbounds float, float* %0, i64 %122
%124 = load float, float* %123
%125 = fadd float %121, %124
%126 = or i64 %108, 4
%127 = getelementptr inbounds float, float* %0, i64 %126
%128 = load float, float* %127
%129 = fadd float %125, %128
%130 = or i64 %108, 5
%131 = getelementptr inbounds float, float* %0, i64 %130
%132 = load float, float* %131
%133 = fadd float %129, %132
%134 = or i64 %108, 6
%135 = getelementptr inbounds float, float* %0, i64 %134
%136 = load float, float* %135
%137 = fadd float %133, %136
%138 = or i64 %108, 7
%139 = getelementptr inbounds float, float* %0, i64 %138
%140 = load float, float* %139
%141 = fadd float %137, %140
%142 = add nuw nsw i64 %108, 8
%143 = add i64 %110, 8
%144 = icmp eq i64 %143, %19
br i1 %144, label %145, label %107
145: ; preds = %107, %12
%146 = phi float [ undef, %12 ], [ %141, %107 ]
%147 = phi i64 [ 0, %12 ], [ %142, %107 ]
%148 = phi float [ 0.000000e+00, %12 ], [ %141, %107 ]
%149 = icmp eq i64 %16, 0
br i1 %149, label %160, label %150
150: ; preds = %145, %150
%151 = phi i64 [ %157, %150 ], [ %147, %145 ]
%152 = phi float [ %156, %150 ], [ %148, %145 ]
%153 = phi i64 [ %158, %150 ], [ 0, %145 ]
%154 = getelementptr inbounds float, float* %0, i64 %151
%155 = load float, float* %154
%156 = fadd float %152, %155
%157 = add nuw nsw i64 %151, 1
%158 = add i64 %153, 1
%159 = icmp eq i64 %158, %16
br i1 %159, label %160, label %150
160: ; preds = %145, %150, %102, %98, %11, %20
%161 = phi float [ 1.000000e+00, %20 ], [ 0.000000e+00, %11 ], [ 1.000000e+00, %98 ], [ 1.000000e+00, %102 ], [ %146, %145 ], [ %156, %150 ]
%162 = icmp sgt i32 %3, 0
br i1 %162, label %163, label %179
163: ; preds = %160, %163
%164 = phi i64 [ %176, %163 ], [ 0, %160 ]
%165 = phi float [ %175, %163 ], [ 1.000000e+00, %160 ]
%166 = phi i32 [ %173, %163 ], [ %3, %160 ]
%167 = getelementptr inbounds float, float* %7, i64 %164
%168 = load float, float* %167
%169 = fmul float %161, %165
%170 = fdiv float %168, %169
%171 = tail call i32 @mint_random_binomial(float %170, i32 %166)
%172 = getelementptr inbounds i32, i32* %1, i64 %164
store i32 %171, i32* %172
%173 = sub nsw i32 %166, %171
%174 = fdiv float %168, %161
%175 = fsub float %165, %174
%176 = add nuw nsw i64 %164, 1
%177 = sext i32 %173 to i64
%178 = icmp slt i64 %176, %177
br i1 %178, label %163, label %179
179: ; preds = %163, %160
ret void
}