float predict_pl(map, a, data, ivars, total_units, samples)
float **map, **a, **data;
int ivars, total_units, samples;
{
  int i, j, winner, ma = ivars + 1;
  float mean, std, rms_error, tmp;
  mean = 0;
  for (i = 1; i <= samples; i++)
    mean += data[i][ma];
  mean = mean / ((float)samples);
  rms_error = 0;
  std = 0;
  for (i = 1; i <= samples; i++) {
    tmp = data[i][ma];
    std += (tmp - mean) * (tmp - mean);
    winner = find_winner(map, data, i, ivars, total_units);
    data[i][ma] = a[winner][1];
    for (j = 1; j <= ivars; j++)
      data[i][ma] += data[i][j] * a[winner][j + 1];
    rms_error += (data[i][ma] - tmp) * (data[i][ma] - tmp);
  }
  std = sqrt(std / ((float)samples));
  rms_error = sqrt(rms_error / ((float)samples));
  printf("mean=%f std=%f err=%f ", mean, std, rms_error);
  rms_error = rms_error / std;
  printf("serr=%f\n", rms_error);
  return (rms_error);
}


x86 assembly (clang, -O0):
<predict_pl>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x38(%rbp)
movl   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jg     11a0 <predict_pl+0x70>
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1167 <predict_pl+0x37>
movss  -0x38(%rbp),%xmm0
cvtsi2ssl -0x24(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3c(%rbp)
movl   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jg     131e <predict_pl+0x1ee>
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x44(%rbp),%xmm0
subss  -0x38(%rbp),%xmm0
movss  -0x44(%rbp),%xmm2
subss  -0x38(%rbp),%xmm2
movss  -0x3c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x3c(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0x18(%rbp),%rsi
mov    -0x28(%rbp),%edx
mov    -0x1c(%rbp),%ecx
mov    -0x20(%rbp),%r8d
mov    $0x0,%al
call   1040 <find_winner@plt>
mov    %eax,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movss  0x4(%rax),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     12ca <predict_pl+0x19a>
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    125f <predict_pl+0x12f>
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x44(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
subss  -0x44(%rbp),%xmm2
movss  -0x40(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    11ca <predict_pl+0x9a>
movss  -0x3c(%rbp),%xmm0
cvtsi2ssl -0x24(%rbp),%xmm1
divss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
call   1050 <sqrt@plt>
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x40(%rbp),%xmm0
cvtsi2ssl -0x24(%rbp),%xmm1
divss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
call   1050 <sqrt@plt>
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x38(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x3c(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
movss  -0x40(%rbp),%xmm2
cvtss2sd %xmm2,%xmm2
lea    0xc80(%rip),%rdi
mov    $0x3,%al
call   1030 <printf@plt>
movss  -0x40(%rbp),%xmm0
divss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x40(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xc71(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movss  -0x40(%rbp),%xmm0
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<predict_pl>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %r9d,%r15d
mov    %r8d,0x1c(%rsp)
mov    %rsi,0x30(%rsp)
mov    %rdi,0x28(%rsp)
test   %r9d,%r9d
jle    1188 <predict_pl+0x58>
mov    %ecx,%ebp
mov    %rdx,%rbx
lea    0x1(%rbp),%eax
mov    %eax,0x18(%rsp)
movslq %eax,%r12
mov    %r15d,%edx
lea    -0x1(%rdx),%rcx
mov    %edx,%eax
and    $0x3,%eax
cmp    $0x3,%rcx
jae    119c <predict_pl+0x6c>
xorps  %xmm2,%xmm2
mov    $0x1,%ecx
test   %rax,%rax
jne    11ee <predict_pl+0xbe>
jmp    1213 <predict_pl+0xe3>
cvtsi2ss %r15d,%xmm1
xorps  %xmm0,%xmm0
xorps  %xmm2,%xmm2
divss  %xmm1,%xmm2
jmp    13b0 <predict_pl+0x280>
and    $0xfffffffc,%edx
xorps  %xmm2,%xmm2
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    0x8(%rbx,%rcx,8),%rsi
mov    0x10(%rbx,%rcx,8),%rdi
addss  (%rsi,%r12,4),%xmm2
addss  (%rdi,%r12,4),%xmm2
mov    0x18(%rbx,%rcx,8),%rsi
addss  (%rsi,%r12,4),%xmm2
mov    0x20(%rbx,%rcx,8),%rsi
addss  (%rsi,%r12,4),%xmm2
add    $0x4,%rcx
cmp    %rcx,%rdx
jne    11b0 <predict_pl+0x80>
add    $0x1,%rcx
test   %rax,%rax
je     1213 <predict_pl+0xe3>
lea    (%rbx,%rcx,8),%rcx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    (%rcx,%rdx,8),%rsi
addss  (%rsi,%r12,4),%xmm2
add    $0x1,%rdx
cmp    %rdx,%rax
jne    1200 <predict_pl+0xd0>
cvtsi2ss %r15d,%xmm1
divss  %xmm1,%xmm2
test   %r15d,%r15d
jle    13ad <predict_pl+0x27d>
movss  %xmm1,0x10(%rsp)
add    $0x1,%r15d
mov    0x18(%rsp),%r13d
add    $0xffffffffffffffff,%r13
mov    %r13,0x20(%rsp)
and    $0xfffffffffffffffe,%r13
neg    %r13
xorps  %xmm3,%xmm3
mov    $0x1,%r14d
xorps  %xmm0,%xmm0
movss  %xmm2,0xc(%rsp)
jmp    129d <predict_pl+0x16d>
nopl   0x0(%rax,%rax,1)
movss  0xc(%rsp),%xmm2
movss  0x4(%rsp),%xmm3
movss  0x8(%rsp),%xmm0
movss  0x14(%rsp),%xmm5
subss  %xmm5,%xmm4
subss  %xmm2,%xmm5
mulss  %xmm5,%xmm5
addss  %xmm5,%xmm0
mulss  %xmm4,%xmm4
addss  %xmm4,%xmm3
add    $0x1,%r14
cmp    %r15,%r14
je     13a5 <predict_pl+0x275>
movss  %xmm0,0x8(%rsp)
movss  %xmm3,0x4(%rsp)
mov    (%rbx,%r14,8),%rax
movss  (%rax,%r12,4),%xmm0
movss  %xmm0,0x14(%rsp)
mov    0x28(%rsp),%rdi
mov    %rbx,%rsi
mov    %r14d,%edx
mov    %ebp,%ecx
mov    0x1c(%rsp),%r8d
xor    %eax,%eax
call   1050 <find_winner@plt>
cltq
mov    0x30(%rsp),%rcx
mov    (%rcx,%rax,8),%rcx
movss  0x4(%rcx),%xmm4
mov    (%rbx,%r14,8),%rax
movss  %xmm4,(%rax,%r12,4)
test   %ebp,%ebp
jle    1260 <predict_pl+0x130>
mov    $0x1,%edx
cmpl   $0x2,0x18(%rsp)
jne    1320 <predict_pl+0x1f0>
movss  0xc(%rsp),%xmm2
movss  0x4(%rsp),%xmm3
movss  0x8(%rsp),%xmm0
movss  0x14(%rsp),%xmm5
jmp    1380 <predict_pl+0x250>
nopw   0x0(%rax,%rax,1)
mov    $0x1,%esi
movss  0xc(%rsp),%xmm2
movss  0x4(%rsp),%xmm3
movss  0x8(%rsp),%xmm0
movss  0x14(%rsp),%xmm5
nopl   (%rax)
movss  (%rax,%rsi,4),%xmm1
mulss  0x4(%rcx,%rsi,4),%xmm1
addss  %xmm4,%xmm1
movss  %xmm1,(%rax,%r12,4)
movss  0x4(%rax,%rsi,4),%xmm4
lea    0x2(%rsi),%rdx
mulss  0x8(%rcx,%rsi,4),%xmm4
addss  %xmm1,%xmm4
movss  %xmm4,(%rax,%r12,4)
lea    (%rsi,%r13,1),%rdi
add    $0x2,%rdi
mov    %rdx,%rsi
cmp    $0x1,%rdi
jne    1340 <predict_pl+0x210>
testb  $0x1,0x20(%rsp)
je     1278 <predict_pl+0x148>
movss  (%rax,%rdx,4),%xmm1
mulss  0x4(%rcx,%rdx,4),%xmm1
addss  %xmm1,%xmm4
movss  %xmm4,(%rax,%r12,4)
jmp    1278 <predict_pl+0x148>
movss  0x10(%rsp),%xmm1
jmp    13b3 <predict_pl+0x283>
xorps  %xmm0,%xmm0
xorps  %xmm3,%xmm3
divss  %xmm1,%xmm0
xorps  %xmm4,%xmm4
ucomiss %xmm4,%xmm0
movss  %xmm2,0xc(%rsp)
jb     13d1 <predict_pl+0x2a1>
sqrtss %xmm0,%xmm0
movss  %xmm0,0x8(%rsp)
jmp    13f7 <predict_pl+0x2c7>
movss  %xmm1,0x10(%rsp)
movss  %xmm3,0x4(%rsp)
call   1040 <sqrtf@plt>
xorps  %xmm4,%xmm4
movss  0x4(%rsp),%xmm3
movss  %xmm0,0x8(%rsp)
movss  0x10(%rsp),%xmm1
divss  %xmm1,%xmm3
ucomiss %xmm4,%xmm3
jb     1409 <predict_pl+0x2d9>
xorps  %xmm2,%xmm2
sqrtss %xmm3,%xmm2
jmp    1414 <predict_pl+0x2e4>
movaps %xmm3,%xmm0
call   1040 <sqrtf@plt>
movaps %xmm0,%xmm2
movss  0xc(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
movss  0x8(%rsp),%xmm1
cvtss2sd %xmm1,%xmm1
movss  %xmm2,0x4(%rsp)
cvtss2sd %xmm2,%xmm2
lea    0xbc7(%rip),%rdi
mov    $0x3,%al
call   1030 <printf@plt>
movss  0x4(%rsp),%xmm0
divss  0x8(%rsp),%xmm0
movss  %xmm0,0x4(%rsp)
cvtss2sd %xmm0,%xmm0
lea    0xbba(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movss  0x4(%rsp),%xmm0
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [23 x i8] c"mean=%f std=%f err=%f \00"
@.str.1 = [9 x i8] c"serr=%f\0A\00"
define float @predict_pl(float** %0, float** %1, float** %2, i32 %3, i32 %4, i32 %5) {
%7 = alloca float**
%8 = alloca float**
%9 = alloca float**
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca float
%18 = alloca float
%19 = alloca float
%20 = alloca float
store float** %0, float*** %7
store float** %1, float*** %8
store float** %2, float*** %9
store i32 %3, i32* %10
store i32 %4, i32* %11
store i32 %5, i32* %12
%21 = load i32, i32* %10
%22 = add nsw i32 %21, 1
store i32 %22, i32* %16
store float 0.000000e+00, float* %17
store i32 1, i32* %13
br label %23
23: ; preds = %39, %6
%24 = load i32, i32* %13
%25 = load i32, i32* %12
%26 = icmp sle i32 %24, %25
br i1 %26, label %27, label %42
27: ; preds = %23
%28 = load float**, float*** %9
%29 = load i32, i32* %13
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds float*, float** %28, i64 %30
%32 = load float*, float** %31
%33 = load i32, i32* %16
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds float, float* %32, i64 %34
%36 = load float, float* %35
%37 = load float, float* %17
%38 = fadd float %37, %36
store float %38, float* %17
br label %39
39: ; preds = %27
%40 = load i32, i32* %13
%41 = add nsw i32 %40, 1
store i32 %41, i32* %13
br label %23
42: ; preds = %23
%43 = load float, float* %17
%44 = load i32, i32* %12
%45 = sitofp i32 %44 to float
%46 = fdiv float %43, %45
store float %46, float* %17
store float 0.000000e+00, float* %19
store float 0.000000e+00, float* %18
store i32 1, i32* %13
br label %47
47: ; preds = %152, %42
%48 = load i32, i32* %13
%49 = load i32, i32* %12
%50 = icmp sle i32 %48, %49
br i1 %50, label %51, label %155
51: ; preds = %47
%52 = load float**, float*** %9
%53 = load i32, i32* %13
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds float*, float** %52, i64 %54
%56 = load float*, float** %55
%57 = load i32, i32* %16
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds float, float* %56, i64 %58
%60 = load float, float* %59
store float %60, float* %20
%61 = load float, float* %20
%62 = load float, float* %17
%63 = fsub float %61, %62
%64 = load float, float* %20
%65 = load float, float* %17
%66 = fsub float %64, %65
%67 = load float, float* %18
%68 = call float @llvm.fmuladd.f32(float %63, float %66, float %67)
store float %68, float* %18
%69 = load float**, float*** %7
%70 = load float**, float*** %9
%71 = load i32, i32* %13
%72 = load i32, i32* %10
%73 = load i32, i32* %11
%74 = call i32 (float**, float**, i32, i32, i32, ...) bitcast (i32 (...)* @find_winner to i32 (float**, float**, i32, i32, i32, ...)*)(float** %69, float** %70, i32 %71, i32 %72, i32 %73)
store i32 %74, i32* %15
%75 = load float**, float*** %8
%76 = load i32, i32* %15
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds float*, float** %75, i64 %77
%79 = load float*, float** %78
%80 = getelementptr inbounds float, float* %79, i64 1
%81 = load float, float* %80
%82 = load float**, float*** %9
%83 = load i32, i32* %13
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds float*, float** %82, i64 %84
%86 = load float*, float** %85
%87 = load i32, i32* %16
%88 = sext i32 %87 to i64
%89 = getelementptr inbounds float, float* %86, i64 %88
store float %81, float* %89
store i32 1, i32* %14
br label %90
90: ; preds = %124, %51
%91 = load i32, i32* %14
%92 = load i32, i32* %10
%93 = icmp sle i32 %91, %92
br i1 %93, label %94, label %127
94: ; preds = %90
%95 = load float**, float*** %9
%96 = load i32, i32* %13
%97 = sext i32 %96 to i64
%98 = getelementptr inbounds float*, float** %95, i64 %97
%99 = load float*, float** %98
%100 = load i32, i32* %14
%101 = sext i32 %100 to i64
%102 = getelementptr inbounds float, float* %99, i64 %101
%103 = load float, float* %102
%104 = load float**, float*** %8
%105 = load i32, i32* %15
%106 = sext i32 %105 to i64
%107 = getelementptr inbounds float*, float** %104, i64 %106
%108 = load float*, float** %107
%109 = load i32, i32* %14
%110 = add nsw i32 %109, 1
%111 = sext i32 %110 to i64
%112 = getelementptr inbounds float, float* %108, i64 %111
%113 = load float, float* %112
%114 = load float**, float*** %9
%115 = load i32, i32* %13
%116 = sext i32 %115 to i64
%117 = getelementptr inbounds float*, float** %114, i64 %116
%118 = load float*, float** %117
%119 = load i32, i32* %16
%120 = sext i32 %119 to i64
%121 = getelementptr inbounds float, float* %118, i64 %120
%122 = load float, float* %121
%123 = call float @llvm.fmuladd.f32(float %103, float %113, float %122)
store float %123, float* %121
br label %124
124: ; preds = %94
%125 = load i32, i32* %14
%126 = add nsw i32 %125, 1
store i32 %126, i32* %14
br label %90
127: ; preds = %90
%128 = load float**, float*** %9
%129 = load i32, i32* %13
%130 = sext i32 %129 to i64
%131 = getelementptr inbounds float*, float** %128, i64 %130
%132 = load float*, float** %131
%133 = load i32, i32* %16
%134 = sext i32 %133 to i64
%135 = getelementptr inbounds float, float* %132, i64 %134
%136 = load float, float* %135
%137 = load float, float* %20
%138 = fsub float %136, %137
%139 = load float**, float*** %9
%140 = load i32, i32* %13
%141 = sext i32 %140 to i64
%142 = getelementptr inbounds float*, float** %139, i64 %141
%143 = load float*, float** %142
%144 = load i32, i32* %16
%145 = sext i32 %144 to i64
%146 = getelementptr inbounds float, float* %143, i64 %145
%147 = load float, float* %146
%148 = load float, float* %20
%149 = fsub float %147, %148
%150 = load float, float* %19
%151 = call float @llvm.fmuladd.f32(float %138, float %149, float %150)
store float %151, float* %19
br label %152
152: ; preds = %127
%153 = load i32, i32* %13
%154 = add nsw i32 %153, 1
store i32 %154, i32* %13
br label %47
155: ; preds = %47
%156 = load float, float* %18
%157 = load i32, i32* %12
%158 = sitofp i32 %157 to float
%159 = fdiv float %156, %158
%160 = fpext float %159 to double
%161 = call double @sqrt(double %160)
%162 = fptrunc double %161 to float
store float %162, float* %18
%163 = load float, float* %19
%164 = load i32, i32* %12
%165 = sitofp i32 %164 to float
%166 = fdiv float %163, %165
%167 = fpext float %166 to double
%168 = call double @sqrt(double %167)
%169 = fptrunc double %168 to float
store float %169, float* %19
%170 = load float, float* %17
%171 = fpext float %170 to double
%172 = load float, float* %18
%173 = fpext float %172 to double
%174 = load float, float* %19
%175 = fpext float %174 to double
%176 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([23 x i8], [23 x i8]* @.str, i64 0, i64 0), double %171, double %173, double %175)
%177 = load float, float* %19
%178 = load float, float* %18
%179 = fdiv float %177, %178
store float %179, float* %19
%180 = load float, float* %19
%181 = fpext float %180 to double
%182 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([9 x i8], [9 x i8]* @.str.1, i64 0, i64 0), double %181)
%183 = load float, float* %19
ret float %183
}

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