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

```asm
<biquad_lowpass>:
  endbr64
  push   %rbp
  mov    %rdi,%rbp
  sub    $0x20,%rsp
  movsd  D69(%rip),%xmm5
  comisd %xmm0,%xmm5
  ja     L19
  sub    $0x8,%rsp
  xor    %edx,%edx
  xor    %r9d,%r9d
  mov    $0x1,%r8d
  push   $0x0
  xor    %ecx,%ecx
  mov    $0x1,%esi
  xor    %eax,%eax
  call   <set_coefficient@plt>
  pop    %rax
  pop    %rdx
  add    $0x20,%rsp
  pop    %rbp
  ret
  nopl   0x0(%rax)
L19:
  pxor   %xmm2,%xmm2
  movapd %xmm0,%xmm3
  comisd %xmm0,%xmm2
  jbe    L28
L73:
  sub    $0x8,%rsp
  xor    %ecx,%ecx
  xor    %esi,%esi
  mov    %rbp,%rdi
  push   $0x0
  xor    %r9d,%r9d
  mov    $0x1,%r8d
  xor    %edx,%edx
  xor    %eax,%eax
  call   <set_coefficient@plt>
  pop    %rcx
  pop    %rsi
  add    $0x20,%rsp
  pop    %rbp
  ret
L28:
  comisd %xmm2,%xmm0
  jbe    L73
  comisd %xmm1,%xmm2
  ja     L76
  mulsd  D68(%rip),%xmm1
  movsd  D38(%rip),%xmm0
  movsd  %xmm3,0x8(%rsp)
  call   <pow@plt>
  pxor   %xmm2,%xmm2
  movsd  0x8(%rsp),%xmm3
  movapd %xmm0,%xmm1
  mulsd  %xmm0,%xmm1
  movsd  D6(%rip),%xmm0
  movapd %xmm0,%xmm4
  divsd  %xmm1,%xmm4
  subsd  %xmm4,%xmm0
  ucomisd %xmm0,%xmm2
  movapd %xmm0,%xmm1
  ja     L69
L21:
  sqrtsd %xmm1,%xmm1
L71:
  movsd  D80(%rip),%xmm0
  subsd  %xmm1,%xmm0
  movsd  D52(%rip),%xmm1
  mulsd  %xmm1,%xmm0
  ucomisd %xmm0,%xmm2
  ja     L6
  movapd %xmm0,%xmm2
  sqrtsd %xmm2,%xmm2
L34:
  lea    0x18(%rsp),%rdi
  mulsd  D58(%rip),%xmm3
  lea    0x10(%rsp),%rsi
  movsd  %xmm2,0x8(%rsp)
  movapd %xmm3,%xmm0
  call   <sincos@plt>
  mov    D52(%rip),%rax
  mov    %rbp,%rdi
  movsd  0x8(%rsp),%xmm2
  movsd  0x18(%rsp),%xmm0
  movsd  0x10(%rsp),%xmm3
  mov    $0x1,%esi
  movq   %rax,%xmm1
  movsd  D69(%rip),%xmm4
  mov    $0x5,%eax
  mulsd  %xmm1,%xmm2
  mulsd  %xmm2,%xmm0
  subsd  %xmm0,%xmm4
  addsd  D69(%rip),%xmm0
  mulsd  %xmm1,%xmm4
  divsd  %xmm0,%xmm4
  addsd  %xmm4,%xmm1
  addsd  %xmm4,%xmm4
  mulsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm1
  mulsd  D26(%rip),%xmm1
  xorpd  D122(%rip),%xmm3
  addsd  %xmm3,%xmm3
  movapd %xmm1,%xmm2
  addsd  %xmm1,%xmm2
  mulsd  D80(%rip),%xmm1
  add    $0x20,%rsp
  pop    %rbp
  movapd %xmm2,%xmm0
  jmp    1090 <set_coefficient@plt>
  nopl   0x0(%rax)
L76:
  movapd %xmm2,%xmm1
  jmp    L21
L69:
  movsd  %xmm3,0x8(%rsp)
  call   <sqrt@plt>
  movsd  0x8(%rsp),%xmm3
  pxor   %xmm2,%xmm2
  movapd %xmm0,%xmm1
  jmp    L71
L6:
  movsd  %xmm3,0x8(%rsp)
  call   <sqrt@plt>
  movsd  0x8(%rsp),%xmm3
  movapd %xmm0,%xmm2
  jmp    L34
```

tool_calls: name: unpack_data
arguments: data_label: D69
type_infered: f64
name: unpack_data
arguments: data_label: D68
type_infered: f64
name: unpack_data
arguments: data_label: D38
type_infered: f64
name: unpack_data
arguments: data_label: D6
type_infered: f64
name: unpack_data
arguments: data_label: D80
type_infered: f64
name: unpack_data
arguments: data_label: D52
type_infered: f64
name: unpack_data
arguments: data_label: D58
type_infered: f64
name: unpack_data
arguments: data_label: D26
type_infered: f64
name: unpack_data
arguments: data_label: D122
type_infered: f32[4]

tool: D69: 1

tool: D68: 0.05

tool: D38: 10

tool: D6: 16

tool: D80: 4

tool: D52: 0.5

tool: D58: 3.141592653589793

tool: D26: 0.25

tool: {"D122": [0.0,-0.0,0.0,0.0]}

assistant: ```c
void biquad_lowpass(struct biquad *bq, double cutoff, double resonance) {
  cutoff = ({
    __typeof__(0.0) _a = (0.0);
    __typeof__(({
      __typeof__(cutoff) _a = (cutoff);
      __typeof__(1.0) _b = (1.0);
      _a < _b ? _a : _b;
    })) _b = (({
      __typeof__(cutoff) _a = (cutoff);
      __typeof__(1.0) _b = (1.0);
      _a < _b ? _a : _b;
    }));
    _a > _b ? _a : _b;
  });
  if (cutoff == 1) {
    set_coefficient(bq, 1, 0, 0, 1, 0, 0);
  } else if (cutoff > 0) {
    resonance = ({
      __typeof__(0.0) _a = (0.0);
      __typeof__(resonance) _b = (resonance);
      _a > _b ? _a : _b;
    });
    double g = pow(10.0, 0.05 * resonance);
    double d = sqrt((4 - sqrt(16 - 16 / (g * g))) / 2);
    double theta = 3.14159265358979323846 * cutoff;
    double sn = 0.5 * d * sin(theta);
    double beta = 0.5 * (1 - sn) / (1 + sn);
    double gamma = (0.5 + beta) * cos(theta);
    double alpha = 0.25 * (0.5 + beta - gamma);
    double b0 = 2 * alpha;
    double b1 = 2 * 2 * alpha;
    double b2 = 2 * alpha;
    double a1 = 2 * -gamma;
    double a2 = 2 * beta;
    set_coefficient(bq, b0, b1, b2, 1, a1, a2);
  } else {
    set_coefficient(bq, 0, 0, 0, 1, 0, 0);
  }
}
```