void lev_durb(double *corr, double *lpc_coef) {
  double k[11], tab[11];
  double err, acc;
  int i, j;
  double *a = tab;
  double *prev_a = lpc_coef;
  double *exch;
  for (i = 0; i < 11; i++) {
    prev_a[i] = 0;
    a[i] = 0;
  };
  err = corr[0];
  for (i = 1; i < 11; i++) {
    prev_a[0] = 1;
    acc = 0;
    for (j = 0; j < i; j++) {
      acc = acc + prev_a[j] * corr[i - j];
    };
    a[i] = k[i] = -acc / err;
    for (j = 1; j < i; j++) {
      a[j] = prev_a[j] + k[i] * prev_a[i - j];
    };
    err = (1 - k[i] * k[i]) * err;
    exch = prev_a;
    prev_a = a;
    a = exch;
  };
}


x86 assembly (clang, -O0):
<lev_durb>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
lea    -0xd0(%rbp),%rax
mov    %rax,-0xf0(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0xf8(%rbp)
movl   $0x0,-0xe4(%rbp)
cmpl   $0xb,-0xe4(%rbp)
jge    1183 <lev_durb+0x83>
mov    -0xf8(%rbp),%rax
movslq -0xe4(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xf0(%rbp),%rax
movslq -0xe4(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xe4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xe4(%rbp)
jmp    1136 <lev_durb+0x36>
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
movl   $0x1,-0xe4(%rbp)
cmpl   $0xb,-0xe4(%rbp)
jge    1397 <lev_durb+0x297>
mov    -0xf8(%rbp),%rax
movsd  0xe47(%rip),%xmm0
movsd  %xmm0,(%rax)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
movl   $0x0,-0xe8(%rbp)
mov    -0xe8(%rbp),%eax
cmp    -0xe4(%rbp),%eax
jge    123b <lev_durb+0x13b>
movsd  -0xe0(%rbp),%xmm1
ff
mov    -0xf8(%rbp),%rax
movslq -0xe8(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0xe4(%rbp),%ecx
sub    -0xe8(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
mov    -0xe8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xe8(%rbp)
jmp    11d2 <lev_durb+0xd2>
movsd  -0xe0(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
divsd  -0xd8(%rbp),%xmm0
ff
movslq -0xe4(%rbp),%rax
movsd  %xmm0,-0x70(%rbp,%rax,8)
mov    -0xf0(%rbp),%rax
movslq -0xe4(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x1,-0xe8(%rbp)
mov    -0xe8(%rbp),%eax
cmp    -0xe4(%rbp),%eax
jge    1308 <lev_durb+0x208>
mov    -0xf8(%rbp),%rax
movslq -0xe8(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movslq -0xe4(%rbp),%rax
movsd  -0x70(%rbp,%rax,8),%xmm0
mov    -0xf8(%rbp),%rax
mov    -0xe4(%rbp),%ecx
sub    -0xe8(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0xf0(%rbp),%rax
movslq -0xe8(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xe8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xe8(%rbp)
jmp    128c <lev_durb+0x18c>
movslq -0xe4(%rbp),%rax
movsd  -0x70(%rbp,%rax,8),%xmm0
movslq -0xe4(%rbp),%rax
movsd  -0x70(%rbp,%rax,8),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm1,%xmm0
movsd  0xcbb(%rip),%xmm1
addsd  %xmm1,%xmm0
mulsd  -0xd8(%rbp),%xmm0
ff
movsd  %xmm0,-0xd8(%rbp)
ff
mov    -0xf8(%rbp),%rax
mov    %rax,-0x100(%rbp)
mov    -0xf0(%rbp),%rax
mov    %rax,-0xf8(%rbp)
mov    -0x100(%rbp),%rax
mov    %rax,-0xf0(%rbp)
mov    -0xe4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xe4(%rbp)
jmp    119d <lev_durb+0x9d>
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<lev_durb>:
push   %r15
push   %r14
push   %rbx
sub    $0x40,%rsp
xorps  %xmm0,%xmm0
movups %xmm0,0x40(%rsi)
movups %xmm0,0x30(%rsi)
movups %xmm0,0x20(%rsi)
movups %xmm0,0x10(%rsi)
movups %xmm0,(%rsi)
movq   $0x0,0x50(%rsi)
movaps %xmm0,-0x40(%rsp)
movaps %xmm0,-0x50(%rsp)
movaps %xmm0,-0x60(%rsp)
movaps %xmm0,-0x70(%rsp)
movaps %xmm0,-0x80(%rsp)
movq   $0x0,-0x30(%rsp)
lea    -0x80(%rsp),%r11
movsd  (%rdi),%xmm0
lea    0x8(%rdi),%r8
mov    $0x1,%r15d
xor    %r10d,%r10d
movabs $0x3ff0000000000000,%r9
xorpd  %xmm1,%xmm1
movapd 0xe8b(%rip),%xmm2
movsd  0xe93(%rip),%xmm3
mov    %rdi,%rcx
jmp    11c4 <lev_durb+0xc4>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  -0x20(%rsp,%r15,8),%xmm4
mulsd  %xmm4,%xmm4
movapd %xmm3,%xmm5
subsd  %xmm4,%xmm5
mulsd  %xmm5,%xmm0
add    $0x1,%r15
add    $0x1,%r10
add    $0x8,%rcx
add    $0x8,%r8
mov    %r14,%rsi
cmp    $0xb,%r15
je     1365 <lev_durb+0x265>
mov    %r11,%r14
mov    %rsi,%r11
mov    %r9,(%rsi)
movsd  (%rdi,%r15,8),%xmm4
addsd  %xmm1,%xmm4
cmp    $0x1,%r15
je     12aa <lev_durb+0x1aa>
lea    -0x1(%r10),%rax
mov    $0x1,%esi
cmp    $0x3,%rax
jb     1263 <lev_durb+0x163>
mov    %r10,%rbx
and    $0xfffffffffffffffc,%rbx
lea    0x20(%r11),%rax
mov    $0xffffffffffffffff,%rsi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  -0x18(%rax),%xmm5
movsd  -0x10(%rax),%xmm6
mulsd  0x8(%rcx,%rsi,8),%xmm5
addsd  %xmm4,%xmm5
mulsd  (%rcx,%rsi,8),%xmm6
addsd  %xmm5,%xmm6
movsd  -0x8(%rax),%xmm5
mulsd  -0x8(%rcx,%rsi,8),%xmm5
addsd  %xmm6,%xmm5
movsd  (%rax),%xmm4
mulsd  -0x10(%rcx,%rsi,8),%xmm4
addsd  %xmm5,%xmm4
add    $0x20,%rax
lea    (%rbx,%rsi,1),%rdx
add    $0xfffffffffffffffc,%rdx
add    $0xfffffffffffffffc,%rsi
cmp    $0xffffffffffffffff,%rdx
jne    1210 <lev_durb+0x110>
neg    %rsi
test   $0x3,%r10b
je     12aa <lev_durb+0x1aa>
mov    %r10d,%eax
and    $0x3,%eax
lea    0x0(,%rsi,8),%rdx
mov    %r8,%rbx
sub    %rdx,%rbx
lea    (%r11,%rsi,8),%rsi
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rsi,%rdx,8),%xmm5
mulsd  (%rbx),%xmm5
addsd  %xmm5,%xmm4
add    $0x1,%rdx
add    $0xfffffffffffffff8,%rbx
cmp    %rdx,%rax
jne    1290 <lev_durb+0x190>
xorpd  %xmm2,%xmm4
divsd  %xmm0,%xmm4
movsd  %xmm4,-0x20(%rsp,%r15,8)
movsd  %xmm4,(%r14,%r15,8)
cmp    $0x2,%r15
jb     1190 <lev_durb+0x90>
mov    $0x1,%esi
cmp    $0x1,%r10
je     1337 <lev_durb+0x237>
mov    %r10,%rax
and    $0xfffffffffffffffe,%rax
neg    %rax
lea    (%r11,%r10,8),%rbx
mov    $0x1,%esi
nopw   0x0(%rax,%rax,1)
movsd  -0x20(%rsp,%r15,8),%xmm4
mulsd  (%rbx),%xmm4
addsd  (%r11,%rsi,8),%xmm4
movsd  %xmm4,(%r14,%rsi,8)
movsd  -0x20(%rsp,%r15,8),%xmm4
mulsd  -0x8(%rbx),%xmm4
addsd  0x8(%r11,%rsi,8),%xmm4
movsd  %xmm4,0x8(%r14,%rsi,8)
add    $0xfffffffffffffff0,%rbx
lea    (%rax,%rsi,1),%rdx
add    $0x2,%rdx
add    $0x2,%rsi
cmp    $0x1,%rdx
jne    12f0 <lev_durb+0x1f0>
test   $0x1,%r10b
je     1190 <lev_durb+0x90>
movsd  -0x20(%rsp,%r15,8),%xmm4
mov    %r15,%rax
sub    %rsi,%rax
mulsd  (%r11,%rax,8),%xmm4
addsd  (%r11,%rsi,8),%xmm4
movsd  %xmm4,(%r14,%rsi,8)
jmp    1190 <lev_durb+0x90>
add    $0x40,%rsp
pop    %rbx
pop    %r14
pop    %r15
ret

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

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