user: 请将这段汇编代码片段转换成等效的C代码：

```asm
<blackscholes_put>:
  endbr64
  sub    $0x38,%rsp
  movsd  %xmm0,0x20(%rsp)
  movsd  %xmm1,0x8(%rsp)
  movsd  %xmm2,(%rsp)
  movsd  %xmm3,0x28(%rsp)
  movsd  %xmm4,0x10(%rsp)
  ucomisd %xmm0,%xmm0
  mov    $0x0,%edi
  jp     L102
  movapd %xmm0,%xmm4
  mov    $0x0,%edi
  comisd D116(%rip),%xmm0
  jb     L102
  movsd  D33(%rip),%xmm0
  comisd %xmm4,%xmm0
  setae  %dil
  movzbl %dil,%edi
L102:
  mov    $0x0,%eax
  call   <assert@plt>
  movsd  0x8(%rsp),%xmm6
  ucomisd %xmm6,%xmm6
  mov    $0x0,%edi
  jp     L45
  mov    $0x0,%edi
  comisd D116(%rip),%xmm6
  jb     L45
  movsd  D33(%rip),%xmm0
  comisd %xmm6,%xmm0
  setae  %dil
  movzbl %dil,%edi
L45:
  mov    $0x0,%eax
  call   <assert@plt>
  movsd  (%rsp),%xmm2
  ucomisd %xmm2,%xmm2
  mov    $0x0,%edi
  jp     L15
  mov    $0x0,%edi
  comisd D83(%rip),%xmm2
  jb     L15
  movsd  D123(%rip),%xmm0
  comisd %xmm2,%xmm0
  setae  %dil
  movzbl %dil,%edi
L15:
  mov    $0x0,%eax
  call   <assert@plt>
  movsd  0x28(%rsp),%xmm7
  ucomisd %xmm7,%xmm7
  mov    $0x0,%edi
  jp     L79
  mov    $0x0,%edi
  comisd D29(%rip),%xmm7
  jb     L79
  movsd  D7(%rip),%xmm0
  comisd %xmm7,%xmm0
  setae  %dil
  movzbl %dil,%edi
L79:
  mov    $0x0,%eax
  call   <assert@plt>
  movsd  0x10(%rsp),%xmm5
  ucomisd %xmm5,%xmm5
  mov    $0x0,%edi
  jp     L2
  mov    $0x0,%edi
  comisd D46(%rip),%xmm5
  jb     L2
  movsd  D127(%rip),%xmm0
  comisd %xmm5,%xmm0
  setae  %dil
  movzbl %dil,%edi
L2:
  mov    $0x0,%eax
  call   <assert@plt>
  pxor   %xmm0,%xmm0
  movsd  (%rsp),%xmm3
  ucomisd %xmm3,%xmm0
  ja     L60
  sqrtsd %xmm3,%xmm3
  movapd %xmm3,%xmm0
L96:
  movapd %xmm0,%xmm6
  mulsd  0x10(%rsp),%xmm6
  movsd  %xmm6,0x18(%rsp)
  movsd  0x20(%rsp),%xmm7
  divsd  0x8(%rsp),%xmm7
  movapd %xmm7,%xmm0
  call   <log@plt>
  movapd %xmm0,%xmm1
  movsd  0x10(%rsp),%xmm0
  mulsd  %xmm0,%xmm0
  mulsd  D95(%rip),%xmm0
  movsd  0x28(%rsp),%xmm5
  addsd  %xmm5,%xmm0
  movsd  (%rsp),%xmm4
  mulsd  %xmm4,%xmm0
  addsd  %xmm1,%xmm0
  divsd  0x18(%rsp),%xmm0
  movsd  %xmm0,(%rsp)
  xorpd  D113(%rip),%xmm5
  movapd %xmm5,%xmm0
  mulsd  %xmm4,%xmm0
  call   <exp@plt>
  movapd %xmm0,%xmm1
  mulsd  0x8(%rsp),%xmm1
  movsd  %xmm1,0x8(%rsp)
  movsd  (%rsp),%xmm0
  subsd  0x18(%rsp),%xmm0
  xorpd  D113(%rip),%xmm0
  mov    $0x1,%eax
  call   <cnd@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movapd %xmm0,%xmm2
  mulsd  0x8(%rsp),%xmm2
  movsd  %xmm2,0x8(%rsp)
  movsd  (%rsp),%xmm0
  xorpd  D113(%rip),%xmm0
  mov    $0x1,%eax
  call   <cnd@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  mulsd  0x20(%rsp),%xmm1
  movsd  0x8(%rsp),%xmm0
  subsd  %xmm1,%xmm0
  add    $0x38,%rsp
  ret
L60:
  movsd  (%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L96
```

tool_calls: name: read_data
arguments: data_label: D116
type_infered: f64[2]
name: read_data
arguments: data_label: D33
type_infered: f64
name: read_data
arguments: data_label: D83
type_infered: byte[16]
name: read_data
arguments: data_label: D123
type_infered: f64
name: read_data
arguments: data_label: D29
type_infered: f64
name: read_data
arguments: data_label: D7
type_infered: f64
name: read_data
arguments: data_label: D46
type_infered: f64[2]
name: read_data
arguments: data_label: D127
type_infered: f64
name: read_data
arguments: data_label: D95
type_infered: byte[8]
name: read_data
arguments: data_label: D113
type_infered: byte[16]

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

tool: D33: 1000000000000

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

tool: D123: 20

tool: D29: 0

tool: D7: 1000

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

tool: D127: 100.99999

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

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

assistant: ```c
extern double blackscholes_put(double S, double X, double T, double r, double v) {
  double vst, d1;
  assert(((S) == (S)) && (S) >= 0.000001 && (S) <= 1000000000000.0);
  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);
  vst = v * sqrt(T);
  d1 = (log(S / X) + (r + ((v) * (v)) / 2) * T) / (vst);
  return X * exp(-r * T) * cnd(-(d1 - vst)) - S * cnd(-d1);
}
```