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

```asm
<SpreadApproximation>:
  endbr64
  push   %rbp
  mov    %edi,%ebp
  push   %rbx
  movq   %xmm0,%rbx
  sub    $0x48,%rsp
  ucomisd %xmm2,%xmm2
  movsd  %xmm1,0x20(%rsp)
  movsd  %xmm2,0x28(%rsp)
  movsd  %xmm3,(%rsp)
  movsd  %xmm4,0x8(%rsp)
  movsd  %xmm5,0x10(%rsp)
  movsd  %xmm6,0x18(%rsp)
  movsd  %xmm7,0x30(%rsp)
  jp     L53
  comisd D5(%rip),%xmm2
  jae    L79
L53:
  xor    %edi,%edi
L120:
  xor    %eax,%eax
  call   <assert@plt>
  movsd  (%rsp),%xmm3
  ucomisd %xmm3,%xmm3
  jp     L95
  comisd D99(%rip),%xmm3
  jae    L129
L95:
  xor    %edi,%edi
L36:
  xor    %eax,%eax
  call   <assert@plt>
  movsd  0x8(%rsp),%xmm4
  ucomisd %xmm4,%xmm4
  jp     L135
  pxor   %xmm6,%xmm6
  xor    %edi,%edi
  comisd %xmm6,%xmm4
  movsd  %xmm6,0x38(%rsp)
  jb     L105
  movsd  D91(%rip),%xmm0
  xor    %edi,%edi
  comisd %xmm4,%xmm0
  setae  %dil
L105:
  xor    %eax,%eax
  call   <assert@plt>
  movsd  0x10(%rsp),%xmm4
  ucomisd %xmm4,%xmm4
  jp     L90
  comisd D73(%rip),%xmm4
  jae    L148
L90:
  xor    %edi,%edi
L71:
  xor    %eax,%eax
  call   <assert@plt>
  movsd  0x18(%rsp),%xmm7
  ucomisd %xmm7,%xmm7
  jp     L7
  comisd D73(%rip),%xmm7
  jae    L158
L7:
  xor    %edi,%edi
L128:
  xor    %eax,%eax
  call   <assert@plt>
  movsd  0x30(%rsp),%xmm1
  movsd  0x10(%rsp),%xmm7
  movsd  0x18(%rsp),%xmm3
  movsd  0x20(%rsp),%xmm2
  addsd  %xmm1,%xmm1
  movsd  0x28(%rsp),%xmm6
  movapd %xmm7,%xmm5
  mulsd  %xmm7,%xmm5
  movapd %xmm3,%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm2,%xmm6
  mulsd  %xmm7,%xmm1
  divsd  %xmm6,%xmm0
  mulsd  %xmm3,%xmm1
  movsd  0x38(%rsp),%xmm3
  mulsd  %xmm2,%xmm1
  divsd  %xmm6,%xmm1
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm5
  subsd  %xmm1,%xmm5
  ucomisd %xmm5,%xmm3
  ja     L140
  sqrtsd %xmm5,%xmm5
L21:
  movq   %rbx,%xmm0
  movsd  (%rsp),%xmm2
  mov    %ebp,%edi
  pxor   %xmm4,%xmm4
  divsd  %xmm6,%xmm0
  movsd  0x8(%rsp),%xmm3
  mov    $0x6,%eax
  movsd  D8(%rip),%xmm1
  movsd  %xmm6,0x10(%rsp)
  call   <gbs@plt>
  pxor   %xmm0,%xmm0
  movsd  0x10(%rsp),%xmm6
  add    $0x48,%rsp
  cvtsi2sd %eax,%xmm0
  pop    %rbx
  pop    %rbp
  mulsd  %xmm6,%xmm0
  ret
  nop
L135:
  pxor   %xmm5,%xmm5
  xor    %edi,%edi
  movsd  %xmm5,0x38(%rsp)
  jmp    L105
  nopl   0x0(%rax)
L129:
  movsd  D42(%rip),%xmm0
  xor    %edi,%edi
  comisd %xmm3,%xmm0
  setae  %dil
  jmp    L36
  nop
L79:
  movsd  D62(%rip),%xmm0
  xor    %edi,%edi
  comisd %xmm2,%xmm0
  setae  %dil
  jmp    L120
  nopw   0x0(%rax,%rax,1)
L158:
  movsd  D43(%rip),%xmm0
  xor    %edi,%edi
  comisd %xmm7,%xmm0
  setae  %dil
  jmp    L128
  nopw   0x0(%rax,%rax,1)
L148:
  movsd  D43(%rip),%xmm0
  xor    %edi,%edi
  comisd %xmm4,%xmm0
  setae  %dil
  jmp    L71
L140:
  movapd %xmm5,%xmm0
  movsd  %xmm6,0x10(%rsp)
  call   <sqrt@plt>
  movsd  0x10(%rsp),%xmm6
  movapd %xmm0,%xmm5
  jmp    L21
```

tool_calls: name: read_data
arguments: label: D5
type: f64[2]
name: read_data
arguments: label: D99
type: byte[16]
name: read_data
arguments: label: D91
type: f64
name: read_data
arguments: label: D73
type: f64[2]
name: read_data
arguments: label: D8
type: f64
name: read_data
arguments: label: D42
type: f64
name: read_data
arguments: label: D62
type: f64
name: read_data
arguments: label: D43
type: f64

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

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

tool: D91: 1000

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

tool: D8: 1

tool: D42: 20

tool: D62: 1000000000000

tool: D43: 100.99999

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;
}
```