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

```asm
<lfr_filter_calculate_f32>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  mov    %esi,%r12d
  push   %rbp
  push   %rbx
  mov    %edx,%ebx
  sub    $0x68,%rsp
  movsd  %xmm0,0x50(%rsp)
  movapd %xmm2,%xmm0
  movsd  %xmm1,0x8(%rsp)
  movsd  %xmm2,0x38(%rsp)
  mov    %rdi,0x40(%rsp)
  mov    %edx,0x4c(%rsp)
  call   <bessel_i0@plt>
  movsd  0x8(%rsp),%xmm1
  movapd %xmm1,%xmm2
  addsd  %xmm1,%xmm2
  mulsd  D29(%rip),%xmm1
  divsd  %xmm0,%xmm2
  movsd  %xmm1,0x30(%rsp)
  movsd  %xmm2,0x28(%rsp)
  test   %ebx,%ebx
  jle    L86
  lea    -0x1(%r12),%edx
  pxor   %xmm7,%xmm7
  mov    %edx,%eax
  shr    $0x1f,%eax
  add    %edx,%eax
  sar    %eax
  cvtsi2sd %eax,%xmm7
  movsd  %xmm7,0x58(%rsp)
  test   %r12d,%r12d
  jle    L86
  lea    -0x2(%r12),%eax
  pxor   %xmm0,%xmm0
  xor    %ebp,%ebp
  xor    %r14d,%r14d
  cvtsi2sd %eax,%xmm0
  movslq %r12d,%r13
  pxor   %xmm3,%xmm3
  movsd  D9(%rip),%xmm5
  movsd  D33(%rip),%xmm6
  movsd  D54(%rip),%xmm7
  divsd  %xmm0,%xmm5
  nopl   0x0(%rax,%rax,1)
L10:
  mov    0x40(%rsp),%rax
  pxor   %xmm4,%xmm4
  movslq %ebp,%rcx
  xor    %ebx,%ebx
  cvtsi2sd %r14d,%xmm4
  mulsd  0x50(%rsp),%xmm4
  addsd  0x58(%rsp),%xmm4
  lea    (%rax,%rcx,4),%r15
  jmp    L35
  nopl   (%rax)
L32:
  movss  %xmm0,(%r15,%rbx,4)
  add    $0x1,%rbx
  cmp    %r13,%rbx
  je     L43
L35:
  pxor   %xmm1,%xmm1
  movapd %xmm5,%xmm2
  movaps %xmm3,%xmm0
  cvtsi2sd %ebx,%xmm1
  subsd  %xmm4,%xmm1
  mulsd  %xmm1,%xmm2
  comisd %xmm2,%xmm6
  jae    L32
  comisd %xmm7,%xmm2
  jae    L32
  mulsd  %xmm2,%xmm2
  movapd %xmm7,%xmm0
  subsd  %xmm2,%xmm0
  pxor   %xmm2,%xmm2
  ucomisd %xmm0,%xmm2
  ja     L5
  sqrtsd %xmm0,%xmm0
L25:
  mulsd  0x38(%rsp),%xmm0
  movsd  %xmm4,0x18(%rsp)
  movsd  %xmm5,0x10(%rsp)
  movsd  %xmm1,0x20(%rsp)
  call   <bessel_i0@plt>
  movsd  0x20(%rsp),%xmm1
  mulsd  0x28(%rsp),%xmm0
  mulsd  0x30(%rsp),%xmm1
  movsd  %xmm0,0x8(%rsp)
  movapd %xmm1,%xmm0
  call   <sinc@plt>
  mulsd  0x8(%rsp),%xmm0
  pxor   %xmm3,%xmm3
  mov    D54(%rip),%rax
  movsd  0x18(%rsp),%xmm4
  movsd  0x10(%rsp),%xmm5
  movq   %rax,%xmm7
  mov    D33(%rip),%rax
  movq   %rax,%xmm6
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%r15,%rbx,4)
  add    $0x1,%rbx
  cmp    %r13,%rbx
  jne    L35
  nopw   0x0(%rax,%rax,1)
L43:
  add    $0x1,%r14d
  add    %r12d,%ebp
  cmp    %r14d,0x4c(%rsp)
  jne    L10
L86:
  add    $0x68,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L5:
  movsd  %xmm1,0x18(%rsp)
  movsd  %xmm4,0x10(%rsp)
  movsd  %xmm5,0x8(%rsp)
  call   <sqrt@plt>
  movsd  0x18(%rsp),%xmm1
  movsd  0x10(%rsp),%xmm4
  movsd  0x8(%rsp),%xmm5
  jmp    L25
```

tool_calls: name: parse_data
arguments: label: D29
type_guessed: byte[8]
name: parse_data
arguments: label: D9
type_guessed: f64
name: parse_data
arguments: label: D33
type_guessed: f32[2]
name: parse_data
arguments: label: D54
type_guessed: f64

tool: {"D29": [0x18,0x2d,0x44,0x54,0xfb,0x21,0x19,0x40]}

tool: D9: 2

tool: {"D33": [0.0,-1.875]}

tool: D54: 1

assistant: ```c
void lfr_filter_calculate_f32(float *data, int nsamp, int nfilt, double offset, double cutoff, double beta) {
  double x, x0, t, y, xscale, yscale;
  int i, j;
  yscale = 2.0 * cutoff / bessel_i0(beta);
  xscale = (8.0 * atan(1.0)) * cutoff;
  for (i = 0; i < nfilt; ++i) {
    x0 = (nsamp - 1) / 2 + offset * i;
    for (j = 0; j < nsamp; ++j) {
      x = j - x0;
      t = x * (2.0 / (nsamp - 2));
      if (t <= -1.0 || t >= 1.0) {
        y = 0.0;
      } else {
        y = yscale * bessel_i0(beta * sqrt(1.0 - t * t)) * sinc(xscale * x);
      }
      data[i * nsamp + j] = (float)y;
    }
  }
}
```