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

```asm
<gbs_put>:
  endbr64
  sub    $0x48,%rsp
  movsd  %xmm0,0x38(%rsp)
  movsd  %xmm1,0x10(%rsp)
  movsd  %xmm2,0x8(%rsp)
  movsd  %xmm3,0x18(%rsp)
  movsd  %xmm4,0x20(%rsp)
  movsd  %xmm5,0x28(%rsp)
  ucomisd %xmm0,%xmm0
  mov    $0x0,%edi
  jp     L78
  movapd %xmm0,%xmm6
  mov    $0x0,%edi
  comisd D52(%rip),%xmm0
  jb     L78
  movsd  D81(%rip),%xmm0
  comisd %xmm6,%xmm0
  setae  %dil
  movzbl %dil,%edi
L78:
  mov    $0x0,%eax
  call   <assert@plt>
  movsd  0x10(%rsp),%xmm4
  ucomisd %xmm4,%xmm4
  mov    $0x0,%edi
  jp     L75
  mov    $0x0,%edi
  comisd D52(%rip),%xmm4
  jb     L75
  movsd  D81(%rip),%xmm0
  comisd %xmm4,%xmm0
  setae  %dil
  movzbl %dil,%edi
L75:
  mov    $0x0,%eax
  call   <assert@plt>
  movsd  0x8(%rsp),%xmm1
  ucomisd %xmm1,%xmm1
  mov    $0x0,%edi
  jp     L123
  mov    $0x0,%edi
  comisd D12(%rip),%xmm1
  jb     L123
  movsd  D96(%rip),%xmm0
  comisd %xmm1,%xmm0
  setae  %dil
  movzbl %dil,%edi
L123:
  mov    $0x0,%eax
  call   <assert@plt>
  movsd  0x18(%rsp),%xmm2
  ucomisd %xmm2,%xmm2
  mov    $0x0,%edi
  jp     L17
  mov    $0x0,%edi
  comisd D76(%rip),%xmm2
  jb     L17
  movsd  D102(%rip),%xmm0
  comisd %xmm2,%xmm0
  setae  %dil
  movzbl %dil,%edi
L17:
  mov    $0x0,%eax
  call   <assert@plt>
  movsd  0x20(%rsp),%xmm3
  ucomisd %xmm3,%xmm3
  mov    $0x0,%edi
  jp     L27
  mov    $0x0,%edi
  comisd D120(%rip),%xmm3
  jb     L27
  movsd  D102(%rip),%xmm0
  comisd %xmm3,%xmm0
  setae  %dil
  movzbl %dil,%edi
L27:
  mov    $0x0,%eax
  call   <assert@plt>
  movsd  0x28(%rsp),%xmm7
  ucomisd %xmm7,%xmm7
  mov    $0x0,%edi
  jp     L33
  mov    $0x0,%edi
  comisd D31(%rip),%xmm7
  jb     L33
  movsd  D136(%rip),%xmm0
  comisd %xmm7,%xmm0
  setae  %dil
  movzbl %dil,%edi
L33:
  mov    $0x0,%eax
  call   <assert@plt>
  pxor   %xmm0,%xmm0
  movsd  0x8(%rsp),%xmm5
  ucomisd %xmm5,%xmm0
  ja     L67
  sqrtsd %xmm5,%xmm5
  movapd %xmm5,%xmm0
L120:
  movapd %xmm0,%xmm7
  mulsd  0x28(%rsp),%xmm7
  movsd  %xmm7,0x30(%rsp)
  movsd  0x38(%rsp),%xmm2
  divsd  0x10(%rsp),%xmm2
  movapd %xmm2,%xmm0
  call   <log@plt>
  movapd %xmm0,%xmm1
  movsd  0x28(%rsp),%xmm0
  mulsd  %xmm0,%xmm0
  mulsd  D62(%rip),%xmm0
  addsd  0x20(%rsp),%xmm0
  movsd  0x8(%rsp),%xmm3
  mulsd  %xmm3,%xmm0
  addsd  %xmm1,%xmm0
  divsd  0x30(%rsp),%xmm0
  movsd  %xmm0,0x28(%rsp)
  movsd  0x18(%rsp),%xmm4
  xorpd  D71(%rip),%xmm4
  movapd %xmm4,%xmm0
  mulsd  %xmm3,%xmm0
  call   <exp@plt>
  movapd %xmm0,%xmm5
  mulsd  0x10(%rsp),%xmm5
  movsd  %xmm5,0x10(%rsp)
  movsd  0x28(%rsp),%xmm0
  subsd  0x30(%rsp),%xmm0
  xorpd  D71(%rip),%xmm0
  mov    $0x1,%eax
  call   <cnd@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movapd %xmm0,%xmm5
  mulsd  0x10(%rsp),%xmm5
  movsd  %xmm5,0x10(%rsp)
  movsd  0x20(%rsp),%xmm6
  subsd  0x18(%rsp),%xmm6
  movapd %xmm6,%xmm0
  mulsd  0x8(%rsp),%xmm0
  call   <exp@plt>
  movapd %xmm0,%xmm7
  mulsd  0x38(%rsp),%xmm7
  movsd  %xmm7,0x8(%rsp)
  movsd  0x28(%rsp),%xmm0
  xorpd  D71(%rip),%xmm0
  mov    $0x1,%eax
  call   <cnd@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  mulsd  0x8(%rsp),%xmm1
  movsd  0x10(%rsp),%xmm0
  subsd  %xmm1,%xmm0
  add    $0x48,%rsp
  ret
L67:
  movsd  0x8(%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L120
```

tool_calls: name: parse_data
arguments: label: D52
type_guessed: f64[2]
name: parse_data
arguments: label: D81
type_guessed: f64
name: parse_data
arguments: label: D12
type_guessed: byte[16]
name: parse_data
arguments: label: D96
type_guessed: f64
name: parse_data
arguments: label: D76
type_guessed: f64
name: parse_data
arguments: label: D102
type_guessed: f64
name: parse_data
arguments: label: D120
type_guessed: byte[16]
name: parse_data
arguments: label: D31
type_guessed: f64[2]
name: parse_data
arguments: label: D136
type_guessed: f64
name: parse_data
arguments: label: D62
type_guessed: byte[8]
name: parse_data
arguments: label: D71
type_guessed: byte[16]

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

tool: D81: 1000000000000

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

tool: D96: 20

tool: D76: 0

tool: D102: 1000

tool: {"D120": [0x00,0x00,0x00,0x00,0x00,0x40,0x8f,0xc0,0x48,0xaf,0xbc,0x9a,0xf2,0xd7,0x7a,0x3e]}

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

tool: D136: 100.99999

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

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

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