user: 请将以下的汇编代码转换成C语言代码:

```asm
<biquad_lowpass>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0xb0,%rsp
  mov    %rdi,-0x98(%rbp)
  movsd  %xmm0,-0xa0(%rbp)
  movsd  %xmm1,-0xa8(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x90(%rbp)
  movsd  -0xa0(%rbp),%xmm0
  movsd  %xmm0,-0x88(%rbp)
  movsd  D69(%rip),%xmm0
  movsd  %xmm0,-0x80(%rbp)
  movsd  -0x80(%rbp),%xmm0
  comisd -0x88(%rbp),%xmm0
  jbe    L1
  movsd  -0x88(%rbp),%xmm0
  jmp    L72
L1:
  movsd  -0x80(%rbp),%xmm0
L72:
  movsd  %xmm0,-0x78(%rbp)
  movsd  -0x90(%rbp),%xmm0
  comisd -0x78(%rbp),%xmm0
  jbe    L69
  movsd  -0x90(%rbp),%xmm0
  jmp    L90
L69:
  movsd  -0x78(%rbp),%xmm0
L90:
  movsd  %xmm0,-0xa0(%rbp)
  movsd  D69(%rip),%xmm0
  ucomisd -0xa0(%rbp),%xmm0
  jp     L87
  movsd  D69(%rip),%xmm0
  ucomisd -0xa0(%rbp),%xmm0
  jne    L87
  mov    -0x98(%rbp),%rax
  sub    $0x8,%rsp
  push   $0x0
  mov    $0x0,%r9d
  mov    $0x1,%r8d
  mov    $0x0,%ecx
  mov    $0x0,%edx
  mov    $0x1,%esi
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <set_coefficient@plt>
  add    $0x10,%rsp
  jmp    L42
L87:
  movsd  -0xa0(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm1,%xmm0
  jbe    L28
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x70(%rbp)
  movsd  -0xa8(%rbp),%xmm0
  movsd  %xmm0,-0x68(%rbp)
  movsd  -0x70(%rbp),%xmm0
  comisd -0x68(%rbp),%xmm0
  jbe    L50
  movsd  -0x70(%rbp),%xmm0
  jmp    L109
L50:
  movsd  -0x68(%rbp),%xmm0
L109:
  movsd  %xmm0,-0xa8(%rbp)
  movsd  -0xa8(%rbp),%xmm1
  movsd  D25(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  mov    D180(%rip),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <pow@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x60(%rbp)
  movsd  -0x60(%rbp),%xmm0
  movapd %xmm0,%xmm2
  mulsd  %xmm0,%xmm2
  movsd  D46(%rip),%xmm0
  movapd %xmm0,%xmm1
  divsd  %xmm2,%xmm1
  movsd  D46(%rip),%xmm0
  subsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movapd %xmm0,%xmm1
  movsd  D67(%rip),%xmm0
  subsd  %xmm1,%xmm0
  movsd  D160(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x58(%rbp)
  movsd  -0xa0(%rbp),%xmm1
  movsd  D75(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  -0x58(%rbp),%xmm1
  movsd  D26(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,-0xb0(%rbp)
  mov    -0x50(%rbp),%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  mulsd  -0xb0(%rbp),%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  D69(%rip),%xmm0
  movapd %xmm0,%xmm1
  subsd  -0x48(%rbp),%xmm1
  movsd  D26(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x48(%rbp),%xmm2
  movsd  D69(%rip),%xmm1
  addsd  %xmm2,%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  movsd  -0x40(%rbp),%xmm1
  movsd  D26(%rip),%xmm0
  addsd  %xmm0,%xmm1
  movsd  %xmm1,-0xb0(%rbp)
  mov    -0x50(%rbp),%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  mulsd  -0xb0(%rbp),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x40(%rbp),%xmm1
  movsd  D26(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  subsd  -0x38(%rbp),%xmm1
  movsd  D142(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0x30(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x30(%rbp),%xmm1
  movsd  D67(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x30(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x38(%rbp),%xmm0
  movq   D181(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  addsd  %xmm0,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x40(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x8(%rbp),%xmm3
  movsd  -0x10(%rbp),%xmm2
  movsd  -0x18(%rbp),%xmm1
  movsd  -0x20(%rbp),%xmm0
  mov    -0x28(%rbp),%rdx
  mov    -0x98(%rbp),%rax
  movapd %xmm3,%xmm4
  movapd %xmm2,%xmm3
  mov    $0x1,%esi
  movapd %xmm1,%xmm2
  movapd %xmm0,%xmm1
  movq   %rdx,%xmm0
  mov    %rax,%rdi
  mov    $0x5,%eax
  call   <set_coefficient@plt>
  jmp    L42
L28:
  mov    -0x98(%rbp),%rax
  sub    $0x8,%rsp
  push   $0x0
  mov    $0x0,%r9d
  mov    $0x1,%r8d
  mov    $0x0,%ecx
  mov    $0x0,%edx
  mov    $0x0,%esi
  mov    %rax,%rdi
  mov    $0x0,%eax
  call   <set_coefficient@plt>
  add    $0x10,%rsp
L42:
  nop
  leave
  ret
```

tool_calls: name: retrieve_data
arguments: data_tag: D69
type: f64
name: retrieve_data
arguments: data_tag: D25
type: f64
name: retrieve_data
arguments: data_tag: D180
type: f64
name: retrieve_data
arguments: data_tag: D46
type: f64
name: retrieve_data
arguments: data_tag: D67
type: f64
name: retrieve_data
arguments: data_tag: D160
type: f64
name: retrieve_data
arguments: data_tag: D75
type: f64
name: retrieve_data
arguments: data_tag: D26
type: f64
name: retrieve_data
arguments: data_tag: D142
type: f64
name: retrieve_data
arguments: data_tag: D181
type: f32[2]

tool: D69: 1

tool: D25: 0.05

tool: D180: 10

tool: D46: 16

tool: D67: 4

tool: D160: 2

tool: D75: 3.141592653589793

tool: D26: 0.5

tool: D142: 0.25

tool: {"D181": [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);
  }
}
```