user: Please show me the C implementation that matches the assembly code listed below:

```asm
<SKP_Silk_LPC_analysis_filter>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x48,%rsp
  mov    %rdi,%rbx
  mov    %rdx,%rbp
  mov    %rcx,%r13
  mov    %fs:0x28,%rax
  mov    %rax,0x38(%rsp)
  xor    %eax,%eax
  mov    %r9d,%r12d
  sar    %r12d
  test   %r12d,%r12d
  jle    L17
L129:
  movzwl (%rsi,%rax,4),%edx
  movswl 0x2(%rsi,%rax,4),%ecx
  shl    $0x10,%ecx
  or     %ecx,%edx
  mov    %edx,0x10(%rsp,%rax,4)
  add    $0x1,%rax
  cmp    %eax,%r12d
  jg     L129
L17:
  test   %r8d,%r8d
  jle    L51
  movslq %r9d,%r9
  lea    -0x2(%rbp,%r9,2),%r14
  mov    %r8d,%eax
  mov    %rax,0x8(%rsp)
  lea    -0x2(%r12),%eax
  lea    0x8(%rbp,%rax,4),%r10
  mov    $0x0,%r11d
  sub    $0x1,%r12d
  movslq %r12d,%r15
  jmp    L32
L107:
  mov    $0x0,%eax
  jmp    L30
L159:
  mov    $0xffff8000,%eax
L35:
  mov    %ax,0x0(%r13,%r11,2)
  movzwl (%rbx,%r11,2),%eax
  mov    %ax,0x0(%rbp)
  add    $0x1,%r11
  cmp    0x8(%rsp),%r11
  je     L51
L32:
  movzwl 0x0(%rbp),%edx
  test   %r12d,%r12d
  jle    L107
  lea    0x10(%rsp),%rdi
  lea    0x4(%rbp),%rcx
  mov    $0x0,%eax
L44:
  mov    (%rdi),%esi
  movzwl -0x2(%rcx),%r8d
  mov    %dx,-0x2(%rcx)
  movswl %dx,%edx
  movswl %si,%r9d
  imul   %r9d,%edx
  add    %eax,%edx
  sar    $0x10,%esi
  movswl %r8w,%eax
  imul   %eax,%esi
  lea    (%rsi,%rdx,1),%eax
  movzwl (%rcx),%edx
  mov    %r8w,(%rcx)
  add    $0x4,%rdi
  add    $0x4,%rcx
  cmp    %r10,%rcx
  jne    L44
L30:
  mov    0x10(%rsp,%r15,4),%ecx
  movzwl (%r14),%esi
  mov    %dx,(%r14)
  movswl %dx,%edx
  movswl %cx,%edi
  imul   %edi,%edx
  add    %eax,%edx
  sar    $0x10,%ecx
  movswl %si,%esi
  imul   %esi,%ecx
  add    %edx,%ecx
  movswl (%rbx,%r11,2),%eax
  shl    $0xc,%eax
  mov    %eax,%edx
  sub    %ecx,%edx
  js     L110
  add    $0x80000000,%ecx
  test   %eax,%ecx
  js     L159
  sar    $0xb,%edx
  add    $0x1,%edx
  mov    %edx,%eax
  sar    %eax
  cmp    $0xffff,%edx
  jle    L35
  mov    $0x7fff,%eax
  jmp    L35
L110:
  add    $0x80000000,%eax
  test   %ecx,%eax
  js     L39
  sar    $0xb,%edx
  add    $0x1,%edx
  mov    %edx,%eax
  sar    %eax
  cmp    $0xffff0000,%edx
  jge    L35
  mov    $0xffff8000,%eax
  jmp    L35
L39:
  mov    $0x7fff,%eax
  jmp    L35
L51:
  mov    0x38(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L43
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L43:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
void SKP_Silk_LPC_analysis_filter(const short *in, const short *B, short *S, short *out, const int len, const int Order) {
  int k, j, idx, Order_half = ((Order) >> (1));
  int Btmp, B_align_Q12[16 >> 1], out32_Q12, out32;
  short SA, SB;
  ;
  for (k = 0; k < Order_half; k++) {
    idx = ((int)((short)(2)) * (int)((short)(k)));
    B_align_Q12[k] = (((int)B[idx]) & 0x0000ffff) | (((int)B[idx + 1]) << (16));
  }
  for (k = 0; k < len; k++) {
    SA = S[0];
    out32_Q12 = 0;
    for (j = 0; j < (Order_half - 1); j++) {
      idx = ((int)((short)(2)) * (int)((short)(j))) + 1;
      Btmp = B_align_Q12[j];
      SB = S[idx];
      S[idx] = SA;
      out32_Q12 = ((out32_Q12) + ((int)((short)(SA))) * (int)((short)(Btmp)));
      out32_Q12 = ((out32_Q12) + ((int)((short)(SB))) * ((Btmp) >> 16));
      SA = S[idx + 1];
      S[idx + 1] = SB;
    }
    Btmp = B_align_Q12[Order_half - 1];
    SB = S[Order - 1];
    S[Order - 1] = SA;
    out32_Q12 = ((out32_Q12) + ((int)((short)(SA))) * (int)((short)(Btmp)));
    out32_Q12 = ((out32_Q12) + ((int)((short)(SB))) * ((Btmp) >> 16));
    out32_Q12 = (((((((int)in[k]) << (12))) - (out32_Q12)) & 0x80000000) == 0 ? ((((((int)in[k]) << (12))) & ((out32_Q12) ^ 0x80000000) & 0x80000000) ? ((int)0x80000000) : ((((int)in[k]) << (12))) - (out32_Q12)) : (((((((int)in[k]) << (12))) ^ 0x80000000) & (out32_Q12)&0x80000000) ? 0x7FFFFFFF : ((((int)in[k]) << (12))) - (out32_Q12)));
    out32 = ((12) == 1 ? ((out32_Q12) >> 1) + ((out32_Q12)&1) : (((out32_Q12) >> ((12) - 1)) + 1) >> 1);
    out[k] = (short)((out32) > 0x7FFF ? 0x7FFF : ((out32) < ((short)0x8000) ? ((short)0x8000) : (out32)));
    S[0] = in[k];
  }
}
```