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

```asm
<VasicekBondOption>:
  endbr64
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x40,%rsp
  mov    %edi,%r12d
  movsd  %xmm0,0x20(%rsp)
  movsd  %xmm1,0x28(%rsp)
  movq   %xmm2,%rbx
  movsd  %xmm3,0x8(%rsp)
  movsd  %xmm4,0x18(%rsp)
  movq   %xmm5,%rbp
  movsd  %xmm6,(%rsp)
  movsd  %xmm7,0x10(%rsp)
  ucomisd %xmm1,%xmm1
  mov    $0x0,%edi
  jp     L36
  mov    $0x0,%edi
  comisd D137(%rip),%xmm1
  jb     L36
  movsd  D20(%rip),%xmm0
  comisd %xmm1,%xmm0
  setae  %dil
  movzbl %dil,%edi
L36:
  mov    $0x0,%eax
  call   <assert@plt>
  movsd  0x8(%rsp),%xmm2
  ucomisd %xmm2,%xmm2
  mov    $0x0,%edi
  jp     L96
  mov    $0x0,%edi
  comisd D39(%rip),%xmm2
  jb     L96
  movsd  D102(%rip),%xmm0
  comisd %xmm2,%xmm0
  setae  %dil
  movzbl %dil,%edi
L96:
  mov    $0x0,%eax
  call   <assert@plt>
  movsd  0x18(%rsp),%xmm6
  ucomisd %xmm6,%xmm6
  mov    $0x0,%edi
  jp     L106
  mov    $0x0,%edi
  comisd D24(%rip),%xmm6
  jb     L106
  movsd  D25(%rip),%xmm0
  comisd %xmm6,%xmm0
  setae  %dil
  movzbl %dil,%edi
L106:
  mov    $0x0,%eax
  call   <assert@plt>
  movsd  0x10(%rsp),%xmm4
  ucomisd %xmm4,%xmm4
  mov    $0x0,%edi
  jp     L0
  mov    $0x0,%edi
  comisd D158(%rip),%xmm4
  jb     L0
  movsd  D72(%rip),%xmm0
  comisd %xmm4,%xmm0
  setae  %dil
  movzbl %dil,%edi
L0:
  mov    $0x0,%eax
  call   <assert@plt>
  movsd  0x28(%rsp),%xmm2
  divsd  0x20(%rsp),%xmm2
  movsd  %xmm2,0x30(%rsp)
  movsd  0x10(%rsp),%xmm5
  movsd  (%rsp),%xmm4
  movq   %rbp,%xmm3
  movsd  0x18(%rsp),%xmm2
  movsd  0x8(%rsp),%xmm1
  pxor   %xmm0,%xmm0
  mov    $0x6,%eax
  call   <VasicekBondPrice@plt>
  pxor   %xmm7,%xmm7
  cvtsi2sd %eax,%xmm7
  movsd  %xmm7,0x38(%rsp)
  movsd  0x10(%rsp),%xmm5
  movsd  (%rsp),%xmm4
  movq   %rbp,%xmm3
  movsd  0x18(%rsp),%xmm2
  movq   %rbx,%xmm1
  pxor   %xmm0,%xmm0
  mov    $0x6,%eax
  call   <VasicekBondPrice@plt>
  pxor   %xmm3,%xmm3
  cvtsi2sd %eax,%xmm3
  movsd  %xmm3,0x18(%rsp)
  movsd  (%rsp),%xmm7
  mulsd  D14(%rip),%xmm7
  movsd  %xmm7,0x28(%rsp)
  mulsd  0x8(%rsp),%xmm7
  movapd %xmm7,%xmm0
  call   <exp@plt>
  movapd %xmm0,%xmm2
  movsd  0x10(%rsp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  D56(%rip),%xmm1
  subsd  %xmm2,%xmm1
  mulsd  %xmm1,%xmm0
  xorpd  D108(%rip),%xmm0
  divsd  0x28(%rsp),%xmm0
  pxor   %xmm5,%xmm5
  ucomisd %xmm0,%xmm5
  ja     L13
  sqrtsd %xmm0,%xmm0
  movsd  %xmm0,0x10(%rsp)
L89:
  movsd  (%rsp),%xmm0
  xorpd  D108(%rip),%xmm0
  movq   %rbx,%xmm1
  subsd  0x8(%rsp),%xmm1
  mulsd  %xmm1,%xmm0
  call   <exp@plt>
  movsd  D56(%rip),%xmm5
  subsd  %xmm0,%xmm5
  movapd %xmm5,%xmm0
  mulsd  0x10(%rsp),%xmm0
  divsd  (%rsp),%xmm0
  movsd  %xmm0,(%rsp)
  movsd  0x30(%rsp),%xmm3
  mulsd  0x38(%rsp),%xmm3
  movsd  0x18(%rsp),%xmm5
  movsd  %xmm3,0x10(%rsp)
  divsd  %xmm3,%xmm5
  movapd %xmm5,%xmm0
  call   <log@plt>
  movapd %xmm0,%xmm1
  movsd  D56(%rip),%xmm4
  movsd  (%rsp),%xmm6
  divsd  %xmm6,%xmm4
  movapd %xmm4,%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm6,%xmm1
  mulsd  D29(%rip),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,0x8(%rsp)
  test   %r12d,%r12d
  je     L54
  mov    $0x1,%eax
  call   <cnd@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movapd %xmm0,%xmm7
  mulsd  0x18(%rsp),%xmm7
  movq   %xmm7,%rbx
  movsd  0x8(%rsp),%xmm0
  subsd  (%rsp),%xmm0
  mov    $0x1,%eax
  call   <cnd@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  mulsd  0x10(%rsp),%xmm1
  movq   %rbx,%xmm0
  subsd  %xmm1,%xmm0
  mulsd  0x20(%rsp),%xmm0
L38:
  add    $0x40,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  ret
L13:
  call   <sqrt@plt>
  movsd  %xmm0,0x10(%rsp)
  jmp    L89
L54:
  movsd  (%rsp),%xmm2
  subsd  0x8(%rsp),%xmm2
  movapd %xmm2,%xmm0
  mov    $0x1,%eax
  call   <cnd@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movapd %xmm0,%xmm3
  mulsd  0x10(%rsp),%xmm3
  movq   %xmm3,%rbx
  movsd  0x8(%rsp),%xmm0
  xorpd  D108(%rip),%xmm0
  mov    $0x1,%eax
  call   <cnd@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  mulsd  0x18(%rsp),%xmm1
  movq   %rbx,%xmm0
  subsd  %xmm1,%xmm0
  mulsd  0x20(%rsp),%xmm0
  jmp    L38
```

tool_calls: name: extract_symbolic_data
arguments: data_tag: D137
type_infered: f64[2]
name: extract_symbolic_data
arguments: data_tag: D20
type_infered: f64
name: extract_symbolic_data
arguments: data_tag: D39
type_infered: byte[16]
name: extract_symbolic_data
arguments: data_tag: D102
type_infered: f64
name: extract_symbolic_data
arguments: data_tag: D24
type_infered: f64
name: extract_symbolic_data
arguments: data_tag: D25
type_infered: f64
name: extract_symbolic_data
arguments: data_tag: D158
type_infered: f64[2]
name: extract_symbolic_data
arguments: data_tag: D72
type_infered: f64
name: extract_symbolic_data
arguments: data_tag: D14
type_infered: byte[8]
name: extract_symbolic_data
arguments: data_tag: D56
type_infered: f64
name: extract_symbolic_data
arguments: data_tag: D108
type_infered: byte[16]
name: extract_symbolic_data
arguments: data_tag: D29
type_infered: byte[8]

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

tool: D20: 1000000000000

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

tool: D102: 20

tool: D24: 0

tool: D25: 1000

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

tool: D72: 100.99999

tool: {"D14": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0xc0]}

tool: D56: 1

tool: {"D108": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00]}

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

assistant: ```c
double VasicekBondOption(int fCall, double F, double X, double tau, double T, double r, double Theta, double kappa, double v) {
  double PtT, Pt_tau, vp, H, result;
  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(((v) == (v)) && (v) >= 0.0000001 && v <= 100.99999);
  X = X / F;
  PtT = VasicekBondPrice(0.0, T, r, Theta, kappa, v);
  Pt_tau = VasicekBondPrice(0.0, tau, r, Theta, kappa, v);
  vp = sqrt(((v) * (v)) * (1.0 - exp(-2.0 * kappa * T)) / (2.0 * kappa)) * (1.0 - exp(-kappa * (tau - T))) / kappa;
  H = 1.0 / vp * log(Pt_tau / (PtT * X)) + vp / 2.0;
  if (fCall)
    result = F * (Pt_tau * cnd(H) - X * PtT * cnd(H - vp));
  else
    result = F * (X * PtT * cnd(-H + vp) - Pt_tau * cnd(-H));
  return result;
}
```