user: What would the C source look like for the assembly code provided here?

```asm
<lfr_filter_calculate_f32>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  mov    %esi,%ebp
  push   %rbx
  mov    %edx,%ebx
  sub    $0x68,%rsp
  movsd  %xmm0,0x40(%rsp)
  movapd %xmm2,%xmm0
  movsd  %xmm1,0x8(%rsp)
  movsd  %xmm2,0x38(%rsp)
  mov    %rdi,0x50(%rsp)
  mov    %edx,0x58(%rsp)
  call   <bessel_i0@plt>
  movsd  0x8(%rsp),%xmm1
  movapd %xmm1,%xmm2
  addsd  %xmm1,%xmm2
  mulsd  D35(%rip),%xmm1
  divsd  %xmm0,%xmm2
  movsd  %xmm1,0x30(%rsp)
  movsd  %xmm2,0x28(%rsp)
  test   %ebx,%ebx
  jle    L96
  lea    -0x1(%rbp),%edx
  pxor   %xmm6,%xmm6
  xor    %r13d,%r13d
  xor    %r12d,%r12d
  mov    %edx,%eax
  mov    %edx,0x5c(%rsp)
  pxor   %xmm5,%xmm5
  shr    $0x1f,%eax
  add    %edx,%eax
  sar    %eax
  cvtsi2sd %eax,%xmm6
  movsd  %xmm6,0x48(%rsp)
  movsd  D12(%rip),%xmm6
  nopl   0x0(%rax,%rax,1)
L77:
  pxor   %xmm3,%xmm3
  cvtsi2sd %r12d,%xmm3
  mulsd  0x40(%rsp),%xmm3
  addsd  0x48(%rsp),%xmm3
  test   %ebp,%ebp
  jle    L24
  lea    -0x2(%rbp),%eax
  pxor   %xmm0,%xmm0
  mov    0x50(%rsp),%rcx
  xor    %ebx,%ebx
  cvtsi2sd %eax,%xmm0
  movslq %r13d,%rax
  movsd  D18(%rip),%xmm4
  mov    0x5c(%rsp),%r15d
  movsd  D48(%rip),%xmm7
  lea    (%rcx,%rax,4),%r14
  divsd  %xmm0,%xmm4
  jmp    L23
  nopl   0x0(%rax)
L61:
  movss  %xmm0,(%r14,%rbx,4)
  lea    0x1(%rbx),%rax
  cmp    %rbx,%r15
  je     L24
L46:
  mov    %rax,%rbx
L23:
  pxor   %xmm1,%xmm1
  movapd %xmm4,%xmm2
  movaps %xmm5,%xmm0
  cvtsi2sd %ebx,%xmm1
  subsd  %xmm3,%xmm1
  mulsd  %xmm1,%xmm2
  comisd %xmm2,%xmm6
  jae    L61
  comisd %xmm7,%xmm2
  pxor   %xmm0,%xmm0
  jae    L61
  mulsd  %xmm2,%xmm2
  movapd %xmm7,%xmm0
  subsd  %xmm2,%xmm0
  pxor   %xmm2,%xmm2
  ucomisd %xmm0,%xmm2
  ja     L65
  sqrtsd %xmm0,%xmm0
L31:
  mulsd  0x38(%rsp),%xmm0
  movsd  %xmm3,0x18(%rsp)
  movsd  %xmm4,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   %xmm5,%xmm5
  mov    D12(%rip),%rax
  movsd  0x18(%rsp),%xmm3
  movsd  0x10(%rsp),%xmm4
  movq   %rax,%xmm6
  mov    D48(%rip),%rax
  movq   %rax,%xmm7
  lea    0x1(%rbx),%rax
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%r14,%rbx,4)
  cmp    %rbx,%r15
  jne    L46
  cs nopw 0x0(%rax,%rax,1)
L24:
  add    $0x1,%r12d
  add    %ebp,%r13d
  cmp    %r12d,0x58(%rsp)
  jne    L77
L96:
  add    $0x68,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L65:
  movsd  %xmm1,0x18(%rsp)
  movsd  %xmm3,0x10(%rsp)
  movsd  %xmm4,0x8(%rsp)
  call   <sqrt@plt>
  movsd  0x18(%rsp),%xmm1
  movsd  0x10(%rsp),%xmm3
  movsd  0x8(%rsp),%xmm4
  jmp    L31
```

tool_calls: name: read_data
arguments: data_tag: D35
data_type: byte[8]
name: read_data
arguments: data_tag: D12
data_type: f32[2]
name: read_data
arguments: data_tag: D18
data_type: f64
name: read_data
arguments: data_tag: D48
data_type: f64

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

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

tool: D18: 2

tool: D48: 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;
    }
  }
}
```