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

```asm
<errf>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x258,%rsp
  mov    %rdx,%r13
  mov    %rcx,%r12
  mov    %fs:0x28,%rax
  mov    %rax,0x248(%rsp)
  xor    %eax,%eax
  mov    (%rdi),%r14
  mov    (%rsi),%r15
  movq   D169(%rip),%xmm0
  movq   %r14,%xmm5
  andpd  %xmm0,%xmm5
  movsd  %xmm5,(%rsp)
  movq   %r15,%xmm6
  andpd  %xmm6,%xmm0
  movapd %xmm0,%xmm6
  movsd  %xmm0,0x8(%rsp)
  movsd  D8(%rip),%xmm0
  comisd %xmm6,%xmm0
  jbe    L34
  movsd  D65(%rip),%xmm0
  comisd %xmm5,%xmm0
  ja     L94
L34:
  movq   $0x0,0x30(%rsp)
  movq   $0x0,0x140(%rsp)
  pxor   %xmm1,%xmm1
  movapd %xmm1,%xmm0
  mov    $0x9,%eax
  movsd  D19(%rip),%xmm3
L156:
  pxor   %xmm2,%xmm2
  cvtsi2sd %eax,%xmm2
  mulsd  %xmm2,%xmm0
  addsd  0x8(%rsp),%xmm0
  mulsd  %xmm1,%xmm2
  movsd  (%rsp),%xmm1
  subsd  %xmm2,%xmm1
  movapd %xmm0,%xmm2
  mulsd  %xmm0,%xmm2
  movapd %xmm1,%xmm4
  mulsd  %xmm1,%xmm4
  addsd  %xmm4,%xmm2
  mulsd  %xmm3,%xmm0
  divsd  %xmm2,%xmm0
  mulsd  %xmm3,%xmm1
  divsd  %xmm2,%xmm1
  sub    $0x1,%eax
  jne    L156
  movsd  D23(%rip),%xmm2
  mulsd  %xmm2,%xmm0
  movq   %xmm0,%rbx
  mulsd  %xmm2,%xmm1
  movq   %xmm1,%rbp
L72:
  pxor   %xmm0,%xmm0
  movq   %r15,%xmm5
  comisd %xmm5,%xmm0
  ja     L101
  pxor   %xmm0,%xmm0
  movq   %r14,%xmm5
  comisd %xmm5,%xmm0
  ja     L4
L121:
  mov    %rbx,0x0(%r13)
  mov    %rbp,(%r12)
  mov    0x248(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L14
  add    $0x258,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L94:
  movsd  D42(%rip),%xmm0
  movapd %xmm6,%xmm1
  divsd  D8(%rip),%xmm1
  movapd %xmm0,%xmm7
  subsd  %xmm1,%xmm7
  movsd  %xmm7,0x10(%rsp)
  movq   %r14,%xmm1
  divsd  D65(%rip),%xmm1
  mulsd  %xmm1,%xmm1
  subsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L135
  sqrtsd %xmm0,%xmm0
L129:
  mulsd  0x10(%rsp),%xmm0
  movq   %xmm0,%rbp
  mulsd  D133(%rip),%xmm0
  movsd  D42(%rip),%xmm1
  movapd %xmm1,%xmm6
  divsd  %xmm0,%xmm6
  movq   %rbp,%xmm0
  mulsd  D109(%rip),%xmm0
  cvttsd2si %xmm0,%ebx
  add    $0x7,%ebx
  movsd  %xmm6,0x10(%rsp)
  movapd %xmm6,%xmm0
  call   <log@plt>
  movapd %xmm0,%xmm1
  mov    $0x1,%eax
  sub    %ebx,%eax
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  mulsd  %xmm1,%xmm0
  call   <exp@plt>
  movapd %xmm0,%xmm3
  movsd  D19(%rip),%xmm5
  divsd  0x10(%rsp),%xmm5
  addsd  0x8(%rsp),%xmm5
  movq   %rbp,%xmm0
  mulsd  D94(%rip),%xmm0
  cvttsd2si %xmm0,%eax
  add    $0xa,%eax
  movslq %eax,%rdx
  movq   $0x0,0x30(%rsp,%rdx,8)
  movq   $0x0,0x140(%rsp,%rdx,8)
  test   %eax,%eax
  jle    L134
  mov    %rdx,%rax
  lea    0x30(%rsp),%rcx
  lea    0x140(%rsp),%rdx
  movsd  D19(%rip),%xmm4
L76:
  pxor   %xmm2,%xmm2
  cvtsi2sd %eax,%xmm2
  movapd %xmm2,%xmm0
  mulsd  (%rcx,%rax,8),%xmm0
  addsd  %xmm5,%xmm0
  mulsd  (%rdx,%rax,8),%xmm2
  movsd  (%rsp),%xmm1
  subsd  %xmm2,%xmm1
  movapd %xmm0,%xmm2
  mulsd  %xmm0,%xmm2
  mulsd  %xmm4,%xmm0
  divsd  %xmm2,%xmm0
  movsd  %xmm0,-0x8(%rcx,%rax,8)
  mulsd  %xmm4,%xmm1
  movapd %xmm1,%xmm0
  divsd  %xmm2,%xmm0
  movsd  %xmm0,-0x8(%rdx,%rax,8)
  sub    $0x1,%rax
  test   %eax,%eax
  jg     L76
L134:
  test   %ebx,%ebx
  jle    L103
  movslq %ebx,%rax
  pxor   %xmm0,%xmm0
  movapd %xmm0,%xmm2
L137:
  movapd %xmm2,%xmm4
  addsd  %xmm3,%xmm4
  movsd  0x28(%rsp,%rax,8),%xmm5
  movsd  0x138(%rsp,%rax,8),%xmm1
  movapd %xmm5,%xmm2
  mulsd  %xmm4,%xmm2
  movapd %xmm1,%xmm6
  mulsd  %xmm0,%xmm6
  subsd  %xmm6,%xmm2
  mulsd  %xmm5,%xmm0
  mulsd  %xmm4,%xmm1
  addsd  %xmm1,%xmm0
  mulsd  0x10(%rsp),%xmm3
  sub    $0x1,%rax
  test   %eax,%eax
  jg     L137
L80:
  movsd  D23(%rip),%xmm1
  mulsd  %xmm1,%xmm2
  movq   %xmm2,%rbx
  mulsd  %xmm1,%xmm0
  movq   %xmm0,%rbp
  jmp    L72
L135:
  call   <sqrt@plt>
  jmp    L129
L103:
  pxor   %xmm0,%xmm0
  movapd %xmm0,%xmm2
  jmp    L80
L101:
  movq   %r15,%xmm1
  mulsd  %xmm5,%xmm1
  movq   %r14,%xmm0
  mulsd  %xmm0,%xmm0
  subsd  %xmm0,%xmm1
  movsd  %xmm1,0x10(%rsp)
  movapd %xmm1,%xmm0
  call   <exp@plt>
  movsd  %xmm0,0x18(%rsp)
  lea    0x28(%rsp),%rdi
  lea    0x20(%rsp),%rsi
  movsd  (%rsp),%xmm7
  addsd  %xmm7,%xmm7
  movapd %xmm7,%xmm0
  mulsd  0x8(%rsp),%xmm0
  call   <sincos@plt>
  mov    0x28(%rsp),%r15
  movsd  0x18(%rsp),%xmm0
  addsd  %xmm0,%xmm0
  mulsd  0x20(%rsp),%xmm0
  movq   %rbx,%xmm6
  subsd  %xmm6,%xmm0
  movq   %xmm0,%rbx
  movsd  0x10(%rsp),%xmm0
  call   <exp@plt>
  mulsd  D5(%rip),%xmm0
  movq   %r15,%xmm7
  mulsd  %xmm7,%xmm0
  movq   %rbp,%xmm5
  subsd  %xmm5,%xmm0
  movq   %xmm0,%rbp
  movq   %r14,%xmm6
  comisd D13(%rip),%xmm6
  jbe    L121
  movapd %xmm0,%xmm5
  xorpd  D49(%rip),%xmm5
  movq   %xmm5,%rbp
  jmp    L121
L4:
  movq   %rbp,%xmm5
  xorpd  D49(%rip),%xmm5
  movq   %xmm5,%rbp
  jmp    L121
L14:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: data_label: D169
datatype: qword
name: read_data
arguments: data_label: D8
datatype: f64
name: read_data
arguments: data_label: D65
datatype: f64
name: read_data
arguments: data_label: D19
datatype: f64
name: read_data
arguments: data_label: D23
datatype: f64
name: read_data
arguments: data_label: D42
datatype: f64
name: read_data
arguments: data_label: D133
datatype: f64
name: read_data
arguments: data_label: D109
datatype: f64
name: read_data
arguments: data_label: D94
datatype: f32[2]
name: read_data
arguments: data_label: D5
datatype: f32[2]
name: read_data
arguments: data_label: D13
datatype: f32[4]
name: read_data
arguments: data_label: D49
datatype: f32[4]

tool: D169: 9223372036854775807

tool: D8: 4.29

tool: D65: 5.33

tool: D19: 0.5

tool: D23: 1.12837916709551

tool: D42: 1

tool: D133: 3.2

tool: D109: 23

tool: {"D94": [0.0,2.828125]}

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

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

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

assistant: ```c
void errf(double *var_xx, double *var_yy, double *var_wx, double *var_wy) {
  int n, nc, nu;
  double cc, h, q, rx[33], ry[33], saux, sx, sy, tn, tx, ty, wx, wy, x, xh, xl, xlim, xx, y, yh, ylim, yy;
  cc = 1.12837916709551;
  xlim = 5.33;
  ylim = 4.29;
  xx = *var_xx;
  yy = *var_yy;
  wx = *var_wx;
  wy = *var_wy;
  x = fabs(xx);
  y = fabs(yy);
  if (y < ylim && x < xlim) {
    q = (1.0 - y / ylim) * sqrt(1.0 - (x / xlim) * (x / xlim));
    h = 1.0 / (3.2 * q);
    nc = 7 + (int)(23.0 * q);
    xl = exp((1 - nc) * log(h));
    xh = y + 0.5 / h;
    yh = x;
    nu = 10 + (int)(21 * q);
    rx[nu] = 0.0;
    ry[nu] = 0.0;
    for (n = nu; n >= 1; n--) {
      tx = xh + (double)(n)*rx[n];
      ty = yh - (double)(n)*ry[n];
      tn = tx * tx;
      rx[n - 1] = (0.5 * tx) / tn;
      ry[n - 1] = (0.5 * ty) / tn;
    }
    sx = 0.0;
    sy = 0.0;
    for (n = nc; n >= 1; n--) {
      saux = sx + xl;
      sx = rx[n - 1] * saux - ry[n - 1] * sy;
      sy = rx[n - 1] * sy + ry[n - 1] * saux;
      xl = h * xl;
    }
    wx = cc * sx;
    wy = cc * sy;
  } else {
    xh = y;
    yh = x;
    rx[0] = 0.0;
    ry[0] = 0.0;
    for (n = 9; n >= 1; n--) {
      tx = xh + (double)(n)*rx[0];
      ty = yh - (double)(n)*ry[0];
      tn = tx * tx + ty * ty;
      rx[0] = (0.5 * tx) / tn;
      ry[0] = (0.5 * ty) / tn;
    }
    wx = cc * rx[0];
    wy = cc * ry[0];
  }
  if (yy < 0.0) {
    wx = (2.0 * exp(y * y - x * x)) * cos((2.0 * x) * y) - wx;
    wy = ((-1.0 * 2.0) * exp(y * y - x * x)) * sin((2.0 * x) * y) - wy;
    if (xx > 0.0)
      wy = -1.0 * wy;
  } else if (xx < 0.0)
    wy = -1.0 * wy;
  *var_xx = xx;
  *var_yy = yy;
  *var_wx = wx;
  *var_wy = wy;
}
```