user: 请显示与下面列出的汇编代码相匹配的C语言实现:

```asm
<_bal_pja_invnorm>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x40,%rsp
  mov    %rdi,-0x38(%rbp)
  mov    %rsi,-0x40(%rbp)
  movsd  D84(%rip),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  D204(%rip),%xmm0
  subsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movq   $0x0,-0x28(%rbp)
  jmp    L80
L48:
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  comisd -0x20(%rbp),%xmm0
  setae  %dl
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,8),%rcx
  mov    -0x38(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rax),%xmm1
  movsd  -0x18(%rbp),%xmm0
  comisd %xmm1,%xmm0
  setae  %al
  and    %edx,%eax
  test   %al,%al
  je     L65
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  D187(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x10(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  D226(%rip),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x8(%rbp),%xmm1
  movsd  D374(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x8(%rbp),%xmm1
  movsd  D131(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x8(%rbp),%xmm1
  movsd  D328(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x8(%rbp),%xmm1
  movsd  D89(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x8(%rbp),%xmm1
  movsd  D223(%rip),%xmm0
  addsd  %xmm1,%xmm0
  mulsd  -0x10(%rbp),%xmm0
  movsd  D266(%rip),%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x8(%rbp),%xmm2
  movsd  D3(%rip),%xmm1
  addsd  %xmm2,%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x8(%rbp),%xmm2
  movsd  D234(%rip),%xmm1
  addsd  %xmm2,%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x8(%rbp),%xmm2
  movsd  D303(%rip),%xmm1
  addsd  %xmm2,%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x8(%rbp),%xmm2
  movsd  D376(%rip),%xmm1
  addsd  %xmm2,%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x8(%rbp),%xmm2
  movsd  D204(%rip),%xmm1
  addsd  %xmm2,%xmm1
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  jmp    L81
L65:
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm1,%xmm0
  seta   %dl
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,8),%rcx
  mov    -0x38(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rax),%xmm1
  movsd  -0x20(%rbp),%xmm0
  comisd %xmm1,%xmm0
  seta   %al
  and    %edx,%eax
  test   %al,%al
  je     L36
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,%edi
  call   <log@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  movsd  D394(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,%edi
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x10(%rbp)
  movsd  D96(%rip),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x10(%rbp),%xmm1
  movsd  D270(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x10(%rbp),%xmm1
  movsd  D174(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x10(%rbp),%xmm1
  movsd  D163(%rip),%xmm0
  addsd  %xmm1,%xmm0
  mulsd  -0x10(%rbp),%xmm0
  movapd %xmm0,%xmm1
  movsd  D16(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x10(%rbp),%xmm1
  movsd  D219(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  D103(%rip),%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x10(%rbp),%xmm2
  movsd  D4(%rip),%xmm1
  addsd  %xmm2,%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x10(%rbp),%xmm2
  movsd  D154(%rip),%xmm1
  addsd  %xmm2,%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x10(%rbp),%xmm2
  movsd  D151(%rip),%xmm1
  addsd  %xmm2,%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x10(%rbp),%xmm2
  movsd  D204(%rip),%xmm1
  addsd  %xmm2,%xmm1
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  movq   %rax,%xmm0
  call   <isnan@plt>
  test   %rax,%rax
  je     L81
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  movsd  D109(%rip),%xmm0
  movsd  %xmm0,(%rax)
  jmp    L81
L36:
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  comisd -0x18(%rbp),%xmm0
  seta   %dl
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,8),%rcx
  mov    -0x38(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rax),%xmm1
  movsd  D204(%rip),%xmm0
  comisd %xmm1,%xmm0
  seta   %al
  and    %edx,%eax
  test   %al,%al
  je     L52
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  D204(%rip),%xmm0
  subsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,%edi
  call   <log@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  movsd  D394(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,%edi
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x10(%rbp)
  movsd  D96(%rip),%xmm0
  mulsd  -0x10(%rbp),%xmm0
  movapd %xmm0,%xmm1
  movsd  D270(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x10(%rbp),%xmm1
  movsd  D174(%rip),%xmm0
  addsd  %xmm1,%xmm0
  mulsd  -0x10(%rbp),%xmm0
  movapd %xmm0,%xmm1
  movsd  D163(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x10(%rbp),%xmm1
  movsd  D16(%rip),%xmm0
  addsd  %xmm1,%xmm0
  mulsd  -0x10(%rbp),%xmm0
  movapd %xmm0,%xmm1
  movsd  D219(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movq   D143(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  D103(%rip),%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x10(%rbp),%xmm2
  movsd  D4(%rip),%xmm1
  addsd  %xmm2,%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x10(%rbp),%xmm2
  movsd  D154(%rip),%xmm1
  addsd  %xmm2,%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x10(%rbp),%xmm2
  movsd  D151(%rip),%xmm1
  addsd  %xmm2,%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x10(%rbp),%xmm2
  movsd  D204(%rip),%xmm1
  addsd  %xmm2,%xmm1
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  movq   %rax,%xmm0
  call   <isnan@plt>
  test   %rax,%rax
  je     L81
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  movsd  D77(%rip),%xmm0
  movsd  %xmm0,(%rax)
  jmp    L81
L52:
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
L81:
  addq   $0x1,-0x28(%rbp)
L80:
  mov    -0x28(%rbp),%rax
  cmp    -0x40(%rbp),%rax
  jb     L48
  nop
  nop
  leave
  ret
```

tool_calls: name: read_data
arguments: label: D84
data_type: f64
name: read_data
arguments: label: D204
data_type: f64
name: read_data
arguments: label: D187
data_type: f32[2]
name: read_data
arguments: label: D226
data_type: byte[8]
name: read_data
arguments: label: D374
data_type: byte[8]
name: read_data
arguments: label: D131
data_type: byte[8]
name: read_data
arguments: label: D328
data_type: byte[8]
name: read_data
arguments: label: D89
data_type: byte[8]
name: read_data
arguments: label: D223
data_type: byte[8]
name: read_data
arguments: label: D266
data_type: byte[8]
name: read_data
arguments: label: D3
data_type: byte[8]
name: read_data
arguments: label: D234
data_type: byte[8]
name: read_data
arguments: label: D303
data_type: byte[8]
name: read_data
arguments: label: D376
data_type: byte[8]
name: read_data
arguments: label: D394
data_type: f32[2]
name: read_data
arguments: label: D96
data_type: byte[8]
name: read_data
arguments: label: D270
data_type: byte[8]
name: read_data
arguments: label: D174
data_type: byte[8]
name: read_data
arguments: label: D163
data_type: byte[8]
name: read_data
arguments: label: D16
data_type: byte[8]
name: read_data
arguments: label: D219
data_type: byte[8]
name: read_data
arguments: label: D103
data_type: byte[8]
name: read_data
arguments: label: D4
data_type: byte[8]
name: read_data
arguments: label: D154
data_type: byte[8]
name: read_data
arguments: label: D151
data_type: byte[8]
name: read_data
arguments: label: D109
data_type: byte[8]
name: read_data
arguments: label: D143
data_type: f32[2]
name: read_data
arguments: label: D77
data_type: f64

tool: D84: 0.02425

tool: D204: 1

tool: {"D187": [0.0,1.75]}

tool: {"D226": [0x25,0x05,0x1e,0xbc,0x31,0xd9,0x43,0xc0]}

tool: {"D374": [0x9b,0x03,0x34,0x70,0x46,0x9e,0x6b,0x40]}

tool: {"D131": [0x99,0x3b,0xc5,0x2d,0xdb,0x3e,0x71,0xc0]}

tool: {"D328": [0x1b,0x40,0x0b,0xb4,0x72,0x4b,0x61,0x40]}

tool: {"D89": [0xcd,0x8b,0xc0,0x34,0x30,0xaa,0x3e,0xc0]}

tool: {"D223": [0x79,0x54,0x57,0x20,0x93,0x0d,0x04,0x40]}

tool: {"D266": [0xc4,0x04,0x30,0xce,0xf0,0x3c,0x4b,0xc0]}

tool: {"D3": [0x77,0x42,0xf0,0x2c,0xbf,0x32,0x64,0x40]}

tool: {"D234": [0xd2,0xd8,0x02,0x0b,0x5e,0x76,0x63,0xc0]}

tool: {"D303": [0xbe,0xe9,0xa7,0xb1,0x48,0xb3,0x50,0x40]}

tool: {"D376": [0x26,0x78,0x14,0x7e,0xb5,0x8f,0x2a,0xc0]}

tool: {"D394": [0.0,-2.0]}

tool: {"D96": [0xe0,0xcf,0x4a,0x92,0x0d,0xe3,0x7f,0xbf]}

tool: {"D270": [0xb8,0x81,0xe8,0xc0,0x24,0xa2,0xd4,0xbf]}

tool: {"D174": [0x58,0x17,0x70,0xc1,0xc0,0x34,0x03,0xc0]}

tool: {"D163": [0x1a,0x3a,0x70,0x2c,0xda,0x65,0x04,0xc0]}

tool: {"D16": [0xc7,0x4d,0xea,0xf4,0xa7,0x7f,0x11,0x40]}

tool: {"D219": [0xa2,0xcf,0x3f,0x1e,0x5c,0x81,0x07,0x40]}

tool: {"D103": [0xd4,0x9f,0xac,0x57,0xd8,0xe2,0x7f,0x3f]}

tool: {"D4": [0x64,0x03,0x59,0x2b,0x4d,0xa3,0xd4,0x3f]}

tool: {"D154": [0x16,0xe6,0x8a,0x7c,0xa2,0x8f,0x03,0x40]}

tool: {"D151": [0x19,0xb1,0x95,0x68,0x07,0x09,0x0e,0x40]}

tool: {"D109": [0x21,0x93,0x8c,0x9c,0x85,0x39,0x43,0xc0]}

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

tool: D77: 8.20954

assistant: ```c
void _bal_pja_invnorm(double *p, size_t num) {
  const double a[] = {-3.969683028665376e+01, 2.209460984245205e+02, -2.759285104469687e+02, 1.383577518672690e+02, -3.066479806614716e+01, 2.506628277459239e+00};
  const double b[] = {-5.447609879822406e+01, 1.615858368580409e+02, -1.556989798598866e+02, 6.680131188771972e+01, -1.328068155288572e+01};
  const double c[] = {-7.784894002430293e-03, -3.223964580411365e-01, -2.400758277161838e+00, -2.549732539343734e+00, 4.374664141464968e+00, 2.938163982698783e+00};
  const double d[] = {7.784695709041462e-03, 3.224671290700398e-01, 2.445134137142996e+00, 3.754408661907416e+00};
  const double p_low = 0.02425, p_high = 1. - p_low;
  double q, r;
  size_t i;
  for (i = 0; i < num; i++) {
    if (p_low <= p[i] & p[i] <= p_high) {
      q = p[i] - .5;
      r = q * q;
      p[i] = ((((((a[0] * r + a[1]) * r + a[2]) * r + a[3]) * r + a[4]) * r + a[5]) * q) / (((((b[0] * r + b[1]) * r + b[2]) * r + b[3]) * r + b[4]) * r + 1.);
    } else if (0. < p[i] & p[i] < p_low) {
      q = sqrt(-2. * log(p[i]));
      p[i] = (((((c[0] * q + c[1]) * q + c[2]) * q + c[3]) * q + c[4]) * q + c[5]) / ((((d[0] * q + d[1]) * q + d[2]) * q + d[3]) * q + 1.);
      if (isnan(p[i])) {
        p[i] = -38.44939;
      }
    } else if (p_high < p[i] & p[i] < 1.) {
      q = sqrt(-2. * log(1. - p[i]));
      p[i] = -(((((c[0] * q + c[1]) * q + c[2]) * q + c[3]) * q + c[4]) * q + c[5]) / ((((d[0] * q + d[1]) * q + d[2]) * q + d[3]) * q + 1.);
      if (isnan(p[i])) {
        p[i] = 8.20954;
      }
    } else {
      p[i] = 0.;
    }
  }
}
```