user: 下面的汇编指令对应的C源代码是什么:

```asm
<xfir_gen_bandstop>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r12
  push   %rbx
  sub    $0x50,%rsp
  mov    %rdi,-0x48(%rbp)
  mov    %rsi,-0x50(%rbp)
  movss  %xmm0,-0x54(%rbp)
  movss  %xmm1,-0x58(%rbp)
  mov    %edx,-0x5c(%rbp)
  movss  %xmm2,-0x60(%rbp)
  pxor   %xmm0,%xmm0
  cvtss2sd -0x54(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  pxor   %xmm0,%xmm0
  cvtss2sd -0x58(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0x50(%rbp),%rax
  sub    $0x1,%rax
  test   %rax,%rax
  js     L6
  pxor   %xmm0,%xmm0
  cvtsi2sd %rax,%xmm0
  jmp    L21
L6:
  mov    %rax,%rdx
  shr    %rdx
  and    $0x1,%eax
  or     %rax,%rdx
  pxor   %xmm0,%xmm0
  cvtsi2sd %rdx,%xmm0
  addsd  %xmm0,%xmm0
L21:
  movsd  D21(%rip),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movq   $0x0,-0x30(%rbp)
  jmp    L48
L43:
  movsd  -0x38(%rbp),%xmm1
  movsd  D30(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x18(%rbp),%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movq   %xmm0,%rbx
  movsd  -0x20(%rbp),%xmm0
  mulsd  -0x18(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movq   %xmm0,%r12
  movsd  -0x28(%rbp),%xmm0
  mulsd  -0x18(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movq   %r12,%xmm1
  subsd  %xmm0,%xmm1
  movq   %rbx,%xmm0
  subsd  %xmm1,%xmm0
  divsd  -0x18(%rbp),%xmm0
  mov    -0x30(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rax)
  movsd  -0x38(%rbp),%xmm1
  movsd  D32(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  addq   $0x1,-0x30(%rbp)
L48:
  mov    -0x30(%rbp),%rax
  cmp    -0x50(%rbp),%rax
  jb     L43
  mov    -0x50(%rbp),%rax
  and    $0x1,%eax
  test   %rax,%rax
  je     L25
  movss  -0x58(%rbp),%xmm0
  subss  -0x54(%rbp),%xmm0
  cvtss2sd %xmm0,%xmm0
  movapd %xmm0,%xmm1
  addsd  %xmm0,%xmm1
  movsd  D32(%rip),%xmm0
  subsd  %xmm1,%xmm0
  mov    -0x50(%rbp),%rax
  sub    $0x1,%rax
  shr    %rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rax)
L25:
  mov    -0x60(%rbp),%esi
  mov    -0x5c(%rbp),%edx
  mov    -0x50(%rbp),%rcx
  mov    -0x48(%rbp),%rax
  movd   %esi,%xmm0
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <window@plt>
  nop
  add    $0x50,%rsp
  pop    %rbx
  pop    %r12
  pop    %rbp
  ret
```

tool_calls: name: disassemble_data
arguments: label: D21
type: byte[8]
name: disassemble_data
arguments: label: D30
type: f64
name: disassemble_data
arguments: label: D32
type: f64

tool: {"D21": [0x00,0x00,0x00,0x00,0x00,0x00,0xe0,0xbf]}

tool: D30: 3.141592653589793

tool: D32: 1

assistant: ```c
void xfir_gen_bandstop(float taps[], size_t n_taps, float low_freq, float high_freq, int win_type, float alpha) {
  const double tlf = 2.0 * low_freq;
  const double thf = 2.0 * high_freq;
  double n = -0.5 * (n_taps - 1);
  size_t i;
  for (i = 0; i < n_taps; i++) {
    const double x = (3.14159265358979323846) * n;
    taps[i] = (sin(x) - (sin(thf * x) - sin(tlf * x))) / x;
    n += 1.0;
  }
  if (n_taps & 1) {
    taps[(n_taps - 1) / 2] = 1.0 - (2.0 * (high_freq - low_freq));
  }
  window(taps, n_taps, win_type, alpha);
}
```