user: What is the c source code of the assembly code below:

```asm
<CriticalValueChooser>:
  endbr64
  push   %r15
  movapd %xmm0,%xmm8
  movq   %xmm4,%r15
  movapd %xmm7,%xmm4
  push   %r14
  movq   %xmm2,%r14
  push   %rbp
  movq   %xmm5,%rbp
  push   %rbx
  movq   %xmm1,%rbx
  sub    $0x38,%rsp
  ucomisd %xmm0,%xmm0
  movsd  %xmm3,(%rsp)
  movapd %xmm6,%xmm3
  jp     L42
  comisd D47(%rip),%xmm0
  jae    L67
L42:
  xor    %edi,%edi
L6:
  movsd  %xmm4,0x18(%rsp)
  movsd  %xmm3,0x10(%rsp)
  movsd  %xmm8,0x8(%rsp)
  call   <assert@plt>
  movsd  (%rsp),%xmm1
  movsd  0x8(%rsp),%xmm8
  movsd  0x10(%rsp),%xmm3
  movsd  0x18(%rsp),%xmm4
  ucomisd %xmm1,%xmm1
  jp     L64
  comisd D9(%rip),%xmm1
  jae    L50
L64:
  xor    %edi,%edi
L5:
  movsd  %xmm4,0x18(%rsp)
  movsd  %xmm3,0x10(%rsp)
  movsd  %xmm8,0x8(%rsp)
  call   <assert@plt>
  movq   %r15,%xmm2
  subsd  (%rsp),%xmm2
  movsd  0x8(%rsp),%xmm8
  movsd  0x60(%rsp),%xmm5
  movsd  0x18(%rsp),%xmm4
  movq   %rbx,%xmm1
  movsd  0x10(%rsp),%xmm3
  movapd %xmm8,%xmm0
  movq   %xmm2,%r15
  call   <gbs_call@plt>
  movq   %rbp,%xmm2
  subsd  (%rsp),%xmm2
  movsd  0x8(%rsp),%xmm8
  movsd  0x60(%rsp),%xmm5
  movsd  0x18(%rsp),%xmm4
  movq   %r14,%xmm1
  movsd  %xmm0,0x28(%rsp)
  movsd  0x10(%rsp),%xmm3
  movapd %xmm8,%xmm0
  movq   %xmm2,%rbp
  call   <gbs_put@plt>
  movsd  0x60(%rsp),%xmm5
  movq   %r15,%xmm2
  movsd  0x8(%rsp),%xmm8
  movsd  0x18(%rsp),%xmm4
  movsd  0x10(%rsp),%xmm3
  movsd  %xmm0,(%rsp)
  movq   %rbx,%xmm1
  movapd %xmm8,%xmm0
  call   <delta_call@plt>
  movsd  0x60(%rsp),%xmm5
  movq   %r14,%xmm1
  movsd  0x8(%rsp),%xmm8
  movsd  0x18(%rsp),%xmm4
  movsd  0x10(%rsp),%xmm3
  movsd  %xmm0,0x20(%rsp)
  movq   %rbp,%xmm2
  movapd %xmm8,%xmm0
  call   <delta_put@plt>
  movsd  0x20(%rsp),%xmm7
  movsd  0x28(%rsp),%xmm6
  subsd  (%rsp),%xmm6
  movapd %xmm7,%xmm1
  subsd  %xmm0,%xmm1
  jmp    L10
  nop
L55:
  divsd  %xmm1,%xmm6
  movq   %r15,%xmm2
  movsd  0x60(%rsp),%xmm5
  movq   %rbx,%xmm1
  movsd  %xmm4,0x10(%rsp)
  movsd  %xmm3,0x8(%rsp)
  subsd  %xmm6,%xmm8
  movapd %xmm8,%xmm0
  movsd  %xmm8,(%rsp)
  call   <gbs_call@plt>
  movsd  (%rsp),%xmm8
  movq   %rbp,%xmm2
  movsd  0x10(%rsp),%xmm4
  movsd  0x8(%rsp),%xmm3
  movsd  0x60(%rsp),%xmm5
  movsd  %xmm0,0x28(%rsp)
  movq   %r14,%xmm1
  movapd %xmm8,%xmm0
  movsd  %xmm4,0x18(%rsp)
  movsd  %xmm3,0x10(%rsp)
  movsd  %xmm8,0x8(%rsp)
  call   <gbs_put@plt>
  movsd  0x60(%rsp),%xmm5
  movq   %r15,%xmm2
  movsd  0x8(%rsp),%xmm8
  movsd  0x18(%rsp),%xmm4
  movsd  0x10(%rsp),%xmm3
  movsd  %xmm0,(%rsp)
  movq   %rbx,%xmm1
  movapd %xmm8,%xmm0
  call   <delta_call@plt>
  movsd  0x60(%rsp),%xmm5
  movq   %r14,%xmm1
  movsd  0x8(%rsp),%xmm8
  movsd  0x18(%rsp),%xmm4
  movsd  0x10(%rsp),%xmm3
  movsd  %xmm0,0x20(%rsp)
  movq   %rbp,%xmm2
  movapd %xmm8,%xmm0
  call   <delta_put@plt>
  movsd  0x20(%rsp),%xmm7
  movsd  0x28(%rsp),%xmm6
  subsd  (%rsp),%xmm6
  subsd  %xmm0,%xmm7
  movapd %xmm7,%xmm1
L10:
  movapd %xmm6,%xmm2
  andpd  D15(%rip),%xmm2
  movsd  0x8(%rsp),%xmm8
  comisd D29(%rip),%xmm2
  movsd  0x10(%rsp),%xmm3
  movsd  0x18(%rsp),%xmm4
  ja     L55
  add    $0x38,%rsp
  movapd %xmm8,%xmm0
  pop    %rbx
  pop    %rbp
  pop    %r14
  pop    %r15
  ret
  nopl   (%rax)
L50:
  movsd  D32(%rip),%xmm0
  xor    %edi,%edi
  comisd %xmm1,%xmm0
  setae  %dil
  jmp    L5
  nop
L67:
  movsd  D49(%rip),%xmm0
  xor    %edi,%edi
  comisd %xmm8,%xmm0
  setae  %dil
  jmp    L6
```

tool_calls: name: disassemble_data
arguments: label: D47
datatype: f64[2]
name: disassemble_data
arguments: label: D9
datatype: byte[16]
name: disassemble_data
arguments: label: D15
datatype: qword[2]
name: disassemble_data
arguments: label: D29
datatype: f64[2]
name: disassemble_data
arguments: label: D32
datatype: f64
name: disassemble_data
arguments: label: D49
datatype: f64

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

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

tool: {"D15": [0x7fffffffffffffff,0x0000000000000000]}

tool: {"D29": [0.001,0.0]}

tool: D32: 20

tool: D49: 1000000000000

assistant: ```c
double CriticalValueChooser(double S, double Xc, double Xp, double T, double Tc, double Tp, double r, double b, double v) {
  double Sv, ci, Pi, dc, dp, yi, di;
  const double epsilon = 0.001;
  assert(((S) == (S)) && (S) >= 0.000001 && (S) <= 1000000000000.0);
  assert(((T) == (T)) && (T) >= (1.0 / (253.0 * 24.0)) && (T) <= (20.0));
  Sv = S;
  ci = gbs_call(Sv, Xc, Tc - T, r, b, v);
  Pi = gbs_put(Sv, Xp, Tp - T, r, b, v);
  dc = delta_call(Sv, Xc, Tc - T, r, b, v);
  dp = delta_put(Sv, Xp, Tp - T, r, b, v);
  yi = ci - Pi;
  di = dc - dp;
  while (fabs(yi) > epsilon) {
    Sv = Sv - (yi) / di;
    ci = gbs_call(Sv, Xc, Tc - T, r, b, v);
    Pi = gbs_put(Sv, Xp, Tp - T, r, b, v);
    dc = delta_call(Sv, Xc, Tc - T, r, b, v);
    dp = delta_put(Sv, Xp, Tp - T, r, b, v);
    yi = ci - Pi;
    di = dc - dp;
  }
  return Sv;
}
```