user: 请显示与下面列出的汇编代码相匹配的C语言实现:

```asm
<filter_fft>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x90,%rsp
  mov    %edi,-0x74(%rbp)
  mov    %esi,-0x78(%rbp)
  mov    %rdx,-0x80(%rbp)
  mov    %rcx,-0x88(%rbp)
  mov    -0x78(%rbp),%eax
  mov    $0x1,%edx
  mov    %eax,%ecx
  shl    %cl,%edx
  mov    %edx,%eax
  cltq
  mov    %rax,-0x20(%rbp)
  mov    -0x20(%rbp),%rax
  sar    %rax
  mov    %rax,-0x18(%rbp)
  movq   $0x0,-0x40(%rbp)
  movq   $0x0,-0x48(%rbp)
  jmp    L10
L99:
  mov    -0x48(%rbp),%rax
  cmp    -0x40(%rbp),%rax
  jge    L97
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  movss  %xmm0,-0x50(%rbp)
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  movss  %xmm0,-0x4c(%rbp)
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x80(%rbp),%rax
  add    %rax,%rdx
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,4),%rcx
  mov    -0x80(%rbp),%rax
  add    %rcx,%rax
  movss  (%rdx),%xmm0
  movss  %xmm0,(%rax)
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x88(%rbp),%rax
  add    %rax,%rdx
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,4),%rcx
  mov    -0x88(%rbp),%rax
  add    %rcx,%rax
  movss  (%rdx),%xmm0
  movss  %xmm0,(%rax)
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movss  -0x50(%rbp),%xmm0
  movss  %xmm0,(%rax)
  mov    -0x40(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  movss  -0x4c(%rbp),%xmm0
  movss  %xmm0,(%rax)
L97:
  mov    -0x18(%rbp),%rax
  mov    %rax,-0x38(%rbp)
  jmp    L103
L3:
  mov    -0x38(%rbp),%rax
  sub    %rax,-0x40(%rbp)
  sarq   -0x38(%rbp)
L103:
  mov    -0x38(%rbp),%rax
  cmp    -0x40(%rbp),%rax
  jle    L3
  mov    -0x38(%rbp),%rax
  add    %rax,-0x40(%rbp)
  addq   $0x1,-0x48(%rbp)
L10:
  mov    -0x20(%rbp),%rax
  sub    $0x1,%rax
  cmp    %rax,-0x48(%rbp)
  jl     L99
  movss  D2(%rip),%xmm0
  movss  %xmm0,-0x6c(%rbp)
  pxor   %xmm0,%xmm0
  movss  %xmm0,-0x68(%rbp)
  movq   $0x1,-0x28(%rbp)
  movq   $0x0,-0x30(%rbp)
  jmp    L6
L44:
  mov    -0x28(%rbp),%rax
  mov    %rax,-0x10(%rbp)
  shlq   -0x28(%rbp)
  movss  D32(%rip),%xmm0
  movss  %xmm0,-0x64(%rbp)
  pxor   %xmm0,%xmm0
  movss  %xmm0,-0x60(%rbp)
  movq   $0x0,-0x40(%rbp)
  jmp    L102
L72:
  mov    -0x40(%rbp),%rax
  mov    %rax,-0x48(%rbp)
  jmp    L21
L93:
  mov    -0x48(%rbp),%rdx
  mov    -0x10(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x8(%rbp)
  mov    -0x8(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mulss  -0x64(%rbp),%xmm0
  mov    -0x8(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  mulss  -0x60(%rbp),%xmm1
  subss  %xmm1,%xmm0
  movss  %xmm0,-0x58(%rbp)
  mov    -0x8(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  movaps %xmm0,%xmm1
  mulss  -0x64(%rbp),%xmm1
  mov    -0x8(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mulss  -0x60(%rbp),%xmm0
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x54(%rbp)
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mov    -0x8(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  subss  -0x58(%rbp),%xmm0
  movss  %xmm0,(%rax)
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mov    -0x8(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  subss  -0x54(%rbp),%xmm0
  movss  %xmm0,(%rax)
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  addss  -0x58(%rbp),%xmm0
  movss  %xmm0,(%rax)
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  addss  -0x54(%rbp),%xmm0
  movss  %xmm0,(%rax)
  mov    -0x28(%rbp),%rax
  add    %rax,-0x48(%rbp)
L21:
  mov    -0x48(%rbp),%rax
  cmp    -0x20(%rbp),%rax
  jl     L93
  movss  -0x64(%rbp),%xmm0
  mulss  -0x6c(%rbp),%xmm0
  movss  -0x60(%rbp),%xmm1
  mulss  -0x68(%rbp),%xmm1
  subss  %xmm1,%xmm0
  movss  %xmm0,-0x5c(%rbp)
  movss  -0x64(%rbp),%xmm0
  movaps %xmm0,%xmm1
  mulss  -0x68(%rbp),%xmm1
  movss  -0x60(%rbp),%xmm0
  mulss  -0x6c(%rbp),%xmm0
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x60(%rbp)
  movss  -0x5c(%rbp),%xmm0
  movss  %xmm0,-0x64(%rbp)
  addq   $0x1,-0x40(%rbp)
L102:
  mov    -0x40(%rbp),%rax
  cmp    -0x10(%rbp),%rax
  jl     L72
  pxor   %xmm1,%xmm1
  cvtss2sd -0x6c(%rbp),%xmm1
  movsd  D6(%rip),%xmm0
  subsd  %xmm1,%xmm0
  movsd  D23(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movd   %xmm0,%eax
  mov    %eax,-0x68(%rbp)
  cmpl   $0x0,-0x74(%rbp)
  je     L50
  movss  -0x68(%rbp),%xmm0
  movss  D53(%rip),%xmm1
  xorps  %xmm1,%xmm0
  movss  %xmm0,-0x68(%rbp)
L50:
  pxor   %xmm1,%xmm1
  cvtss2sd -0x6c(%rbp),%xmm1
  movsd  D6(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  D23(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movd   %xmm0,%eax
  mov    %eax,-0x6c(%rbp)
  addq   $0x1,-0x30(%rbp)
L6:
  mov    -0x78(%rbp),%eax
  cltq
  cmp    %rax,-0x30(%rbp)
  jl     L44
  cmpl   $0x0,-0x74(%rbp)
  je     L117
  movq   $0x0,-0x48(%rbp)
  jmp    L89
L96:
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2ssq -0x20(%rbp),%xmm1
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  divss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2ssq -0x20(%rbp),%xmm1
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  divss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  addq   $0x1,-0x48(%rbp)
L89:
  mov    -0x48(%rbp),%rax
  cmp    -0x20(%rbp),%rax
  jl     L96
L117:
  nop
  leave
  ret
```

tool_calls: name: disassemble_data
arguments: data_label: D2
type_infered: byte[4]
name: disassemble_data
arguments: data_label: D32
type_infered: f32
name: disassemble_data
arguments: data_label: D6
type_infered: f32[2]
name: disassemble_data
arguments: data_label: D23
type_infered: f64
name: disassemble_data
arguments: data_label: D53
type_infered: byte[4]

tool: {"D2": [0x00,0x00,0x80,0xbf]}

tool: D32: 1

tool: {"D6": [0.0,1.875]}

tool: D23: 2

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

assistant: ```c
void filter_fft(int dir, int m, float *real, float *imag) {
  long n, i, i1, j, k, i2, l, l1, l2;
  float c1, c2, tx, ty, t1, t2, u1, u2, z;
  n = 1 << m;
  i2 = n >> 1;
  j = 0;
  for (i = 0; i < n - 1; i++) {
    if (i < j) {
      tx = real[i];
      ty = imag[i];
      real[i] = real[j];
      imag[i] = imag[j];
      real[j] = tx;
      imag[j] = ty;
    }
    k = i2;
    while (k <= j) {
      j -= k;
      k >>= 1;
    }
    j += k;
  }
  c1 = -1.0;
  c2 = 0.0;
  l2 = 1;
  for (l = 0; l < m; l++) {
    l1 = l2;
    l2 <<= 1;
    u1 = 1.0;
    u2 = 0.0;
    for (j = 0; j < l1; j++) {
      for (i = j; i < n; i += l2) {
        i1 = i + l1;
        t1 = u1 * real[i1] - u2 * imag[i1];
        t2 = u1 * imag[i1] + u2 * real[i1];
        real[i1] = real[i] - t1;
        imag[i1] = imag[i] - t2;
        real[i] += t1;
        imag[i] += t2;
      }
      z = u1 * c1 - u2 * c2;
      u2 = u1 * c2 + u2 * c1;
      u1 = z;
    }
    c2 = sqrt((1.0 - c1) / 2.0);
    if (dir) {
      c2 = -c2;
    }
    c1 = sqrt((1.0 + c1) / 2.0);
  }
  if (dir) {
    for (i = 0; i < n; i++) {
      real[i] /= n;
      imag[i] /= n;
    }
  }
}
```