#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 */
void *calloc(int, int);
double drand();
int /*<<< orphan*/ free(double *);
double rand_double(int /*<<< orphan*/, int);
int /*<<< orphan*/ rowMeans(double *, int, int, double *);
int /*<<< orphan*/ strcmp(char *, char *);
int DataImputation(double *G, int **missing, int nrow, int ncol, char *method) {
  int i, j, geno, na = 0, count = 0;
  double r;
  double *AF = calloc(nrow, sizeof(double));
  for (i = 0; i < ncol * nrow; i++)
    if (G[i] == 9)
      na++;
  *missing = calloc(2 * na, sizeof(int));
  if (!strcmp(method, "naive")) {
    rowMeans(G, nrow, ncol, AF);
    for (i = 0; i < nrow; i++) {
      for (j = 0; j < ncol; j++) {
        geno = 0;
        if (G[i * ncol + j] == 9) {
          r = drand();
          if (r < AF[i] / 2)
            geno++;
          r = rand_double(0, 1);
          if (r < AF[i] / 2)
            geno++;
          G[i * ncol + j] = geno;
          (*missing)[count * 2] = i;
          (*missing)[count * 2 + 1] = j;
          count++;
        }
      }
    }
  }
  free(AF);
  return na;
}


x86 assembly (clang, -O0):
<DataImputation>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
movl   $0x0,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x14(%rbp),%edi
mov    $0x8,%esi
call   1060 <calloc@plt>
mov    %rax,-0x48(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x18(%rbp),%ecx
imul   -0x14(%rbp),%ecx
cmp    %ecx,%eax
jge    11f3 <DataImputation+0x93>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xe39(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    11e0 <DataImputation+0x80>
jp     11e0 <DataImputation+0x80>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11e5 <DataImputation+0x85>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11a0 <DataImputation+0x40>
mov    -0x30(%rbp),%edi
shl    $0x1,%edi
mov    $0x4,%esi
call   1060 <calloc@plt>
mov    %rax,%rcx
mov    -0x10(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x20(%rbp),%rdi
lea    0xdf8(%rip),%rsi
call   1070 <strcmp@plt>
cmp    $0x0,%eax
jne    1391 <DataImputation+0x231>
mov    -0x8(%rbp),%rdi
mov    -0x14(%rbp),%esi
mov    -0x18(%rbp),%edx
mov    -0x48(%rbp),%rcx
call   1080 <rowMeans@plt>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    138c <DataImputation+0x22c>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    1379 <DataImputation+0x219>
movl   $0x0,-0x2c(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xd7c(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    1366 <DataImputation+0x206>
jp     1366 <DataImputation+0x206>
mov    $0x0,%al
call   1040 <drand@plt>
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm1
mov    -0x48(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xd4e(%rip),%xmm2
divsd  %xmm2,%xmm0
ucomisd %xmm1,%xmm0
jbe    12d1 <DataImputation+0x171>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
xor    %edi,%edi
mov    $0x1,%esi
call   1050 <rand_double@plt>
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm1
mov    -0x48(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xd0c(%rip),%xmm2
divsd  %xmm2,%xmm0
ucomisd %xmm1,%xmm0
jbe    1313 <DataImputation+0x1b3>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
cvtsi2sdl -0x2c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%edx
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
mov    -0x34(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x28(%rbp),%edx
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
mov    -0x34(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    136b <DataImputation+0x20b>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1253 <DataImputation+0xf3>
jmp    137e <DataImputation+0x21e>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1240 <DataImputation+0xe0>
jmp    1391 <DataImputation+0x231>
mov    -0x48(%rbp),%rdi
call   1030 <free@plt>
mov    -0x30(%rbp),%eax
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<DataImputation>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %r8,%r15
mov    %ecx,%ebx
mov    %edx,%r14d
mov    %rsi,0x10(%rsp)
mov    %rdi,%rbp
mov    %edx,%edi
mov    $0x8,%esi
call   1060 <calloc@plt>
mov    %rax,0x8(%rsp)
mov    %ebx,0x4(%rsp)
imul   %r14d,%ebx
test   %ebx,%ebx
jle    11ab <DataImputation+0x4b>
mov    %ebx,%eax
cmp    $0x4,%ebx
jae    11b2 <DataImputation+0x52>
xor    %ecx,%ecx
xor    %esi,%esi
jmp    128f <DataImputation+0x12f>
xor    %esi,%esi
jmp    12bb <DataImputation+0x15b>
mov    %eax,%ecx
and    $0xfffffffc,%ecx
lea    -0x4(%rcx),%rsi
mov    %rsi,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
test   %rsi,%rsi
je     1418 <DataImputation+0x2b8>
mov    %rdx,%rdi
and    $0xfffffffffffffffe,%rdi
pxor   %xmm0,%xmm0
xor    %esi,%esi
movapd 0xe1c(%rip),%xmm2
pxor   %xmm1,%xmm1
nopl   0x0(%rax,%rax,1)
movupd 0x0(%rbp,%rsi,8),%xmm3
movupd 0x10(%rbp,%rsi,8),%xmm4
movupd 0x20(%rbp,%rsi,8),%xmm5
movupd 0x30(%rbp,%rsi,8),%xmm6
cmpeqpd %xmm2,%xmm3
shufps $0xe8,%xmm3,%xmm3
psubd  %xmm3,%xmm0
cmpeqpd %xmm2,%xmm4
shufps $0xe8,%xmm4,%xmm4
psubd  %xmm4,%xmm1
cmpeqpd %xmm2,%xmm5
shufps $0xe8,%xmm5,%xmm5
psubd  %xmm5,%xmm0
cmpeqpd %xmm2,%xmm6
shufps $0xe8,%xmm6,%xmm6
psubd  %xmm6,%xmm1
add    $0x8,%rsi
add    $0xfffffffffffffffe,%rdi
jne    11f0 <DataImputation+0x90>
test   $0x1,%dl
je     1279 <DataImputation+0x119>
movupd 0x0(%rbp,%rsi,8),%xmm2
movupd 0x10(%rbp,%rsi,8),%xmm3
movapd 0xda1(%rip),%xmm4
cmpeqpd %xmm4,%xmm2
shufps $0xe8,%xmm2,%xmm2
psubd  %xmm2,%xmm0
cmpeqpd %xmm4,%xmm3
shufps $0xe8,%xmm3,%xmm3
psubd  %xmm3,%xmm1
paddd  %xmm1,%xmm0
pshufd $0x55,%xmm0,%xmm1
paddd  %xmm0,%xmm1
movd   %xmm1,%esi
cmp    %rax,%rcx
je     12bb <DataImputation+0x15b>
movsd  0xd79(%rip),%xmm0
nopw   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rcx,8),%xmm1
cmpeqsd %xmm0,%xmm1
movq   %xmm1,%rdx
sub    %edx,%esi
add    $0x1,%rcx
cmp    %rcx,%rax
jne    12a0 <DataImputation+0x140>
mov    %rsi,0x18(%rsp)
lea    (%rsi,%rsi,1),%edi
mov    $0x4,%esi
call   1060 <calloc@plt>
mov    0x10(%rsp),%rcx
mov    %rax,(%rcx)
lea    0xd44(%rip),%rsi
mov    %r15,%rdi
call   1070 <strcmp@plt>
test   %eax,%eax
je     1306 <DataImputation+0x1a6>
mov    0x8(%rsp),%rdi
call   1030 <free@plt>
mov    0x18(%rsp),%rax
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %rbp,%rdi
mov    %r14d,%esi
mov    0x4(%rsp),%edx
mov    0x8(%rsp),%rcx
call   1080 <rowMeans@plt>
test   %r14d,%r14d
jle    12e8 <DataImputation+0x188>
mov    0x4(%rsp),%eax
movslq %eax,%rcx
mov    %r14d,%edx
mov    %rdx,0x20(%rsp)
mov    %eax,%eax
mov    %rax,0x30(%rsp)
shl    $0x3,%rcx
mov    %rcx,0x28(%rsp)
xor    %r12d,%r12d
xor    %r15d,%r15d
jmp    1360 <DataImputation+0x200>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%r12
add    0x28(%rsp),%rbp
cmp    0x20(%rsp),%r12
je     12e8 <DataImputation+0x188>
cmpl   $0x0,0x4(%rsp)
jle    1350 <DataImputation+0x1f0>
xor    %r13d,%r13d
jmp    137b <DataImputation+0x21b>
nopl   0x0(%rax)
add    $0x1,%r13
cmp    %r13,0x30(%rsp)
je     1350 <DataImputation+0x1f0>
movsd  0x0(%rbp,%r13,8),%xmm0
ucomisd 0xc86(%rip),%xmm0
jne    1370 <DataImputation+0x210>
jp     1370 <DataImputation+0x210>
xor    %eax,%eax
call   1040 <drand@plt>
mov    %rbp,%r14
mov    0x8(%rsp),%rbp
movsd  0x0(%rbp,%r12,8),%xmm1
movsd  0xc6c(%rip),%xmm2
mulsd  %xmm2,%xmm1
xor    %ebx,%ebx
ucomisd %xmm0,%xmm1
seta   %bl
xor    %edi,%edi
mov    $0x1,%esi
call   1050 <rand_double@plt>
movsd  0x0(%rbp,%r12,8),%xmm1
mov    %r14,%rbp
mulsd  0xc41(%rip),%xmm1
ucomisd %xmm0,%xmm1
lea    0x1(%rbx),%eax
cmovbe %ebx,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%r14,%r13,8)
mov    0x10(%rsp),%rax
mov    (%rax),%rax
lea    (%r15,%r15,1),%ecx
movslq %ecx,%rcx
mov    %r12d,(%rax,%rcx,4)
lea    (%r15,%r15,1),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    %r13d,(%rax,%rcx,4)
add    $0x1,%r15d
jmp    1370 <DataImputation+0x210>
pxor   %xmm0,%xmm0
xor    %esi,%esi
pxor   %xmm1,%xmm1
test   $0x1,%dl
jne    124b <DataImputation+0xeb>
jmp    1279 <DataImputation+0x119>

LLVM IR (-O0):
@.str = [6 x i8] c"naive\00"
define i32 @DataImputation(double* %0, i32** %1, i32 %2, i32 %3, i8* %4) {
%6 = alloca double*
%7 = alloca i32**
%8 = alloca i32
%9 = alloca i32
%10 = alloca i8*
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca double
%17 = alloca double*
store double* %0, double** %6
store i32** %1, i32*** %7
store i32 %2, i32* %8
store i32 %3, i32* %9
store i8* %4, i8** %10
store i32 0, i32* %14
store i32 0, i32* %15
%18 = load i32, i32* %8
%19 = call i8* @calloc(i32 %18, i32 8)
%20 = bitcast i8* %19 to double*
store double* %20, double** %17
store i32 0, i32* %11
br label %21
21: ; preds = %38, %5
%22 = load i32, i32* %11
%23 = load i32, i32* %9
%24 = load i32, i32* %8
%25 = mul nsw i32 %23, %24
%26 = icmp slt i32 %22, %25
br i1 %26, label %27, label %41
27: ; preds = %21
%28 = load double*, double** %6
%29 = load i32, i32* %11
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds double, double* %28, i64 %30
%32 = load double, double* %31
%33 = fcmp oeq double %32, 9.000000e+00
br i1 %33, label %34, label %37
34: ; preds = %27
%35 = load i32, i32* %14
%36 = add nsw i32 %35, 1
store i32 %36, i32* %14
br label %37
37: ; preds = %34, %27
br label %38
38: ; preds = %37
%39 = load i32, i32* %11
%40 = add nsw i32 %39, 1
store i32 %40, i32* %11
br label %21
41: ; preds = %21
%42 = load i32, i32* %14
%43 = mul nsw i32 2, %42
%44 = call i8* @calloc(i32 %43, i32 4)
%45 = bitcast i8* %44 to i32*
%46 = load i32**, i32*** %7
store i32* %45, i32** %46
%47 = load i8*, i8** %10
%48 = call i32 @strcmp(i8* %47, i8* getelementptr inbounds ([6 x i8], [6 x i8]* @.str, i64 0, i64 0))
%49 = icmp ne i32 %48, 0
br i1 %49, label %139, label %50
50: ; preds = %41
%51 = load double*, double** %6
%52 = load i32, i32* %8
%53 = load i32, i32* %9
%54 = load double*, double** %17
%55 = call i32 @rowMeans(double* %51, i32 %52, i32 %53, double* %54)
store i32 0, i32* %11
br label %56
56: ; preds = %135, %50
%57 = load i32, i32* %11
%58 = load i32, i32* %8
%59 = icmp slt i32 %57, %58
br i1 %59, label %60, label %138
60: ; preds = %56
store i32 0, i32* %12
br label %61
61: ; preds = %131, %60
%62 = load i32, i32* %12
%63 = load i32, i32* %9
%64 = icmp slt i32 %62, %63
br i1 %64, label %65, label %134
65: ; preds = %61
store i32 0, i32* %13
%66 = load double*, double** %6
%67 = load i32, i32* %11
%68 = load i32, i32* %9
%69 = mul nsw i32 %67, %68
%70 = load i32, i32* %12
%71 = add nsw i32 %69, %70
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds double, double* %66, i64 %72
%74 = load double, double* %73
%75 = fcmp oeq double %74, 9.000000e+00
br i1 %75, label %76, label %130
76: ; preds = %65
%77 = call double (...) @drand()
store double %77, double* %16
%78 = load double, double* %16
%79 = load double*, double** %17
%80 = load i32, i32* %11
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds double, double* %79, i64 %81
%83 = load double, double* %82
%84 = fdiv double %83, 2.000000e+00
%85 = fcmp olt double %78, %84
br i1 %85, label %86, label %89
86: ; preds = %76
%87 = load i32, i32* %13
%88 = add nsw i32 %87, 1
store i32 %88, i32* %13
br label %89
89: ; preds = %86, %76
%90 = call double @rand_double(i32 0, i32 1)
store double %90, double* %16
%91 = load double, double* %16
%92 = load double*, double** %17
%93 = load i32, i32* %11
%94 = sext i32 %93 to i64
%95 = getelementptr inbounds double, double* %92, i64 %94
%96 = load double, double* %95
%97 = fdiv double %96, 2.000000e+00
%98 = fcmp olt double %91, %97
br i1 %98, label %99, label %102
99: ; preds = %89
%100 = load i32, i32* %13
%101 = add nsw i32 %100, 1
store i32 %101, i32* %13
br label %102
102: ; preds = %99, %89
%103 = load i32, i32* %13
%104 = sitofp i32 %103 to double
%105 = load double*, double** %6
%106 = load i32, i32* %11
%107 = load i32, i32* %9
%108 = mul nsw i32 %106, %107
%109 = load i32, i32* %12
%110 = add nsw i32 %108, %109
%111 = sext i32 %110 to i64
%112 = getelementptr inbounds double, double* %105, i64 %111
store double %104, double* %112
%113 = load i32, i32* %11
%114 = load i32**, i32*** %7
%115 = load i32*, i32** %114
%116 = load i32, i32* %15
%117 = mul nsw i32 %116, 2
%118 = sext i32 %117 to i64
%119 = getelementptr inbounds i32, i32* %115, i64 %118
store i32 %113, i32* %119
%120 = load i32, i32* %12
%121 = load i32**, i32*** %7
%122 = load i32*, i32** %121
%123 = load i32, i32* %15
%124 = mul nsw i32 %123, 2
%125 = add nsw i32 %124, 1
%126 = sext i32 %125 to i64
%127 = getelementptr inbounds i32, i32* %122, i64 %126
store i32 %120, i32* %127
%128 = load i32, i32* %15
%129 = add nsw i32 %128, 1
store i32 %129, i32* %15
br label %130
130: ; preds = %102, %65
br label %131
131: ; preds = %130
%132 = load i32, i32* %12
%133 = add nsw i32 %132, 1
store i32 %133, i32* %12
br label %61
134: ; preds = %61
br label %135
135: ; preds = %134
%136 = load i32, i32* %11
%137 = add nsw i32 %136, 1
store i32 %137, i32* %11
br label %56
138: ; preds = %56
br label %139
139: ; preds = %138, %41
%140 = load double*, double** %17
%141 = call i32 @free(double* %140)
%142 = load i32, i32* %14
ret i32 %142
}

LLVM IR (-O2):
@.str = [6 x i8] c"naive\00"
define i32 @DataImputation(double* %0, i32** nocapture %1, i32 %2, i32 %3, i8* nocapture %4) {
%6 = tail call i8* @calloc(i32 %2, i32 8)
%7 = bitcast i8* %6 to double*
%8 = mul i32 %3, %2
%9 = icmp sgt i32 %8, 0
br i1 %9, label %10, label %94
10: ; preds = %5
%11 = zext i32 %8 to i64
%12 = icmp ult i32 %8, 4
br i1 %12, label %81, label %13
13: ; preds = %10
%14 = and i64 %11, 4294967292
%15 = add nsw i64 %14, -4
%16 = lshr exact i64 %15, 2
%17 = add nuw nsw i64 %16, 1
%18 = and i64 %17, 1
%19 = icmp eq i64 %15, 0
br i1 %19, label %55, label %20
20: ; preds = %13
%21 = and i64 %17, 9223372036854775806
br label %22
22: ; preds = %22, %20
%23 = phi i64 [ 0, %20 ], [ %52, %22 ]
%24 = phi <2 x i32> [ zeroinitializer, %20 ], [ %50, %22 ]
%25 = phi <2 x i32> [ zeroinitializer, %20 ], [ %51, %22 ]
%26 = phi i64 [ 0, %20 ], [ %53, %22 ]
%27 = getelementptr inbounds double, double* %0, i64 %23
%28 = bitcast double* %27 to <2 x double>*
%29 = load <2 x double>, <2 x double>* %28
%30 = getelementptr inbounds double, double* %27, i64 2
%31 = bitcast double* %30 to <2 x double>*
%32 = load <2 x double>, <2 x double>* %31
%33 = fcmp oeq <2 x double> %29, <double 9.000000e+00, double 9.000000e+00>
%34 = fcmp oeq <2 x double> %32, <double 9.000000e+00, double 9.000000e+00>
%35 = zext <2 x i1> %33 to <2 x i32>
%36 = zext <2 x i1> %34 to <2 x i32>
%37 = add <2 x i32> %24, %35
%38 = add <2 x i32> %25, %36
%39 = or i64 %23, 4
%40 = getelementptr inbounds double, double* %0, i64 %39
%41 = bitcast double* %40 to <2 x double>*
%42 = load <2 x double>, <2 x double>* %41
%43 = getelementptr inbounds double, double* %40, i64 2
%44 = bitcast double* %43 to <2 x double>*
%45 = load <2 x double>, <2 x double>* %44
%46 = fcmp oeq <2 x double> %42, <double 9.000000e+00, double 9.000000e+00>
%47 = fcmp oeq <2 x double> %45, <double 9.000000e+00, double 9.000000e+00>
%48 = zext <2 x i1> %46 to <2 x i32>
%49 = zext <2 x i1> %47 to <2 x i32>
%50 = add <2 x i32> %37, %48
%51 = add <2 x i32> %38, %49
%52 = add nuw i64 %23, 8
%53 = add i64 %26, 2
%54 = icmp eq i64 %53, %21
br i1 %54, label %55, label %22
55: ; preds = %22, %13
%56 = phi <2 x i32> [ undef, %13 ], [ %50, %22 ]
%57 = phi <2 x i32> [ undef, %13 ], [ %51, %22 ]
%58 = phi i64 [ 0, %13 ], [ %52, %22 ]
%59 = phi <2 x i32> [ zeroinitializer, %13 ], [ %50, %22 ]
%60 = phi <2 x i32> [ zeroinitializer, %13 ], [ %51, %22 ]
%61 = icmp eq i64 %18, 0
br i1 %61, label %75, label %62
62: ; preds = %55
%63 = getelementptr inbounds double, double* %0, i64 %58
%64 = bitcast double* %63 to <2 x double>*
%65 = load <2 x double>, <2 x double>* %64
%66 = getelementptr inbounds double, double* %63, i64 2
%67 = bitcast double* %66 to <2 x double>*
%68 = load <2 x double>, <2 x double>* %67
%69 = fcmp oeq <2 x double> %65, <double 9.000000e+00, double 9.000000e+00>
%70 = fcmp oeq <2 x double> %68, <double 9.000000e+00, double 9.000000e+00>
%71 = zext <2 x i1> %69 to <2 x i32>
%72 = zext <2 x i1> %70 to <2 x i32>
%73 = add <2 x i32> %59, %71
%74 = add <2 x i32> %60, %72
br label %75
75: ; preds = %55, %62
%76 = phi <2 x i32> [ %56, %55 ], [ %73, %62 ]
%77 = phi <2 x i32> [ %57, %55 ], [ %74, %62 ]
%78 = add <2 x i32> %77, %76
%79 = call i32 @llvm.vector.reduce.add.v2i32(<2 x i32> %78)
%80 = icmp eq i64 %14, %11
br i1 %80, label %94, label %81
81: ; preds = %10, %75
%82 = phi i64 [ 0, %10 ], [ %14, %75 ]
%83 = phi i32 [ 0, %10 ], [ %79, %75 ]
br label %84
84: ; preds = %81, %84
%85 = phi i64 [ %92, %84 ], [ %82, %81 ]
%86 = phi i32 [ %91, %84 ], [ %83, %81 ]
%87 = getelementptr inbounds double, double* %0, i64 %85
%88 = load double, double* %87
%89 = fcmp oeq double %88, 9.000000e+00
%90 = zext i1 %89 to i32
%91 = add nuw nsw i32 %86, %90
%92 = add nuw nsw i64 %85, 1
%93 = icmp eq i64 %92, %11
br i1 %93, label %94, label %84
94: ; preds = %84, %75, %5
%95 = phi i32 [ 0, %5 ], [ %79, %75 ], [ %91, %84 ]
%96 = shl nuw nsw i32 %95, 1
%97 = tail call i8* @calloc(i32 %96, i32 4)
%98 = bitcast i32** %1 to i8**
store i8* %97, i8** %98
%99 = tail call i32 @strcmp(i8* dereferenceable(1) %4, i8* dereferenceable(6) getelementptr inbounds ([6 x i8], [6 x i8]* @.str, i64 0, i64 0))
%100 = icmp eq i32 %99, 0
br i1 %100, label %101, label %153
101: ; preds = %94
%102 = tail call i32 @rowMeans(double* %0, i32 %2, i32 %3, double* %7)
%103 = icmp sgt i32 %2, 0
br i1 %103, label %104, label %153
104: ; preds = %101
%105 = icmp sgt i32 %3, 0
%106 = sext i32 %3 to i64
%107 = zext i32 %2 to i64
%108 = zext i32 %3 to i64
br label %109
109: ; preds = %104, %149
%110 = phi i64 [ 0, %104 ], [ %151, %149 ]
%111 = phi i32 [ 0, %104 ], [ %150, %149 ]
br i1 %105, label %112, label %149
112: ; preds = %109
%113 = mul nsw i64 %110, %106
%114 = getelementptr inbounds double, double* %7, i64 %110
%115 = trunc i64 %110 to i32
br label %116
116: ; preds = %112, %145
%117 = phi i64 [ 0, %112 ], [ %147, %145 ]
%118 = phi i32 [ %111, %112 ], [ %146, %145 ]
%119 = add nsw i64 %117, %113
%120 = getelementptr inbounds double, double* %0, i64 %119
%121 = load double, double* %120
%122 = fcmp oeq double %121, 9.000000e+00
br i1 %122, label %123, label %145
123: ; preds = %116
%124 = tail call double (...) @drand()
%125 = load double, double* %114
%126 = fmul double %125, 5.000000e-01
%127 = fcmp olt double %124, %126
%128 = zext i1 %127 to i32
%129 = tail call double @rand_double(i32 0, i32 1)
%130 = load double, double* %114
%131 = fmul double %130, 5.000000e-01
%132 = fcmp olt double %129, %131
%133 = select i1 %127, i32 2, i32 1
%134 = select i1 %132, i32 %133, i32 %128
%135 = sitofp i32 %134 to double
store double %135, double* %120
%136 = load i32*, i32** %1
%137 = shl nsw i32 %118, 1
%138 = sext i32 %137 to i64
%139 = getelementptr inbounds i32, i32* %136, i64 %138
store i32 %115, i32* %139
%140 = or i32 %137, 1
%141 = sext i32 %140 to i64
%142 = getelementptr inbounds i32, i32* %136, i64 %141
%143 = trunc i64 %117 to i32
store i32 %143, i32* %142
%144 = add nsw i32 %118, 1
br label %145
145: ; preds = %116, %123
%146 = phi i32 [ %144, %123 ], [ %118, %116 ]
%147 = add nuw nsw i64 %117, 1
%148 = icmp eq i64 %147, %108
br i1 %148, label %149, label %116
149: ; preds = %145, %109
%150 = phi i32 [ %111, %109 ], [ %146, %145 ]
%151 = add nuw nsw i64 %110, 1
%152 = icmp eq i64 %151, %107
br i1 %152, label %153, label %109
153: ; preds = %149, %101, %94
%154 = tail call i32 @free(double* %7)
ret i32 %95
}