#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 fabsf(float);
void levinson_durbin(float R[], float lpcs[], int order) {
  float a[order + 1][order + 1];
  float sum, e, k;
  int i, j;
  e = 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 = -1.0 * (R[i] + sum) / e;
    if (fabsf(k) > 1.0)
      k = 0.0;
    a[i][i] = k;
    for (j = 1; j <= i - 1; j++)
      a[i][j] = a[i - 1][j] + k * a[i - 1][i - j];
    e *= (1 - k * k);
  }
  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    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %rax,-0x58(%rbp)
mov    %rax,%rcx
mov    %rsp,%rdx
mov    %rdx,-0x20(%rbp)
mov    %rax,%rdx
imul   %rdx,%rdx
lea    0xf(,%rdx,4),%rsi
and    $0xfffffffffffffff0,%rsi
mov    %rsp,%rdx
sub    %rsi,%rdx
mov    %rdx,-0x50(%rbp)
mov    %rdx,%rsp
mov    %rcx,-0x28(%rbp)
mov    %rax,-0x30(%rbp)
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x38(%rbp)
movl   $0x1,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     134e <levinson_durbin+0x23e>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x34(%rbp)
movl   $0x1,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     11f6 <levinson_durbin+0xe6>
mov    -0x50(%rbp),%rax
mov    -0x58(%rbp),%rdx
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
imul   %rdx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
movslq -0x44(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    -0x44(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x34(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    118e <levinson_durbin+0x7e>
mov    -0x8(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x34(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xdfc(%rip),%xmm0
mulsd  %xmm1,%xmm0
movss  -0x38(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x3c(%rbp),%xmm0
call   1030 <fabsf@plt>
movsd  0xdc8(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1252 <levinson_durbin+0x142>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3c(%rbp)
mov    -0x50(%rbp),%rax
mov    -0x58(%rbp),%rdx
movss  -0x3c(%rbp),%xmm0
movslq -0x40(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
movslq -0x40(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movl   $0x1,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     130f <levinson_durbin+0x1ff>
mov    -0x50(%rbp),%rax
mov    -0x58(%rbp),%rdx
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rsi
imul   %rdx,%rsi
shl    $0x2,%rsi
mov    %rax,%rcx
add    %rsi,%rcx
movslq -0x44(%rbp),%rsi
movss  (%rcx,%rsi,4),%xmm1
movss  -0x3c(%rbp),%xmm0
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rsi
imul   %rdx,%rsi
shl    $0x2,%rsi
mov    %rax,%rcx
add    %rsi,%rcx
mov    -0x40(%rbp),%esi
sub    -0x44(%rbp),%esi
movslq %esi,%rsi
movss  (%rcx,%rsi,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movslq -0x40(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
movslq -0x44(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    127e <levinson_durbin+0x16e>
movss  -0x3c(%rbp),%xmm0
movss  -0x3c(%rbp),%xmm1
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm1,%xmm0
movss  0xcce(%rip),%xmm1
addss  %xmm1,%xmm0
mulss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1173 <levinson_durbin+0x63>
movl   $0x1,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     139c <levinson_durbin+0x28c>
mov    -0x50(%rbp),%rax
mov    -0x58(%rbp),%rdx
movslq -0x14(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1355 <levinson_durbin+0x245>
mov    -0x10(%rbp),%rax
movss  0xc58(%rip),%xmm0
movss  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
mov    %rax,%rsp
mov    %rbp,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<levinson_durbin>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x98,%rsp
mov    %rsi,%rbx
lea    0x1(%rdx),%r11d
mov    %r11,%rax
imul   %r11,%rax
mov    %rsp,%r10
lea    0xf(,%rax,4),%rax
and    $0xfffffffffffffff0,%rax
sub    %rax,%r10
mov    %r10,%rsp
test   %edx,%edx
jle    1599 <levinson_durbin+0x479>
mov    %rdx,-0x78(%rbp)
mov    %rbx,-0x80(%rbp)
movss  (%rdi),%xmm5
lea    0x0(,%r11,4),%rax
lea    0x4(,%r11,4),%rcx
mov    %rcx,-0x70(%rbp)
lea    0x1(%r11),%rcx
mov    %rcx,-0x98(%rbp)
lea    0x4(%r10,%r11,4),%r12
lea    0x4(%r10),%rcx
mov    %rcx,-0x88(%rbp)
lea    0x10(%r10),%r8
lea    0x4(%rdi),%r9
lea    0x14(%r10),%rcx
mov    %rcx,-0x58(%rbp)
lea    -0xc(%r10),%rcx
mov    %rcx,-0x50(%rbp)
mov    $0x1,%r13d
xor    %ebx,%ebx
mov    %r12,-0x90(%rbp)
mov    %r10,-0x68(%rbp)
mov    %r10,%r15
mov    %rdi,%r14
mov    %r11,-0xa8(%rbp)
mov    %rdi,-0xa0(%rbp)
mov    %r10,-0x38(%rbp)
mov    %rax,-0x60(%rbp)
jmp    1253 <levinson_durbin+0x133>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mulss  %xmm6,%xmm6
movss  0xe1c(%rip),%xmm0
subss  %xmm6,%xmm0
mulss  %xmm0,%xmm5
add    $0x1,%r13
mov    -0x30(%rbp),%rbx
add    $0x1,%rbx
add    $0x4,%r14
mov    -0x60(%rbp),%rax
mov    -0xc0(%rbp),%r8
add    %rax,%r8
mov    -0xb8(%rbp),%r9
add    $0x4,%r9
add    %rax,%r15
add    %rax,-0x58(%rbp)
mov    -0x70(%rbp),%rcx
add    %rcx,-0x50(%rbp)
add    %rcx,-0x68(%rbp)
add    %rax,%r12
cmp    %r11,%r13
mov    -0xa0(%rbp),%rdi
je     1572 <levinson_durbin+0x452>
lea    -0x1(%rbx),%r10
xorps  %xmm0,%xmm0
cmp    $0x2,%r13
jb     133a <levinson_durbin+0x21a>
cmp    $0x3,%r10
jae    1280 <levinson_durbin+0x160>
xorps  %xmm0,%xmm0
mov    $0x1,%eax
test   $0x3,%bl
jne    12fc <levinson_durbin+0x1dc>
jmp    133a <levinson_durbin+0x21a>
mov    %rbx,%rcx
and    $0xfffffffffffffffc,%rcx
xorps  %xmm0,%xmm0
mov    $0xffffffffffffffff,%rax
mov    %r8,%rdx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  -0xc(%rdx),%xmm1
movss  -0x8(%rdx),%xmm2
mulss  0x4(%r14,%rax,4),%xmm1
addss  %xmm0,%xmm1
mulss  (%r14,%rax,4),%xmm2
addss  %xmm1,%xmm2
movss  -0x4(%rdx),%xmm1
mulss  -0x4(%r14,%rax,4),%xmm1
addss  %xmm2,%xmm1
movss  (%rdx),%xmm0
mulss  -0x8(%r14,%rax,4),%xmm0
addss  %xmm1,%xmm0
add    $0x10,%rdx
lea    (%rcx,%rax,1),%rsi
add    $0xfffffffffffffffc,%rsi
add    $0xfffffffffffffffc,%rax
cmp    $0xffffffffffffffff,%rsi
jne    12a0 <levinson_durbin+0x180>
neg    %rax
test   $0x3,%bl
je     133a <levinson_durbin+0x21a>
mov    %ebx,%ecx
and    $0x3,%ecx
lea    0x0(,%rax,4),%rsi
mov    %r9,%rdx
sub    %rsi,%rdx
lea    (%r15,%rax,4),%rax
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nop
movss  (%rax,%rsi,4),%xmm1
mulss  (%rdx),%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rsi
add    $0xfffffffffffffffc,%rdx
cmp    %rsi,%rcx
jne    1320 <levinson_durbin+0x200>
mov    %r10,-0xb0(%rbp)
mov    %rbx,-0x30(%rbp)
mov    %r9,-0xb8(%rbp)
mov    %r8,-0xc0(%rbp)
addss  (%rdi,%r13,4),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  0xc9c(%rip),%xmm0
movss  %xmm5,-0x44(%rbp)
xorps  %xmm1,%xmm1
cvtss2sd %xmm5,%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x40(%rbp)
call   1030 <fabsf@plt>
xorps  %xmm6,%xmm6
ucomisd 0xc83(%rip),%xmm0
ja     1394 <levinson_durbin+0x274>
movss  -0x40(%rbp),%xmm6
mov    %r13,%rax
mov    -0xa8(%rbp),%r11
imul   %r11,%rax
mov    -0x38(%rbp),%r10
lea    (%r10,%rax,4),%r9
movss  %xmm6,(%r9,%r13,4)
cmp    $0x2,%r13
movss  -0x44(%rbp),%xmm5
jb     11f0 <levinson_durbin+0xd0>
mov    $0x1,%esi
mov    -0x30(%rbp),%rax
cmp    $0x8,%rax
jae    1460 <levinson_durbin+0x340>
mov    %rsi,%rax
test   $0x1,%r13b
jne    1405 <levinson_durbin+0x2e5>
lea    -0x1(%r13),%rax
imul   %r11,%rax
lea    (%r10,%rax,4),%rax
mov    %r13,%rcx
sub    %rsi,%rcx
movss  (%rax,%rcx,4),%xmm0
mulss  %xmm6,%xmm0
addss  (%rax,%rsi,4),%xmm0
movss  %xmm0,(%r9,%rsi,4)
lea    0x1(%rsi),%rax
cmp    %rsi,-0x30(%rbp)
je     11f0 <levinson_durbin+0xd0>
lea    0x0(,%rax,4),%rdx
mov    -0x68(%rbp),%rcx
sub    %rdx,%rcx
xchg   %ax,%ax
movss  0x4(%rcx),%xmm0
mulss  %xmm6,%xmm0
addss  (%r15,%rax,4),%xmm0
movss  %xmm0,-0x4(%r12,%rax,4)
movss  (%rcx),%xmm0
mulss  %xmm6,%xmm0
addss  0x4(%r15,%rax,4),%xmm0
movss  %xmm0,(%r12,%rax,4)
add    $0x2,%rax
add    $0xfffffffffffffff8,%rcx
cmp    %rax,%r13
jne    1420 <levinson_durbin+0x300>
jmp    11f0 <levinson_durbin+0xd0>
nopl   (%rax)
mov    -0x70(%rbp),%rdi
imul   %rax,%rdi
mov    -0xb0(%rbp),%rax
mov    $0x4,%ecx
mul    %rcx
lea    (%r10,%rdi,1),%rcx
seto   %r8b
mov    %rcx,%rdx
sub    %rax,%rdx
cmp    %rcx,%rdx
ja     13d2 <levinson_durbin+0x2b2>
test   %r8b,%r8b
jne    13d2 <levinson_durbin+0x2b2>
mov    -0x30(%rbp),%rcx
mov    %rcx,%rax
imul   %r11,%rax
mov    %rcx,%rbx
and    $0xfffffffffffffff8,%rbx
mov    -0x98(%rbp),%rcx
add    %rax,%rcx
mov    %rbx,%rdx
neg    %rdx
lea    (%r10,%rcx,4),%r8
mov    -0x90(%rbp),%rcx
add    %rdi,%rcx
mov    %rcx,-0x40(%rbp)
lea    (%r10,%rax,4),%rax
add    $0x4,%rax
add    -0x88(%rbp),%rdi
cmp    %r8,%rdi
seta   %cl
cmp    %rax,-0x40(%rbp)
seta   %al
and    %cl,%al
test   %al,%al
jne    1569 <levinson_durbin+0x449>
mov    -0x38(%rbp),%r10
jne    13d2 <levinson_durbin+0x2b2>
mov    %rbx,%r8
lea    0x1(%rbx),%rsi
movaps %xmm6,%xmm0
shufps $0x0,%xmm6,%xmm0
xor    %edi,%edi
mov    -0x58(%rbp),%rax
mov    -0x60(%rbp),%rcx
mov    -0x50(%rbp),%rbx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movups -0x10(%rax),%xmm1
movups (%rax),%xmm2
movups -0x10(%rbx,%rdi,4),%xmm3
movups (%rbx,%rdi,4),%xmm4
shufps $0x1b,%xmm4,%xmm4
shufps $0x1b,%xmm3,%xmm3
mulps  %xmm0,%xmm4
addps  %xmm1,%xmm4
mulps  %xmm0,%xmm3
addps  %xmm2,%xmm3
movups %xmm4,-0x10(%rax,%rcx,1)
movups %xmm3,(%rax,%rcx,1)
add    $0x20,%rax
add    $0xfffffffffffffff8,%rdi
cmp    %rdi,%rdx
jne    1520 <levinson_durbin+0x400>
cmp    %r8,-0x30(%rbp)
je     11f0 <levinson_durbin+0xd0>
jmp    13d2 <levinson_durbin+0x2b2>
mov    -0x38(%rbp),%r10
jmp    13d2 <levinson_durbin+0x2b2>
mov    -0x78(%rbp),%rax
test   %eax,%eax
mov    -0x80(%rbp),%rbx
jle    1599 <levinson_durbin+0x479>
mov    %eax,%edx
imul   %rdx,%r11
lea    0x4(%rbx),%rdi
lea    (%r10,%r11,4),%rsi
add    $0x4,%rsi
shl    $0x2,%rdx
call   1040 <memcpy@plt>
movl   $0x3f800000,(%rbx)
lea    -0x28(%rbp),%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 i8*
%8 = alloca i64
%9 = alloca i64
%10 = alloca float
%11 = alloca float
%12 = alloca float
%13 = alloca i32
%14 = alloca i32
store float* %0, float** %4
store float* %1, float** %5
store i32 %2, i32* %6
%15 = load i32, i32* %6
%16 = add nsw i32 %15, 1
%17 = zext i32 %16 to i64
%18 = load i32, i32* %6
%19 = add nsw i32 %18, 1
%20 = zext i32 %19 to i64
%21 = call i8* @llvm.stacksave()
store i8* %21, i8** %7
%22 = mul nuw i64 %17, %20
%23 = alloca float, i64 %22
store i64 %17, i64* %8
store i64 %20, i64* %9
%24 = load float*, float** %4
%25 = getelementptr inbounds float, float* %24, i64 0
%26 = load float, float* %25
store float %26, float* %11
store i32 1, i32* %13
br label %27
27: ; preds = %131, %3
%28 = load i32, i32* %13
%29 = load i32, i32* %6
%30 = icmp sle i32 %28, %29
br i1 %30, label %31, label %134
31: ; preds = %27
store float 0.000000e+00, float* %10
store i32 1, i32* %14
br label %32
32: ; preds = %56, %31
%33 = load i32, i32* %14
%34 = load i32, i32* %13
%35 = sub nsw i32 %34, 1
%36 = icmp sle i32 %33, %35
br i1 %36, label %37, label %59
37: ; preds = %32
%38 = load i32, i32* %13
%39 = sub nsw i32 %38, 1
%40 = sext i32 %39 to i64
%41 = mul nsw i64 %40, %20
%42 = getelementptr inbounds float, float* %23, i64 %41
%43 = load i32, i32* %14
%44 = sext i32 %43 to i64
%45 = getelementptr inbounds float, float* %42, i64 %44
%46 = load float, float* %45
%47 = load float*, float** %4
%48 = load i32, i32* %13
%49 = load i32, i32* %14
%50 = sub nsw i32 %48, %49
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds float, float* %47, i64 %51
%53 = load float, float* %52
%54 = load float, float* %10
%55 = call float @llvm.fmuladd.f32(float %46, float %53, float %54)
store float %55, float* %10
br label %56
56: ; preds = %37
%57 = load i32, i32* %14
%58 = add nsw i32 %57, 1
store i32 %58, i32* %14
br label %32
59: ; preds = %32
%60 = load float*, float** %4
%61 = load i32, i32* %13
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds float, float* %60, i64 %62
%64 = load float, float* %63
%65 = load float, float* %10
%66 = fadd float %64, %65
%67 = fpext float %66 to double
%68 = fmul double -1.000000e+00, %67
%69 = load float, float* %11
%70 = fpext float %69 to double
%71 = fdiv double %68, %70
%72 = fptrunc double %71 to float
store float %72, float* %12
%73 = load float, float* %12
%74 = call double @fabsf(float %73)
%75 = fcmp ogt double %74, 1.000000e+00
br i1 %75, label %76, label %77
76: ; preds = %59
store float 0.000000e+00, float* %12
br label %77
77: ; preds = %76, %59
%78 = load float, float* %12
%79 = load i32, i32* %13
%80 = sext i32 %79 to i64
%81 = mul nsw i64 %80, %20
%82 = getelementptr inbounds float, float* %23, i64 %81
%83 = load i32, i32* %13
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds float, float* %82, i64 %84
store float %78, float* %85
store i32 1, i32* %14
br label %86
86: ; preds = %121, %77
%87 = load i32, i32* %14
%88 = load i32, i32* %13
%89 = sub nsw i32 %88, 1
%90 = icmp sle i32 %87, %89
br i1 %90, label %91, label %124
91: ; preds = %86
%92 = load i32, i32* %13
%93 = sub nsw i32 %92, 1
%94 = sext i32 %93 to i64
%95 = mul nsw i64 %94, %20
%96 = getelementptr inbounds float, float* %23, i64 %95
%97 = load i32, i32* %14
%98 = sext i32 %97 to i64
%99 = getelementptr inbounds float, float* %96, i64 %98
%100 = load float, float* %99
%101 = load float, float* %12
%102 = load i32, i32* %13
%103 = sub nsw i32 %102, 1
%104 = sext i32 %103 to i64
%105 = mul nsw i64 %104, %20
%106 = getelementptr inbounds float, float* %23, i64 %105
%107 = load i32, i32* %13
%108 = load i32, i32* %14
%109 = sub nsw i32 %107, %108
%110 = sext i32 %109 to i64
%111 = getelementptr inbounds float, float* %106, i64 %110
%112 = load float, float* %111
%113 = call float @llvm.fmuladd.f32(float %101, float %112, float %100)
%114 = load i32, i32* %13
%115 = sext i32 %114 to i64
%116 = mul nsw i64 %115, %20
%117 = getelementptr inbounds float, float* %23, i64 %116
%118 = load i32, i32* %14
%119 = sext i32 %118 to i64
%120 = getelementptr inbounds float, float* %117, i64 %119
store float %113, float* %120
br label %121
121: ; preds = %91
%122 = load i32, i32* %14
%123 = add nsw i32 %122, 1
store i32 %123, i32* %14
br label %86
124: ; preds = %86
%125 = load float, float* %12
%126 = load float, float* %12
%127 = fneg float %125
%128 = call float @llvm.fmuladd.f32(float %127, float %126, float 1.000000e+00)
%129 = load float, float* %11
%130 = fmul float %129, %128
store float %130, float* %11
br label %131
131: ; preds = %124
%132 = load i32, i32* %13
%133 = add nsw i32 %132, 1
store i32 %133, i32* %13
br label %27
134: ; preds = %27
store i32 1, i32* %13
br label %135
135: ; preds = %152, %134
%136 = load i32, i32* %13
%137 = load i32, i32* %6
%138 = icmp sle i32 %136, %137
br i1 %138, label %139, label %155
139: ; preds = %135
%140 = load i32, i32* %6
%141 = sext i32 %140 to i64
%142 = mul nsw i64 %141, %20
%143 = getelementptr inbounds float, float* %23, i64 %142
%144 = load i32, i32* %13
%145 = sext i32 %144 to i64
%146 = getelementptr inbounds float, float* %143, i64 %145
%147 = load float, float* %146
%148 = load float*, float** %5
%149 = load i32, i32* %13
%150 = sext i32 %149 to i64
%151 = getelementptr inbounds float, float* %148, i64 %150
store float %147, float* %151
br label %152
152: ; preds = %139
%153 = load i32, i32* %13
%154 = add nsw i32 %153, 1
store i32 %154, i32* %13
br label %135
155: ; preds = %135
%156 = load float*, float** %5
%157 = getelementptr inbounds float, float* %156, i64 0
store float 1.000000e+00, float* %157
%158 = load i8*, i8** %7
ret void
}

LLVM IR (-O2):
define void @levinson_durbin(float* nocapture %0, float* nocapture writeonly %1, i32 %2) {
%4 = add i32 %2, 1
%5 = zext i32 %4 to i64
%6 = mul nuw i64 %5, %5
%7 = alloca float, i64 %6
%8 = bitcast float* %7 to i8*
%9 = icmp slt i32 %2, 1
br i1 %9, label %226, label %10
10: ; preds = %3
%11 = load float, float* %0
%12 = shl nuw nsw i64 %5, 2
%13 = add nuw nsw i64 %12, 4
%14 = add nuw nsw i64 %5, 1
%15 = getelementptr float, float* %7, i64 %14
%16 = bitcast float* %15 to i8*
%17 = shl nuw nsw i64 %5, 2
%18 = add nuw nsw i64 %17, 4
%19 = getelementptr inbounds float, float* %7, i64 1
%20 = bitcast float* %19 to i8*
br label %21
21: ; preds = %10, %219
%22 = phi i64 [ 0, %10 ], [ %225, %219 ]
%23 = phi i64 [ 1, %10 ], [ %223, %219 ]
%24 = phi float [ %11, %10 ], [ %222, %219 ]
%25 = mul i64 %22, %5
%26 = add i64 %14, %25
%27 = getelementptr float, float* %7, i64 %26
%28 = bitcast float* %27 to i8*
%29 = mul i64 %18, %22
%30 = getelementptr i8, i8* %16, i64 %29
%31 = add i64 %25, 1
%32 = getelementptr float, float* %7, i64 %31
%33 = bitcast float* %32 to i8*
%34 = getelementptr i8, i8* %20, i64 %29
%35 = add i64 %22, -1
%36 = mul i64 %13, %22
%37 = getelementptr i8, i8* %8, i64 %36
%38 = add nsw i64 %23, -1
%39 = icmp ugt i64 %23, 1
br i1 %39, label %40, label %111
40: ; preds = %21
%41 = add i64 %22, -1
%42 = mul nuw nsw i64 %38, %5
%43 = getelementptr inbounds float, float* %7, i64 %42
%44 = and i64 %22, 3
%45 = icmp ult i64 %41, 3
br i1 %45, label %93, label %46
46: ; preds = %40
%47 = and i64 %22, -4
br label %59
48: ; preds = %219
br i1 %9, label %226, label %49
49: ; preds = %48
%50 = zext i32 %2 to i64
%51 = mul nuw nsw i64 %5, %50
%52 = getelementptr float, float* %1, i64 1
%53 = bitcast float* %52 to i8*
%54 = add nuw nsw i64 %51, 1
%55 = getelementptr float, float* %7, i64 %54
%56 = bitcast float* %55 to i8*
%57 = zext i32 %2 to i64
%58 = shl nuw nsw i64 %57, 2
br label %226
59: ; preds = %59, %46
%60 = phi i64 [ 1, %46 ], [ %90, %59 ]
%61 = phi float [ 0.000000e+00, %46 ], [ %89, %59 ]
%62 = phi i64 [ 0, %46 ], [ %91, %59 ]
%63 = getelementptr inbounds float, float* %43, i64 %60
%64 = load float, float* %63
%65 = sub nsw i64 %23, %60
%66 = getelementptr inbounds float, float* %0, i64 %65
%67 = load float, float* %66
%68 = tail call float @llvm.fmuladd.f32(float %64, float %67, float %61)
%69 = add nuw nsw i64 %60, 1
%70 = getelementptr inbounds float, float* %43, i64 %69
%71 = load float, float* %70
%72 = sub nsw i64 %23, %69
%73 = getelementptr inbounds float, float* %0, i64 %72
%74 = load float, float* %73
%75 = tail call float @llvm.fmuladd.f32(float %71, float %74, float %68)
%76 = add nuw nsw i64 %60, 2
%77 = getelementptr inbounds float, float* %43, i64 %76
%78 = load float, float* %77
%79 = sub nsw i64 %23, %76
%80 = getelementptr inbounds float, float* %0, i64 %79
%81 = load float, float* %80
%82 = tail call float @llvm.fmuladd.f32(float %78, float %81, float %75)
%83 = add nuw nsw i64 %60, 3
%84 = getelementptr inbounds float, float* %43, i64 %83
%85 = load float, float* %84
%86 = sub nsw i64 %23, %83
%87 = getelementptr inbounds float, float* %0, i64 %86
%88 = load float, float* %87
%89 = tail call float @llvm.fmuladd.f32(float %85, float %88, float %82)
%90 = add nuw nsw i64 %60, 4
%91 = add nuw i64 %62, 4
%92 = icmp eq i64 %91, %47
br i1 %92, label %93, label %59
93: ; preds = %59, %40
%94 = phi float [ undef, %40 ], [ %89, %59 ]
%95 = phi i64 [ 1, %40 ], [ %90, %59 ]
%96 = phi float [ 0.000000e+00, %40 ], [ %89, %59 ]
%97 = icmp eq i64 %44, 0
br i1 %97, label %111, label %98
98: ; preds = %93, %98
%99 = phi i64 [ %108, %98 ], [ %95, %93 ]
%100 = phi float [ %107, %98 ], [ %96, %93 ]
%101 = phi i64 [ %109, %98 ], [ 0, %93 ]
%102 = getelementptr inbounds float, float* %43, i64 %99
%103 = load float, float* %102
%104 = sub nsw i64 %23, %99
%105 = getelementptr inbounds float, float* %0, i64 %104
%106 = load float, float* %105
%107 = tail call float @llvm.fmuladd.f32(float %103, float %106, float %100)
%108 = add nuw nsw i64 %99, 1
%109 = add i64 %101, 1
%110 = icmp eq i64 %109, %44
br i1 %110, label %111, label %98
111: ; preds = %93, %98, %21
%112 = phi float [ 0.000000e+00, %21 ], [ %94, %93 ], [ %107, %98 ]
%113 = getelementptr inbounds float, float* %0, i64 %23
%114 = load float, float* %113
%115 = fadd float %112, %114
%116 = fpext float %115 to double
%117 = fneg double %116
%118 = fpext float %24 to double
%119 = fdiv double %117, %118
%120 = fptrunc double %119 to float
%121 = tail call double @fabsf(float %120)
%122 = fcmp ogt double %121, 1.000000e+00
%123 = select i1 %122, float 0.000000e+00, float %120
%124 = mul nuw nsw i64 %23, %5
%125 = getelementptr inbounds float, float* %7, i64 %124
%126 = getelementptr inbounds float, float* %125, i64 %23
store float %123, float* %126
br i1 %39, label %127, label %219
127: ; preds = %111
%128 = mul nuw nsw i64 %38, %5
%129 = getelementptr inbounds float, float* %7, i64 %128
%130 = icmp ult i64 %22, 8
br i1 %130, label %184, label %131
131: ; preds = %127
%132 = call { i64, i1 } @llvm.umul.with.overflow.i64(i64 %35, i64 4)
%133 = extractvalue { i64, i1 } %132, 0
%134 = extractvalue { i64, i1 } %132, 1
%135 = sub i64 0, %133
%136 = getelementptr i8, i8* %37, i64 %135
%137 = icmp ugt i8* %136, %37
%138 = or i1 %137, %134
br i1 %138, label %184, label %139
139: ; preds = %131
%140 = icmp ugt i8* %34, %28
%141 = icmp ugt i8* %30, %33
%142 = and i1 %140, %141
%143 = icmp ugt i8* %34, %28
%144 = icmp ugt i8* %30, %33
%145 = and i1 %143, %144
%146 = or i1 %142, %145
br i1 %146, label %184, label %147
147: ; preds = %139
%148 = and i64 %22, -8
%149 = or i64 %148, 1
%150 = insertelement <4 x float> poison, float %123, i64 0
%151 = shufflevector <4 x float> %150, <4 x float> poison, <4 x i32> zeroinitializer
%152 = insertelement <4 x float> poison, float %123, i64 0
%153 = shufflevector <4 x float> %152, <4 x float> poison, <4 x i32> zeroinitializer
br label %154
154: ; preds = %154, %147
%155 = phi i64 [ 0, %147 ], [ %180, %154 ]
%156 = or i64 %155, 1
%157 = getelementptr inbounds float, float* %129, i64 %156
%158 = bitcast float* %157 to <4 x float>*
%159 = load <4 x float>, <4 x float>* %158
%160 = getelementptr inbounds float, float* %157, i64 4
%161 = bitcast float* %160 to <4 x float>*
%162 = load <4 x float>, <4 x float>* %161
%163 = sub nsw i64 %23, %156
%164 = getelementptr inbounds float, float* %129, i64 %163
%165 = getelementptr inbounds float, float* %164, i64 -3
%166 = bitcast float* %165 to <4 x float>*
%167 = load <4 x float>, <4 x float>* %166
%168 = shufflevector <4 x float> %167, <4 x float> poison, <4 x i32> <i32 3, i32 2, i32 1, i32 0>
%169 = getelementptr inbounds float, float* %164, i64 -4
%170 = getelementptr inbounds float, float* %169, i64 -3
%171 = bitcast float* %170 to <4 x float>*
%172 = load <4 x float>, <4 x float>* %171
%173 = shufflevector <4 x float> %172, <4 x float> poison, <4 x i32> <i32 3, i32 2, i32 1, i32 0>
%174 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %151, <4 x float> %168, <4 x float> %159)
%175 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %153, <4 x float> %173, <4 x float> %162)
%176 = getelementptr inbounds float, float* %125, i64 %156
%177 = bitcast float* %176 to <4 x float>*
store <4 x float> %174, <4 x float>* %177
%178 = getelementptr inbounds float, float* %176, i64 4
%179 = bitcast float* %178 to <4 x float>*
store <4 x float> %175, <4 x float>* %179
%180 = add nuw i64 %155, 8
%181 = icmp eq i64 %180, %148
br i1 %181, label %182, label %154
182: ; preds = %154
%183 = icmp eq i64 %22, %148
br i1 %183, label %219, label %184
184: ; preds = %139, %131, %127, %182
%185 = phi i64 [ 1, %139 ], [ 1, %131 ], [ 1, %127 ], [ %149, %182 ]
%186 = and i64 %23, 1
%187 = icmp eq i64 %186, 0
br i1 %187, label %188, label %197
188: ; preds = %184
%189 = getelementptr inbounds float, float* %129, i64 %185
%190 = load float, float* %189
%191 = sub nsw i64 %23, %185
%192 = getelementptr inbounds float, float* %129, i64 %191
%193 = load float, float* %192
%194 = tail call float @llvm.fmuladd.f32(float %123, float %193, float %190)
%195 = getelementptr inbounds float, float* %125, i64 %185
store float %194, float* %195
%196 = add nuw nsw i64 %185, 1
br label %197
197: ; preds = %188, %184
%198 = phi i64 [ %185, %184 ], [ %196, %188 ]
%199 = icmp eq i64 %22, %185
br i1 %199, label %219, label %200
200: ; preds = %197, %200
%201 = phi i64 [ %217, %200 ], [ %198, %197 ]
%202 = getelementptr inbounds float, float* %129, i64 %201
%203 = load float, float* %202
%204 = sub nsw i64 %23, %201
%205 = getelementptr inbounds float, float* %129, i64 %204
%206 = load float, float* %205
%207 = tail call float @llvm.fmuladd.f32(float %123, float %206, float %203)
%208 = getelementptr inbounds float, float* %125, i64 %201
store float %207, float* %208
%209 = add nuw nsw i64 %201, 1
%210 = getelementptr inbounds float, float* %129, i64 %209
%211 = load float, float* %210
%212 = sub nsw i64 %23, %209
%213 = getelementptr inbounds float, float* %129, i64 %212
%214 = load float, float* %213
%215 = tail call float @llvm.fmuladd.f32(float %123, float %214, float %211)
%216 = getelementptr inbounds float, float* %125, i64 %209
store float %215, float* %216
%217 = add nuw nsw i64 %201, 2
%218 = icmp eq i64 %217, %23
br i1 %218, label %219, label %200
219: ; preds = %197, %200, %182, %111
%220 = fneg float %123
%221 = tail call float @llvm.fmuladd.f32(float %220, float %123, float 1.000000e+00)
%222 = fmul float %24, %221
%223 = add nuw nsw i64 %23, 1
%224 = icmp eq i64 %223, %5
%225 = add i64 %22, 1
br i1 %224, label %48, label %21
226: ; preds = %3, %49, %48
store float 1.000000e+00, float* %1
ret void
}