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

```asm
<PointNormal>:
  endbr64
  sub    $0x18,%rsp
  movsd  D50(%rip),%xmm4
  movapd %xmm0,%xmm2
  movapd %xmm0,%xmm1
  comisd %xmm0,%xmm4
  ja     L5
  movsd  D137(%rip),%xmm1
  subsd  %xmm0,%xmm1
L5:
  movsd  D65(%rip),%xmm3
  movsd  D29(%rip),%xmm0
  movsd  %xmm2,0x8(%rsp)
  comisd %xmm1,%xmm3
  jbe    L1
L22:
  add    $0x18,%rsp
  ret
  nopl   0x0(%rax,%rax,1)
L1:
  mulsd  %xmm1,%xmm1
  movsd  D137(%rip),%xmm0
  divsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%edi
  call   <log@plt>
  mov    %eax,%edi
  call   <sqrt@plt>
  movsd  D12(%rip),%xmm1
  movsd  D50(%rip),%xmm5
  movapd %xmm0,%xmm3
  movsd  D18(%rip),%xmm0
  movsd  0x8(%rsp),%xmm2
  mulsd  %xmm3,%xmm1
  mulsd  %xmm3,%xmm0
  comisd %xmm2,%xmm5
  addsd  D144(%rip),%xmm1
  subsd  D165(%rip),%xmm0
  mulsd  %xmm3,%xmm1
  mulsd  %xmm3,%xmm0
  addsd  D6(%rip),%xmm1
  subsd  D117(%rip),%xmm0
  mulsd  %xmm3,%xmm1
  mulsd  %xmm3,%xmm0
  addsd  D59(%rip),%xmm1
  subsd  D137(%rip),%xmm0
  mulsd  %xmm3,%xmm1
  mulsd  %xmm3,%xmm0
  addsd  D163(%rip),%xmm1
  subsd  D143(%rip),%xmm0
  divsd  %xmm1,%xmm0
  addsd  %xmm3,%xmm0
  jbe    L22
  xorpd  D131(%rip),%xmm0
  add    $0x18,%rsp
  ret
```

tool_calls: name: read_data
arguments: data_tag: D50
type_guessed: f64
name: read_data
arguments: data_tag: D137
type_guessed: f64
name: read_data
arguments: data_tag: D65
type_guessed: byte[8]
name: read_data
arguments: data_tag: D29
type_guessed: f32[2]
name: read_data
arguments: data_tag: D12
type_guessed: byte[8]
name: read_data
arguments: data_tag: D18
type_guessed: byte[8]
name: read_data
arguments: data_tag: D144
type_guessed: byte[8]
name: read_data
arguments: data_tag: D165
type_guessed: byte[8]
name: read_data
arguments: data_tag: D6
type_guessed: byte[8]
name: read_data
arguments: data_tag: D117
type_guessed: byte[8]
name: read_data
arguments: data_tag: D59
type_guessed: byte[8]
name: read_data
arguments: data_tag: D163
type_guessed: byte[8]
name: read_data
arguments: data_tag: D143
type_guessed: byte[8]
name: read_data
arguments: data_tag: D131
type_guessed: f32[4]

tool: D50: 0.5

tool: D137: 1

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

tool: {"D29": [0.0,-6.11029052734375]}

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

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

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

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

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

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

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

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

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

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

assistant: ```c
double PointNormal(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)
    return (-9999);
  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);
}
```