user: Please show me the C implementation that matches the assembly code listed below:

```asm
<PointNormal>:
  endbr64
  movsd  D61(%rip),%xmm3
  movapd %xmm0,%xmm1
  comisd %xmm0,%xmm3
  ja     L28
  movsd  D120(%rip),%xmm5
  movsd  D57(%rip),%xmm4
  movapd %xmm5,%xmm2
  subsd  %xmm0,%xmm2
  movsd  D84(%rip),%xmm0
  comisd %xmm2,%xmm4
  jbe    L59
  comisd %xmm1,%xmm3
  ja     L54
L0:
  ret
  nopw   0x0(%rax,%rax,1)
L28:
  movsd  D57(%rip),%xmm2
  movsd  D159(%rip),%xmm0
  comisd %xmm1,%xmm2
  ja     L0
  movsd  D120(%rip),%xmm5
  movapd %xmm1,%xmm2
L59:
  mulsd  %xmm2,%xmm2
  movapd %xmm5,%xmm0
  sub    $0x18,%rsp
  movsd  %xmm1,0x8(%rsp)
  movsd  %xmm5,(%rsp)
  divsd  %xmm2,%xmm0
  cvttsd2si %xmm0,%edi
  call   <log@plt>
  mov    %eax,%edi
  call   <sqrt@plt>
  movsd  (%rsp),%xmm5
  movsd  D187(%rip),%xmm2
  movapd %xmm0,%xmm4
  movsd  D87(%rip),%xmm0
  mov    D61(%rip),%rax
  mulsd  %xmm4,%xmm2
  movsd  0x8(%rsp),%xmm1
  mulsd  %xmm4,%xmm0
  addsd  D74(%rip),%xmm2
  movq   %rax,%xmm3
  subsd  D49(%rip),%xmm0
  comisd %xmm1,%xmm3
  mulsd  %xmm4,%xmm2
  mulsd  %xmm4,%xmm0
  addsd  D6(%rip),%xmm2
  subsd  D21(%rip),%xmm0
  mulsd  %xmm4,%xmm2
  mulsd  %xmm4,%xmm0
  addsd  D11(%rip),%xmm2
  mulsd  %xmm4,%xmm2
  addsd  D4(%rip),%xmm2
  subsd  %xmm5,%xmm0
  mulsd  %xmm4,%xmm0
  subsd  D168(%rip),%xmm0
  divsd  %xmm2,%xmm0
  addsd  %xmm4,%xmm0
  ja     L47
L6:
  add    $0x18,%rsp
  ret
  nopw   0x0(%rax,%rax,1)
L54:
  xorpd  D188(%rip),%xmm0
  ret
  nopl   0x0(%rax)
L47:
  xorpd  D188(%rip),%xmm0
  jmp    L6
```

tool_calls: name: read_data
arguments: label: D61
data_type: f64
name: read_data
arguments: label: D120
data_type: f64
name: read_data
arguments: label: D57
data_type: byte[8]
name: read_data
arguments: label: D84
data_type: f32[2]
name: read_data
arguments: label: D159
data_type: f32[2]
name: read_data
arguments: label: D187
data_type: byte[8]
name: read_data
arguments: label: D87
data_type: byte[8]
name: read_data
arguments: label: D74
data_type: byte[8]
name: read_data
arguments: label: D49
data_type: byte[8]
name: read_data
arguments: label: D6
data_type: byte[8]
name: read_data
arguments: label: D21
data_type: byte[8]
name: read_data
arguments: label: D11
data_type: byte[8]
name: read_data
arguments: label: D4
data_type: byte[8]
name: read_data
arguments: label: D168
data_type: byte[8]
name: read_data
arguments: label: D188
data_type: f32[4]

tool: D61: 0.5

tool: D120: 1

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

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

tool: {"D159": [0.0,-4.4755859375]}

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

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

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

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

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

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

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

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

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

tool: {"D188": [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)
    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);
}
```