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

```asm
<_bal_pja_invnorm>:
  endbr64
  test   %rsi,%rsi
  je     L114
  push   %r12
  mov    %rsi,%r12
  push   %rbp
  xor    %ebp,%ebp
  push   %rbx
  mov    %rdi,%rbx
  jmp    L88
  xchg   %ax,%ax
L126:
  movapd %xmm0,%xmm1
  subsd  D99(%rip),%xmm1
  movsd  D265(%rip),%xmm0
  movapd %xmm1,%xmm2
  mulsd  %xmm1,%xmm2
  mulsd  %xmm2,%xmm0
  addsd  D30(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  subsd  D157(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D221(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  subsd  D403(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D45(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D291(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  addsd  D85(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  subsd  D294(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  addsd  D131(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  subsd  D132(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  addsd  D322(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rbx)
L102:
  add    $0x1,%rbp
  add    $0x8,%rbx
  cmp    %rbp,%r12
  je     L62
L88:
  movsd  (%rbx),%xmm0
  comisd D92(%rip),%xmm0
  jb     L3
  movsd  D385(%rip),%xmm1
  comisd %xmm0,%xmm1
  jae    L126
L3:
  pxor   %xmm3,%xmm3
  comisd %xmm3,%xmm0
  jbe    L134
  movsd  D92(%rip),%xmm1
  comisd %xmm0,%xmm1
  jbe    L134
  cvttsd2si %xmm0,%edi
  call   <log@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  mulsd  D348(%rip),%xmm0
  cvttsd2si %xmm0,%edi
  call   <sqrt@plt>
  movsd  D189(%rip),%xmm1
  movapd %xmm0,%xmm2
  movsd  D63(%rip),%xmm0
  mulsd  %xmm2,%xmm1
  mulsd  %xmm2,%xmm0
  addsd  D161(%rip),%xmm1
  subsd  D141(%rip),%xmm0
  mulsd  %xmm2,%xmm1
  mulsd  %xmm2,%xmm0
  addsd  D154(%rip),%xmm1
  subsd  D289(%rip),%xmm0
  mulsd  %xmm2,%xmm1
  mulsd  %xmm2,%xmm0
  addsd  D327(%rip),%xmm1
  subsd  D444(%rip),%xmm0
  mulsd  %xmm2,%xmm1
  mulsd  %xmm2,%xmm0
  addsd  D322(%rip),%xmm1
  addsd  D224(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D252(%rip),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rbx)
  call   <isnan@plt>
  test   %rax,%rax
  je     L102
  mov    D328(%rip),%rax
  add    $0x1,%rbp
  add    $0x8,%rbx
  mov    %rax,-0x8(%rbx)
  cmp    %rbp,%r12
  jne    L88
L62:
  pop    %rbx
  pop    %rbp
  pop    %r12
  ret
  nopw   0x0(%rax,%rax,1)
L134:
  comisd D385(%rip),%xmm0
  jbe    L29
  movsd  D322(%rip),%xmm4
  comisd %xmm0,%xmm4
  jbe    L29
  movapd %xmm4,%xmm1
  subsd  %xmm0,%xmm1
  cvttsd2si %xmm1,%edi
  call   <log@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  mulsd  D348(%rip),%xmm0
  cvttsd2si %xmm0,%edi
  call   <sqrt@plt>
  movsd  D189(%rip),%xmm1
  movapd %xmm0,%xmm2
  movsd  D63(%rip),%xmm0
  mulsd  %xmm2,%xmm1
  mulsd  %xmm2,%xmm0
  addsd  D161(%rip),%xmm1
  subsd  D141(%rip),%xmm0
  mulsd  %xmm2,%xmm1
  mulsd  %xmm2,%xmm0
  addsd  D154(%rip),%xmm1
  subsd  D289(%rip),%xmm0
  mulsd  %xmm2,%xmm1
  mulsd  %xmm2,%xmm0
  addsd  D327(%rip),%xmm1
  subsd  D444(%rip),%xmm0
  mulsd  %xmm2,%xmm1
  mulsd  %xmm2,%xmm0
  addsd  D322(%rip),%xmm1
  addsd  D224(%rip),%xmm0
  mulsd  %xmm2,%xmm0
  addsd  D252(%rip),%xmm0
  xorpd  D281(%rip),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rbx)
  call   <isnan@plt>
  test   %rax,%rax
  je     L102
  mov    D186(%rip),%rax
  mov    %rax,(%rbx)
  jmp    L102
  nop
L29:
  movq   $0x0,(%rbx)
  jmp    L102
  nopl   0x0(%rax)
L114:
  ret
```

tool_calls: name: retrieve_data
arguments: data_label: D99
type_guessed: f32[2]
name: retrieve_data
arguments: data_label: D265
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D30
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D157
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D221
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D403
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D45
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D291
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D85
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D294
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D131
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D132
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D322
type_guessed: f64
name: retrieve_data
arguments: data_label: D92
type_guessed: f64[2]
name: retrieve_data
arguments: data_label: D385
type_guessed: byte[16]
name: retrieve_data
arguments: data_label: D348
type_guessed: f32[2]
name: retrieve_data
arguments: data_label: D189
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D63
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D161
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D141
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D154
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D289
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D327
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D444
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D224
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D252
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D328
type_guessed: byte[8]
name: retrieve_data
arguments: data_label: D281
type_guessed: f32[4]
name: retrieve_data
arguments: data_label: D186
type_guessed: f64

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

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

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

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

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

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

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

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

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

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

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

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

tool: D322: 1

tool: {"D92": [0.02425,0.97575]}

tool: {"D385": [0xd3,0x4d,0x62,0x10,0x58,0x39,0xef,0x3f,0x00,0x00,0x00,0x00,0x00,0x00,0xe0,0x3f]}

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

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

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

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

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

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

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

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

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

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

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

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

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

tool: D186: 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.;
    }
  }
}
```