user: 请从给定的汇编指令推导出C源代码：

```asm
<quadratic_spline>:
  endbr64
  sub    $0x58,%rsp
  movapd %xmm0,%xmm8
  movapd %xmm2,%xmm0
  movsd  %xmm2,0x8(%rsp)
  movsd  %xmm3,0x10(%rsp)
  movsd  %xmm4,0x18(%rsp)
  movsd  %xmm5,0x20(%rsp)
  movsd  %xmm6,0x30(%rsp)
  movsd  %xmm7,0x48(%rsp)
  jmp    L11
  nopw   0x0(%rax,%rax,1)
L48:
  movsd  0x10(%rsp),%xmm7
  movsd  0x8(%rsp),%xmm2
  movapd %xmm8,%xmm0
  movsd  D14(%rip),%xmm5
  movsd  D14(%rip),%xmm4
  movapd %xmm7,%xmm3
  movsd  (%rsp),%xmm6
  mulsd  %xmm7,%xmm5
  addsd  %xmm1,%xmm3
  movsd  0x28(%rsp),%xmm7
  mulsd  %xmm2,%xmm4
  addsd  0x20(%rsp),%xmm5
  addsd  %xmm8,%xmm2
  addsd  0x18(%rsp),%xmm4
  mulsd  D33(%rip),%xmm5
  mulsd  D33(%rip),%xmm4
  mulsd  D138(%rip),%xmm3
  mulsd  D138(%rip),%xmm2
  call   <quadratic_spline>
L2:
  movsd  (%rsp),%xmm0
  subsd  0x30(%rsp),%xmm0
  call   <fabs@plt>
  test   %eax,%eax
  jle    L43
L44:
  movsd  0x20(%rsp),%xmm4
  movsd  D14(%rip),%xmm0
  movsd  0x18(%rsp),%xmm3
  movsd  0x30(%rsp),%xmm1
  mulsd  %xmm4,%xmm0
  movsd  0x48(%rsp),%xmm2
  movsd  (%rsp),%xmm8
  addsd  0x10(%rsp),%xmm0
  mulsd  D33(%rip),%xmm0
  addsd  %xmm3,%xmm1
  mulsd  D138(%rip),%xmm1
  addsd  %xmm4,%xmm2
  mulsd  D138(%rip),%xmm2
  movsd  %xmm0,0x10(%rsp)
  movsd  D14(%rip),%xmm0
  movsd  %xmm1,0x18(%rsp)
  movsd  0x28(%rsp),%xmm1
  mulsd  %xmm3,%xmm0
  movsd  %xmm2,0x20(%rsp)
  addsd  0x8(%rsp),%xmm0
  mulsd  D33(%rip),%xmm0
  movsd  %xmm0,0x8(%rsp)
L11:
  addsd  0x18(%rsp),%xmm0
  movsd  D138(%rip),%xmm6
  movapd %xmm8,%xmm5
  movsd  %xmm1,0x40(%rsp)
  movsd  %xmm8,0x38(%rsp)
  mulsd  %xmm0,%xmm6
  movsd  0x10(%rsp),%xmm0
  addsd  0x20(%rsp),%xmm0
  mulsd  D138(%rip),%xmm0
  subsd  %xmm6,%xmm5
  movsd  %xmm6,(%rsp)
  movsd  %xmm0,0x28(%rsp)
  movapd %xmm5,%xmm0
  call   <fabs@plt>
  movsd  0x38(%rsp),%xmm8
  movsd  0x40(%rsp),%xmm1
  test   %eax,%eax
  jg     L48
  movapd %xmm1,%xmm7
  subsd  0x28(%rsp),%xmm7
  movsd  %xmm8,0x40(%rsp)
  movsd  %xmm1,0x38(%rsp)
  movapd %xmm7,%xmm0
  call   <fabs@plt>
  movsd  0x38(%rsp),%xmm1
  movsd  0x40(%rsp),%xmm8
  test   %eax,%eax
  jg     L48
  movsd  0x28(%rsp),%xmm1
  movsd  (%rsp),%xmm0
  xor    %edi,%edi
  call   <splinerel@plt>
  jmp    L2
  nopl   0x0(%rax,%rax,1)
L43:
  movsd  0x28(%rsp),%xmm0
  subsd  0x48(%rsp),%xmm0
  call   <fabs@plt>
  test   %eax,%eax
  jg     L44
  movsd  0x48(%rsp),%xmm1
  movsd  0x30(%rsp),%xmm0
  xor    %edi,%edi
  add    $0x58,%rsp
  jmp    1070 <splinerel@plt>
```

tool_calls: name: unpack_data
arguments: data_label: D14
data_type: f64
name: unpack_data
arguments: data_label: D33
data_type: byte[8]
name: unpack_data
arguments: data_label: D138
data_type: byte[8]

tool: D14: 3

tool: {"D33": [0x00,0x00,0x00,0x00,0x00,0x00,0xd0,0x3f]}

tool: {"D138": [0x00,0x00,0x00,0x00,0x00,0x00,0xe0,0x3f]}

assistant: ```c
void quadratic_spline(double a_1, double b_1, double a_2, double b_2, double a_3, double b_3, double a_4, double b_4) {
  double x_1, y_1, x_4, y_4;
  double x_mid, y_mid;
  x_1 = a_1;
  y_1 = b_1;
  x_4 = a_4;
  y_4 = b_4;
  x_mid = (a_2 + a_3) / 2.0;
  y_mid = (b_2 + b_3) / 2.0;
  if ((fabs(x_1 - x_mid) < 1) && (fabs(y_1 - y_mid) < 1)) {
    splinerel(x_mid, y_mid, 0);
  } else {
    quadratic_spline(x_1, y_1, ((x_1 + a_2) / 2.0), ((y_1 + b_2) / 2.0), ((3.0 * a_2 + a_3) / 4.0), ((3.0 * b_2 + b_3) / 4.0), x_mid, y_mid);
  }
  if ((fabs(x_mid - x_4) < 1) && (fabs(y_mid - y_4) < 1)) {
    splinerel(x_4, y_4, 0);
  } else {
    quadratic_spline(x_mid, y_mid, ((a_2 + 3.0 * a_3) / 4.0), ((b_2 + 3.0 * b_3) / 4.0), ((a_3 + x_4) / 2.0), ((b_3 + y_4) / 2.0), x_4, y_4);
  }
}
```