user: Please show me the C implementation that matches the assembly code listed below:

```asm
<SpreadApproximation>:
  endbr64
  push   %rbp
  push   %rbx
  sub    $0x48,%rsp
  mov    %edi,%ebp
  movq   %xmm0,%rbx
  movsd  %xmm1,0x8(%rsp)
  movsd  %xmm2,0x30(%rsp)
  movsd  %xmm3,0x10(%rsp)
  movsd  %xmm4,0x18(%rsp)
  movsd  %xmm5,0x20(%rsp)
  movsd  %xmm6,0x28(%rsp)
  movsd  %xmm7,0x38(%rsp)
  ucomisd %xmm2,%xmm2
  mov    $0x0,%edi
  jp     L76
  mov    $0x0,%edi
  comisd D85(%rip),%xmm2
  jb     L76
  movsd  D67(%rip),%xmm0
  comisd %xmm2,%xmm0
  setae  %dil
  movzbl %dil,%edi
L76:
  mov    $0x0,%eax
  call   <assert@plt>
  movsd  0x10(%rsp),%xmm3
  ucomisd %xmm3,%xmm3
  mov    $0x0,%edi
  jp     L59
  mov    $0x0,%edi
  comisd D35(%rip),%xmm3
  jb     L59
  movsd  D31(%rip),%xmm0
  comisd %xmm3,%xmm0
  setae  %dil
  movzbl %dil,%edi
L59:
  mov    $0x0,%eax
  call   <assert@plt>
  movsd  0x18(%rsp),%xmm4
  ucomisd %xmm4,%xmm4
  mov    $0x0,%edi
  jp     L52
  mov    $0x0,%edi
  comisd D42(%rip),%xmm4
  jb     L52
  movsd  D10(%rip),%xmm0
  comisd %xmm4,%xmm0
  setae  %dil
  movzbl %dil,%edi
L52:
  mov    $0x0,%eax
  call   <assert@plt>
  movsd  0x20(%rsp),%xmm7
  ucomisd %xmm7,%xmm7
  mov    $0x0,%edi
  jp     L31
  mov    $0x0,%edi
  comisd D75(%rip),%xmm7
  jb     L31
  movsd  D51(%rip),%xmm0
  comisd %xmm7,%xmm0
  setae  %dil
  movzbl %dil,%edi
L31:
  mov    $0x0,%eax
  call   <assert@plt>
  movsd  0x28(%rsp),%xmm4
  ucomisd %xmm4,%xmm4
  mov    $0x0,%edi
  jp     L48
  mov    $0x0,%edi
  comisd D75(%rip),%xmm4
  jb     L48
  movsd  D51(%rip),%xmm0
  comisd %xmm4,%xmm0
  setae  %dil
  movzbl %dil,%edi
L48:
  mov    $0x0,%eax
  call   <assert@plt>
  movsd  0x30(%rsp),%xmm6
  movsd  0x8(%rsp),%xmm2
  addsd  %xmm2,%xmm6
  movsd  %xmm6,0x8(%rsp)
  movsd  0x28(%rsp),%xmm3
  movapd %xmm3,%xmm0
  mulsd  %xmm2,%xmm0
  divsd  %xmm6,%xmm0
  movsd  0x20(%rsp),%xmm7
  movapd %xmm7,%xmm5
  mulsd  %xmm7,%xmm5
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm5
  movsd  0x38(%rsp),%xmm1
  addsd  %xmm1,%xmm1
  mulsd  %xmm7,%xmm1
  mulsd  %xmm3,%xmm1
  mulsd  %xmm2,%xmm1
  divsd  %xmm6,%xmm1
  subsd  %xmm1,%xmm5
  pxor   %xmm0,%xmm0
  ucomisd %xmm5,%xmm0
  ja     L107
  sqrtsd %xmm5,%xmm5
L94:
  movq   %rbx,%xmm0
  divsd  0x8(%rsp),%xmm0
  pxor   %xmm4,%xmm4
  movsd  0x18(%rsp),%xmm3
  movsd  0x10(%rsp),%xmm2
  movsd  D37(%rip),%xmm1
  mov    %ebp,%edi
  mov    $0x6,%eax
  call   <gbs@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  mulsd  0x8(%rsp),%xmm0
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  ret
L107:
  movapd %xmm5,%xmm0
  call   <sqrt@plt>
  movapd %xmm0,%xmm5
  jmp    L94
```

tool_calls: name: disassemble_data
arguments: data_label: D85
type_guessed: f64[2]
name: disassemble_data
arguments: data_label: D67
type_guessed: f64
name: disassemble_data
arguments: data_label: D35
type_guessed: byte[16]
name: disassemble_data
arguments: data_label: D31
type_guessed: f64
name: disassemble_data
arguments: data_label: D42
type_guessed: f64[2]
name: disassemble_data
arguments: data_label: D10
type_guessed: f64
name: disassemble_data
arguments: data_label: D75
type_guessed: f64[2]
name: disassemble_data
arguments: data_label: D51
type_guessed: f64
name: disassemble_data
arguments: data_label: D37
type_guessed: f64

tool: {"D85": [1e-6,1000000000000.0]}

tool: D67: 1000000000000

tool: {"D35": [0xce,0xd2,0x29,0x9c,0x17,0x96,0x25,0x3f,0x00,0x00,0x00,0x00,0x00,0x00,0x34,0x40]}

tool: D31: 20

tool: {"D42": [0.0,1000.0]}

tool: D10: 1000

tool: {"D75": [1e-7,100.99999]}

tool: D51: 100.99999

tool: D37: 1

assistant: ```c
double SpreadApproximation(int fCall, double f1, double f2, double X, double T, double r, double v1, double v2, double Rho) {
  double f2X, v, F;
  assert(((X) == (X)) && (X) >= 0.000001 && (X) <= 1000000000000.0);
  assert(((T) == (T)) && (T) >= (1.0 / (253.0 * 24.0)) && (T) <= (20.0));
  assert(((r) == (r)) && (r) >= 0.0 && (r) <= 1000.0);
  assert(((v1) == (v1)) && (v1) >= 0.0000001 && v1 <= 100.99999);
  assert(((v2) == (v2)) && (v2) >= 0.0000001 && v2 <= 100.99999);
  f2X = f2 + X;
  v = sqrt(((v1) * (v1)) + ((v2 * f2 / f2X) * (v2 * f2 / f2X)) - 2.0 * Rho * v1 * v2 * f2 / f2X);
  F = f1 / f2X;
  return gbs(fCall, F, 1.0, T, r, 0.0, v) * f2X;
}
```