#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 */
int /*<<< orphan*/ exit(int);
int fabs(float);
int /*<<< orphan*/ printf(char *, int);
void check_constant_column(float *C, int N, int K) {
  double mean, cov;
  int i, k, count;
  for (k = 0; k < K; k++) {
    mean = 0;
    count = 0;
    for (i = 0; i < N; i++) {
      if (fabs(C[i * K + k]) != 9) {
        mean += C[i * K + k];
        count++;
      }
    }
    if (count)
      mean /= count;
    else {
      printf("Error: SNP '%d' contains only missing data.\n\n", k + 1);
      exit(1);
    }
    cov = 0;
    for (i = 0; i < N; i++) {
      if (fabs(C[i * K + k]) != 9)
        cov += (C[i * K + k] - mean) * (C[i * K + k] - mean);
    }
    cov /= (count - 1);
    if (!cov) {
      printf("Warning: SNP %d is constant "
             "among individuals.\n\n",
             k + 1);
    }
  }
}


x86 assembly (clang, -O0):
<check_constant_column>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1310 <check_constant_column+0x1e0>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x18(%rbp)
movl   $0x0,-0x2c(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11db <check_constant_column+0xab>
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
call   1040 <fabs@plt>
cmp    $0x9,%eax
je     11c8 <check_constant_column+0x98>
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11cd <check_constant_column+0x9d>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    116b <check_constant_column+0x3b>
cmpl   $0x0,-0x2c(%rbp)
je     11fd <check_constant_column+0xcd>
cvtsi2sdl -0x2c(%rbp),%xmm1
movsd  -0x18(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
jmp    1219 <check_constant_column+0xe9>
mov    -0x28(%rbp),%esi
add    $0x1,%esi
lea    0xdf6(%rip),%rdi
call   1030 <printf@plt>
mov    $0x1,%edi
call   1050 <exit@plt>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12bb <check_constant_column+0x18b>
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
call   1040 <fabs@plt>
cmp    $0x9,%eax
je     12a8 <check_constant_column+0x178>
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
subsd  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm2
subsd  -0x18(%rbp),%xmm2
movsd  -0x20(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    12ad <check_constant_column+0x17d>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1228 <check_constant_column+0xf8>
mov    -0x2c(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
movsd  -0x20(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    12fd <check_constant_column+0x1cd>
jp     12fd <check_constant_column+0x1cd>
mov    -0x28(%rbp),%esi
add    $0x1,%esi
lea    0xd36(%rip),%rdi
call   1030 <printf@plt>
jmp    1302 <check_constant_column+0x1d2>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1149 <check_constant_column+0x19>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<check_constant_column>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
test   %edx,%edx
jle    140c <check_constant_column+0x2dc>
mov    %rdi,%r15
mov    %edx,%r14d
mov    %esi,%r13d
mov    %r13,0x28(%rsp)
and    $0xfffffffe,%r13d
lea    0x0(,%r14,8),%rax
mov    %rax,0x40(%rsp)
xor    %ebp,%ebp
mov    %rdi,0x20(%rsp)
mov    %rdi,0x30(%rsp)
mov    %esi,0x1c(%rsp)
jmp    1193 <check_constant_column+0x63>
nopw   0x0(%rax,%rax,1)
addq   $0x4,0x20(%rsp)
cmp    %r14,%rbp
mov    0x1c(%rsp),%esi
je     140c <check_constant_column+0x2dc>
test   %esi,%esi
jle    141b <check_constant_column+0x2eb>
mov    %rbp,0x38(%rsp)
xorps  %xmm0,%xmm0
cmp    $0x1,%esi
jne    12a0 <check_constant_column+0x170>
movsd  %xmm0,0x8(%rsp)
xor    %ebx,%ebx
xor    %r12d,%r12d
testb  $0x1,0x28(%rsp)
mov    0x30(%rsp),%r15
mov    0x38(%rsp),%rbp
je     11fd <check_constant_column+0xcd>
imul   %r14,%rbx
add    %rbp,%rbx
movss  (%r15,%rbx,4),%xmm0
call   1040 <fabs@plt>
cmp    $0x9,%eax
je     11fd <check_constant_column+0xcd>
movss  (%r15,%rbx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0x8(%rsp),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsp)
add    $0x1,%r12d
test   %r12d,%r12d
mov    0x1c(%rsp),%eax
je     141b <check_constant_column+0x2eb>
test   %eax,%eax
jle    1330 <check_constant_column+0x200>
xorps  %xmm0,%xmm0
cvtsi2sd %r12d,%xmm0
movsd  0x8(%rsp),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsp)
xorpd  %xmm0,%xmm0
cmp    $0x1,%eax
jne    1380 <check_constant_column+0x250>
movsd  %xmm0,0x10(%rsp)
xor    %ebx,%ebx
testb  $0x1,0x28(%rsp)
mov    0x30(%rsp),%r15
mov    0x38(%rsp),%rbp
je     133a <check_constant_column+0x20a>
imul   %r14,%rbx
add    %rbp,%rbx
movss  (%r15,%rbx,4),%xmm0
call   1040 <fabs@plt>
cmp    $0x9,%eax
je     133a <check_constant_column+0x20a>
movss  (%r15,%rbx,4),%xmm0
cvtss2sd %xmm0,%xmm0
subsd  0x8(%rsp),%xmm0
mulsd  %xmm0,%xmm0
movsd  0x10(%rsp),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x10(%rsp)
jmp    133a <check_constant_column+0x20a>
nopl   0x0(%rax,%rax,1)
movsd  %xmm0,0x8(%rsp)
xor    %ebx,%ebx
mov    0x20(%rsp),%r15
xor    %r12d,%r12d
mov    0x40(%rsp),%rbp
jmp    12d0 <check_constant_column+0x1a0>
nopw   0x0(%rax,%rax,1)
add    $0x2,%rbx
add    %rbp,%r15
cmp    %rbx,%r13
je     11b7 <check_constant_column+0x87>
movss  (%r15),%xmm0
call   1040 <fabs@plt>
cmp    $0x9,%eax
je     12fc <check_constant_column+0x1cc>
movss  (%r15),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0x8(%rsp),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsp)
add    $0x1,%r12d
movss  (%r15,%r14,4),%xmm0
call   1040 <fabs@plt>
cmp    $0x9,%eax
je     12c0 <check_constant_column+0x190>
movss  (%r15,%r14,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0x8(%rsp),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsp)
add    $0x1,%r12d
jmp    12c0 <check_constant_column+0x190>
nopl   0x0(%rax)
xorpd  %xmm0,%xmm0
movsd  %xmm0,0x10(%rsp)
add    $0xffffffff,%r12d
xorps  %xmm0,%xmm0
cvtsi2sd %r12d,%xmm0
movsd  0x10(%rsp),%xmm1
divsd  %xmm0,%xmm1
add    $0x1,%rbp
ucomisd 0xca4(%rip),%xmm1
jne    1180 <check_constant_column+0x50>
jp     1180 <check_constant_column+0x50>
lea    0xcc7(%rip),%rdi
mov    %ebp,%esi
call   1030 <printf@plt>
jmp    1180 <check_constant_column+0x50>
nopl   0x0(%rax,%rax,1)
movsd  %xmm0,0x10(%rsp)
mov    0x20(%rsp),%rbp
xor    %ebx,%ebx
jmp    13a2 <check_constant_column+0x272>
nop
add    $0x2,%rbx
add    0x40(%rsp),%rbp
cmp    %rbx,%r13
je     123f <check_constant_column+0x10f>
movss  0x0(%rbp),%xmm0
call   1040 <fabs@plt>
cmp    $0x9,%eax
je     13d4 <check_constant_column+0x2a4>
movss  0x0(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
subsd  0x8(%rsp),%xmm0
mulsd  %xmm0,%xmm0
movsd  0x10(%rsp),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x10(%rsp)
movss  0x0(%rbp,%r14,4),%xmm0
call   1040 <fabs@plt>
cmp    $0x9,%eax
je     1390 <check_constant_column+0x260>
movss  0x0(%rbp,%r14,4),%xmm0
cvtss2sd %xmm0,%xmm0
subsd  0x8(%rsp),%xmm0
mulsd  %xmm0,%xmm0
movsd  0x10(%rsp),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x10(%rsp)
jmp    1390 <check_constant_column+0x260>
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
add    $0x1,%ebp
lea    0xbe3(%rip),%rdi
mov    %ebp,%esi
call   1030 <printf@plt>
mov    $0x1,%edi
call   1050 <exit@plt>

LLVM IR (-O0):
@.str = [46 x i8] c"Error: SNP '%d' contains only missing data.\0A\0A\00"
@.str.1 = [49 x i8] c"Warning: SNP %d is constant among individuals.\0A\0A\00"
define void @check_constant_column(float* %0, i32 %1, i32 %2) {
%4 = alloca float*
%5 = alloca i32
%6 = alloca i32
%7 = alloca double
%8 = alloca double
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
store float* %0, float** %4
store i32 %1, i32* %5
store i32 %2, i32* %6
store i32 0, i32* %10
br label %12
12: ; preds = %126, %3
%13 = load i32, i32* %10
%14 = load i32, i32* %6
%15 = icmp slt i32 %13, %14
br i1 %15, label %16, label %129
16: ; preds = %12
store double 0.000000e+00, double* %7
store i32 0, i32* %11
store i32 0, i32* %9
br label %17
17: ; preds = %49, %16
%18 = load i32, i32* %9
%19 = load i32, i32* %5
%20 = icmp slt i32 %18, %19
br i1 %20, label %21, label %52
21: ; preds = %17
%22 = load float*, float** %4
%23 = load i32, i32* %9
%24 = load i32, i32* %6
%25 = mul nsw i32 %23, %24
%26 = load i32, i32* %10
%27 = add nsw i32 %25, %26
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds float, float* %22, i64 %28
%30 = load float, float* %29
%31 = call i32 @fabs(float %30)
%32 = icmp ne i32 %31, 9
br i1 %32, label %33, label %48
33: ; preds = %21
%34 = load float*, float** %4
%35 = load i32, i32* %9
%36 = load i32, i32* %6
%37 = mul nsw i32 %35, %36
%38 = load i32, i32* %10
%39 = add nsw i32 %37, %38
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds float, float* %34, i64 %40
%42 = load float, float* %41
%43 = fpext float %42 to double
%44 = load double, double* %7
%45 = fadd double %44, %43
store double %45, double* %7
%46 = load i32, i32* %11
%47 = add nsw i32 %46, 1
store i32 %47, i32* %11
br label %48
48: ; preds = %33, %21
br label %49
49: ; preds = %48
%50 = load i32, i32* %9
%51 = add nsw i32 %50, 1
store i32 %51, i32* %9
br label %17
52: ; preds = %17
%53 = load i32, i32* %11
%54 = icmp ne i32 %53, 0
br i1 %54, label %55, label %60
55: ; preds = %52
%56 = load i32, i32* %11
%57 = sitofp i32 %56 to double
%58 = load double, double* %7
%59 = fdiv double %58, %57
store double %59, double* %7
br label %65
60: ; preds = %52
%61 = load i32, i32* %10
%62 = add nsw i32 %61, 1
%63 = call i32 @printf(i8* getelementptr inbounds ([46 x i8], [46 x i8]* @.str, i64 0, i64 0), i32 %62)
%64 = call i32 @exit(i32 1)
unreachable
65: ; preds = %55
store double 0.000000e+00, double* %8
store i32 0, i32* %9
br label %66
66: ; preds = %110, %65
%67 = load i32, i32* %9
%68 = load i32, i32* %5
%69 = icmp slt i32 %67, %68
br i1 %69, label %70, label %113
70: ; preds = %66
%71 = load float*, float** %4
%72 = load i32, i32* %9
%73 = load i32, i32* %6
%74 = mul nsw i32 %72, %73
%75 = load i32, i32* %10
%76 = add nsw i32 %74, %75
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds float, float* %71, i64 %77
%79 = load float, float* %78
%80 = call i32 @fabs(float %79)
%81 = icmp ne i32 %80, 9
br i1 %81, label %82, label %109
82: ; preds = %70
%83 = load float*, float** %4
%84 = load i32, i32* %9
%85 = load i32, i32* %6
%86 = mul nsw i32 %84, %85
%87 = load i32, i32* %10
%88 = add nsw i32 %86, %87
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds float, float* %83, i64 %89
%91 = load float, float* %90
%92 = fpext float %91 to double
%93 = load double, double* %7
%94 = fsub double %92, %93
%95 = load float*, float** %4
%96 = load i32, i32* %9
%97 = load i32, i32* %6
%98 = mul nsw i32 %96, %97
%99 = load i32, i32* %10
%100 = add nsw i32 %98, %99
%101 = sext i32 %100 to i64
%102 = getelementptr inbounds float, float* %95, i64 %101
%103 = load float, float* %102
%104 = fpext float %103 to double
%105 = load double, double* %7
%106 = fsub double %104, %105
%107 = load double, double* %8
%108 = call double @llvm.fmuladd.f64(double %94, double %106, double %107)
store double %108, double* %8
br label %109
109: ; preds = %82, %70
br label %110
110: ; preds = %109
%111 = load i32, i32* %9
%112 = add nsw i32 %111, 1
store i32 %112, i32* %9
br label %66
113: ; preds = %66
%114 = load i32, i32* %11
%115 = sub nsw i32 %114, 1
%116 = sitofp i32 %115 to double
%117 = load double, double* %8
%118 = fdiv double %117, %116
store double %118, double* %8
%119 = load double, double* %8
%120 = fcmp une double %119, 0.000000e+00
br i1 %120, label %125, label %121
121: ; preds = %113
%122 = load i32, i32* %10
%123 = add nsw i32 %122, 1
%124 = call i32 @printf(i8* getelementptr inbounds ([49 x i8], [49 x i8]* @.str.1, i64 0, i64 0), i32 %123)
br label %125
125: ; preds = %121, %113
br label %126
126: ; preds = %125
%127 = load i32, i32* %10
%128 = add nsw i32 %127, 1
store i32 %128, i32* %10
br label %12
129: ; preds = %12
ret void
}

LLVM IR (-O2):
@.str = [46 x i8] c"Error: SNP '%d' contains only missing data.\0A\0A\00"
@.str.1 = [49 x i8] c"Warning: SNP %d is constant among individuals.\0A\0A\00"
define void @check_constant_column(float* nocapture %0, i32 %1, i32 %2) {
%4 = icmp sgt i32 %2, 0
br i1 %4, label %5, label %151
5: ; preds = %3
%6 = icmp sgt i32 %1, 0
%7 = zext i32 %2 to i64
%8 = zext i32 %2 to i64
%9 = zext i32 %1 to i64
%10 = and i64 %9, 1
%11 = icmp eq i32 %1, 1
%12 = and i64 %9, 4294967294
%13 = icmp eq i64 %10, 0
%14 = and i64 %9, 1
%15 = icmp eq i32 %1, 1
%16 = and i64 %9, 4294967294
%17 = icmp eq i64 %14, 0
br label %18
18: ; preds = %5, %149
%19 = phi i64 [ 0, %5 ], [ %145, %149 ]
br i1 %6, label %20, label %84
20: ; preds = %18
br i1 %11, label %58, label %21
21: ; preds = %20, %52
%22 = phi i64 [ %55, %52 ], [ 0, %20 ]
%23 = phi i32 [ %54, %52 ], [ 0, %20 ]
%24 = phi double [ %53, %52 ], [ 0.000000e+00, %20 ]
%25 = phi i64 [ %56, %52 ], [ 0, %20 ]
%26 = mul nsw i64 %22, %7
%27 = add nuw nsw i64 %26, %19
%28 = getelementptr inbounds float, float* %0, i64 %27
%29 = load float, float* %28
%30 = tail call i32 @fabs(float %29)
%31 = icmp eq i32 %30, 9
br i1 %31, label %37, label %32
32: ; preds = %21
%33 = load float, float* %28
%34 = fpext float %33 to double
%35 = fadd double %24, %34
%36 = add nsw i32 %23, 1
br label %37
37: ; preds = %21, %32
%38 = phi double [ %35, %32 ], [ %24, %21 ]
%39 = phi i32 [ %36, %32 ], [ %23, %21 ]
%40 = or i64 %22, 1
%41 = mul nsw i64 %40, %7
%42 = add nuw nsw i64 %41, %19
%43 = getelementptr inbounds float, float* %0, i64 %42
%44 = load float, float* %43
%45 = tail call i32 @fabs(float %44)
%46 = icmp eq i32 %45, 9
br i1 %46, label %52, label %47
47: ; preds = %37
%48 = load float, float* %43
%49 = fpext float %48 to double
%50 = fadd double %38, %49
%51 = add nsw i32 %39, 1
br label %52
52: ; preds = %47, %37
%53 = phi double [ %50, %47 ], [ %38, %37 ]
%54 = phi i32 [ %51, %47 ], [ %39, %37 ]
%55 = add nuw nsw i64 %22, 2
%56 = add i64 %25, 2
%57 = icmp eq i64 %56, %12
br i1 %57, label %58, label %21
58: ; preds = %52, %20
%59 = phi double [ undef, %20 ], [ %53, %52 ]
%60 = phi i32 [ undef, %20 ], [ %54, %52 ]
%61 = phi i64 [ 0, %20 ], [ %55, %52 ]
%62 = phi i32 [ 0, %20 ], [ %54, %52 ]
%63 = phi double [ 0.000000e+00, %20 ], [ %53, %52 ]
br i1 %13, label %76, label %64
64: ; preds = %58
%65 = mul nsw i64 %61, %7
%66 = add nuw nsw i64 %65, %19
%67 = getelementptr inbounds float, float* %0, i64 %66
%68 = load float, float* %67
%69 = tail call i32 @fabs(float %68)
%70 = icmp eq i32 %69, 9
br i1 %70, label %76, label %71
71: ; preds = %64
%72 = load float, float* %67
%73 = fpext float %72 to double
%74 = fadd double %63, %73
%75 = add nsw i32 %62, 1
br label %76
76: ; preds = %64, %71, %58
%77 = phi double [ %59, %58 ], [ %74, %71 ], [ %63, %64 ]
%78 = phi i32 [ %60, %58 ], [ %75, %71 ], [ %62, %64 ]
%79 = icmp eq i32 %78, 0
br i1 %79, label %84, label %80
80: ; preds = %76
%81 = sitofp i32 %78 to double
%82 = fdiv double %77, %81
br i1 %6, label %83, label %139
83: ; preds = %80
br i1 %15, label %123, label %89
84: ; preds = %18, %76
%85 = trunc i64 %19 to i32
%86 = add nuw nsw i32 %85, 1
%87 = tail call i32 @printf(i8* getelementptr inbounds ([46 x i8], [46 x i8]* @.str, i64 0, i64 0), i32 %86)
%88 = tail call i32 @exit(i32 1)
unreachable
89: ; preds = %83, %118
%90 = phi i64 [ %120, %118 ], [ 0, %83 ]
%91 = phi double [ %119, %118 ], [ 0.000000e+00, %83 ]
%92 = phi i64 [ %121, %118 ], [ 0, %83 ]
%93 = mul nsw i64 %90, %7
%94 = add nuw nsw i64 %93, %19
%95 = getelementptr inbounds float, float* %0, i64 %94
%96 = load float, float* %95
%97 = tail call i32 @fabs(float %96)
%98 = icmp eq i32 %97, 9
br i1 %98, label %104, label %99
99: ; preds = %89
%100 = load float, float* %95
%101 = fpext float %100 to double
%102 = fsub double %101, %82
%103 = tail call double @llvm.fmuladd.f64(double %102, double %102, double %91)
br label %104
104: ; preds = %89, %99
%105 = phi double [ %103, %99 ], [ %91, %89 ]
%106 = or i64 %90, 1
%107 = mul nsw i64 %106, %7
%108 = add nuw nsw i64 %107, %19
%109 = getelementptr inbounds float, float* %0, i64 %108
%110 = load float, float* %109
%111 = tail call i32 @fabs(float %110)
%112 = icmp eq i32 %111, 9
br i1 %112, label %118, label %113
113: ; preds = %104
%114 = load float, float* %109
%115 = fpext float %114 to double
%116 = fsub double %115, %82
%117 = tail call double @llvm.fmuladd.f64(double %116, double %116, double %105)
br label %118
118: ; preds = %113, %104
%119 = phi double [ %117, %113 ], [ %105, %104 ]
%120 = add nuw nsw i64 %90, 2
%121 = add i64 %92, 2
%122 = icmp eq i64 %121, %16
br i1 %122, label %123, label %89
123: ; preds = %118, %83
%124 = phi double [ undef, %83 ], [ %119, %118 ]
%125 = phi i64 [ 0, %83 ], [ %120, %118 ]
%126 = phi double [ 0.000000e+00, %83 ], [ %119, %118 ]
br i1 %17, label %139, label %127
127: ; preds = %123
%128 = mul nsw i64 %125, %7
%129 = add nuw nsw i64 %128, %19
%130 = getelementptr inbounds float, float* %0, i64 %129
%131 = load float, float* %130
%132 = tail call i32 @fabs(float %131)
%133 = icmp eq i32 %132, 9
br i1 %133, label %139, label %134
134: ; preds = %127
%135 = load float, float* %130
%136 = fpext float %135 to double
%137 = fsub double %136, %82
%138 = tail call double @llvm.fmuladd.f64(double %137, double %137, double %126)
br label %139
139: ; preds = %123, %134, %127, %80
%140 = phi double [ 0.000000e+00, %80 ], [ %124, %123 ], [ %138, %134 ], [ %126, %127 ]
%141 = add nsw i32 %78, -1
%142 = sitofp i32 %141 to double
%143 = fdiv double %140, %142
%144 = fcmp une double %143, 0.000000e+00
%145 = add nuw nsw i64 %19, 1
br i1 %144, label %149, label %146
146: ; preds = %139
%147 = trunc i64 %145 to i32
%148 = tail call i32 @printf(i8* getelementptr inbounds ([49 x i8], [49 x i8]* @.str.1, i64 0, i64 0), i32 %147)
br label %149
149: ; preds = %139, %146
%150 = icmp eq i64 %145, %8
br i1 %150, label %151, label %18
151: ; preds = %149, %3
ret void
}