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

```asm
<SKP_Silk_burg_modified>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x268,%rsp
  mov    %rdi,0xd8(%rsp)
  lea    0xfc(%rsp),%rdi
  mov    %rsi,0xe0(%rsp)
  lea    0xf8(%rsp),%rsi
  mov    %rdx,0xe8(%rsp)
  mov    %rcx,%rdx
  mov    %rcx,0xc0(%rsp)
  mov    %r8d,0xc8(%rsp)
  mov    %r9d,0x80(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x258(%rsp)
  mov    %r9d,%eax
  imul   %r8d,%eax
  mov    %eax,%ecx
  xor    %eax,%eax
  call   <SKP_Silk_sum_sqr_shift@plt>
  mov    0xf8(%rsp),%eax
  cmp    $0x8,%eax
  jle    L104
  movl   $0x8,0xf8(%rsp)
  pxor   %xmm0,%xmm0
  sub    $0x8,%eax
  movl   $0x8,0xb0(%rsp)
  mov    %eax,%ecx
  movaps %xmm0,0x100(%rsp)
  shll   %cl,0xfc(%rsp)
  movaps %xmm0,0x110(%rsp)
  movaps %xmm0,0x120(%rsp)
  movaps %xmm0,0x130(%rsp)
L17:
  mov    0x80(%rsp),%edi
  test   %edi,%edi
  jle    L15
  movslq 0xc8(%rsp),%rbp
  mov    0xc0(%rsp),%r14
  xor    %ebx,%ebx
  mov    0x2a8(%rsp),%r13d
  mov    %rbp,%r15
  add    %rbp,%rbp
  mov    %rbp,0x8(%rsp)
  nopl   0x0(%rax,%rax,1)
L191:
  mov    $0x1,%r12d
  lea    0x100(%rsp),%rbp
  test   %r13d,%r13d
  jle    L301
  nopl   0x0(%rax,%rax,1)
L361:
  mov    %r15d,%edx
  lea    (%r14,%r12,2),%rsi
  mov    %r14,%rdi
  xor    %eax,%eax
  sub    %r12d,%edx
  call   <SKP_Silk_inner_prod16_aligned_64@plt>
  mov    0xf8(%rsp),%ecx
  sar    %cl,%eax
  add    %eax,-0x4(%rbp,%r12,4)
  add    $0x1,%r12
  cmp    %r12d,%r13d
  jge    L361
L301:
  add    $0x1,%ebx
  add    0x8(%rsp),%r14
  cmp    %ebx,0x80(%rsp)
  jne    L191
L66:
  mov    0xf8(%rsp),%eax
  mov    %eax,0xb0(%rsp)
L15:
  movslq 0x2a0(%rsp),%rax
  mov    0x2a8(%rsp),%esi
  movdqa 0x100(%rsp),%xmm1
  movdqa 0x110(%rsp),%xmm2
  mov    %rax,%rbx
  mov    %rax,0xd0(%rsp)
  movslq 0xfc(%rsp),%rax
  movdqa 0x120(%rsp),%xmm3
  movaps %xmm1,0x140(%rsp)
  movdqa 0x130(%rsp),%xmm4
  mov    %rax,%rdx
  imul   %rbx,%rax
  movaps %xmm2,0x150(%rsp)
  movaps %xmm3,0x160(%rsp)
  movaps %xmm4,0x170(%rsp)
  sar    $0x20,%rax
  lea    0x1(%rdx,%rax,1),%r14d
  mov    %eax,%ecx
  mov    %r14d,0x28(%rsp)
  mov    %r14d,0x1c0(%rsp)
  mov    %r14d,0x210(%rsp)
  test   %esi,%esi
  jle    L32
  lea    0x1c0(%rsp),%rax
  movslq %r14d,%r15
  lea    0x180(%rsp),%r13
  movq   $0xfffffffffffffffc,0xa8(%rsp)
  mov    %rax,0x60(%rsp)
  lea    0x1c8(%rsp),%rax
  mov    %rax,0xb8(%rsp)
  mov    0x2a8(%rsp),%eax
  movl   $0xffffffff,0xa4(%rsp)
  sub    $0x1,%eax
  mov    %r13,0x38(%rsp)
  mov    %eax,0xcc(%rsp)
  movslq 0xc8(%rsp),%rax
  movq   $0x0,0x70(%rsp)
  movq   $0x1,0x20(%rsp)
  add    %rax,%rax
  mov    %rax,0x88(%rsp)
  lea    0x100(%rsp),%rax
  mov    %rax,0x48(%rsp)
  lea    0x140(%rsp),%rax
  mov    %rax,0x40(%rsp)
  lea    0x210(%rsp),%rax
  movq   $0x0,0x58(%rsp)
  mov    %rax,0x30(%rsp)
  nopl   0x0(%rax)
L56:
  mov    0x58(%rsp),%eax
  cmpl   $0xffffffff,0xb0(%rsp)
  mov    %eax,0xb4(%rsp)
  jl     L198
  mov    0x80(%rsp),%eax
  test   %eax,%eax
  jle    L402
  mov    0xc8(%rsp),%eax
  sub    0xb4(%rsp),%eax
  mov    %r15d,0x6c(%rsp)
  mov    0xb0(%rsp),%ebx
  cltq
  movl   $0x0,0x68(%rsp)
  lea    (%rax,%rax,1),%rdi
  mov    $0x10,%eax
  sub    %ebx,%eax
  mov    %rdi,0x90(%rsp)
  mov    %eax,0x98(%rsp)
  mov    $0x8,%eax
  sub    %ebx,%eax
  mov    0x70(%rsp),%rbx
  mov    %eax,0xa0(%rsp)
  mov    0xc0(%rsp),%rax
  add    %rax,%rbx
  add    %rdi,%rax
  mov    %rax,%r14
  mov    0xa4(%rsp),%eax
  mov    %rbx,0x50(%rsp)
  mov    %rax,0x18(%rsp)
  nopl   0x0(%rax,%rax,1)
L269:
  mov    0x70(%rsp),%rbx
  mov    %r14,%rax
  sub    0x90(%rsp),%rax
  movswl -0x2(%r14),%r9d
  movswl (%rax,%rbx,1),%r8d
  mov    0x98(%rsp),%ebx
  mov    %r9d,%r10d
  mov    %ebx,%ecx
  mov    %r8d,%eax
  mov    %r8d,%r11d
  shl    $0x8,%r10d
  shl    %cl,%eax
  shl    $0x8,%r11d
  neg    %eax
  mov    %eax,%edx
  mov    %r9d,%eax
  shl    %cl,%eax
  neg    %eax
  cmpq   $0x0,0x58(%rsp)
  je     L404
  mov    %edx,%ebx
  mov    %r8d,0x78(%rsp)
  movzwl %dx,%edi
  mov    0x40(%rsp),%r13
  sar    $0x10,%ebx
  mov    0x38(%rsp),%r12
  mov    0x48(%rsp),%r15
  movzwl %ax,%ebp
  mov    %ebx,0x8(%rsp)
  mov    %eax,%ebx
  xor    %edx,%edx
  sar    $0x10,%ebx
  mov    %r9d,0x84(%rsp)
  mov    %edi,%r9d
  mov    %ebx,0x14(%rsp)
  mov    0x50(%rsp),%rbx
  nopl   0x0(%rax,%rax,1)
L51:
  movswl -0x2(%rbx),%esi
  mov    0x8(%rsp),%ecx
  sub    $0x2,%rbx
  mov    0x14(%rsp),%edi
  mov    %esi,%eax
  imul   %esi,%ecx
  imul   %r9d,%eax
  sar    $0x10,%eax
  add    %ecx,%eax
  movswl (%r14,%rdx,2),%ecx
  add    %eax,(%r15,%rdx,4)
  mov    %ecx,%eax
  imul   %ecx,%edi
  imul   %ebp,%eax
  sar    $0x10,%eax
  add    %edi,%eax
  add    %eax,0x0(%r13,%rdx,4)
  mov    (%r12,%rdx,4),%eax
  mov    %esi,%edi
  mov    %eax,%r8d
  movzwl %ax,%eax
  imul   %eax,%edi
  sar    $0x10,%r8d
  imul   %ecx,%eax
  imul   %r8d,%esi
  imul   %r8d,%ecx
  sar    $0x10,%edi
  sar    $0x10,%eax
  add    %esi,%edi
  add    %ecx,%eax
  add    %edi,%r11d
  add    %eax,%r10d
  mov    %rdx,%rax
  add    $0x1,%rdx
  cmp    %rax,0x18(%rsp)
  jne    L51
  mov    0x78(%rsp),%r8d
  mov    0x84(%rsp),%r9d
L404:
  mov    0xa0(%rsp),%eax
  neg    %r11d
  neg    %r10d
  xor    %edx,%edx
  mov    0x20(%rsp),%r15
  mov    0x6c(%rsp),%ebx
  mov    %eax,%ecx
  mov    0x28(%rsp),%ebp
  mov    0x30(%rsp),%r12
  shl    %cl,%r11d
  shl    %cl,%r10d
  mov    0x60(%rsp),%r13
  mov    0x50(%rsp),%rcx
  mov    %r11d,%edi
  mov    %r10d,%esi
  movzwl %r11w,%r11d
  movzwl %r10w,%r10d
  sar    $0x10,%edi
  sar    $0x10,%esi
  jmp    L256
  nopl   0x0(%rax,%rax,1)
L354:
  mov    0x0(%r13,%rdx,4),%ebx
  mov    (%r12,%rdx,4),%ebp
  movswl (%rcx),%r8d
  movswl -0x2(%r14,%rdx,2),%r9d
L256:
  mov    %r8d,%eax
  imul   %edi,%r8d
  sub    $0x2,%rcx
  imul   %r11d,%eax
  sar    $0x10,%eax
  add    %r8d,%eax
  add    %ebx,%eax
  mov    %eax,0x0(%r13,%rdx,4)
  mov    %r9d,%eax
  imul   %esi,%r9d
  imul   %r10d,%eax
  sar    $0x10,%eax
  add    %r9d,%eax
  add    %ebp,%eax
  mov    %eax,(%r12,%rdx,4)
  add    $0x1,%rdx
  cmp    %rdx,%r15
  jne    L354
  mov    0x210(%rsp),%ebx
  mov    0x88(%rsp),%rdi
  addl   $0x1,0x68(%rsp)
  mov    0x68(%rsp),%eax
  mov    %ebx,0x28(%rsp)
  mov    0x1c0(%rsp),%ebx
  add    %rdi,%r14
  add    %rdi,0x50(%rsp)
  mov    %ebx,0x6c(%rsp)
  cmp    %eax,0x80(%rsp)
  jne    L269
L186:
  movslq %ebx,%r15
L402:
  mov    0x58(%rsp),%rax
  mov    0x48(%rsp),%rbx
  mov    0x28(%rsp),%r14d
  mov    (%rbx,%rax,4),%ebx
  add    %r15d,%r14d
  mov    %ebx,0x8(%rsp)
  mov    0x40(%rsp),%rbx
  mov    (%rbx,%rax,4),%r12d
  mov    0x20(%rsp),%ebx
  mov    %ebx,0x68(%rsp)
  test   %rax,%rax
  je     L49
  mov    0xa4(%rsp),%eax
  mov    $0x1,%r13d
  movl   $0x0,0x14(%rsp)
  mov    %r15d,0x50(%rsp)
  mov    %r13,%rbp
  mov    0xa8(%rsp),%rbx
  add    $0x2,%rax
  mov    0x60(%rsp),%r13
  mov    %rax,0x18(%rsp)
  lea    0x214(%rsp),%rax
  mov    %rax,0x28(%rsp)
  nopl   0x0(%rax)
L28:
  mov    0x38(%rsp),%rax
  mov    -0x4(%rax,%rbp,4),%r15d
  mov    %r15d,%edi
  neg    %edi
  cmovs  %r15d,%edi
  xor    %eax,%eax
  call   <SKP_Silk_CLZ32@plt>
  mov    $0x8,%ecx
  mov    0x30(%rsp),%rdx
  sub    $0x1,%eax
  cmp    %ecx,%eax
  cmovg  %ecx,%eax
  mov    %eax,%ecx
  shl    %cl,%r15d
  mov    $0x8,%ecx
  sub    %eax,%ecx
  mov    0x40(%rsp),%rax
  movslq %r15d,%rdi
  movslq (%rax,%rbx,1),%rax
  imul   %rdi,%rax
  sar    $0x20,%rax
  shl    %cl,%eax
  add    %eax,0x8(%rsp)
  mov    0x48(%rsp),%rax
  movslq (%rax,%rbx,1),%rax
  imul   %rdi,%rax
  sar    $0x20,%rax
  shl    %cl,%eax
  add    %eax,%r12d
  mov    0x28(%rsp),%rax
  movslq (%rax,%rbx,1),%rax
  sub    $0x4,%rbx
  imul   %rdi,%rax
  sar    $0x20,%rax
  shl    %cl,%eax
  add    %eax,0x14(%rsp)
  mov    0x0(%r13,%rbp,4),%eax
  add    (%rdx,%rbp,4),%eax
  add    $0x1,%rbp
  cltq
  imul   %rdi,%rax
  sar    $0x20,%rax
  shl    %cl,%eax
  add    %eax,%r14d
  cmp    %rbp,0x18(%rsp)
  jne    L28
  mov    0x20(%rsp),%rax
  mov    0x60(%rsp),%rbx
  mov    0x8(%rsp),%edi
  movslq 0x50(%rsp),%r15
  mov    %edi,(%rbx,%rax,4)
  mov    0x14(%rsp),%edi
  mov    0x30(%rsp),%rbx
  add    %r12d,%edi
  mov    %r12d,(%rbx,%rax,4)
  mov    %edi,%ebx
  neg    %ebx
  lea    (%rbx,%rbx,1),%edi
  mov    %edi,%eax
  neg    %eax
  cmovs  %edi,%eax
  cmp    %r14d,%eax
  jge    L196
  mov    $0x1f,%edx
  mov    %r14d,%esi
  xor    %eax,%eax
  call   <SKP_DIV32_varQ@plt>
  mov    0x68(%rsp),%r10d
  mov    %eax,%r11d
  sar    %r10d
L11:
  mov    0x38(%rsp),%rsi
  movslq %r11d,%rcx
  xor    %r8d,%r8d
  mov    %rsi,%rdi
  add    0xa8(%rsp),%rsi
  xchg   %ax,%ax
L309:
  movslq (%rsi),%rdx
  movslq (%rdi),%rax
  add    $0x1,%r8d
  add    $0x4,%rdi
  sub    $0x4,%rsi
  mov    %rdx,%r9
  imul   %rcx,%rdx
  sar    $0x20,%rdx
  lea    (%rax,%rdx,2),%edx
  imul   %rcx,%rax
  mov    %edx,-0x4(%rdi)
  sar    $0x20,%rax
  lea    (%r9,%rax,2),%eax
  mov    %eax,0x4(%rsi)
  cmp    %r10d,%r8d
  jl     L309
L364:
  mov    0x38(%rsp),%rdi
  mov    0x58(%rsp),%rbx
  mov    %r11d,%eax
  sar    $0x7,%eax
  mov    0x60(%rsp),%rdx
  mov    %eax,(%rdi,%rbx,4)
  mov    0x30(%rsp),%rax
  mov    0xb8(%rsp),%rdi
  lea    (%rax,%rbx,4),%rsi
  jmp    L436
  nopw   0x0(%rax,%rax,1)
L50:
  movslq (%rdx),%r15
  mov    (%rsi),%r12d
  sub    $0x4,%rsi
L436:
  movslq %r12d,%rax
  add    $0x4,%rdx
  imul   %rcx,%rax
  sar    $0x20,%rax
  lea    (%r15,%rax,2),%eax
  imul   %rcx,%r15
  mov    %eax,-0x4(%rdx)
  sar    $0x20,%r15
  lea    (%r12,%r15,2),%eax
  mov    %eax,0x4(%rsi)
  cmp    %rdi,%rdx
  jne    L50
  mov    0xf8(%rsp),%eax
  addq   $0x1,0x58(%rsp)
  addl   $0x1,0xa4(%rsp)
  movslq 0x1c0(%rsp),%r15
  mov    %eax,0xb0(%rsp)
  mov    0xa4(%rsp),%eax
  addq   $0x4,0xb8(%rsp)
  addq   $0x4,0xa8(%rsp)
  cmp    %eax,0xcc(%rsp)
  je     L222
  mov    0x210(%rsp),%eax
  addq   $0x1,0x20(%rsp)
  addq   $0x2,0x70(%rsp)
  mov    %eax,0x28(%rsp)
  jmp    L56
L198:
  mov    0x80(%rsp),%edx
  test   %edx,%edx
  jle    L402
  mov    0xc8(%rsp),%eax
  sub    0xb4(%rsp),%eax
  mov    %r15d,0x6c(%rsp)
  cltq
  movl   $0x0,0x68(%rsp)
  lea    (%rax,%rax,1),%rdi
  mov    0xb0(%rsp),%eax
  mov    %rdi,0x78(%rsp)
  mov    %eax,%ebx
  not    %eax
  neg    %ebx
  mov    %eax,%r11d
  mov    0xc0(%rsp),%rax
  mov    %ebx,0x84(%rsp)
  mov    0x70(%rsp),%rbx
  mov    %r11d,0x90(%rsp)
  add    %rax,%rbx
  add    %rdi,%rax
  mov    %rax,%r10
  mov    0xa4(%rsp),%eax
  mov    %rbx,0x50(%rsp)
  mov    %rax,0x98(%rsp)
L349:
  mov    0x70(%rsp),%rbx
  mov    %r10,%rax
  sub    0x78(%rsp),%rax
  movswl -0x2(%r10),%edx
  movswl (%rax,%rbx,1),%eax
  mov    0x84(%rsp),%ebx
  mov    %edx,%r14d
  mov    %edx,%edi
  mov    %ebx,%ecx
  mov    %eax,%r12d
  mov    %eax,%r8d
  shl    $0x11,%edi
  shl    %cl,%r12d
  shl    %cl,%r14d
  neg    %r12d
  neg    %r14d
  shl    $0x11,%r8d
  cmpq   $0x0,0x58(%rsp)
  je     L366
  mov    0x50(%rsp),%rbp
  mov    0x40(%rsp),%r11
  mov    %eax,0x8(%rsp)
  xor    %ecx,%ecx
  mov    0x38(%rsp),%r13
  mov    0x48(%rsp),%r15
  mov    0x98(%rsp),%rax
  nopw   0x0(%rax,%rax,1)
L64:
  movswl -0x2(%rbp),%ebx
  movswl (%r10,%rcx,2),%r9d
  sub    $0x2,%rbp
  mov    %ebx,%esi
  imul   %r12d,%esi
  add    %esi,(%r15,%rcx,4)
  mov    %r9d,%esi
  imul   %r14d,%esi
  add    %esi,(%r11,%rcx,4)
  mov    0x0(%r13,%rcx,4),%esi
  sar    $0x6,%esi
  add    $0x1,%esi
  sar    %esi
  imul   %esi,%ebx
  imul   %esi,%r9d
  mov    %rcx,%rsi
  add    $0x1,%rcx
  add    %ebx,%r8d
  add    %r9d,%edi
  cmp    %rsi,%rax
  jne    L64
  mov    0x8(%rsp),%eax
L366:
  neg    %r8d
  neg    %edi
  mov    0x50(%rsp),%r9
  mov    0x28(%rsp),%r13d
  mov    %r8d,%ebx
  mov    %edi,%r12d
  mov    %r8d,0x18(%rsp)
  mov    0x30(%rsp),%r14
  sar    $0x10,%ebx
  sar    $0x10,%r12d
  mov    0x90(%rsp),%ecx
  xor    %esi,%esi
  mov    %ebx,0x8(%rsp)
  mov    0x60(%rsp),%r15
  movzwl %r8w,%ebp
  movzwl %di,%ebx
  mov    %r12d,0x14(%rsp)
  mov    0x6c(%rsp),%r12d
  jmp    L31
  nopw   0x0(%rax,%rax,1)
L195:
  mov    (%r15,%rsi,4),%r12d
  mov    (%r14,%rsi,4),%r13d
  movswl (%r9),%eax
  movswl -0x2(%r10,%rsi,2),%edx
L31:
  shl    %cl,%eax
  shl    %cl,%edx
  sub    $0x2,%r9
  movswl %ax,%r11d
  sar    $0xf,%eax
  mov    %r11d,%r8d
  imul   0x8(%rsp),%r11d
  add    $0x1,%eax
  imul   %ebp,%r8d
  sar    %eax
  imul   0x18(%rsp),%eax
  sar    $0x10,%r8d
  add    %r11d,%r8d
  add    %r12d,%r8d
  add    %eax,%r8d
  mov    %r8d,(%r15,%rsi,4)
  movswl %dx,%r8d
  sar    $0xf,%edx
  mov    %r8d,%eax
  imul   0x14(%rsp),%r8d
  add    $0x1,%edx
  imul   %ebx,%eax
  sar    %edx
  imul   %edi,%edx
  sar    $0x10,%eax
  add    %r8d,%eax
  add    %r13d,%eax
  add    %edx,%eax
  mov    %eax,(%r14,%rsi,4)
  add    $0x1,%rsi
  cmp    %rsi,0x20(%rsp)
  jne    L195
  mov    0x210(%rsp),%ebx
  mov    0x88(%rsp),%rdi
  addl   $0x1,0x68(%rsp)
  mov    0x68(%rsp),%eax
  mov    %ebx,0x28(%rsp)
  mov    0x1c0(%rsp),%ebx
  add    %rdi,%r10
  add    %rdi,0x50(%rsp)
  mov    %ebx,0x6c(%rsp)
  cmp    %eax,0x80(%rsp)
  jne    L349
  jmp    L186
L49:
  mov    0x20(%rsp),%rax
  mov    0x60(%rsp),%rbx
  mov    %r12d,%edi
  mov    0x8(%rsp),%ecx
  neg    %edi
  add    %edi,%edi
  mov    %ecx,(%rbx,%rax,4)
  mov    0x30(%rsp),%rbx
  mov    %r12d,(%rbx,%rax,4)
  mov    %edi,%eax
  neg    %eax
  cmovs  %edi,%eax
  cmp    %r14d,%eax
  jge    L399
  xor    %eax,%eax
  mov    $0x1f,%edx
  mov    %r14d,%esi
  call   <SKP_DIV32_varQ@plt>
  mov    0x68(%rsp),%r10d
  mov    %eax,%r11d
  sar    %r10d
  jne    L11
  movslq %eax,%rcx
  jmp    L364
L222:
  mov    0x38(%rsp),%r13
  mov    %r15d,%r14d
L242:
  mov    0x2a8(%rsp),%r11d
  mov    0x60(%rsp),%r10
  xor    %edi,%edi
  mov    $0x10000,%eax
  mov    0xe8(%rsp),%r9
  nopl   0x0(%rax)
L252:
  mov    0x0(%r13,%rdi,4),%edx
  mov    0x4(%r10,%rdi,4),%ebp
  sar    $0x7,%edx
  movzwl %bp,%ecx
  mov    %ebp,%esi
  add    $0x1,%edx
  sar    $0x10,%esi
  mov    %edx,%r8d
  sar    %r8d
  movswl %r8w,%ebx
  imul   %ebx,%ecx
  imul   %ebx,%esi
  sar    $0x10,%ecx
  add    %esi,%ecx
  mov    %edx,%esi
  sar    $0x11,%edx
  sar    $0x10,%esi
  imul   %ebx,%edx
  add    %r14d,%ecx
  add    $0x1,%esi
  sar    %esi
  imul   %esi,%ebp
  imul   %r8d,%esi
  lea    0x0(%rbp,%rcx,1),%r14d
  movzwl %r8w,%ecx
  neg    %r8d
  imul   %ebx,%ecx
  mov    %r8d,(%r9,%rdi,4)
  add    $0x1,%rdi
  sar    $0x10,%ecx
  add    %ecx,%edx
  add    %edx,%eax
  add    %esi,%eax
  cmp    %r11,%rdi
  jne    L252
  movslq 0xfc(%rsp),%rdx
  mov    %eax,%edi
  neg    %eax
  imul   0xd0(%rsp),%rdx
  neg    %edi
  sar    $0xf,%eax
  movswl %di,%edi
  add    $0x1,%eax
  sar    %eax
  sar    $0x20,%rdx
  movzwl %dx,%esi
  mov    %edx,%ecx
  sar    $0x10,%edx
  imul   %edi,%esi
  imul   %edi,%edx
  sar    $0x10,%esi
  add    %esi,%edx
  lea    (%rdx,%r14,1),%ebx
  mov    %ebx,0x28(%rsp)
L27:
  imul   %eax,%ecx
  mov    0x28(%rsp),%eax
  mov    0xd8(%rsp),%rbx
  add    %ecx,%eax
  mov    %eax,(%rbx)
  mov    0xb0(%rsp),%eax
  mov    0xe0(%rsp),%rbx
  neg    %eax
  mov    %eax,(%rbx)
  mov    0x258(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L379
  add    $0x268,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L196:
  mov    0xf8(%rsp),%eax
  mov    0x38(%rsp),%r13
  mov    %r15d,%r14d
  mov    %eax,0xb0(%rsp)
L192:
  movslq 0xb4(%rsp),%rax
  mov    0x2a8(%rsp),%edx
  xor    %esi,%esi
  sub    %eax,%edx
  lea    0x0(%r13,%rax,4),%rdi
  movslq %edx,%rdx
  shl    $0x2,%rdx
  call   <memset@plt>
  jmp    L242
L104:
  mov    0xfc(%rsp),%edi
  xor    %eax,%eax
  call   <SKP_Silk_CLZ32@plt>
  mov    0xf8(%rsp),%esi
  mov    %eax,%r8d
  mov    $0x3,%eax
  sub    %r8d,%eax
  test   %eax,%eax
  jle    L331
  mov    $0x8,%edx
  sub    %esi,%edx
  cmp    %eax,%edx
  cmovle %edx,%eax
  mov    0xfc(%rsp),%edx
  mov    %eax,%ecx
  sar    %cl,%edx
L333:
  pxor   %xmm0,%xmm0
  add    %esi,%eax
  mov    %edx,0xfc(%rsp)
  mov    %eax,0xb0(%rsp)
  mov    %eax,0xf8(%rsp)
  movaps %xmm0,0x100(%rsp)
  movaps %xmm0,0x110(%rsp)
  movaps %xmm0,0x120(%rsp)
  movaps %xmm0,0x130(%rsp)
  test   %eax,%eax
  jg     L17
  mov    0x80(%rsp),%r8d
  test   %r8d,%r8d
  jle    L15
  movslq 0xc8(%rsp),%rbp
  mov    0xc0(%rsp),%r14
  xor    %ebx,%ebx
  mov    0x2a8(%rsp),%r13d
  mov    %rbp,%r15
  add    %rbp,%rbp
  mov    %rbp,0x8(%rsp)
L322:
  mov    $0x1,%r12d
  lea    0x100(%rsp),%rbp
  test   %r13d,%r13d
  jle    L179
  cs nopw 0x0(%rax,%rax,1)
L108:
  mov    %r15d,%edx
  lea    (%r14,%r12,2),%rsi
  mov    %r14,%rdi
  xor    %eax,%eax
  sub    %r12d,%edx
  call   <SKP_Silk_inner_prod_aligned@plt>
  mov    0xf8(%rsp),%ecx
  neg    %ecx
  shl    %cl,%eax
  add    %eax,-0x4(%rbp,%r12,4)
  add    $0x1,%r12
  cmp    %r12d,%r13d
  jge    L108
L179:
  add    $0x1,%ebx
  add    0x8(%rsp),%r14
  cmp    %ebx,0x80(%rsp)
  jne    L322
  jmp    L66
L331:
  mov    $0xfffffff0,%edx
  sub    %esi,%edx
  cmp    %eax,%edx
  cmovge %edx,%eax
  mov    0xfc(%rsp),%edx
  mov    %eax,%ecx
  neg    %ecx
  shl    %cl,%edx
  jmp    L333
L32:
  mov    $0xffffffff,%eax
  jmp    L27
L399:
  mov    0x38(%rsp),%r13
  mov    %r15d,%r14d
  jmp    L192
L379:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
void SKP_Silk_burg_modified(int *res_nrg, int *res_nrg_Q, int A_Q16[], const short x[], const int subfr_length, const int nb_subfr, const int WhiteNoiseFrac_Q32, const int D) {
  int k, n, s, lz, rshifts, rshifts_extra;
  int C0, num, nrg, rc_Q31, Atmp_QA, Atmp1, tmp1, tmp2, x1, x2;
  const short *x_ptr;
  int C_first_row[16];
  int C_last_row[16];
  int Af_QA[16];
  int CAf[16 + 1];
  int CAb[16 + 1];
  ;
  ;
  SKP_Silk_sum_sqr_shift(&C0, &rshifts, x, nb_subfr * subfr_length);
  if (rshifts > (32 - 24)) {
    C0 = ((C0) << (rshifts - (32 - 24)));
    ;
    rshifts = (32 - 24);
  } else {
    lz = SKP_Silk_CLZ32(C0) - 1;
    rshifts_extra = 2 - lz;
    if (rshifts_extra > 0) {
      rshifts_extra = (((rshifts_extra) < ((32 - 24) - rshifts)) ? (rshifts_extra) : ((32 - 24) - rshifts));
      C0 = ((C0) >> (rshifts_extra));
    } else {
      rshifts_extra = (((rshifts_extra) > (-16 - rshifts)) ? (rshifts_extra) : (-16 - rshifts));
      C0 = ((C0) << (-rshifts_extra));
    }
    rshifts += rshifts_extra;
  }
  memset((C_first_row), (0), (16 * sizeof(int)));
  if (rshifts > 0) {
    for (s = 0; s < nb_subfr; s++) {
      x_ptr = x + s * subfr_length;
      for (n = 1; n < D + 1; n++) {
        C_first_row[n - 1] += (int)((SKP_Silk_inner_prod16_aligned_64(x_ptr, x_ptr + n, subfr_length - n)) >> (rshifts));
      }
    }
  } else {
    for (s = 0; s < nb_subfr; s++) {
      x_ptr = x + s * subfr_length;
      for (n = 1; n < D + 1; n++) {
        C_first_row[n - 1] += ((SKP_Silk_inner_prod_aligned(x_ptr, x_ptr + n, subfr_length - n)) << (-rshifts));
      }
    }
  }
  memcpy((C_last_row), (C_first_row), (16 * sizeof(int)));
  CAb[0] = CAf[0] = C0 + (int)((((long long)((WhiteNoiseFrac_Q32)) * ((C0)))) >> (32)) + 1;
  for (n = 0; n < D; n++) {
    if (rshifts > -2) {
      for (s = 0; s < nb_subfr; s++) {
        x_ptr = x + s * subfr_length;
        x1 = -(((int)x_ptr[n]) << (16 - rshifts));
        x2 = -(((int)x_ptr[subfr_length - n - 1]) << (16 - rshifts));
        tmp1 = (((int)x_ptr[n]) << (24 - 16));
        tmp2 = (((int)x_ptr[subfr_length - n - 1]) << (24 - 16));
        for (k = 0; k < n; k++) {
          C_first_row[k] = ((C_first_row[k]) + ((((x1) >> 16) * (int)((short)(x_ptr[n - k - 1]))) + ((((x1)&0x0000FFFF) * (int)((short)(x_ptr[n - k - 1]))) >> 16)));
          C_last_row[k] = ((C_last_row[k]) + ((((x2) >> 16) * (int)((short)(x_ptr[subfr_length - n + k]))) + ((((x2)&0x0000FFFF) * (int)((short)(x_ptr[subfr_length - n + k]))) >> 16)));
          Atmp_QA = Af_QA[k];
          tmp1 = ((tmp1) + ((((Atmp_QA) >> 16) * (int)((short)(x_ptr[n - k - 1]))) + ((((Atmp_QA)&0x0000FFFF) * (int)((short)(x_ptr[n - k - 1]))) >> 16)));
          tmp2 = ((tmp2) + ((((Atmp_QA) >> 16) * (int)((short)(x_ptr[subfr_length - n + k]))) + ((((Atmp_QA)&0x0000FFFF) * (int)((short)(x_ptr[subfr_length - n + k]))) >> 16)));
        }
        tmp1 = ((-tmp1) << (32 - 24 - rshifts));
        tmp2 = ((-tmp2) << (32 - 24 - rshifts));
        for (k = 0; k <= n; k++) {
          CAf[k] = ((CAf[k]) + ((((tmp1) >> 16) * (int)((short)(x_ptr[n - k]))) + ((((tmp1)&0x0000FFFF) * (int)((short)(x_ptr[n - k]))) >> 16)));
          CAb[k] = ((CAb[k]) + ((((tmp2) >> 16) * (int)((short)(x_ptr[subfr_length - n + k - 1]))) + ((((tmp2)&0x0000FFFF) * (int)((short)(x_ptr[subfr_length - n + k - 1]))) >> 16)));
        }
      }
    } else {
      for (s = 0; s < nb_subfr; s++) {
        x_ptr = x + s * subfr_length;
        x1 = -(((int)x_ptr[n]) << (-rshifts));
        x2 = -(((int)x_ptr[subfr_length - n - 1]) << (-rshifts));
        tmp1 = (((int)x_ptr[n]) << (17));
        tmp2 = (((int)x_ptr[subfr_length - n - 1]) << (17));
        for (k = 0; k < n; k++) {
          C_first_row[k] = (((C_first_row[k])) + (((x1) * (x_ptr[n - k - 1]))));
          C_last_row[k] = (((C_last_row[k])) + (((x2) * (x_ptr[subfr_length - n + k]))));
          Atmp1 = ((24 - 17) == 1 ? ((Af_QA[k]) >> 1) + ((Af_QA[k]) & 1) : (((Af_QA[k]) >> ((24 - 17) - 1)) + 1) >> 1);
          tmp1 = (((tmp1)) + (((x_ptr[n - k - 1]) * (Atmp1))));
          tmp2 = (((tmp2)) + (((x_ptr[subfr_length - n + k]) * (Atmp1))));
        }
        tmp1 = -tmp1;
        tmp2 = -tmp2;
        for (k = 0; k <= n; k++) {
          CAf[k] = ((((((CAf[k])) + (((((tmp1)) >> 16) * (int)((short)(((((int)x_ptr[n - k]) << (-rshifts - 1)))))) + (((((tmp1)) & 0x0000FFFF) * (int)((short)(((((int)x_ptr[n - k]) << (-rshifts - 1)))))) >> 16))))) + ((((tmp1)) * (((16) == 1 ? ((((((int)x_ptr[n - k]) << (-rshifts - 1)))) >> 1) + ((((((int)x_ptr[n - k]) << (-rshifts - 1)))) & 1) : (((((((int)x_ptr[n - k]) << (-rshifts - 1)))) >> ((16) - 1)) + 1) >> 1)))));
          CAb[k] = ((((((CAb[k])) + (((((tmp2)) >> 16) * (int)((short)(((((int)x_ptr[subfr_length - n + k - 1]) << (-rshifts - 1)))))) + (((((tmp2)) & 0x0000FFFF) * (int)((short)(((((int)x_ptr[subfr_length - n + k - 1]) << (-rshifts - 1)))))) >> 16))))) + ((((tmp2)) * (((16) == 1 ? ((((((int)x_ptr[subfr_length - n + k - 1]) << (-rshifts - 1)))) >> 1) + ((((((int)x_ptr[subfr_length - n + k - 1]) << (-rshifts - 1)))) & 1) : (((((((int)x_ptr[subfr_length - n + k - 1]) << (-rshifts - 1)))) >> ((16) - 1)) + 1) >> 1)))));
        }
      }
    }
    tmp1 = C_first_row[n];
    tmp2 = C_last_row[n];
    num = 0;
    nrg = ((CAb[0]) + (CAf[0]));
    for (k = 0; k < n; k++) {
      Atmp_QA = Af_QA[k];
      lz = SKP_Silk_CLZ32((((Atmp_QA) > 0) ? (Atmp_QA) : -(Atmp_QA))) - 1;
      lz = (((32 - 24) < (lz)) ? (32 - 24) : (lz));
      Atmp1 = ((Atmp_QA) << (lz));
      tmp1 = (((tmp1)) + (((((int)((((long long)((C_last_row[n - k - 1])) * ((Atmp1)))) >> (32)))) << ((32 - 24 - lz)))));
      tmp2 = (((tmp2)) + (((((int)((((long long)((C_first_row[n - k - 1])) * ((Atmp1)))) >> (32)))) << ((32 - 24 - lz)))));
      num = (((num)) + (((((int)((((long long)((CAb[n - k])) * ((Atmp1)))) >> (32)))) << ((32 - 24 - lz)))));
      nrg = (((nrg)) + (((((int)((((long long)((((CAb[k + 1]) + (CAf[k + 1])))) * ((Atmp1)))) >> (32)))) << ((32 - 24 - lz)))));
    }
    CAf[n + 1] = tmp1;
    CAb[n + 1] = tmp2;
    num = ((num) + (tmp2));
    num = ((-num) << (1));
    if ((((num) > 0) ? (num) : -(num)) < nrg) {
      rc_Q31 = SKP_DIV32_varQ(num, nrg, 31);
    } else {
      memset((&Af_QA[n]), (0), ((D - n) * sizeof(int)));
      ;
      break;
    }
    for (k = 0; k < (n + 1) >> 1; k++) {
      tmp1 = Af_QA[k];
      tmp2 = Af_QA[n - k - 1];
      Af_QA[k] = (((tmp1)) + (((((int)((((long long)((tmp2)) * ((rc_Q31)))) >> (32)))) << ((1)))));
      Af_QA[n - k - 1] = (((tmp2)) + (((((int)((((long long)((tmp1)) * ((rc_Q31)))) >> (32)))) << ((1)))));
    }
    Af_QA[n] = ((rc_Q31) >> (31 - 24));
    for (k = 0; k <= n + 1; k++) {
      tmp1 = CAf[k];
      tmp2 = CAb[n - k + 1];
      CAf[k] = (((tmp1)) + (((((int)((((long long)((tmp2)) * ((rc_Q31)))) >> (32)))) << ((1)))));
      CAb[n - k + 1] = (((tmp2)) + (((((int)((((long long)((tmp1)) * ((rc_Q31)))) >> (32)))) << ((1)))));
    }
  }
  nrg = CAf[0];
  tmp1 = 1 << 16;
  for (k = 0; k < D; k++) {
    Atmp1 = ((24 - 16) == 1 ? ((Af_QA[k]) >> 1) + ((Af_QA[k]) & 1) : (((Af_QA[k]) >> ((24 - 16) - 1)) + 1) >> 1);
    nrg = ((((((nrg)) + (((((CAf[k + 1])) >> 16) * (int)((short)((Atmp1)))) + (((((CAf[k + 1])) & 0x0000FFFF) * (int)((short)((Atmp1)))) >> 16))))) + ((((CAf[k + 1])) * (((16) == 1 ? (((Atmp1)) >> 1) + (((Atmp1)) & 1) : ((((Atmp1)) >> ((16) - 1)) + 1) >> 1)))));
    tmp1 = ((((((tmp1)) + (((((Atmp1)) >> 16) * (int)((short)((Atmp1)))) + (((((Atmp1)) & 0x0000FFFF) * (int)((short)((Atmp1)))) >> 16))))) + ((((Atmp1)) * (((16) == 1 ? (((Atmp1)) >> 1) + (((Atmp1)) & 1) : ((((Atmp1)) >> ((16) - 1)) + 1) >> 1)))));
    A_Q16[k] = -Atmp1;
  }
  *res_nrg = ((((((nrg)) + ((((((int)((((long long)((WhiteNoiseFrac_Q32)) * ((C0)))) >> (32)))) >> 16) * (int)((short)((-tmp1)))) + ((((((int)((((long long)((WhiteNoiseFrac_Q32)) * ((C0)))) >> (32)))) & 0x0000FFFF) * (int)((short)((-tmp1)))) >> 16))))) + (((((int)((((long long)((WhiteNoiseFrac_Q32)) * ((C0)))) >> (32)))) * (((16) == 1 ? (((-tmp1)) >> 1) + (((-tmp1)) & 1) : ((((-tmp1)) >> ((16) - 1)) + 1) >> 1)))));
  *res_nrg_Q = -rshifts;
}
```