float silk_schur_FLP(float refl_coef[], const float auto_corr[], int order) {
  int k, n;
  double C[24 + 1][2];
  double Ctmp1, Ctmp2, rc_tmp;
  ;
  k = 0;
  do {
    C[k][0] = C[k][1] = auto_corr[k];
  } while (++k <= order);
  for (k = 0; k < order; k++) {
    rc_tmp = -C[k + 1][0] / (((C[0][1]) > (1e-9f)) ? (C[0][1]) : (1e-9f));
    refl_coef[k] = (float)rc_tmp;
    for (n = 0; n < order - k; n++) {
      Ctmp1 = C[n + k + 1][0];
      Ctmp2 = C[n][1];
      C[n + k + 1][0] = Ctmp1 + Ctmp2 * rc_tmp;
      C[n][1] = Ctmp2 + Ctmp1 * rc_tmp;
    }
  }
  return (float)C[0][1];
}


x86 assembly (clang, -O0):
<silk_schur_FLP>:
push   %rbp
mov    %rsp,%rbp
sub    $0x160,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movslq -0x18(%rbp),%rcx
lea    -0x1b0(%rbp),%rax
shl    $0x4,%rcx
add    %rcx,%rax
movsd  %xmm0,0x8(%rax)
movslq -0x18(%rbp),%rcx
lea    -0x1b0(%rbp),%rax
shl    $0x4,%rcx
add    %rcx,%rax
movsd  %xmm0,(%rax)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
cmp    -0x14(%rbp),%eax
jle    111d <silk_schur_FLP+0x1d>
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1326 <silk_schur_FLP+0x226>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
movslq %eax,%rcx
lea    -0x1b0(%rbp),%rax
shl    $0x4,%rcx
add    %rcx,%rax
movsd  (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x1d0(%rbp)
ff
movsd  -0x1a8(%rbp),%xmm0
ff
movsd  0xe36(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    11e9 <silk_schur_FLP+0xe9>
movsd  -0x1a8(%rbp),%xmm0
ff
movsd  %xmm0,-0x1d8(%rbp)
ff
jmp    11fe <silk_schur_FLP+0xfe>
movsd  0xe0f(%rip),%xmm0
movsd  %xmm0,-0x1d8(%rbp)
ff
jmp    11fe <silk_schur_FLP+0xfe>
movsd  -0x1d0(%rbp),%xmm0
ff
movsd  -0x1d8(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x1c8(%rbp)
ff
movsd  -0x1c8(%rbp),%xmm0
ff
cvtsd2ss %xmm0,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    -0x18(%rbp),%ecx
cmp    %ecx,%eax
jge    1313 <silk_schur_FLP+0x213>
mov    -0x1c(%rbp),%eax
add    -0x18(%rbp),%eax
add    $0x1,%eax
movslq %eax,%rcx
lea    -0x1b0(%rbp),%rax
shl    $0x4,%rcx
add    %rcx,%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x1b8(%rbp)
ff
movslq -0x1c(%rbp),%rcx
lea    -0x1b0(%rbp),%rax
shl    $0x4,%rcx
add    %rcx,%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x1c0(%rbp)
ff
movsd  -0x1b8(%rbp),%xmm1
ff
movsd  -0x1c0(%rbp),%xmm0
ff
movsd  -0x1c8(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x1c(%rbp),%eax
add    -0x18(%rbp),%eax
add    $0x1,%eax
movslq %eax,%rcx
lea    -0x1b0(%rbp),%rax
shl    $0x4,%rcx
add    %rcx,%rax
movsd  %xmm0,(%rax)
movsd  -0x1c0(%rbp),%xmm1
ff
movsd  -0x1b8(%rbp),%xmm0
ff
movsd  -0x1c8(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movslq -0x1c(%rbp),%rcx
lea    -0x1b0(%rbp),%rax
shl    $0x4,%rcx
add    %rcx,%rax
movsd  %xmm0,0x8(%rax)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    123a <silk_schur_FLP+0x13a>
jmp    1318 <silk_schur_FLP+0x218>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1174 <silk_schur_FLP+0x74>
movsd  -0x1a8(%rbp),%xmm0
ff
cvtsd2ss %xmm0,%xmm0
add    $0x160,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<silk_schur_FLP>:
sub    $0x118,%rsp
xor    %r9d,%r9d
test   %edx,%edx
mov    $0x0,%eax
cmovg  %edx,%eax
lea    0x1(%rax),%r8d
cmp    $0x3,%eax
jb     1171 <silk_schur_FLP+0x71>
mov    %r8d,%r9d
and    $0xfffffffc,%r9d
lea    -0x60(%rsp),%rcx
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
cvtps2pd (%rsi,%rax,4),%xmm0
cvtps2pd 0x8(%rsi,%rax,4),%xmm1
movapd %xmm0,%xmm2
unpcklpd %xmm0,%xmm2
unpckhpd %xmm0,%xmm0
movapd %xmm0,-0x10(%rcx)
movapd %xmm2,-0x20(%rcx)
movaps %xmm1,%xmm0
movlhps %xmm1,%xmm0
movhlps %xmm1,%xmm1
movaps %xmm1,0x10(%rcx)
movaps %xmm0,(%rcx)
add    $0x4,%rax
add    $0x40,%rcx
cmp    %rax,%r9
jne    1130 <silk_schur_FLP+0x30>
cmp    %r8,%r9
je     11a0 <silk_schur_FLP+0xa0>
mov    %r9,%rax
shl    $0x4,%rax
add    %rsp,%rax
add    $0xffffffffffffff88,%rax
nop
movss  (%rsi,%r9,4),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,(%rax)
movsd  %xmm0,-0x8(%rax)
add    $0x1,%r9
add    $0x10,%rax
cmp    %r9,%r8
jne    1180 <silk_schur_FLP+0x80>
test   %edx,%edx
jle    128d <silk_schur_FLP+0x18d>
mov    %edx,%r9d
movsd  -0x78(%rsp),%xmm0
lea    -0x60(%rsp),%r8
xor    %esi,%esi
movapd 0xe40(%rip),%xmm1
movsd  0xe48(%rip),%xmm2
xor    %r10d,%r10d
mov    %edx,%r11d
jmp    11e1 <silk_schur_FLP+0xe1>
add    $0xffffffff,%r11d
add    $0x10,%rsi
cmp    %r9,%r10
je     1287 <silk_schur_FLP+0x187>
mov    %r10,%rcx
add    $0x1,%r10
mov    %r10,%rax
shl    $0x4,%rax
movsd  -0x80(%rsp,%rax,1),%xmm3
xorpd  %xmm1,%xmm3
movapd %xmm0,%xmm4
maxsd  %xmm2,%xmm4
divsd  %xmm4,%xmm3
xorps  %xmm4,%xmm4
cvtsd2ss %xmm3,%xmm4
movss  %xmm4,(%rdi,%rcx,4)
cmp    %edx,%ecx
jge    11d0 <silk_schur_FLP+0xd0>
movsd  -0x80(%rsp,%rax,1),%xmm4
movapd %xmm0,%xmm5
mulsd  %xmm3,%xmm5
addsd  %xmm4,%xmm5
movsd  %xmm5,-0x80(%rsp,%rax,1)
mulsd  %xmm3,%xmm4
addsd  %xmm4,%xmm0
cmp    $0x1,%r11d
je     11d0 <silk_schur_FLP+0xd0>
mov    %r11d,%eax
add    $0xffffffffffffffff,%rax
mov    %r8,%rcx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  -0x8(%rcx),%xmm4
movsd  (%rcx,%rsi,1),%xmm5
movapd %xmm4,%xmm6
mulsd  %xmm3,%xmm6
addsd  %xmm5,%xmm6
movsd  %xmm6,(%rcx,%rsi,1)
mulsd  %xmm3,%xmm5
addsd  %xmm4,%xmm5
movsd  %xmm5,-0x8(%rcx)
add    $0x10,%rcx
add    $0xffffffffffffffff,%rax
jne    1250 <silk_schur_FLP+0x150>
jmp    11d0 <silk_schur_FLP+0xd0>
movsd  %xmm0,-0x78(%rsp)
movsd  -0x78(%rsp),%xmm0
cvtsd2ss %xmm0,%xmm0
add    $0x118,%rsp
ret

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

LLVM IR (-O2):
define float @silk_schur_FLP(float* nocapture writeonly %0, float* nocapture %1, i32 %2) {
%4 = alloca [25 x [2 x double]]
%5 = bitcast [25 x [2 x double]]* %4 to i8*
%6 = call i32 @llvm.smax.i32(i32 %2, i32 0)
%7 = add nuw i32 %6, 1
%8 = zext i32 %7 to i64
%9 = icmp ult i32 %6, 3
br i1 %9, label %35, label %10
10: ; preds = %3
%11 = and i64 %8, 4294967292
br label %12
12: ; preds = %12, %10
%13 = phi i64 [ 0, %10 ], [ %31, %12 ]
%14 = or i64 %13, 2
%15 = getelementptr inbounds float, float* %1, i64 %13
%16 = bitcast float* %15 to <2 x float>*
%17 = load <2 x float>, <2 x float>* %16
%18 = getelementptr inbounds float, float* %15, i64 2
%19 = bitcast float* %18 to <2 x float>*
%20 = load <2 x float>, <2 x float>* %19
%21 = fpext <2 x float> %17 to <2 x double>
%22 = fpext <2 x float> %20 to <2 x double>
%23 = getelementptr inbounds [25 x [2 x double]], [25 x [2 x double]]* %4, i64 0, i64 %13, i64 0
%24 = getelementptr inbounds [25 x [2 x double]], [25 x [2 x double]]* %4, i64 0, i64 %14, i64 0
%25 = bitcast double* %23 to <4 x double>*
%26 = bitcast double* %24 to <4 x double>*
%27 = shufflevector <2 x double> %21, <2 x double> poison, <4 x i32> <i32 0, i32 1, i32 0, i32 1>
%28 = shufflevector <4 x double> %27, <4 x double> poison, <4 x i32> <i32 0, i32 2, i32 1, i32 3>
store <4 x double> %28, <4 x double>* %25
%29 = shufflevector <2 x double> %22, <2 x double> poison, <4 x i32> <i32 0, i32 1, i32 0, i32 1>
%30 = shufflevector <4 x double> %29, <4 x double> poison, <4 x i32> <i32 0, i32 2, i32 1, i32 3>
store <4 x double> %30, <4 x double>* %26
%31 = add nuw i64 %13, 4
%32 = icmp eq i64 %31, %11
br i1 %32, label %33, label %12
33: ; preds = %12
%34 = icmp eq i64 %11, %8
br i1 %34, label %46, label %35
35: ; preds = %3, %33
%36 = phi i64 [ 0, %3 ], [ %11, %33 ]
br label %37
37: ; preds = %35, %37
%38 = phi i64 [ %44, %37 ], [ %36, %35 ]
%39 = getelementptr inbounds float, float* %1, i64 %38
%40 = load float, float* %39
%41 = fpext float %40 to double
%42 = getelementptr inbounds [25 x [2 x double]], [25 x [2 x double]]* %4, i64 0, i64 %38, i64 1
store double %41, double* %42
%43 = getelementptr inbounds [25 x [2 x double]], [25 x [2 x double]]* %4, i64 0, i64 %38, i64 0
store double %41, double* %43
%44 = add nuw nsw i64 %38, 1
%45 = icmp eq i64 %44, %8
br i1 %45, label %46, label %37
46: ; preds = %37, %33
%47 = icmp sgt i32 %2, 0
br i1 %47, label %48, label %91
48: ; preds = %46
%49 = getelementptr inbounds [25 x [2 x double]], [25 x [2 x double]]* %4, i64 0, i64 0, i64 1
%50 = zext i32 %2 to i64
%51 = load double, double* %49
br label %56
52: ; preds = %78, %71, %56
%53 = phi double [ %76, %71 ], [ %57, %56 ], [ %76, %78 ]
%54 = add i32 %59, -1
%55 = icmp eq i64 %60, %50
br i1 %55, label %90, label %56
56: ; preds = %48, %52
%57 = phi double [ %51, %48 ], [ %53, %52 ]
%58 = phi i64 [ 0, %48 ], [ %60, %52 ]
%59 = phi i32 [ %2, %48 ], [ %54, %52 ]
%60 = add nuw nsw i64 %58, 1
%61 = getelementptr inbounds [25 x [2 x double]], [25 x [2 x double]]* %4, i64 0, i64 %60, i64 0
%62 = load double, double* %61
%63 = fneg double %62
%64 = fcmp ogt double %57, 0x3E112E0BE0000000
%65 = select i1 %64, double %57, double 0x3E112E0BE0000000
%66 = fdiv double %63, %65
%67 = fptrunc double %66 to float
%68 = getelementptr inbounds float, float* %0, i64 %58
store float %67, float* %68
%69 = trunc i64 %58 to i32
%70 = icmp slt i32 %69, %2
br i1 %70, label %71, label %52
71: ; preds = %56
%72 = zext i32 %59 to i64
%73 = getelementptr inbounds [25 x [2 x double]], [25 x [2 x double]]* %4, i64 0, i64 %60, i64 0
%74 = load double, double* %73
%75 = tail call double @llvm.fmuladd.f64(double %57, double %66, double %74)
store double %75, double* %73
%76 = tail call double @llvm.fmuladd.f64(double %74, double %66, double %57)
%77 = icmp eq i32 %59, 1
br i1 %77, label %52, label %78
78: ; preds = %71, %78
%79 = phi i64 [ %88, %78 ], [ 1, %71 ]
%80 = getelementptr inbounds [25 x [2 x double]], [25 x [2 x double]]* %4, i64 0, i64 %79, i64 1
%81 = load double, double* %80
%82 = add nuw nsw i64 %60, %79
%83 = getelementptr inbounds [25 x [2 x double]], [25 x [2 x double]]* %4, i64 0, i64 %82, i64 0
%84 = load double, double* %83
%85 = getelementptr inbounds [25 x [2 x double]], [25 x [2 x double]]* %4, i64 0, i64 %79, i64 1
%86 = tail call double @llvm.fmuladd.f64(double %81, double %66, double %84)
store double %86, double* %83
%87 = tail call double @llvm.fmuladd.f64(double %84, double %66, double %81)
store double %87, double* %85
%88 = add nuw nsw i64 %79, 1
%89 = icmp eq i64 %88, %72
br i1 %89, label %52, label %78
90: ; preds = %52
store double %53, double* %49
br label %91
91: ; preds = %90, %46
%92 = getelementptr inbounds [25 x [2 x double]], [25 x [2 x double]]* %4, i64 0, i64 0, i64 1
%93 = load double, double* %92
%94 = fptrunc double %93 to float
ret float %94
}