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

```asm
<QuantileNormal>:
  endbr64
  sub    $0x18,%rsp
  movapd %xmm0,%xmm3
  movsd  %xmm0,0x8(%rsp)
  movsd  D148(%rip),%xmm0
  comisd %xmm3,%xmm0
  ja     L23
  movsd  D57(%rip),%xmm1
  subsd  %xmm3,%xmm1
  movsd  D152(%rip),%xmm0
  movsd  D40(%rip),%xmm2
  comisd %xmm1,%xmm2
  jbe    L47
L19:
  movsd  D148(%rip),%xmm1
  comisd 0x8(%rsp),%xmm1
  jbe    L59
  jmp    L56
L23:
  movsd  D152(%rip),%xmm0
  movsd  D40(%rip),%xmm1
  comisd 0x8(%rsp),%xmm1
  jbe    L3
L56:
  xorpd  D127(%rip),%xmm0
L59:
  add    $0x18,%rsp
  ret
L3:
  movsd  0x8(%rsp),%xmm1
L47:
  mulsd  %xmm1,%xmm1
  movsd  D57(%rip),%xmm0
  divsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%edi
  call   <log@plt>
  mov    %eax,%edi
  call   <sqrt@plt>
  movapd %xmm0,%xmm2
  mulsd  D141(%rip),%xmm0
  subsd  D2(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  subsd  D116(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  subsd  D57(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  subsd  D115(%rip),%xmm0
  movapd %xmm2,%xmm1
  mulsd  D146(%rip),%xmm1
  addsd  D145(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  addsd  D30(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  addsd  D171(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  addsd  D102(%rip),%xmm1
  divsd  %xmm1,%xmm0
  addsd  %xmm2,%xmm0
  jmp    L19
```

tool_calls: name: parse_data
arguments: data_label: D148
type_guessed: f64
name: parse_data
arguments: data_label: D57
type_guessed: f64
name: parse_data
arguments: data_label: D152
type_guessed: f32[2]
name: parse_data
arguments: data_label: D40
type_guessed: byte[8]
name: parse_data
arguments: data_label: D127
type_guessed: f32[4]
name: parse_data
arguments: data_label: D141
type_guessed: byte[8]
name: parse_data
arguments: data_label: D2
type_guessed: byte[8]
name: parse_data
arguments: data_label: D116
type_guessed: byte[8]
name: parse_data
arguments: data_label: D115
type_guessed: byte[8]
name: parse_data
arguments: data_label: D146
type_guessed: byte[8]
name: parse_data
arguments: data_label: D145
type_guessed: byte[8]
name: parse_data
arguments: data_label: D30
type_guessed: byte[8]
name: parse_data
arguments: data_label: D171
type_guessed: byte[8]
name: parse_data
arguments: data_label: D102
type_guessed: byte[8]

tool: D148: 0.5

tool: D57: 1

tool: {"D152": [0.0,4.4755859375]}

tool: {"D40": [0x23,0x42,0x92,0x0c,0xa1,0x9c,0xc7,0x3b]}

tool: {"D127": [0.0,-0.0,0.0,0.0]}

tool: {"D141": [0x62,0x1a,0xed,0xc3,0xae,0xc8,0x07,0xbf]}

tool: {"D2": [0x30,0xa8,0x87,0x73,0xcc,0xe9,0x94,0x3f]}

tool: {"D116": [0x46,0xf1,0x67,0x5c,0x4b,0xe7,0xd5,0x3f]}

tool: {"D115": [0x37,0xf2,0x4e,0xc6,0x74,0x9f,0xd4,0x3f]}

tool: {"D146": [0x88,0x45,0x07,0xda,0xc3,0x96,0x6f,0x3f]}

tool: {"D145": [0x68,0x5f,0x64,0x3e,0x73,0x81,0xba,0x3f]}

tool: {"D30": [0x54,0xf0,0x34,0xb0,0xcc,0xfe,0xe0,0x3f]}

tool: {"D171": [0x2f,0xbd,0x89,0x04,0xa9,0xd5,0xe2,0x3f]}

tool: {"D102": [0x54,0xf7,0xcc,0x9d,0xe6,0x6e,0xb9,0x3f]}

assistant: ```c
double QuantileNormal(double prob) {
  double a0 = -.322232431088, a1 = -1, a2 = -.342242088547, a3 = -.0204231210245;
  double a4 = -.453642210148e-4, b0 = .0993484626060, b1 = .588581570495;
  double b2 = .531103462366, b3 = .103537752850, b4 = .0038560700634;
  double y, z = 0, p = prob, p1;
  p1 = (p < 0.5 ? p : 1 - p);
  if (p1 < 1e-20)
    z = 999;
  else {
    y = sqrt(log(1 / (p1 * p1)));
    z = y + ((((y * a4 + a3) * y + a2) * y + a1) * y + a0) / ((((y * b4 + b3) * y + b2) * y + b1) * y + b0);
  }
  return (p < 0.5 ? -z : z);
}
```