user: Translate the below assembly code into C programming language:

```asm
<quadratic_spline>:
  endbr64
  push   %rbp
  push   %rbx
  sub    $0x48,%rsp
  movq   %xmm0,%rbx
  movq   %xmm1,%rbp
  movsd  %xmm2,0x18(%rsp)
  movsd  %xmm3,0x20(%rsp)
  movsd  %xmm4,0x28(%rsp)
  movsd  %xmm5,0x30(%rsp)
  movsd  %xmm6,0x10(%rsp)
  movsd  %xmm7,0x38(%rsp)
  movapd %xmm2,%xmm1
  addsd  %xmm4,%xmm1
  movsd  D45(%rip),%xmm0
  movapd %xmm1,%xmm6
  mulsd  %xmm0,%xmm6
  movsd  %xmm6,(%rsp)
  addsd  %xmm5,%xmm3
  movapd %xmm3,%xmm1
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,0x8(%rsp)
  movq   %rbx,%xmm0
  subsd  %xmm6,%xmm0
  call   <fabs@plt>
  test   %eax,%eax
  jle    L40
L50:
  movsd  D65(%rip),%xmm0
  movsd  0x20(%rsp),%xmm7
  movapd %xmm7,%xmm5
  mulsd  %xmm0,%xmm5
  addsd  0x30(%rsp),%xmm5
  movsd  D48(%rip),%xmm4
  movsd  0x18(%rsp),%xmm6
  mulsd  %xmm6,%xmm0
  addsd  0x28(%rsp),%xmm0
  movq   %rbp,%xmm3
  addsd  %xmm7,%xmm3
  movsd  D45(%rip),%xmm1
  movq   %rbx,%xmm2
  addsd  %xmm6,%xmm2
  movsd  0x8(%rsp),%xmm7
  movsd  (%rsp),%xmm6
  mulsd  %xmm4,%xmm5
  mulsd  %xmm0,%xmm4
  mulsd  %xmm1,%xmm3
  mulsd  %xmm1,%xmm2
  movq   %rbp,%xmm1
  movq   %rbx,%xmm0
  call   <quadratic_spline>
L56:
  movsd  (%rsp),%xmm0
  subsd  0x10(%rsp),%xmm0
  call   <fabs@plt>
  test   %eax,%eax
  jle    L6
L13:
  movsd  0x38(%rsp),%xmm7
  movapd %xmm7,%xmm5
  movsd  0x30(%rsp),%xmm3
  addsd  %xmm3,%xmm5
  movsd  D45(%rip),%xmm1
  movsd  0x10(%rsp),%xmm6
  movapd %xmm6,%xmm4
  movsd  0x28(%rsp),%xmm8
  addsd  %xmm8,%xmm4
  movsd  D65(%rip),%xmm2
  mulsd  %xmm2,%xmm3
  addsd  0x20(%rsp),%xmm3
  movsd  D48(%rip),%xmm0
  mulsd  %xmm8,%xmm2
  addsd  0x18(%rsp),%xmm2
  mulsd  %xmm1,%xmm5
  mulsd  %xmm1,%xmm4
  mulsd  %xmm0,%xmm3
  mulsd  %xmm0,%xmm2
  movsd  0x8(%rsp),%xmm1
  movsd  (%rsp),%xmm0
  call   <quadratic_spline>
L33:
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  ret
L40:
  movq   %rbp,%xmm0
  subsd  0x8(%rsp),%xmm0
  call   <fabs@plt>
  test   %eax,%eax
  jg     L50
  mov    $0x0,%edi
  movsd  0x8(%rsp),%xmm1
  movsd  (%rsp),%xmm0
  call   <splinerel@plt>
  jmp    L56
L6:
  movsd  0x8(%rsp),%xmm0
  subsd  0x38(%rsp),%xmm0
  call   <fabs@plt>
  test   %eax,%eax
  jg     L13
  mov    $0x0,%edi
  movsd  0x38(%rsp),%xmm1
  movsd  0x10(%rsp),%xmm0
  call   <splinerel@plt>
  jmp    L33
```

tool_calls: name: unpack_data
arguments: label: D45
data_type: byte[8]
name: unpack_data
arguments: label: D65
data_type: f64
name: unpack_data
arguments: label: D48
data_type: byte[8]

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

tool: D65: 3

tool: {"D48": [0x00,0x00,0x00,0x00,0x00,0x00,0xd0,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);
  }
}
```