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

```asm
<SKP_Silk_LPC_synthesis_order16>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x10,%rsp
  movdqu (%rsi),%xmm2
  movdqu 0x10(%rsi),%xmm3
  mov    %r8,-0x58(%rsp)
  movdqa %xmm2,%xmm0
  movdqa %xmm2,%xmm1
  punpcklwd %xmm3,%xmm2
  punpcklwd %xmm3,%xmm0
  punpckhwd %xmm3,%xmm1
  pxor   %xmm3,%xmm3
  movdqa %xmm0,%xmm4
  punpcklwd %xmm1,%xmm2
  punpcklwd %xmm1,%xmm0
  punpckhwd %xmm1,%xmm4
  punpckhwd %xmm4,%xmm0
  punpcklwd %xmm4,%xmm2
  movdqa %xmm3,%xmm4
  pcmpgtw %xmm0,%xmm4
  movdqa %xmm0,%xmm1
  movdqa %xmm2,%xmm5
  punpcklwd %xmm3,%xmm5
  punpckhwd %xmm3,%xmm2
  punpcklwd %xmm4,%xmm1
  punpckhwd %xmm4,%xmm0
  pslld  $0x10,%xmm1
  pslld  $0x10,%xmm0
  por    %xmm5,%xmm1
  por    %xmm2,%xmm0
  test   %r9d,%r9d
  jle    L89
  mov    0x3c(%rcx),%ebx
  mov    %edx,%r14d
  mov    %rcx,%rdx
  movd   %xmm1,%ecx
  mov    %rdi,%r15
  pshufd $0x55,%xmm1,%xmm2
  movswl %cx,%edi
  sar    $0x10,%ecx
  mov    %r9d,%eax
  mov    0x2c(%rdx),%r13d
  mov    0x28(%rdx),%r12d
  mov    %ecx,-0x48(%rsp)
  movd   %xmm2,%ecx
  movdqa %xmm1,%xmm2
  mov    %ebx,-0x60(%rsp)
  mov    0x24(%rdx),%ebp
  movswl %cx,%ebx
  punpckhdq %xmm1,%xmm2
  sar    $0x10,%ecx
  mov    %edi,-0x4c(%rsp)
  mov    %ecx,-0x40(%rsp)
  movd   %xmm2,%ecx
  pshufd $0xff,%xmm1,%xmm1
  sub    $0x1,%eax
  movswl %cx,%edi
  sar    $0x10,%ecx
  mov    %ebx,-0x44(%rsp)
  mov    0x1c(%rdx),%r11d
  mov    %ecx,-0x38(%rsp)
  movd   %xmm1,%ecx
  pshufd $0x55,%xmm0,%xmm1
  mov    0x18(%rdx),%r10d
  movswl %cx,%esi
  sar    $0x10,%ecx
  mov    %edi,-0x3c(%rsp)
  mov    0x14(%rdx),%r9d
  mov    %ecx,-0x30(%rsp)
  movd   %xmm0,%ecx
  lea    0x2(%r15,%rax,2),%rax
  movswl %cx,%ebx
  sar    $0x10,%ecx
  mov    %esi,-0x34(%rsp)
  mov    %ecx,-0x28(%rsp)
  movd   %xmm1,%ecx
  movdqa %xmm0,%xmm1
  movswl %cx,%edi
  punpckhdq %xmm0,%xmm1
  sar    $0x10,%ecx
  pshufd $0xff,%xmm0,%xmm0
  mov    %ecx,-0x20(%rsp)
  movd   %xmm1,%ecx
  movswl %cx,%esi
  sar    $0x10,%ecx
  mov    %ebx,-0x2c(%rsp)
  mov    %ecx,-0x18(%rsp)
  movd   %xmm0,%ecx
  movswl %cx,%ebx
  mov    %edi,-0x24(%rsp)
  mov    0x38(%rdx),%edi
  sar    $0x10,%ecx
  mov    %ebx,-0x14(%rsp)
  mov    %r14d,%ebx
  sar    $0x10,%ebx
  mov    %esi,-0x1c(%rsp)
  mov    0x34(%rdx),%esi
  mov    %ebx,-0xc(%rsp)
  mov    0x30(%rdx),%ebx
  mov    %ecx,-0x10(%rsp)
  mov    %edi,-0x5c(%rsp)
  mov    %esi,-0x68(%rsp)
  mov    %ebx,-0x64(%rsp)
  mov    0x20(%rdx),%ebx
  mov    0x10(%rdx),%r8d
  mov    0x4(%rdx),%ecx
  mov    %r15,-0x70(%rsp)
  mov    0xc(%rdx),%edi
  mov    0x8(%rdx),%esi
  mov    %rax,-0x8(%rsp)
  mov    %ecx,0x4(%rsp)
  mov    (%rdx),%ecx
  mov    -0x60(%rsp),%r15d
  mov    %ecx,-0x50(%rsp)
  movzwl %r14w,%ecx
  mov    -0x5c(%rsp),%r14d
  mov    %ecx,(%rsp)
  mov    0x4(%rsp),%ecx
  jmp    L24
  xchg   %ax,%ax
L109:
  test   %eax,%r14d
  js     L45
  mov    %r15d,%r14d
  sar    $0x9,%r14d
  add    $0x1,%r14d
  mov    %r14d,%eax
  sar    %eax
  cmp    $0xffff,%r14d
  jle    L116
  mov    $0x7fff,%eax
L116:
  mov    -0x58(%rsp),%r14
  mov    %ax,(%r14)
  mov    $0x7ffffff0,%eax
  cmp    $0x7ffffff,%r15d
  jg     L105
L44:
  mov    $0xf8000000,%eax
  cmp    %eax,%r15d
  cmovge %r15d,%eax
  shl    $0x4,%eax
L105:
  addq   $0x2,-0x70(%rsp)
  mov    -0x70(%rsp),%r15
  addq   $0x2,-0x58(%rsp)
  mov    %ecx,-0x50(%rsp)
  cmp    %r15,-0x8(%rsp)
  je     L6
  mov    -0x68(%rsp),%r14d
  mov    %esi,%ecx
  mov    %edi,%esi
  mov    %r8d,%edi
  mov    %r9d,%r8d
  mov    %r10d,%r9d
  mov    %r11d,%r10d
  mov    %ebx,%r11d
  mov    %ebp,%ebx
  mov    %r12d,%ebp
  mov    %r13d,%r12d
  mov    -0x64(%rsp),%r13d
  mov    %r14d,-0x64(%rsp)
  mov    -0x5c(%rsp),%r14d
  mov    %eax,%r15d
  mov    %r14d,-0x68(%rsp)
  mov    -0x60(%rsp),%r14d
  mov    %eax,-0x60(%rsp)
  mov    %r14d,-0x5c(%rsp)
L24:
  mov    -0x68(%rsp),%eax
  mov    %r14d,0x34(%rdx)
  mov    %r15d,0x38(%rdx)
  mov    %eax,0x30(%rdx)
  mov    -0x64(%rsp),%eax
  mov    %r13d,0x28(%rdx)
  mov    %eax,0x2c(%rdx)
  movzwl %r14w,%eax
  imul   -0x48(%rsp),%eax
  sar    $0x10,%r14d
  imul   -0x48(%rsp),%r14d
  mov    %r12d,0x24(%rdx)
  mov    %ebp,0x20(%rdx)
  sar    $0x10,%eax
  mov    %ebx,0x1c(%rdx)
  add    %r14d,%eax
  movzwl %r15w,%r14d
  sar    $0x10,%r15d
  mov    %r11d,0x18(%rdx)
  imul   -0x4c(%rsp),%r14d
  mov    %r10d,0x14(%rdx)
  imul   -0x4c(%rsp),%r15d
  mov    %r9d,0x10(%rdx)
  mov    %r8d,0xc(%rdx)
  sar    $0x10,%r14d
  mov    %edi,0x8(%rdx)
  add    %r15d,%r14d
  mov    -0x68(%rsp),%r15d
  mov    %esi,0x4(%rdx)
  add    %r14d,%eax
  mov    %ecx,(%rdx)
  movzwl %r15w,%r14d
  imul   -0x44(%rsp),%r14d
  sar    $0x10,%r15d
  imul   -0x44(%rsp),%r15d
  sar    $0x10,%r14d
  add    %r15d,%r14d
  mov    -0x64(%rsp),%r15d
  add    %r14d,%eax
  movzwl %r15w,%r14d
  imul   -0x40(%rsp),%r14d
  sar    $0x10,%r15d
  imul   -0x40(%rsp),%r15d
  sar    $0x10,%r14d
  add    %r15d,%r14d
  mov    %r13d,%r15d
  add    %r14d,%eax
  movzwl %r13w,%r14d
  sar    $0x10,%r15d
  imul   -0x3c(%rsp),%r14d
  imul   -0x3c(%rsp),%r15d
  sar    $0x10,%r14d
  add    %r15d,%r14d
  mov    %r12d,%r15d
  add    %r14d,%eax
  movzwl %r12w,%r14d
  sar    $0x10,%r15d
  imul   -0x38(%rsp),%r14d
  sar    $0x10,%r14d
  imul   -0x38(%rsp),%r15d
  add    %r15d,%r14d
  mov    %ebp,%r15d
  add    %r14d,%eax
  movzwl %bp,%r14d
  sar    $0x10,%r15d
  imul   -0x34(%rsp),%r14d
  imul   -0x34(%rsp),%r15d
  sar    $0x10,%r14d
  add    %r15d,%r14d
  mov    %ebx,%r15d
  add    %r14d,%eax
  movzwl %bx,%r14d
  sar    $0x10,%r15d
  imul   -0x30(%rsp),%r14d
  imul   -0x30(%rsp),%r15d
  sar    $0x10,%r14d
  add    %r15d,%r14d
  mov    %r11d,%r15d
  add    %r14d,%eax
  movzwl %r11w,%r14d
  sar    $0x10,%r15d
  imul   -0x2c(%rsp),%r14d
  imul   -0x2c(%rsp),%r15d
  sar    $0x10,%r14d
  add    %r15d,%r14d
  mov    %r10d,%r15d
  add    %r14d,%eax
  movzwl %r10w,%r14d
  sar    $0x10,%r15d
  imul   -0x28(%rsp),%r14d
  imul   -0x28(%rsp),%r15d
  sar    $0x10,%r14d
  add    %r15d,%r14d
  mov    %r9d,%r15d
  add    %r14d,%eax
  movzwl %r9w,%r14d
  sar    $0x10,%r15d
  imul   -0x24(%rsp),%r14d
  imul   -0x24(%rsp),%r15d
  sar    $0x10,%r14d
  add    %r15d,%r14d
  mov    %r8d,%r15d
  add    %r14d,%eax
  movzwl %r8w,%r14d
  sar    $0x10,%r15d
  imul   -0x20(%rsp),%r14d
  imul   -0x20(%rsp),%r15d
  sar    $0x10,%r14d
  add    %r15d,%r14d
  mov    %edi,%r15d
  add    %r14d,%eax
  movzwl %di,%r14d
  sar    $0x10,%r15d
  imul   -0x1c(%rsp),%r14d
  imul   -0x1c(%rsp),%r15d
  sar    $0x10,%r14d
  add    %r15d,%r14d
  mov    %esi,%r15d
  add    %r14d,%eax
  movzwl %si,%r14d
  sar    $0x10,%r15d
  imul   -0x18(%rsp),%r14d
  imul   -0x18(%rsp),%r15d
  sar    $0x10,%r14d
  add    %r15d,%r14d
  mov    %ecx,%r15d
  add    %r14d,%eax
  movzwl %cx,%r14d
  sar    $0x10,%r15d
  imul   -0x14(%rsp),%r14d
  imul   -0x14(%rsp),%r15d
  sar    $0x10,%r14d
  add    %r15d,%r14d
  mov    -0x50(%rsp),%r15d
  add    %r14d,%eax
  movzwl %r15w,%r14d
  imul   -0x10(%rsp),%r14d
  sar    $0x10,%r15d
  imul   -0x10(%rsp),%r15d
  sar    $0x10,%r14d
  add    %r15d,%r14d
  mov    (%rsp),%r15d
  add    %r14d,%eax
  mov    -0x70(%rsp),%r14
  movswl (%r14),%r14d
  imul   %r14d,%r15d
  sar    $0x10,%r15d
  imul   -0xc(%rsp),%r14d
  add    %r15d,%r14d
  mov    %r14d,%r15d
  add    %eax,%r15d
  jns    L109
  or     %eax,%r14d
  js     L91
  mov    -0x58(%rsp),%rax
  mov    $0x7fff,%r14d
  mov    %r14w,(%rax)
  mov    $0x7ffffff0,%eax
  jmp    L105
  nopl   0x0(%rax)
L45:
  mov    -0x58(%rsp),%rax
  mov    $0xffff8000,%r15d
  mov    %r15w,(%rax)
  mov    $0x80000000,%eax
  jmp    L105
  nopl   0x0(%rax)
L91:
  mov    %r15d,%r14d
  sar    $0x9,%r14d
  add    $0x1,%r14d
  mov    %r14d,%eax
  sar    %eax
  cmp    $0xffff0000,%r14d
  jge    L116
  mov    -0x58(%rsp),%rax
  mov    $0xffff8000,%r14d
  mov    %r14w,(%rax)
  jmp    L44
  nopl   0x0(%rax)
L6:
  mov    %eax,0x3c(%rdx)
L89:
  add    $0x10,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void SKP_Silk_LPC_synthesis_order16(const short *in, const short *A_Q12, const int Gain_Q26, int *S, short *out, const int len) {
  int k;
  int SA, SB, out32_Q10, out32;
  int Atmp, A_align_Q12[8];
  for (k = 0; k < 8; k++) {
    A_align_Q12[k] = (((int)A_Q12[2 * k]) & 0x0000ffff) | (((int)A_Q12[2 * k + 1]) << (16));
  }
  for (k = 0; k < len; k++) {
    SA = S[15];
    Atmp = A_align_Q12[0];
    SB = S[14];
    S[14] = SA;
    out32_Q10 = ((((SA) >> 16) * (int)((short)(Atmp))) + ((((SA)&0x0000FFFF) * (int)((short)(Atmp))) >> 16));
    out32_Q10 = ((int)((unsigned int)((out32_Q10)) + (unsigned int)(((((SB)) >> 16) * (((Atmp)) >> 16) + (((((SB)) & 0x0000FFFF) * (((Atmp)) >> 16)) >> 16)))));
    SA = S[13];
    S[13] = SB;
    Atmp = A_align_Q12[1];
    SB = S[12];
    S[12] = SA;
    out32_Q10 = ((int)((unsigned int)((out32_Q10)) + (unsigned int)((((((SA)) >> 16) * (int)((short)((Atmp)))) + (((((SA)) & 0x0000FFFF) * (int)((short)((Atmp)))) >> 16)))));
    out32_Q10 = ((int)((unsigned int)((out32_Q10)) + (unsigned int)(((((SB)) >> 16) * (((Atmp)) >> 16) + (((((SB)) & 0x0000FFFF) * (((Atmp)) >> 16)) >> 16)))));
    SA = S[11];
    S[11] = SB;
    Atmp = A_align_Q12[2];
    SB = S[10];
    S[10] = SA;
    out32_Q10 = ((int)((unsigned int)((out32_Q10)) + (unsigned int)((((((SA)) >> 16) * (int)((short)((Atmp)))) + (((((SA)) & 0x0000FFFF) * (int)((short)((Atmp)))) >> 16)))));
    out32_Q10 = ((int)((unsigned int)((out32_Q10)) + (unsigned int)(((((SB)) >> 16) * (((Atmp)) >> 16) + (((((SB)) & 0x0000FFFF) * (((Atmp)) >> 16)) >> 16)))));
    SA = S[9];
    S[9] = SB;
    Atmp = A_align_Q12[3];
    SB = S[8];
    S[8] = SA;
    out32_Q10 = ((int)((unsigned int)((out32_Q10)) + (unsigned int)((((((SA)) >> 16) * (int)((short)((Atmp)))) + (((((SA)) & 0x0000FFFF) * (int)((short)((Atmp)))) >> 16)))));
    out32_Q10 = ((int)((unsigned int)((out32_Q10)) + (unsigned int)(((((SB)) >> 16) * (((Atmp)) >> 16) + (((((SB)) & 0x0000FFFF) * (((Atmp)) >> 16)) >> 16)))));
    SA = S[7];
    S[7] = SB;
    Atmp = A_align_Q12[4];
    SB = S[6];
    S[6] = SA;
    out32_Q10 = ((int)((unsigned int)((out32_Q10)) + (unsigned int)((((((SA)) >> 16) * (int)((short)((Atmp)))) + (((((SA)) & 0x0000FFFF) * (int)((short)((Atmp)))) >> 16)))));
    out32_Q10 = ((int)((unsigned int)((out32_Q10)) + (unsigned int)(((((SB)) >> 16) * (((Atmp)) >> 16) + (((((SB)) & 0x0000FFFF) * (((Atmp)) >> 16)) >> 16)))));
    SA = S[5];
    S[5] = SB;
    Atmp = A_align_Q12[5];
    SB = S[4];
    S[4] = SA;
    out32_Q10 = ((int)((unsigned int)((out32_Q10)) + (unsigned int)((((((SA)) >> 16) * (int)((short)((Atmp)))) + (((((SA)) & 0x0000FFFF) * (int)((short)((Atmp)))) >> 16)))));
    out32_Q10 = ((int)((unsigned int)((out32_Q10)) + (unsigned int)(((((SB)) >> 16) * (((Atmp)) >> 16) + (((((SB)) & 0x0000FFFF) * (((Atmp)) >> 16)) >> 16)))));
    SA = S[3];
    S[3] = SB;
    Atmp = A_align_Q12[6];
    SB = S[2];
    S[2] = SA;
    out32_Q10 = ((int)((unsigned int)((out32_Q10)) + (unsigned int)((((((SA)) >> 16) * (int)((short)((Atmp)))) + (((((SA)) & 0x0000FFFF) * (int)((short)((Atmp)))) >> 16)))));
    out32_Q10 = ((int)((unsigned int)((out32_Q10)) + (unsigned int)(((((SB)) >> 16) * (((Atmp)) >> 16) + (((((SB)) & 0x0000FFFF) * (((Atmp)) >> 16)) >> 16)))));
    SA = S[1];
    S[1] = SB;
    Atmp = A_align_Q12[7];
    SB = S[0];
    S[0] = SA;
    out32_Q10 = ((int)((unsigned int)((out32_Q10)) + (unsigned int)((((((SA)) >> 16) * (int)((short)((Atmp)))) + (((((SA)) & 0x0000FFFF) * (int)((short)((Atmp)))) >> 16)))));
    out32_Q10 = ((int)((unsigned int)((out32_Q10)) + (unsigned int)(((((SB)) >> 16) * (((Atmp)) >> 16) + (((((SB)) & 0x0000FFFF) * (((Atmp)) >> 16)) >> 16)))));
    out32_Q10 = ((((out32_Q10) + (((((Gain_Q26) >> 16) * (int)((short)(in[k]))) + ((((Gain_Q26)&0x0000FFFF) * (int)((short)(in[k]))) >> 16)))) & 0x80000000) == 0 ? ((((out32_Q10) & (((((Gain_Q26) >> 16) * (int)((short)(in[k]))) + ((((Gain_Q26)&0x0000FFFF) * (int)((short)(in[k]))) >> 16)))) & 0x80000000) != 0 ? ((int)0x80000000) : (out32_Q10) + (((((Gain_Q26) >> 16) * (int)((short)(in[k]))) + ((((Gain_Q26)&0x0000FFFF) * (int)((short)(in[k]))) >> 16)))) : ((((out32_Q10) | (((((Gain_Q26) >> 16) * (int)((short)(in[k]))) + ((((Gain_Q26)&0x0000FFFF) * (int)((short)(in[k]))) >> 16)))) & 0x80000000) == 0 ? 0x7FFFFFFF : (out32_Q10) + (((((Gain_Q26) >> 16) * (int)((short)(in[k]))) + ((((Gain_Q26)&0x0000FFFF) * (int)((short)(in[k]))) >> 16)))));
    out32 = ((10) == 1 ? ((out32_Q10) >> 1) + ((out32_Q10)&1) : (((out32_Q10) >> ((10) - 1)) + 1) >> 1);
    out[k] = (short)((out32) > 0x7FFF ? 0x7FFF : ((out32) < ((short)0x8000) ? ((short)0x8000) : (out32)));
    S[15] = (((((((((int)0x80000000)) >> ((4)))) > (((0x7FFFFFFF) >> ((4)))) ? (((out32_Q10)) > (((((int)0x80000000)) >> ((4)))) ? (((((int)0x80000000)) >> ((4)))) : (((out32_Q10)) < (((0x7FFFFFFF) >> ((4)))) ? (((0x7FFFFFFF) >> ((4)))) : ((out32_Q10)))) : (((out32_Q10)) > (((0x7FFFFFFF) >> ((4)))) ? (((0x7FFFFFFF) >> ((4)))) : (((out32_Q10)) < (((((int)0x80000000)) >> ((4)))) ? (((((int)0x80000000)) >> ((4)))) : ((out32_Q10)))))) << ((4))));
  }
}
```