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

```asm
<SKP_Silk_LPC_fit>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x40,%rsp
  mov    %rdi,-0x28(%rbp)
  mov    %rsi,-0x30(%rbp)
  mov    %edx,-0x34(%rbp)
  mov    %ecx,-0x38(%rbp)
  movl   $0x0,-0x14(%rbp)
  mov    $0x18,%eax
  sub    -0x34(%rbp),%eax
  mov    %eax,-0xc(%rbp)
L61:
  movl   $0x80000000,-0x10(%rbp)
  movl   $0x0,-0x18(%rbp)
  jmp    L74
L3:
  mov    -0x18(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,%edx
  neg    %edx
  cmovns %edx,%eax
  mov    %eax,-0x4(%rbp)
  mov    -0x4(%rbp),%eax
  cmp    -0x10(%rbp),%eax
  jle    L16
  mov    -0x4(%rbp),%eax
  mov    %eax,-0x10(%rbp)
  mov    -0x18(%rbp),%eax
  mov    %eax,-0x14(%rbp)
L16:
  addl   $0x1,-0x18(%rbp)
L74:
  mov    -0x18(%rbp),%eax
  cmp    -0x38(%rbp),%eax
  jl     L3
  mov    -0xc(%rbp),%eax
  mov    %eax,%ecx
  sarl   %cl,-0x10(%rbp)
  cmpl   $0x7ffe,-0x10(%rbp)
  jle    L77
  mov    -0x10(%rbp),%eax
  sub    $0x7fff,%eax
  imul   $0x3fef,%eax,%eax
  mov    -0x14(%rbp),%edx
  add    $0x1,%edx
  imul   -0x10(%rbp),%edx
  mov    %edx,%esi
  sar    $0x2,%esi
  cltd
  idiv   %esi
  mov    %eax,%edx
  mov    $0xffbe,%eax
  sub    %edx,%eax
  mov    %eax,-0x8(%rbp)
  mov    -0x8(%rbp),%edx
  mov    -0x38(%rbp),%ecx
  mov    -0x30(%rbp),%rax
  mov    %ecx,%esi
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <SKP_Silk_bwexpander_32@plt>
  jmp    L61
L77:
  nop
  movl   $0x0,-0x18(%rbp)
  jmp    L88
L22:
  cmpl   $0x1,-0xc(%rbp)
  jne    L32
  mov    -0x18(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  sar    %eax
  mov    %eax,%ecx
  mov    -0x18(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  and    $0x1,%eax
  add    %ecx,%eax
  mov    %eax,%edx
  jmp    L72
L32:
  mov    -0x18(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%edx
  mov    -0xc(%rbp),%eax
  sub    $0x1,%eax
  mov    %eax,%ecx
  sar    %cl,%edx
  mov    %edx,%eax
  add    $0x1,%eax
  sar    %eax
  mov    %eax,%edx
L72:
  mov    -0x18(%rbp),%eax
  cltq
  lea    (%rax,%rax,1),%rcx
  mov    -0x28(%rbp),%rax
  add    %rcx,%rax
  mov    %dx,(%rax)
  addl   $0x1,-0x18(%rbp)
L88:
  mov    -0x18(%rbp),%eax
  cmp    -0x38(%rbp),%eax
  jl     L22
  nop
  nop
  leave
  ret
```

assistant: ```c
void SKP_Silk_LPC_fit(short *a_QQ, int *a_Q24, const int QQ, const int L) {
  int i, rshift, idx = 0;
  int maxabs, absval, sc_Q16;
  rshift = 24 - QQ;
  while (1) {
    maxabs = ((int)0x80000000);
    for (i = 0; i < L; i++) {
      absval = (((a_Q24[i]) > 0) ? (a_Q24[i]) : -(a_Q24[i]));
      if (absval > maxabs) {
        maxabs = absval;
        idx = i;
      }
    }
    maxabs = ((maxabs) >> (rshift));
    if (maxabs >= 0x7FFF) {
      sc_Q16 = 65470 - ((int)((((65470 >> 2) * (maxabs - 0x7FFF))) / (((((maxabs) * (idx + 1))) >> (2)))));
      SKP_Silk_bwexpander_32(a_Q24, L, sc_Q16);
    } else {
      break;
    }
  };
  ;
  for (i = 0; i < L; i++) {
    a_QQ[i] = (short)((rshift) == 1 ? ((a_Q24[i]) >> 1) + ((a_Q24[i]) & 1) : (((a_Q24[i]) >> ((rshift)-1)) + 1) >> 1);
  }
}
```