user: 请从给定的汇编指令推导出C源代码：

```asm
<lsptopc>:
  endbr64
  push   %r14
  mov    $0x18,%edx
  push   %r13
  push   %r12
  push   %rbp
  mov    %rsi,%rbp
  xor    %esi,%esi
  push   %rbx
  xor    %ebx,%ebx
  sub    $0x180,%rsp
  movdqu (%rdi),%xmm2
  movdqu 0x10(%rdi),%xmm3
  mov    %fs:0x28,%rax
  mov    %rax,0x178(%rsp)
  xor    %eax,%eax
  mov    0x20(%rdi),%rax
  lea    0x50(%rsp),%r14
  lea    0x90(%rsp),%rdi
  movaps %xmm2,0x150(%rsp)
  lea    0x150(%rsp),%r12
  lea    0x70(%rsp),%r13
  movaps %xmm3,0x160(%rsp)
  mov    %rax,0x170(%rsp)
  call   <memset@plt>
  lea    0xb0(%rsp),%rdi
  mov    $0x18,%edx
  xor    %esi,%esi
  call   <memset@plt>
  lea    0xd0(%rsp),%rdi
  mov    $0x18,%edx
  xor    %esi,%esi
  call   <memset@plt>
  lea    0xf0(%rsp),%rdi
  mov    $0x18,%edx
  xor    %esi,%esi
  call   <memset@plt>
  lea    0x110(%rsp),%rdi
  mov    $0x18,%edx
  xor    %esi,%esi
  call   <memset@plt>
  lea    0x130(%rsp),%rdi
  mov    $0x18,%edx
  xor    %esi,%esi
  call   <memset@plt>
L30:
  movss  D58(%rip),%xmm0
  mulss  (%r12,%rbx,2),%xmm0
  call   <cos@plt>
  pxor   %xmm0,%xmm0
  cvtsi2ss %rax,%xmm0
  mulss  D41(%rip),%xmm0
  movss  %xmm0,(%r14,%rbx,1)
  movss  D58(%rip),%xmm0
  mulss  0x4(%r12,%rbx,2),%xmm0
  call   <cos@plt>
  pxor   %xmm0,%xmm0
  cvtsi2ss %rax,%xmm0
  mulss  D41(%rip),%xmm0
  movss  %xmm0,0x0(%r13,%rbx,1)
  add    $0x4,%rbx
  cmp    $0x14,%rbx
  jne    L30
  movss  0x74(%rsp),%xmm7
  movss  0x70(%rsp),%xmm6
  xor    %eax,%eax
  movss  0xb0(%rsp),%xmm13
  movss  0x50(%rsp),%xmm5
  movss  0xd0(%rsp),%xmm3
  movss  %xmm7,0x34(%rsp)
  movss  0xb8(%rsp),%xmm7
  movss  0x110(%rsp),%xmm12
  movss  %xmm6,0x2c(%rsp)
  movss  0xb4(%rsp),%xmm6
  movss  0x130(%rsp),%xmm2
  movss  %xmm7,0x18(%rsp)
  movss  0x118(%rsp),%xmm7
  movss  0xd4(%rsp),%xmm1
  movss  %xmm5,0x28(%rsp)
  movss  0x54(%rsp),%xmm5
  movss  0x134(%rsp),%xmm0
  movss  %xmm6,0x14(%rsp)
  movss  0x58(%rsp),%xmm6
  movss  %xmm7,0x10(%rsp)
  movss  0x5c(%rsp),%xmm7
  movss  %xmm5,0x30(%rsp)
  movss  0x114(%rsp),%xmm5
  movss  %xmm6,0x38(%rsp)
  movss  0x78(%rsp),%xmm6
  movss  %xmm7,0x40(%rsp)
  movss  0xbc(%rsp),%xmm7
  movss  %xmm5,0xc(%rsp)
  movss  0xd8(%rsp),%xmm5
  movss  %xmm6,0x3c(%rsp)
  movss  0x138(%rsp),%xmm6
  movss  %xmm7,0x8(%rsp)
  movss  0xdc(%rsp),%xmm8
  movss  0x7c(%rsp),%xmm4
  movss  0x13c(%rsp),%xmm10
  movss  0x60(%rsp),%xmm7
  movss  0xe0(%rsp),%xmm14
  movss  %xmm4,0x44(%rsp)
  movss  0xc0(%rsp),%xmm4
  movss  0x11c(%rsp),%xmm11
  movss  %xmm7,0x48(%rsp)
  movss  0x80(%rsp),%xmm7
  movss  D19(%rip),%xmm9
  movss  %xmm4,0x20(%rsp)
  movss  0x120(%rsp),%xmm4
  movss  %xmm7,0x4c(%rsp)
  movss  0x140(%rsp),%xmm7
  movss  %xmm4,0x1c(%rsp)
  pxor   %xmm4,%xmm4
  movss  %xmm7,0x24(%rsp)
  movaps %xmm6,%xmm7
  movaps %xmm8,%xmm6
  movaps %xmm10,%xmm8
  movaps %xmm14,%xmm10
  jmp    L44
  nopw   0x0(%rax,%rax,1)
L16:
  movss  0x1c(%rsp),%xmm10
  movaps %xmm11,%xmm8
  movss  %xmm6,0x1c(%rsp)
  movaps %xmm5,%xmm11
  movss  0x8(%rsp),%xmm6
  movss  0x18(%rsp),%xmm5
  add    $0x1,%rax
  movss  %xmm0,0x8(%rsp)
  movss  %xmm10,0x24(%rsp)
  movss  0xc(%rsp),%xmm0
  movss  0x20(%rsp),%xmm10
  movss  %xmm2,0x18(%rsp)
  movaps %xmm12,%xmm2
  movaps %xmm15,%xmm12
  movss  %xmm7,0x20(%rsp)
  movss  0x10(%rsp),%xmm7
  movss  %xmm3,0xc(%rsp)
  movaps %xmm13,%xmm3
  movaps %xmm14,%xmm13
  movss  %xmm1,0x10(%rsp)
  movss  0x14(%rsp),%xmm1
  movss  %xmm4,0x14(%rsp)
  movaps %xmm9,%xmm4
  pxor   %xmm9,%xmm9
L44:
  movaps %xmm9,%xmm14
  movaps %xmm9,%xmm15
  addss  %xmm4,%xmm14
  subss  %xmm4,%xmm15
  movss  0x28(%rsp),%xmm4
  mulss  %xmm13,%xmm4
  addss  %xmm14,%xmm4
  addss  %xmm3,%xmm4
  movss  0x2c(%rsp),%xmm3
  mulss  %xmm12,%xmm3
  addss  %xmm15,%xmm3
  addss  %xmm2,%xmm3
  movss  0x14(%rsp),%xmm2
  mulss  0x30(%rsp),%xmm2
  addss  %xmm4,%xmm2
  addss  %xmm1,%xmm2
  movss  0xc(%rsp),%xmm1
  mulss  0x34(%rsp),%xmm1
  addss  %xmm3,%xmm1
  addss  %xmm0,%xmm1
  movss  0x18(%rsp),%xmm0
  mulss  0x38(%rsp),%xmm0
  addss  %xmm2,%xmm0
  addss  %xmm5,%xmm0
  movss  0x10(%rsp),%xmm5
  mulss  0x3c(%rsp),%xmm5
  addss  %xmm1,%xmm5
  addss  %xmm7,%xmm5
  movss  0x8(%rsp),%xmm7
  mulss  0x40(%rsp),%xmm7
  addss  %xmm0,%xmm7
  addss  %xmm6,%xmm7
  movss  0x44(%rsp),%xmm6
  mulss  %xmm11,%xmm6
  addss  %xmm5,%xmm6
  addss  %xmm8,%xmm6
  test   %rax,%rax
  je     L16
  movss  0x20(%rsp),%xmm8
  mulss  0x48(%rsp),%xmm8
  addss  %xmm7,%xmm8
  addss  %xmm10,%xmm8
  movss  0x1c(%rsp),%xmm10
  mulss  0x4c(%rsp),%xmm10
  addss  %xmm6,%xmm10
  addss  0x24(%rsp),%xmm10
  addss  %xmm10,%xmm8
  mulss  D59(%rip),%xmm8
  movss  %xmm8,-0x4(%rbp,%rax,4)
  cmp    $0xa,%rax
  jne    L16
  movups 0x18(%rbp),%xmm0
  movaps D26(%rip),%xmm1
  movq   D26(%rip),%xmm7
  shufps $0x1b,%xmm0,%xmm0
  xorps  %xmm1,%xmm0
  shufps $0x1b,%xmm0,%xmm0
  movups %xmm0,0x1c(%rbp)
  movups 0x8(%rbp),%xmm0
  shufps $0x1b,%xmm0,%xmm0
  xorps  %xmm1,%xmm0
  shufps $0x1b,%xmm0,%xmm0
  movups %xmm0,0xc(%rbp)
  movq   0x0(%rbp),%xmm0
  movl   $0x3f800000,0x0(%rbp)
  shufps $0xe1,%xmm0,%xmm0
  xorps  %xmm7,%xmm0
  shufps $0xe1,%xmm0,%xmm0
  movlps %xmm0,0x4(%rbp)
  mov    0x178(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L31
  add    $0x180,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
L31:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: extract_symbolic_data
arguments: label: D58
type: byte[4]
name: extract_symbolic_data
arguments: label: D41
type: byte[4]
name: extract_symbolic_data
arguments: label: D19
type: byte[4]
name: extract_symbolic_data
arguments: label: D59
type: byte[4]
name: extract_symbolic_data
arguments: label: D26
type: byte[16]

tool: {"D58": [0xdb,0x0f,0xc9,0x40]}

tool: {"D41": [0x00,0x00,0x00,0xc0]}

tool: {"D19": [0x00,0x00,0x80,0x3f]}

tool: {"D59": [0x00,0x00,0x00,0xbf]}

tool: {"D26": [0x00,0x00,0x00,0x80,0x00,0x00,0x00,0x80,0x00,0x00,0x00,0x80,0x00,0x00,0x00,0x80]}

assistant: ```c
void lsptopc(float *f, float *pc) {
  int i, j, k, lspflag;
  float freq[10], p[10 / 2], q[10 / 2];
  float a[10 / 2 + 1], a1[10 / 2 + 1], a2[10 / 2 + 1];
  float b[10 / 2 + 1], b1[10 / 2 + 1], b2[10 / 2 + 1];
  float xx, xf;
  lspflag = 0;
  for (i = 1; i < 10; i++) {
    if (f[i] <= f[i - 1])
      lspflag = 1;
  }
  for (j = 0; j < 10; j++)
    freq[j] = f[j];
  memset(a, 0, sizeof(a[0]) * (10 / 2 + 1));
  memset(a1, 0, sizeof(a[0]) * (10 / 2 + 1));
  memset(a2, 0, sizeof(a[0]) * (10 / 2 + 1));
  memset(b, 0, sizeof(a[0]) * (10 / 2 + 1));
  memset(b1, 0, sizeof(a[0]) * (10 / 2 + 1));
  memset(b2, 0, sizeof(a[0]) * (10 / 2 + 1));
  for (i = 0; i < (10 / 2); i++) {
    p[i] = -2.f * (float)cos(2.f * 3.14159265358979323846f * freq[2 * i]);
    q[i] = -2.f * (float)cos(2.f * 3.14159265358979323846f * freq[2 * i + 1]);
  }
  xf = 0.0f;
  for (k = 0; k < 10 + 1; k++) {
    xx = 0.0;
    if (k == 0)
      xx = 1.0f;
    a[0] = xx + xf;
    b[0] = xx - xf;
    xf = xx;
    for (i = 0; i < (10 / 2); i++) {
      a[i + 1] = a[i] + p[i] * a1[i] + a2[i];
      b[i + 1] = b[i] + q[i] * b1[i] + b2[i];
      a2[i] = a1[i];
      a1[i] = a[i];
      b2[i] = b1[i];
      b1[i] = b[i];
    }
    if (k != 0)
      pc[k - 1] = -.5f * (a[10 / 2] + b[10 / 2]);
  }
  for (i = 10 - 1; i >= 0; i--)
    pc[i + 1] = -pc[i];
  pc[0] = 1.0f;
}
```