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

```asm
<reals_>:
  endbr64
  push   %r15
  mov    %ecx,%r8d
  mov    %edx,%eax
  push   %r14
  push   %r13
  mov    %rsi,%r13
  mov    %edx,%esi
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x48,%rsp
  neg    %ecx
  cmovs  %r8d,%ecx
  neg    %eax
  cmovs  %edx,%eax
  mov    %r8d,%edx
  imul   %eax,%edx
  test   %edx,%edx
  je     L44
  pxor   %xmm0,%xmm0
  mov    %ecx,%ebx
  mov    %r8d,0x10(%rsp)
  mov    %rdi,%r12
  cvtsi2sd %eax,%xmm0
  movsd  D93(%rip),%xmm1
  movsd  D10(%rip),%xmm4
  mov    %ecx,0xc(%rsp)
  imul   %eax,%ebx
  divsd  %xmm0,%xmm1
  add    $0x2,%ebx
  mov    %ebx,%r9d
  sar    %r9d
  mov    %r9d,0x8(%rsp)
  movapd %xmm1,%xmm0
  movsd  %xmm1,0x28(%rsp)
  mulsd  %xmm4,%xmm0
  mulsd  D79(%rip),%xmm0
  call   <sin@plt>
  movsd  D79(%rip),%xmm5
  mulsd  0x28(%rsp),%xmm5
  movapd %xmm0,%xmm6
  addsd  %xmm0,%xmm6
  mulsd  %xmm0,%xmm6
  movapd %xmm5,%xmm0
  movsd  %xmm6,(%rsp)
  call   <sin@plt>
  mov    0x10(%rsp),%r8d
  movsd  (%rsp),%xmm6
  mov    D10(%rip),%rax
  mov    0x8(%rsp),%r9d
  movapd %xmm0,%xmm5
  test   %r8d,%r8d
  mov    0xc(%rsp),%ecx
  movq   %rax,%xmm4
  js     L79
  xorpd  D105(%rip),%xmm5
  movsd  D70(%rip),%xmm7
  movslq %ebx,%rbx
L56:
  movslq %ecx,%rdx
  shl    $0x3,%rbx
  xor    %ebp,%ebp
  pxor   %xmm3,%xmm3
  shl    $0x3,%rdx
  xor    %r14d,%r14d
  mov    $0x20,%r15d
  mov    $0x1,%eax
  jmp    L43
  nopl   0x0(%rax,%rax,1)
L65:
  movapd %xmm6,%xmm0
  movapd %xmm5,%xmm1
  movapd %xmm6,%xmm2
  mulsd  %xmm3,%xmm1
  mulsd  %xmm7,%xmm0
  mulsd  %xmm3,%xmm2
  addsd  %xmm1,%xmm0
  movapd %xmm5,%xmm1
  mulsd  %xmm7,%xmm1
  subsd  %xmm0,%xmm7
  subsd  %xmm2,%xmm1
  addsd  %xmm1,%xmm3
L38:
  add    %ecx,%eax
  add    %rdx,%rbp
  sub    %rdx,%rbx
  cmp    %eax,%r9d
  jl     L23
L43:
  movsd  (%r12,%rbp,1),%xmm2
  movsd  -0x10(%r12,%rbx,1),%xmm0
  add    $0x1,%r14d
  movsd  -0x10(%r13,%rbx,1),%xmm9
  movapd %xmm2,%xmm8
  subsd  %xmm0,%xmm2
  addsd  %xmm0,%xmm8
  movsd  0x0(%r13,%rbp,1),%xmm0
  movapd %xmm0,%xmm1
  movapd %xmm2,%xmm10
  subsd  %xmm9,%xmm0
  addsd  %xmm9,%xmm1
  mulsd  %xmm7,%xmm2
  mulsd  %xmm3,%xmm10
  movapd %xmm1,%xmm9
  mulsd  %xmm3,%xmm1
  mulsd  %xmm7,%xmm9
  subsd  %xmm2,%xmm1
  addsd  %xmm10,%xmm9
  movapd %xmm1,%xmm2
  subsd  %xmm0,%xmm2
  addsd  %xmm1,%xmm0
  mulsd  %xmm4,%xmm2
  mulsd  %xmm4,%xmm0
  movsd  %xmm2,-0x10(%r13,%rbx,1)
  movsd  %xmm0,0x0(%r13,%rbp,1)
  movapd %xmm8,%xmm0
  subsd  %xmm9,%xmm0
  mulsd  %xmm4,%xmm0
  movsd  %xmm0,-0x10(%r12,%rbx,1)
  movapd %xmm8,%xmm0
  addsd  %xmm9,%xmm0
  mulsd  %xmm4,%xmm0
  movsd  %xmm0,(%r12,%rbp,1)
  cmp    %r15d,%r14d
  jne    L65
  pxor   %xmm0,%xmm0
  lea    0x38(%rsp),%rdi
  lea    0x30(%rsp),%rsi
  mov    %r8d,0x24(%rsp)
  cvtsi2sd %r14d,%xmm0
  mulsd  0x28(%rsp),%xmm0
  mov    %rdx,0x18(%rsp)
  lea    0x20(%r14),%r15d
  mulsd  D79(%rip),%xmm0
  mov    %eax,0x20(%rsp)
  mov    %r9d,0xc(%rsp)
  mov    %ecx,0x8(%rsp)
  movsd  %xmm6,0x10(%rsp)
  movsd  %xmm5,(%rsp)
  call   <sincos@plt>
  mov    0x24(%rsp),%r8d
  movsd  (%rsp),%xmm5
  mov    D10(%rip),%rdi
  movsd  0x38(%rsp),%xmm3
  test   %r8d,%r8d
  movsd  0x30(%rsp),%xmm7
  mov    0x8(%rsp),%ecx
  mov    0xc(%rsp),%r9d
  movsd  0x10(%rsp),%xmm6
  movq   %rdi,%xmm4
  mov    0x20(%rsp),%eax
  mov    0x18(%rsp),%rdx
  jle    L38
  add    %ecx,%eax
  xorpd  D105(%rip),%xmm7
  add    %rdx,%rbp
  sub    %rdx,%rbx
  cmp    %eax,%r9d
  jge    L43
L23:
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L79:
  movsd  (%r12),%xmm0
  movslq %ebx,%rbx
  movsd  D113(%rip),%xmm7
  lea    0x0(,%rbx,8),%rax
  movsd  %xmm0,-0x10(%r12,%rax,1)
  movsd  0x0(%r13),%xmm0
  movsd  %xmm0,-0x10(%r13,%rax,1)
  jmp    L56
  nopw   0x0(%rax,%rax,1)
L44:
  add    $0x48,%rsp
  mov    %r8d,%edx
  lea    D66(%rip),%rdi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1080 <printf@plt>
```

tool_calls: name: unpack_data
arguments: data_tag: D93
data_type: f32[2]
name: unpack_data
arguments: data_tag: D10
data_type: f64
name: unpack_data
arguments: data_tag: D79
data_type: f64
name: unpack_data
arguments: data_tag: D105
data_type: f32[4]
name: unpack_data
arguments: data_tag: D70
data_type: f32[2]
name: unpack_data
arguments: data_tag: D113
data_type: f64
name: unpack_data
arguments: data_tag: D66
data_type: string

tool: {"D93": [0.0,-2.25]}

tool: D10: 0.5

tool: D79: 0.7853981633974483

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

tool: {"D70": [0.0,-1.875]}

tool: D113: 1

tool: D66: 
error - zero in reals parameters : %d : %d 

assistant: ```c
void reals_(double *a, double *b, int n, int isn) {
  int inc, j, k, lim, mm, ml, nf, nk, nh;
  double aa, ab, ba, bb, cd, cn, dr, em, rad, re, sd, sn;
  double xx;
  a--;
  b--;
  inc = abs(isn);
  nf = abs(n);
  if (nf * isn == 0) {
    printf("\nerror - zero in reals parameters : %d : %d ", n, isn);
    return;
  }
  nk = (nf * inc) + 2;
  nh = nk / 2;
  rad = 0.785398163397448278900;
  dr = -4.0 / (double)(nf);
  xx = sin((double)0.5 * dr * rad);
  cd = 2.0 * xx * xx;
  sd = sin(dr * rad);
  lim = 32;
  mm = lim;
  ml = 0;
  sn = 0.0;
  if (isn < 0) {
    cn = 1.0;
    a[nk - 1] = a[1];
    b[nk - 1] = b[1];
  } else {
    cn = -1.0;
    sd = -sd;
  }
  for (j = 1; j <= nh; j += inc) {
    k = nk - j;
    aa = a[j] + a[k];
    ab = a[j] - a[k];
    ba = b[j] + b[k];
    bb = b[j] - b[k];
    re = (cn * ba) + (sn * ab);
    em = (sn * ba) - (cn * ab);
    b[k] = (double)((em - bb) * 0.5);
    b[j] = (double)((em + bb) * 0.5);
    a[k] = (double)((aa - re) * 0.5);
    a[j] = (double)((aa + re) * 0.5);
    ml++;
    if (ml != mm) {
      aa = cn - ((cd * cn) + (sd * sn));
      sn = ((sd * cn) - (cd * sn)) + sn;
      cn = aa;
    } else {
      mm += lim;
      sn = ((double)ml) * dr * rad;
      cn = cos(sn);
      if (isn > 0)
        cn = -cn;
      sn = sin(sn);
    }
  }
  return;
}
```