#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 */
scalar_t__ RAND_MAX;
double random();
void simuHald(int gen, double mu, int *mut) {
  int wild = 1, newWild, newMut, i, j;
  *mut = 0;
  for (i = 1; i <= gen; i++) {
    newMut = newWild = 0;
    for (j = 1; j <= wild; j++) {
      if (random() / ((double)RAND_MAX + 1) < mu)
        newMut++;
      else
        newWild++;
    }
    mut[0] = 2 * (*mut) + newMut;
    wild = wild + newWild;
  }
}


x86 assembly (clang, -O0):
<simuHald>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %edi,-0x4(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
movl   $0x1,-0x1c(%rbp)
mov    -0x18(%rbp),%rax
movl   $0x0,(%rax)
movl   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     11f1 <simuHald+0xe1>
movl   $0x0,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     11c8 <simuHald+0xb8>
mov    $0x0,%al
call   1030 <random@plt>
movaps %xmm0,%xmm1
mov    0x2e6e(%rip),%rax
cvtsi2sdq (%rax),%xmm0
movsd  0xe79(%rip),%xmm2
addsd  %xmm2,%xmm0
divsd  %xmm0,%xmm1
movsd  -0x10(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    11ac <simuHald+0x9c>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11b5 <simuHald+0xa5>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    11ba <simuHald+0xaa>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    115d <simuHald+0x4d>
mov    -0x18(%rbp),%rax
mov    (%rax),%ecx
shl    $0x1,%ecx
add    -0x24(%rbp),%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x1c(%rbp),%eax
add    -0x20(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    113c <simuHald+0x2c>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<simuHald>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %rsi,0x8(%rsp)
movl   $0x0,(%rsi)
mov    %edi,0x4(%rsp)
test   %edi,%edi
jle    11cd <simuHald+0xbd>
mov    $0x1,%r12d
mov    0x2ea6(%rip),%rbx
mov    $0x1,%r13d
movsd  %xmm0,0x10(%rsp)
xchg   %ax,%ax
cmp    $0x2,%r13d
mov    $0x1,%r15d
cmovae %r13d,%r15d
xor    %r14d,%r14d
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
xor    %eax,%eax
call   1030 <random@plt>
movsd  0xe81(%rip),%xmm2
xorps  %xmm1,%xmm1
cvtsi2sdq (%rbx),%xmm1
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
xor    %eax,%eax
xor    %ecx,%ecx
movsd  0x10(%rsp),%xmm1
ucomisd %xmm0,%xmm1
seta   %al
setbe  %cl
add    %ecx,%ebp
add    %eax,%r14d
add    $0xffffffff,%r15d
jne    1170 <simuHald+0x60>
mov    0x8(%rsp),%rcx
mov    (%rcx),%eax
lea    (%r14,%rax,2),%eax
mov    %eax,(%rcx)
add    %ebp,%r13d
lea    0x1(%r12),%eax
cmp    0x4(%rsp),%r12d
mov    %eax,%r12d
jne    1150 <simuHald+0x40>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@RAND_MAX = global i64 0
define void @simuHald(i32 %0, double %1, i32* %2) {
%4 = alloca i32
%5 = alloca double
%6 = alloca i32*
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
store i32 %0, i32* %4
store double %1, double* %5
store i32* %2, i32** %6
store i32 1, i32* %7
%12 = load i32*, i32** %6
store i32 0, i32* %12
store i32 1, i32* %10
br label %13
13: ; preds = %51, %3
%14 = load i32, i32* %10
%15 = load i32, i32* %4
%16 = icmp sle i32 %14, %15
br i1 %16, label %17, label %54
17: ; preds = %13
store i32 0, i32* %8
store i32 0, i32* %9
store i32 1, i32* %11
br label %18
18: ; preds = %37, %17
%19 = load i32, i32* %11
%20 = load i32, i32* %7
%21 = icmp sle i32 %19, %20
br i1 %21, label %22, label %40
22: ; preds = %18
%23 = call double (...) @random()
%24 = load i64, i64* @RAND_MAX
%25 = sitofp i64 %24 to double
%26 = fadd double %25, 1.000000e+00
%27 = fdiv double %23, %26
%28 = load double, double* %5
%29 = fcmp olt double %27, %28
br i1 %29, label %30, label %33
30: ; preds = %22
%31 = load i32, i32* %9
%32 = add nsw i32 %31, 1
store i32 %32, i32* %9
br label %36
33: ; preds = %22
%34 = load i32, i32* %8
%35 = add nsw i32 %34, 1
store i32 %35, i32* %8
br label %36
36: ; preds = %33, %30
br label %37
37: ; preds = %36
%38 = load i32, i32* %11
%39 = add nsw i32 %38, 1
store i32 %39, i32* %11
br label %18
40: ; preds = %18
%41 = load i32*, i32** %6
%42 = load i32, i32* %41
%43 = mul nsw i32 2, %42
%44 = load i32, i32* %9
%45 = add nsw i32 %43, %44
%46 = load i32*, i32** %6
%47 = getelementptr inbounds i32, i32* %46, i64 0
store i32 %45, i32* %47
%48 = load i32, i32* %7
%49 = load i32, i32* %8
%50 = add nsw i32 %48, %49
store i32 %50, i32* %7
br label %51
51: ; preds = %40
%52 = load i32, i32* %10
%53 = add nsw i32 %52, 1
store i32 %53, i32* %10
br label %13
54: ; preds = %13
ret void
}

LLVM IR (-O2):
@RAND_MAX = global i64 0
define void @simuHald(i32 %0, double %1, i32* nocapture %2) {
store i32 0, i32* %2
%4 = icmp slt i32 %0, 1
br i1 %4, label %33, label %5
5: ; preds = %3, %26
%6 = phi i32 [ %31, %26 ], [ 1, %3 ]
%7 = phi i32 [ %30, %26 ], [ 1, %3 ]
%8 = call i32 @llvm.umax.i32(i32 %7, i32 1)
br label %9
9: ; preds = %5, %9
%10 = phi i32 [ %24, %9 ], [ 1, %5 ]
%11 = phi i32 [ %23, %9 ], [ 0, %5 ]
%12 = phi i32 [ %21, %9 ], [ 0, %5 ]
%13 = tail call double (...) @random()
%14 = load i64, i64* @RAND_MAX
%15 = sitofp i64 %14 to double
%16 = fadd double %15, 1.000000e+00
%17 = fdiv double %13, %16
%18 = fcmp olt double %17, %1
%19 = xor i1 %18, true
%20 = zext i1 %19 to i32
%21 = add nuw nsw i32 %12, %20
%22 = zext i1 %18 to i32
%23 = add nuw nsw i32 %11, %22
%24 = add nuw i32 %10, 1
%25 = icmp eq i32 %10, %8
br i1 %25, label %26, label %9
26: ; preds = %9
%27 = load i32, i32* %2
%28 = shl nsw i32 %27, 1
%29 = add nsw i32 %28, %23
store i32 %29, i32* %2
%30 = add nuw nsw i32 %21, %7
%31 = add nuw i32 %6, 1
%32 = icmp eq i32 %6, %0
br i1 %32, label %33, label %5
33: ; preds = %26, %3
ret void
}