user: What is the c source code of the assembly code below?

```asm
<SKP_Silk_LPC_analysis_filter>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  add    $0xffffffffffffff80,%rsp
  mov    %rdi,-0x58(%rbp)
  mov    %rsi,-0x60(%rbp)
  mov    %rdx,-0x68(%rbp)
  mov    %rcx,-0x70(%rbp)
  mov    %r8d,-0x74(%rbp)
  mov    %r9d,-0x78(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  mov    -0x78(%rbp),%eax
  sar    %eax
  mov    %eax,-0x40(%rbp)
  movl   $0x0,-0x4c(%rbp)
  jmp    L87
L71:
  mov    -0x4c(%rbp),%eax
  cwtl
  add    %eax,%eax
  mov    %eax,-0x34(%rbp)
  mov    -0x34(%rbp),%eax
  cltq
  lea    (%rax,%rax,1),%rdx
  mov    -0x60(%rbp),%rax
  add    %rdx,%rax
  movzwl (%rax),%eax
  cwtl
  movzwl %ax,%edx
  mov    -0x34(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    (%rax,%rax,1),%rcx
  mov    -0x60(%rbp),%rax
  add    %rcx,%rax
  movzwl (%rax),%eax
  cwtl
  shl    $0x10,%eax
  or     %eax,%edx
  mov    -0x4c(%rbp),%eax
  cltq
  mov    %edx,-0x30(%rbp,%rax,4)
  addl   $0x1,-0x4c(%rbp)
L87:
  mov    -0x4c(%rbp),%eax
  cmp    -0x40(%rbp),%eax
  jl     L71
  movl   $0x0,-0x4c(%rbp)
  jmp    L21
L59:
  mov    -0x68(%rbp),%rax
  movzwl (%rax),%eax
  mov    %ax,-0x50(%rbp)
  movl   $0x0,-0x44(%rbp)
  movl   $0x0,-0x48(%rbp)
  jmp    L115
L51:
  mov    -0x48(%rbp),%eax
  cwtl
  add    %eax,%eax
  add    $0x1,%eax
  mov    %eax,-0x34(%rbp)
  mov    -0x48(%rbp),%eax
  cltq
  mov    -0x30(%rbp,%rax,4),%eax
  mov    %eax,-0x3c(%rbp)
  mov    -0x34(%rbp),%eax
  cltq
  lea    (%rax,%rax,1),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  movzwl (%rax),%eax
  mov    %ax,-0x4e(%rbp)
  mov    -0x34(%rbp),%eax
  cltq
  lea    (%rax,%rax,1),%rdx
  mov    -0x68(%rbp),%rax
  add    %rax,%rdx
  movzwl -0x50(%rbp),%eax
  mov    %ax,(%rdx)
  movswl -0x50(%rbp),%edx
  mov    -0x3c(%rbp),%eax
  cwtl
  imul   %edx,%eax
  add    %eax,-0x44(%rbp)
  movswl -0x4e(%rbp),%eax
  mov    -0x3c(%rbp),%edx
  sar    $0x10,%edx
  imul   %edx,%eax
  add    %eax,-0x44(%rbp)
  mov    -0x34(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    (%rax,%rax,1),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  movzwl (%rax),%eax
  mov    %ax,-0x50(%rbp)
  mov    -0x34(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    (%rax,%rax,1),%rdx
  mov    -0x68(%rbp),%rax
  add    %rax,%rdx
  movzwl -0x4e(%rbp),%eax
  mov    %ax,(%rdx)
  addl   $0x1,-0x48(%rbp)
L115:
  mov    -0x40(%rbp),%eax
  sub    $0x1,%eax
  cmp    %eax,-0x48(%rbp)
  jl     L51
  mov    -0x40(%rbp),%eax
  sub    $0x1,%eax
  cltq
  mov    -0x30(%rbp,%rax,4),%eax
  mov    %eax,-0x3c(%rbp)
  mov    -0x78(%rbp),%eax
  cltq
  add    %rax,%rax
  lea    -0x2(%rax),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  movzwl (%rax),%eax
  mov    %ax,-0x4e(%rbp)
  mov    -0x78(%rbp),%eax
  cltq
  add    %rax,%rax
  lea    -0x2(%rax),%rdx
  mov    -0x68(%rbp),%rax
  add    %rax,%rdx
  movzwl -0x50(%rbp),%eax
  mov    %ax,(%rdx)
  movswl -0x50(%rbp),%edx
  mov    -0x3c(%rbp),%eax
  cwtl
  imul   %edx,%eax
  add    %eax,-0x44(%rbp)
  movswl -0x4e(%rbp),%eax
  mov    -0x3c(%rbp),%edx
  sar    $0x10,%edx
  imul   %edx,%eax
  add    %eax,-0x44(%rbp)
  mov    -0x4c(%rbp),%eax
  cltq
  lea    (%rax,%rax,1),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movzwl (%rax),%eax
  cwtl
  shl    $0xc,%eax
  sub    -0x44(%rbp),%eax
  test   %eax,%eax
  js     L132
  mov    -0x4c(%rbp),%eax
  cltq
  lea    (%rax,%rax,1),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movzwl (%rax),%eax
  cwtl
  shl    $0xc,%eax
  mov    %eax,%edx
  mov    -0x44(%rbp),%eax
  xor    $0x80000000,%eax
  and    %edx,%eax
  test   %eax,%eax
  js     L138
  mov    -0x4c(%rbp),%eax
  cltq
  lea    (%rax,%rax,1),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movzwl (%rax),%eax
  cwtl
  shl    $0xc,%eax
  sub    -0x44(%rbp),%eax
  jmp    L60
L138:
  mov    $0x80000000,%eax
  jmp    L60
L132:
  mov    -0x4c(%rbp),%eax
  cltq
  lea    (%rax,%rax,1),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movzwl (%rax),%eax
  cwtl
  shl    $0xc,%eax
  xor    $0x80000000,%eax
  mov    %eax,%edx
  mov    -0x44(%rbp),%eax
  and    %edx,%eax
  test   %eax,%eax
  js     L169
  mov    -0x4c(%rbp),%eax
  cltq
  lea    (%rax,%rax,1),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movzwl (%rax),%eax
  cwtl
  shl    $0xc,%eax
  sub    -0x44(%rbp),%eax
  jmp    L60
L169:
  mov    $0x7fffffff,%eax
L60:
  mov    %eax,-0x44(%rbp)
  mov    -0x44(%rbp),%eax
  sar    $0xb,%eax
  add    $0x1,%eax
  sar    %eax
  mov    %eax,-0x38(%rbp)
  cmpl   $0x7fff,-0x38(%rbp)
  jg     L19
  cmpl   $0xffff8000,-0x38(%rbp)
  jl     L106
  mov    -0x38(%rbp),%eax
  jmp    L26
L106:
  mov    $0xffff8000,%eax
  jmp    L26
L19:
  mov    $0x7fff,%eax
L26:
  mov    -0x4c(%rbp),%edx
  movslq %edx,%rdx
  lea    (%rdx,%rdx,1),%rcx
  mov    -0x70(%rbp),%rdx
  add    %rcx,%rdx
  mov    %ax,(%rdx)
  mov    -0x4c(%rbp),%eax
  cltq
  lea    (%rax,%rax,1),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movzwl (%rax),%edx
  mov    -0x68(%rbp),%rax
  mov    %dx,(%rax)
  addl   $0x1,-0x4c(%rbp)
L21:
  mov    -0x4c(%rbp),%eax
  cmp    -0x74(%rbp),%eax
  jl     L59
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L109
  call   <__stack_chk_fail@plt>
L109:
  leave
  ret
```

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];
  }
}
```