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

```asm
<zrati_>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xa8,%rsp
  mov    %rdi,%rbx
  mov    %rsi,%rbp
  mov    %rdx,%r12
  mov    %rdx,0x40(%rsp)
  mov    %rcx,0x38(%rsp)
  mov    %r8,%r13
  mov    %r9,%r14
  mov    0xe0(%rsp),%r15
  mov    %fs:0x28,%rax
  mov    %rax,0x98(%rsp)
  xor    %eax,%eax
  call   <zabs_@plt>
  mov    %eax,%edx
  pxor   %xmm2,%xmm2
  cvtsi2sd %eax,%xmm2
  pxor   %xmm0,%xmm0
  cvtsd2ss (%r12),%xmm0
  cvttss2si %xmm0,%eax
  mov    0x38(%rsp),%rcx
  add    (%rcx),%eax
  sub    $0x1,%eax
  pxor   %xmm0,%xmm0
  cvtsi2ss %edx,%xmm0
  cvttss2si %xmm0,%edx
  lea    0x1(%rdx),%ecx
  pxor   %xmm1,%xmm1
  cvtsi2ss %ecx,%xmm1
  pxor   %xmm0,%xmm0
  cvtsi2ss %eax,%xmm0
  comiss %xmm0,%xmm1
  jb     L153
  pxor   %xmm0,%xmm0
  cvtss2sd %xmm1,%xmm0
L156:
  movsd  D16(%rip),%xmm4
  movapd %xmm4,%xmm3
  divsd  %xmm2,%xmm3
  movsd  %xmm3,0x90(%rsp)
  movsd  (%rbx),%xmm1
  addsd  %xmm1,%xmm1
  mulsd  %xmm3,%xmm1
  movapd %xmm1,%xmm7
  mulsd  %xmm3,%xmm7
  movsd  %xmm7,0x18(%rsp)
  movsd  0x0(%rbp),%xmm1
  addsd  %xmm1,%xmm1
  movq   D119(%rip),%xmm2
  movapd %xmm3,%xmm5
  xorpd  %xmm2,%xmm5
  mulsd  %xmm5,%xmm1
  movapd %xmm1,%xmm6
  mulsd  %xmm3,%xmm6
  movsd  %xmm6,0x20(%rsp)
  movapd %xmm7,%xmm1
  mulsd  %xmm0,%xmm1
  mulsd  %xmm6,%xmm0
  movapd %xmm1,%xmm3
  xorpd  %xmm2,%xmm3
  movsd  %xmm3,0x78(%rsp)
  xorpd  %xmm0,%xmm2
  movsd  %xmm2,0x60(%rsp)
  movsd  %xmm4,0x70(%rsp)
  movq   $0x0,0x58(%rsp)
  addsd  %xmm7,%xmm1
  movsd  %xmm1,0x80(%rsp)
  addsd  %xmm6,%xmm0
  movsd  %xmm0,0x68(%rsp)
  sub    %edx,%eax
  lea    -0x1(%rax),%ecx
  test   %ecx,%ecx
  mov    $0x0,%eax
  cmovle %ecx,%eax
  mov    %eax,0x4c(%rsp)
  lea    0x60(%rsp),%rsi
  lea    0x78(%rsp),%rdi
  mov    $0x0,%eax
  call   <zabs_@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,0x8(%rsp)
  lea    0x58(%rsp),%rsi
  lea    0x70(%rsp),%rdi
  mov    $0x0,%eax
  call   <zabs_@plt>
  pxor   %xmm6,%xmm6
  cvtsi2sd %eax,%xmm6
  movsd  %xmm6,0x10(%rsp)
  movsd  0x8(%rsp),%xmm7
  movapd %xmm7,%xmm0
  addsd  %xmm7,%xmm0
  mulsd  (%r15),%xmm6
  divsd  %xmm6,%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L152
  sqrtsd %xmm0,%xmm0
  movsd  %xmm0,0x30(%rsp)
L35:
  movsd  D16(%rip),%xmm0
  divsd  0x10(%rsp),%xmm0
  movapd %xmm0,%xmm1
  mulsd  0x70(%rsp),%xmm1
  movsd  %xmm1,0x70(%rsp)
  movapd %xmm0,%xmm1
  mulsd  0x58(%rsp),%xmm1
  movsd  %xmm1,0x58(%rsp)
  movapd %xmm0,%xmm1
  mulsd  0x78(%rsp),%xmm1
  movsd  %xmm1,0x78(%rsp)
  movapd %xmm0,%xmm1
  mulsd  0x60(%rsp),%xmm1
  movsd  %xmm1,0x60(%rsp)
  mulsd  0x8(%rsp),%xmm0
  movsd  %xmm0,0x8(%rsp)
  movsd  0x30(%rsp),%xmm3
  movsd  %xmm3,0x28(%rsp)
  mov    $0x1,%r12d
  mov    $0x1,%ebx
  jmp    L122
L153:
  cvtss2sd %xmm0,%xmm0
  jmp    L156
L152:
  call   <sqrt@plt>
  movsd  %xmm0,0x30(%rsp)
  jmp    L35
L91:
  add    $0x2,%ebx
  sub    0x4c(%rsp),%ebx
  pxor   %xmm5,%xmm5
  cvtsi2ss %ebx,%xmm5
  cvtss2sd %xmm5,%xmm5
  movsd  %xmm5,0x80(%rsp)
  movq   $0x0,0x68(%rsp)
  mov    0x38(%rsp),%rax
  mov    (%rax),%edx
  lea    -0x1(%rdx),%eax
  pxor   %xmm7,%xmm7
  cvtsi2ss %eax,%xmm7
  cvtss2sd %xmm7,%xmm7
  mov    0x40(%rsp),%rax
  addsd  (%rax),%xmm7
  movsd  D16(%rip),%xmm3
  divsd  %xmm6,%xmm3
  movsd  %xmm3,0x70(%rsp)
  movq   $0x0,0x58(%rsp)
  movq   $0x0,0x78(%rsp)
  movq   $0x0,0x60(%rsp)
  test   %ebx,%ebx
  jle    L87
  add    $0x1,%ebx
  pxor   %xmm6,%xmm6
  movapd %xmm6,%xmm9
  movapd %xmm6,%xmm4
  mov    $0x1,%eax
  movsd  D16(%rip),%xmm8
  jmp    L77
L129:
  call   <sqrt@plt>
  jmp    L198
L150:
  call   <sqrt@plt>
L130:
  mulsd  0x30(%rsp),%xmm0
  movsd  %xmm0,0x28(%rsp)
  mov    $0x2,%r12d
L117:
  movsd  0x10(%rsp),%xmm7
  movsd  %xmm7,0x8(%rsp)
  mov    %ebp,%ebx
L122:
  lea    0x1(%rbx),%ebp
  movsd  0x78(%rsp),%xmm3
  movsd  %xmm3,0x90(%rsp)
  movsd  0x60(%rsp),%xmm2
  movsd  %xmm2,0x88(%rsp)
  movsd  0x80(%rsp),%xmm1
  movsd  0x68(%rsp),%xmm0
  movapd %xmm3,%xmm5
  mulsd  %xmm1,%xmm5
  movapd %xmm2,%xmm4
  mulsd  %xmm0,%xmm4
  subsd  %xmm4,%xmm5
  movsd  0x70(%rsp),%xmm4
  subsd  %xmm5,%xmm4
  movsd  %xmm4,0x78(%rsp)
  movapd %xmm2,%xmm5
  mulsd  %xmm1,%xmm5
  movapd %xmm3,%xmm4
  mulsd  %xmm0,%xmm4
  addsd  %xmm4,%xmm5
  movsd  0x58(%rsp),%xmm4
  subsd  %xmm5,%xmm4
  movsd  %xmm4,0x60(%rsp)
  movsd  %xmm3,0x70(%rsp)
  movsd  %xmm2,0x58(%rsp)
  addsd  0x18(%rsp),%xmm1
  movsd  %xmm1,0x80(%rsp)
  addsd  0x20(%rsp),%xmm0
  movsd  %xmm0,0x68(%rsp)
  lea    0x60(%rsp),%rsi
  lea    0x78(%rsp),%rdi
  mov    $0x0,%eax
  call   <zabs_@plt>
  pxor   %xmm6,%xmm6
  cvtsi2sd %eax,%xmm6
  movsd  %xmm6,0x10(%rsp)
  movsd  0x28(%rsp),%xmm7
  comisd 0x8(%rsp),%xmm7
  jae    L117
  cmp    $0x2,%r12d
  je     L91
  lea    0x68(%rsp),%rsi
  lea    0x80(%rsp),%rdi
  mov    $0x0,%eax
  call   <zabs_@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  mulsd  D44(%rip),%xmm0
  movsd  %xmm0,0x28(%rsp)
  mulsd  %xmm0,%xmm0
  subsd  D16(%rip),%xmm0
  pxor   %xmm3,%xmm3
  ucomisd %xmm0,%xmm3
  ja     L129
  sqrtsd %xmm0,%xmm0
L198:
  movsd  0x28(%rsp),%xmm1
  addsd  %xmm0,%xmm1
  movsd  0x10(%rsp),%xmm0
  divsd  0x8(%rsp),%xmm0
  movapd %xmm0,%xmm2
  cmplesd %xmm1,%xmm2
  andpd  %xmm2,%xmm0
  andnpd %xmm1,%xmm2
  orpd   %xmm2,%xmm0
  movapd %xmm0,%xmm1
  mulsd  %xmm0,%xmm1
  subsd  D16(%rip),%xmm1
  divsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L150
  sqrtsd %xmm0,%xmm0
  jmp    L130
L56:
  movapd %xmm0,%xmm4
  movapd %xmm1,%xmm3
L77:
  movapd %xmm5,%xmm0
  addsd  %xmm7,%xmm0
  movsd  0x18(%rsp),%xmm2
  mulsd  %xmm0,%xmm2
  mulsd  0x20(%rsp),%xmm0
  movapd %xmm3,%xmm1
  mulsd  %xmm2,%xmm1
  movapd %xmm4,%xmm10
  mulsd  %xmm0,%xmm10
  subsd  %xmm10,%xmm1
  addsd  %xmm9,%xmm1
  mulsd  %xmm3,%xmm0
  mulsd  %xmm4,%xmm2
  addsd  %xmm2,%xmm0
  addsd  %xmm6,%xmm0
  subsd  %xmm8,%xmm5
  add    $0x1,%eax
  movapd %xmm3,%xmm9
  movapd %xmm4,%xmm6
  cmp    %ebx,%eax
  jne    L56
  movsd  %xmm3,0x90(%rsp)
  movsd  %xmm4,0x88(%rsp)
  movsd  %xmm1,0x70(%rsp)
  movsd  %xmm0,0x58(%rsp)
  movsd  %xmm3,0x78(%rsp)
  movsd  %xmm4,0x60(%rsp)
  movsd  %xmm5,0x80(%rsp)
  ucomisd D57(%rip),%xmm1
  jp     L199
  jne    L199
  ucomisd D57(%rip),%xmm0
  jp     L199
  jne    L199
L69:
  movsd  (%r15),%xmm0
  movsd  %xmm0,0x70(%rsp)
  movsd  %xmm0,0x58(%rsp)
L199:
  movslq %edx,%rdx
  lea    -0x8(,%rdx,8),%rax
  lea    0x58(%rsp),%rcx
  lea    0x70(%rsp),%rdx
  lea    0x60(%rsp),%rsi
  lea    0x78(%rsp),%rdi
  lea    (%r14,%rax,1),%r9
  lea    0x0(%r13,%rax,1),%r8
  mov    $0x0,%eax
  call   <zdiv_@plt>
  mov    0x38(%rsp),%rax
  mov    (%rax),%ebp
  cmp    $0x1,%ebp
  je     L14
  lea    -0x1(%rbp),%eax
  pxor   %xmm0,%xmm0
  cvtsi2ss %eax,%xmm0
  cvtss2sd %xmm0,%xmm0
  movsd  %xmm0,0x80(%rsp)
  movq   $0x0,0x68(%rsp)
  mov    0x40(%rsp),%rax
  movsd  (%rax),%xmm0
  movsd  0x18(%rsp),%xmm7
  mulsd  %xmm0,%xmm7
  movsd  %xmm7,0x8(%rsp)
  mulsd  0x20(%rsp),%xmm0
  movsd  %xmm0,0x10(%rsp)
  jle    L14
  lea    -0x2(%rbp),%ebx
  movslq %ebx,%rbx
  lea    0x88(%rsp),%r12
  jmp    L20
L268:
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
L231:
  movsd  D16(%rip),%xmm1
  divsd  %xmm0,%xmm1
  movapd %xmm1,%xmm0
  mulsd  0x90(%rsp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,0x0(%r13,%rbx,8)
  movapd %xmm1,%xmm0
  xorpd  D119(%rip),%xmm0
  mulsd  0x88(%rsp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,(%r14,%rbx,8)
  movsd  0x80(%rsp),%xmm0
  subsd  D16(%rip),%xmm0
  movsd  %xmm0,0x80(%rsp)
  sub    $0x1,%rbx
  mov    %ebp,%eax
  sub    %ebx,%eax
  cmp    %ebp,%eax
  jg     L14
L20:
  movsd  0x80(%rsp),%xmm0
  movsd  0x68(%rsp),%xmm2
  movsd  0x18(%rsp),%xmm7
  movapd %xmm7,%xmm1
  mulsd  %xmm0,%xmm1
  movsd  0x20(%rsp),%xmm6
  movapd %xmm6,%xmm3
  mulsd  %xmm2,%xmm3
  subsd  %xmm3,%xmm1
  addsd  0x8(%rsp),%xmm1
  addsd  0x8(%r13,%rbx,8),%xmm1
  movsd  %xmm1,0x90(%rsp)
  mulsd  %xmm6,%xmm0
  mulsd  %xmm7,%xmm2
  addsd  %xmm2,%xmm0
  addsd  0x10(%rsp),%xmm0
  addsd  0x8(%r14,%rbx,8),%xmm0
  movsd  %xmm0,0x88(%rsp)
  lea    0x90(%rsp),%rdi
  mov    %r12,%rsi
  mov    $0x0,%eax
  call   <zabs_@plt>
  test   %eax,%eax
  jne    L268
  movsd  (%r15),%xmm0
  movsd  %xmm0,0x90(%rsp)
  movsd  %xmm0,0x88(%rsp)
  mulsd  D100(%rip),%xmm0
  jmp    L231
L14:
  mov    0x98(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L110
  mov    $0x0,%eax
  add    $0xa8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L87:
  ucomisd D57(%rip),%xmm3
  jp     L199
  je     L69
  jmp    L199
L110:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: retrieve_data
arguments: data_tag: D16
data_type: f64
name: retrieve_data
arguments: data_tag: D119
data_type: f32[4]
name: retrieve_data
arguments: data_tag: D44
data_type: f32[2]
name: retrieve_data
arguments: data_tag: D57
data_type: f32[2]
name: retrieve_data
arguments: data_tag: D100
data_type: f64

tool: D16: 1

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

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

tool: {"D57": [0.0,0.0]}

tool: D100: 1.4142135623730951

assistant: ```c
int zrati_(double *zr, double *zi, double *fnu, int *n, double *cyr, double *cyi, double *tol) {
  static double czeror = 0.;
  static double czeroi = 0.;
  static double coner = 1.;
  static double conei = 0.;
  static double rt2 = 1.41421356237309505;
  int i__1;
  double d__1;
  int i__, k;
  double ak;
  int id, kk;
  double az, ap1, ap2, p1i, p2i, t1i, p1r, p2r, t1r, arg, rak, rho;
  int inu;
  double pti, tti, rzi, ptr, ttr, rzr, rap1, flam, dfnu, fdnu;
  int magz;
  int idnu;
  double fnup;
  double test, test1, amagz;
  int itime;
  double cdfnui, cdfnur;
  --cyi;
  --cyr;
  az = zabs_(zr, zi);
  inu = (int)((float)(*fnu));
  idnu = inu + *n - 1;
  magz = (int)((float)az);
  amagz = (double)((float)(magz + 1));
  fdnu = (double)((float)idnu);
  fnup = ((amagz) >= (fdnu) ? (amagz) : (fdnu));
  id = idnu - magz - 1;
  itime = 1;
  k = 1;
  ptr = 1. / az;
  rzr = ptr * (*zr + *zr) * ptr;
  rzi = -ptr * (*zi + *zi) * ptr;
  t1r = rzr * fnup;
  t1i = rzi * fnup;
  p2r = -t1r;
  p2i = -t1i;
  p1r = coner;
  p1i = conei;
  t1r += rzr;
  t1i += rzi;
  if (id > 0) {
    id = 0;
  }
  ap2 = zabs_(&p2r, &p2i);
  ap1 = zabs_(&p1r, &p1i);
  arg = (ap2 + ap2) / (ap1 * *tol);
  test1 = sqrt(arg);
  test = test1;
  rap1 = 1. / ap1;
  p1r *= rap1;
  p1i *= rap1;
  p2r *= rap1;
  p2i *= rap1;
  ap2 *= rap1;
L10:
  ++k;
  ap1 = ap2;
  ptr = p2r;
  pti = p2i;
  p2r = p1r - (t1r * ptr - t1i * pti);
  p2i = p1i - (t1r * pti + t1i * ptr);
  p1r = ptr;
  p1i = pti;
  t1r += rzr;
  t1i += rzi;
  ap2 = zabs_(&p2r, &p2i);
  if (ap1 <= test) {
    goto L10;
  }
  if (itime == 2) {
    goto L20;
  }
  ak = zabs_(&t1r, &t1i) * .5;
  flam = ak + sqrt(ak * ak - 1.);
  d__1 = ap2 / ap1;
  rho = ((d__1) <= (flam) ? (d__1) : (flam));
  test = test1 * sqrt(rho / (rho * rho - 1.));
  itime = 2;
  goto L10;
L20:
  kk = k + 1 - id;
  ak = (double)((float)kk);
  t1r = ak;
  t1i = czeroi;
  dfnu = *fnu + (double)((float)(*n - 1));
  p1r = 1. / ap2;
  p1i = czeroi;
  p2r = czeror;
  p2i = czeroi;
  i__1 = kk;
  for (i__ = 1; i__ <= i__1; ++i__) {
    ptr = p1r;
    pti = p1i;
    rap1 = dfnu + t1r;
    ttr = rzr * rap1;
    tti = rzi * rap1;
    p1r = ptr * ttr - pti * tti + p2r;
    p1i = ptr * tti + pti * ttr + p2i;
    p2r = ptr;
    p2i = pti;
    t1r -= coner;
  }
  if (p1r != czeror || p1i != czeroi) {
    goto L40;
  }
  p1r = *tol;
  p1i = *tol;
L40:
  zdiv_(&p2r, &p2i, &p1r, &p1i, &cyr[*n], &cyi[*n]);
  if (*n == 1) {
    return 0;
  }
  k = *n - 1;
  ak = (double)((float)k);
  t1r = ak;
  t1i = czeroi;
  cdfnur = *fnu * rzr;
  cdfnui = *fnu * rzi;
  i__1 = *n;
  for (i__ = 2; i__ <= i__1; ++i__) {
    ptr = cdfnur + (t1r * rzr - t1i * rzi) + cyr[k + 1];
    pti = cdfnui + (t1r * rzi + t1i * rzr) + cyi[k + 1];
    ak = zabs_(&ptr, &pti);
    if (ak != czeror) {
      goto L50;
    }
    ptr = *tol;
    pti = *tol;
    ak = *tol * rt2;
  L50:
    rak = coner / ak;
    cyr[k] = rak * ptr * rak;
    cyi[k] = -rak * pti * rak;
    t1r -= coner;
    --k;
  }
  return 0;
}
```