#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__ PairwiseIdentity(char *, char *);
int sre_random();
float AlignmentIdentityBySampling(char **aseq, int L, int N, int nsample) {
  int x, i, j;
  float sum;
  if (N < 2)
    return 1.0;
  sum = 0.;
  for (x = 0; x < nsample; x++) {
    i = ((int)(sre_random() * (N)));
    do {
      j = ((int)(sre_random() * (N)));
    } while (j == i);
    sum += PairwiseIdentity(aseq[i], aseq[j]);
  }
  return sum / (float)nsample;
}


x86 assembly (clang, -O0):
<AlignmentIdentityBySampling>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %edx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
cmpl   $0x2,-0x18(%rbp)
jge    1151 <AlignmentIdentityBySampling+0x31>
movss  0xeb9(%rip),%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    11e1 <AlignmentIdentityBySampling+0xc1>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x2c(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11ce <AlignmentIdentityBySampling+0xae>
mov    $0x0,%al
call   1040 <sre_random@plt>
imul   -0x18(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    $0x0,%al
call   1040 <sre_random@plt>
imul   -0x18(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
je     117a <AlignmentIdentityBySampling+0x5a>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rdi
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rsi
call   1030 <PairwiseIdentity@plt>
cvtsi2ss %rax,%xmm0
addss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1160 <AlignmentIdentityBySampling+0x40>
movss  -0x2c(%rbp),%xmm0
cvtsi2ssl -0x1c(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<AlignmentIdentityBySampling>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x10,%rsp
cmp    $0x2,%edx
jge    113e <AlignmentIdentityBySampling+0x1e>
movss  0xec7(%rip),%xmm0
jmp    11bc <AlignmentIdentityBySampling+0x9c>
mov    %ecx,%r14d
test   %ecx,%ecx
jle    11ad <AlignmentIdentityBySampling+0x8d>
mov    %edx,%ebp
mov    %rdi,%r15
xorps  %xmm0,%xmm0
xor    %r12d,%r12d
movss  %xmm0,0xc(%rsp)
xor    %eax,%eax
call   1040 <sre_random@plt>
mov    %eax,%ebx
imul   %ebp,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
xor    %eax,%eax
call   1040 <sre_random@plt>
imul   %ebp,%eax
cmp    %ebx,%eax
je     1170 <AlignmentIdentityBySampling+0x50>
movslq %ebx,%rcx
mov    (%r15,%rcx,8),%rdi
cltq
mov    (%r15,%rax,8),%rsi
call   1030 <PairwiseIdentity@plt>
xorps  %xmm1,%xmm1
cvtsi2ss %rax,%xmm1
movss  0xc(%rsp),%xmm0
addss  %xmm1,%xmm0
add    $0x1,%r12d
cmp    %r14d,%r12d
jne    1150 <AlignmentIdentityBySampling+0x30>
jmp    11b0 <AlignmentIdentityBySampling+0x90>
xorps  %xmm0,%xmm0
xorps  %xmm1,%xmm1
cvtsi2ss %r14d,%xmm1
divss  %xmm1,%xmm0
add    $0x10,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define float @AlignmentIdentityBySampling(i8** nocapture %0, i32 %1, i32 %2, i32 %3) {
%5 = icmp slt i32 %2, 2
br i1 %5, label %33, label %6
6: ; preds = %4
%7 = icmp sgt i32 %3, 0
br i1 %7, label %8, label %29
8: ; preds = %6, %17
%9 = phi float [ %26, %17 ], [ 0.000000e+00, %6 ]
%10 = phi i32 [ %27, %17 ], [ 0, %6 ]
%11 = tail call i32 (...) @sre_random()
%12 = mul nsw i32 %11, %2
br label %13
13: ; preds = %13, %8
%14 = tail call i32 (...) @sre_random()
%15 = mul nsw i32 %14, %2
%16 = icmp eq i32 %15, %12
br i1 %16, label %13, label %17
17: ; preds = %13
%18 = sext i32 %12 to i64
%19 = getelementptr inbounds i8*, i8** %0, i64 %18
%20 = load i8*, i8** %19
%21 = sext i32 %15 to i64
%22 = getelementptr inbounds i8*, i8** %0, i64 %21
%23 = load i8*, i8** %22
%24 = tail call i64 @PairwiseIdentity(i8* %20, i8* %23)
%25 = sitofp i64 %24 to float
%26 = fadd float %9, %25
%27 = add nuw nsw i32 %10, 1
%28 = icmp eq i32 %27, %3
br i1 %28, label %29, label %8
29: ; preds = %17, %6
%30 = phi float [ 0.000000e+00, %6 ], [ %26, %17 ]
%31 = sitofp i32 %3 to float
%32 = fdiv float %30, %31
br label %33
33: ; preds = %4, %29
%34 = phi float [ %32, %29 ], [ 1.000000e+00, %4 ]
ret float %34
}