#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 */
double fabs(float);
void levinson_durbin(float R[], float lpcs[], int order) {
  float E[20 + 1];
  float k[20 + 1];
  float a[20 + 1][20 + 1];
  float sum;
  int i, j;
  E[0] = R[0];
  for (i = 1; i <= order; i++) {
    sum = 0.0;
    for (j = 1; j <= i - 1; j++)
      sum += a[i - 1][j] * R[i - j];
    k[i] = -1.0 * (R[i] + sum) / E[i - 1];
    if (fabs(k[i]) > 1.0)
      k[i] = 0.0;
    a[i][i] = k[i];
    for (j = 1; j <= i - 1; j++)
      a[i][j] = a[i - 1][j] + k[i] * a[i - 1][i - j];
    E[i] = (1 - k[i] * k[i]) * E[i - 1];
  }
  for (i = 1; i <= order; i++)
    lpcs[i] = a[order][i];
  lpcs[0] = 1.0;
}


x86 assembly (clang, -O0):
<levinson_durbin>:
push   %rbp
mov    %rsp,%rbp
sub    $0x7d0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x70(%rbp)
movl   $0x1,-0x7c8(%rbp)
mov    -0x7c8(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     13d3 <levinson_durbin+0x2c3>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x7c4(%rbp)
ff
movl   $0x1,-0x7cc(%rbp)
mov    -0x7cc(%rbp),%eax
mov    -0x7c8(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     11e2 <levinson_durbin+0xd2>
mov    -0x7c8(%rbp),%eax
sub    $0x1,%eax
movslq %eax,%rcx
lea    -0x7c0(%rbp),%rax
imul   $0x54,%rcx,%rcx
add    %rcx,%rax
movslq -0x7cc(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x7c8(%rbp),%ecx
sub    -0x7cc(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x7c4(%rbp),%xmm1
ff
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x7c4(%rbp)
ff
mov    -0x7cc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7cc(%rbp)
jmp    1161 <levinson_durbin+0x51>
mov    -0x8(%rbp),%rax
movslq -0x7c8(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x7c4(%rbp),%xmm0
ff
cvtss2sd %xmm0,%xmm1
movsd  0xe0a(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x7c8(%rbp),%eax
sub    $0x1,%eax
cltq
movss  -0x70(%rbp,%rax,4),%xmm1
cvtss2sd %xmm1,%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movslq -0x7c8(%rbp),%rax
movss  %xmm0,-0xd0(%rbp,%rax,4)
ff ff
movslq -0x7c8(%rbp),%rax
movss  -0xd0(%rbp,%rax,4),%xmm0
ff ff
call   1030 <fabs@plt>
movsd  0xdb4(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1271 <levinson_durbin+0x161>
movslq -0x7c8(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,-0xd0(%rbp,%rax,4)
ff ff
movslq -0x7c8(%rbp),%rax
movss  -0xd0(%rbp,%rax,4),%xmm0
ff ff
movslq -0x7c8(%rbp),%rcx
lea    -0x7c0(%rbp),%rax
imul   $0x54,%rcx,%rcx
add    %rcx,%rax
movslq -0x7c8(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movl   $0x1,-0x7cc(%rbp)
mov    -0x7cc(%rbp),%eax
mov    -0x7c8(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     1364 <levinson_durbin+0x254>
mov    -0x7c8(%rbp),%eax
sub    $0x1,%eax
movslq %eax,%rcx
lea    -0x7c0(%rbp),%rax
imul   $0x54,%rcx,%rcx
add    %rcx,%rax
movslq -0x7cc(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
movslq -0x7c8(%rbp),%rax
movss  -0xd0(%rbp,%rax,4),%xmm0
ff ff
mov    -0x7c8(%rbp),%eax
sub    $0x1,%eax
movslq %eax,%rcx
lea    -0x7c0(%rbp),%rax
imul   $0x54,%rcx,%rcx
add    %rcx,%rax
mov    -0x7c8(%rbp),%ecx
sub    -0x7cc(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movslq -0x7c8(%rbp),%rcx
lea    -0x7c0(%rbp),%rax
imul   $0x54,%rcx,%rcx
add    %rcx,%rax
movslq -0x7cc(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x7cc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7cc(%rbp)
jmp    12ac <levinson_durbin+0x19c>
movslq -0x7c8(%rbp),%rax
movss  -0xd0(%rbp,%rax,4),%xmm0
ff ff
movslq -0x7c8(%rbp),%rax
movss  -0xd0(%rbp,%rax,4),%xmm1
ff ff
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm1,%xmm0
movss  0xc63(%rip),%xmm1
addss  %xmm1,%xmm0
mov    -0x7c8(%rbp),%eax
sub    $0x1,%eax
cltq
mulss  -0x70(%rbp,%rax,4),%xmm0
movslq -0x7c8(%rbp),%rax
movss  %xmm0,-0x70(%rbp,%rax,4)
mov    -0x7c8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7c8(%rbp)
jmp    113d <levinson_durbin+0x2d>
movl   $0x1,-0x7c8(%rbp)
mov    -0x7c8(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     142e <levinson_durbin+0x31e>
movslq -0x14(%rbp),%rcx
lea    -0x7c0(%rbp),%rax
imul   $0x54,%rcx,%rcx
add    %rcx,%rax
movslq -0x7c8(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x7c8(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x7c8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7c8(%rbp)
jmp    13dd <levinson_durbin+0x2cd>
mov    -0x10(%rbp),%rax
movss  0xbc6(%rip),%xmm0
movss  %xmm0,(%rax)
add    $0x7d0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<levinson_durbin>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x728,%rsp
mov    %rsi,%rbx
test   %edx,%edx
jle    139f <levinson_durbin+0x27f>
mov    %rbx,0x20(%rsp)
movss  (%rdi),%xmm3
mov    %rdx,0x18(%rsp)
lea    0x1(%rdx),%eax
mov    %rax,0x28(%rsp)
lea    0x44(%rsp),%rbx
lea    0x9c(%rsp),%rbp
mov    $0x1,%r14d
lea    0x40(%rsp),%r15
xor    %r8d,%r8d
mov    %rdi,%r9
jmp    11b7 <levinson_durbin+0x97>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mulss  %xmm1,%xmm1
movss  0xe8c(%rip),%xmm0
subss  %xmm1,%xmm0
mulss  %xmm0,%xmm3
add    $0x1,%r14
add    $0x1,%r8
add    $0x4,%r9
add    $0x54,%rbx
add    $0x58,%r15
add    $0x54,%rbp
cmp    0x28(%rsp),%r14
je     1376 <levinson_durbin+0x256>
mov    %r8,%r12
and    $0xfffffffffffffffe,%r12
lea    -0x1(%r14),%rax
mov    %rax,0x10(%rsp)
xorps  %xmm0,%xmm0
cmp    $0x2,%r14
jb     126c <levinson_durbin+0x14c>
cmp    $0x1,%r8
jne    11f0 <levinson_durbin+0xd0>
xorps  %xmm0,%xmm0
mov    $0x1,%eax
jmp    1245 <levinson_durbin+0x125>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %r12,%rcx
neg    %rcx
xorps  %xmm1,%xmm1
mov    $0x1,%eax
mov    %r9,%rdx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movss  -0x4(%rbx,%rax,4),%xmm2
movss  (%rbx,%rax,4),%xmm0
mulss  (%rdx),%xmm2
addss  %xmm1,%xmm2
mulss  -0x4(%rdx),%xmm0
addss  %xmm2,%xmm0
add    $0xfffffffffffffff8,%rdx
lea    (%rcx,%rax,1),%rsi
add    $0x2,%rsi
add    $0x2,%rax
movaps %xmm0,%xmm1
cmp    $0x1,%rsi
jne    1210 <levinson_durbin+0xf0>
test   $0x1,%r8b
je     126c <levinson_durbin+0x14c>
imul   $0x54,0x10(%rsp),%rcx
add    %rsp,%rcx
add    $0x40,%rcx
movss  (%rcx,%rax,4),%xmm1
mov    %r14,%rcx
sub    %rax,%rcx
mulss  (%rdi,%rcx,4),%xmm1
addss  %xmm1,%xmm0
mov    %r9,0x30(%rsp)
mov    %r8,0x38(%rsp)
mov    %rdi,%r13
addss  (%rdi,%r14,4),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  0xd76(%rip),%xmm0
movss  %xmm3,0xc(%rsp)
xorps  %xmm1,%xmm1
cvtss2sd %xmm3,%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x8(%rsp)
call   1030 <fabs@plt>
ucomisd 0xd5e(%rip),%xmm0
xorps  %xmm1,%xmm1
ja     12bd <levinson_durbin+0x19d>
movss  0x8(%rsp),%xmm1
imul   $0x54,%r14,%rax
add    %rsp,%rax
add    $0x40,%rax
movss  %xmm1,(%rax,%r14,4)
cmp    $0x2,%r14
mov    %r13,%rdi
movss  0xc(%rsp),%xmm3
mov    0x38(%rsp),%r8
mov    0x30(%rsp),%r9
jb     1180 <levinson_durbin+0x60>
mov    $0x1,%ecx
cmp    $0x1,%r8
je     1341 <levinson_durbin+0x221>
mov    $0xffffffffffffffff,%rcx
mov    %rbp,%rdx
movss  0x4(%r15,%rcx,4),%xmm0
mulss  %xmm1,%xmm0
addss  -0x58(%rdx),%xmm0
movss  %xmm0,-0x4(%rdx)
movss  (%r15,%rcx,4),%xmm0
mulss  %xmm1,%xmm0
addss  -0x54(%rdx),%xmm0
movss  %xmm0,(%rdx)
add    $0x8,%rdx
lea    (%r12,%rcx,1),%rsi
add    $0xfffffffffffffffe,%rsi
add    $0xfffffffffffffffe,%rcx
cmp    $0xffffffffffffffff,%rsi
jne    1300 <levinson_durbin+0x1e0>
neg    %rcx
test   $0x1,%r8b
je     1180 <levinson_durbin+0x60>
imul   $0x54,0x10(%rsp),%rdx
add    %rsp,%rdx
add    $0x40,%rdx
mov    %r14,%rsi
sub    %rcx,%rsi
movss  (%rdx,%rsi,4),%xmm0
mulss  %xmm1,%xmm0
addss  (%rdx,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    1180 <levinson_durbin+0x60>
mov    0x18(%rsp),%rax
test   %eax,%eax
mov    0x20(%rsp),%rbx
jle    139f <levinson_durbin+0x27f>
mov    %eax,%edx
lea    0x4(%rbx),%rdi
imul   $0x54,%rdx,%rax
lea    (%rsp,%rax,1),%rsi
add    $0x44,%rsi
shl    $0x2,%rdx
call   1040 <memcpy@plt>
movl   $0x3f800000,(%rbx)
add    $0x728,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @levinson_durbin(float* nocapture %0, float* nocapture writeonly %1, i32 %2) {
%4 = alloca [21 x [21 x float]]
%5 = bitcast [21 x [21 x float]]* %4 to i8*
%6 = icmp slt i32 %2, 1
br i1 %6, label %121, label %7
7: ; preds = %3
%8 = load float, float* %0
%9 = add nuw i32 %2, 1
%10 = zext i32 %9 to i64
br label %11
11: ; preds = %7, %114
%12 = phi i64 [ 0, %7 ], [ %120, %114 ]
%13 = phi float [ %8, %7 ], [ %117, %114 ]
%14 = phi i64 [ 1, %7 ], [ %118, %114 ]
%15 = add nsw i64 %14, -1
%16 = icmp ugt i64 %14, 1
br i1 %16, label %17, label %63
17: ; preds = %11
%18 = and i64 %12, 1
%19 = icmp eq i64 %12, 1
br i1 %19, label %51, label %20
20: ; preds = %17
%21 = and i64 %12, -2
br label %31
22: ; preds = %114
br i1 %6, label %121, label %23
23: ; preds = %22
%24 = zext i32 %2 to i64
%25 = getelementptr float, float* %1, i64 1
%26 = bitcast float* %25 to i8*
%27 = getelementptr [21 x [21 x float]], [21 x [21 x float]]* %4, i64 0, i64 %24, i64 1
%28 = bitcast float* %27 to i8*
%29 = zext i32 %2 to i64
%30 = shl nuw nsw i64 %29, 2
br label %121
31: ; preds = %31, %20
%32 = phi i64 [ 1, %20 ], [ %48, %31 ]
%33 = phi float [ 0.000000e+00, %20 ], [ %47, %31 ]
%34 = phi i64 [ 0, %20 ], [ %49, %31 ]
%35 = getelementptr inbounds [21 x [21 x float]], [21 x [21 x float]]* %4, i64 0, i64 %15, i64 %32
%36 = load float, float* %35
%37 = sub nsw i64 %14, %32
%38 = getelementptr inbounds float, float* %0, i64 %37
%39 = load float, float* %38
%40 = tail call float @llvm.fmuladd.f32(float %36, float %39, float %33)
%41 = add nuw nsw i64 %32, 1
%42 = getelementptr inbounds [21 x [21 x float]], [21 x [21 x float]]* %4, i64 0, i64 %15, i64 %41
%43 = load float, float* %42
%44 = sub nsw i64 %14, %41
%45 = getelementptr inbounds float, float* %0, i64 %44
%46 = load float, float* %45
%47 = tail call float @llvm.fmuladd.f32(float %43, float %46, float %40)
%48 = add nuw nsw i64 %32, 2
%49 = add i64 %34, 2
%50 = icmp eq i64 %49, %21
br i1 %50, label %51, label %31
51: ; preds = %31, %17
%52 = phi float [ undef, %17 ], [ %47, %31 ]
%53 = phi i64 [ 1, %17 ], [ %48, %31 ]
%54 = phi float [ 0.000000e+00, %17 ], [ %47, %31 ]
%55 = icmp eq i64 %18, 0
br i1 %55, label %63, label %56
56: ; preds = %51
%57 = getelementptr inbounds [21 x [21 x float]], [21 x [21 x float]]* %4, i64 0, i64 %15, i64 %53
%58 = load float, float* %57
%59 = sub nsw i64 %14, %53
%60 = getelementptr inbounds float, float* %0, i64 %59
%61 = load float, float* %60
%62 = tail call float @llvm.fmuladd.f32(float %58, float %61, float %54)
br label %63
63: ; preds = %56, %51, %11
%64 = phi float [ 0.000000e+00, %11 ], [ %52, %51 ], [ %62, %56 ]
%65 = getelementptr inbounds float, float* %0, i64 %14
%66 = load float, float* %65
%67 = fadd float %64, %66
%68 = fpext float %67 to double
%69 = fneg double %68
%70 = fpext float %13 to double
%71 = fdiv double %69, %70
%72 = fptrunc double %71 to float
%73 = tail call double @fabs(float %72)
%74 = fcmp ogt double %73, 1.000000e+00
%75 = select i1 %74, float 0.000000e+00, float %72
%76 = getelementptr inbounds [21 x [21 x float]], [21 x [21 x float]]* %4, i64 0, i64 %14, i64 %14
store float %75, float* %76
br i1 %16, label %77, label %114
77: ; preds = %63
%78 = and i64 %12, 1
%79 = icmp eq i64 %12, 1
br i1 %79, label %103, label %80
80: ; preds = %77
%81 = and i64 %12, -2
br label %82
82: ; preds = %82, %80
%83 = phi i64 [ 1, %80 ], [ %100, %82 ]
%84 = phi i64 [ 0, %80 ], [ %101, %82 ]
%85 = getelementptr inbounds [21 x [21 x float]], [21 x [21 x float]]* %4, i64 0, i64 %15, i64 %83
%86 = load float, float* %85
%87 = sub nsw i64 %14, %83
%88 = getelementptr inbounds [21 x [21 x float]], [21 x [21 x float]]* %4, i64 0, i64 %15, i64 %87
%89 = load float, float* %88
%90 = tail call float @llvm.fmuladd.f32(float %75, float %89, float %86)
%91 = getelementptr inbounds [21 x [21 x float]], [21 x [21 x float]]* %4, i64 0, i64 %14, i64 %83
store float %90, float* %91
%92 = add nuw nsw i64 %83, 1
%93 = getelementptr inbounds [21 x [21 x float]], [21 x [21 x float]]* %4, i64 0, i64 %15, i64 %92
%94 = load float, float* %93
%95 = sub nsw i64 %14, %92
%96 = getelementptr inbounds [21 x [21 x float]], [21 x [21 x float]]* %4, i64 0, i64 %15, i64 %95
%97 = load float, float* %96
%98 = tail call float @llvm.fmuladd.f32(float %75, float %97, float %94)
%99 = getelementptr inbounds [21 x [21 x float]], [21 x [21 x float]]* %4, i64 0, i64 %14, i64 %92
store float %98, float* %99
%100 = add nuw nsw i64 %83, 2
%101 = add i64 %84, 2
%102 = icmp eq i64 %101, %81
br i1 %102, label %103, label %82
103: ; preds = %82, %77
%104 = phi i64 [ 1, %77 ], [ %100, %82 ]
%105 = icmp eq i64 %78, 0
br i1 %105, label %114, label %106
106: ; preds = %103
%107 = getelementptr inbounds [21 x [21 x float]], [21 x [21 x float]]* %4, i64 0, i64 %15, i64 %104
%108 = load float, float* %107
%109 = sub nsw i64 %14, %104
%110 = getelementptr inbounds [21 x [21 x float]], [21 x [21 x float]]* %4, i64 0, i64 %15, i64 %109
%111 = load float, float* %110
%112 = tail call float @llvm.fmuladd.f32(float %75, float %111, float %108)
%113 = getelementptr inbounds [21 x [21 x float]], [21 x [21 x float]]* %4, i64 0, i64 %14, i64 %104
store float %112, float* %113
br label %114
114: ; preds = %106, %103, %63
%115 = fneg float %75
%116 = tail call float @llvm.fmuladd.f32(float %115, float %75, float 1.000000e+00)
%117 = fmul float %13, %116
%118 = add nuw nsw i64 %14, 1
%119 = icmp eq i64 %118, %10
%120 = add i64 %12, 1
br i1 %119, label %22, label %11
121: ; preds = %3, %23, %22
store float 1.000000e+00, float* %1
ret void
}