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

```asm
<realtr>:
  endbr64
  push   %r14
  pxor   %xmm0,%xmm0
  mov    %ecx,%r14d
  push   %r13
  cvtsi2sd %rdx,%xmm0
  mov    %rdx,%r13
  push   %r12
  mov    %rsi,%r12
  sar    %r13
  push   %rbp
  mov    %rdi,%rbp
  push   %rbx
  mov    %rdx,%rbx
  sub    $0x10,%rsp
  movsd  D38(%rip),%xmm1
  divsd  %xmm0,%xmm1
  movapd %xmm1,%xmm0
  movsd  %xmm1,0x8(%rsp)
  call   <sin@plt>
  movsd  0x8(%rsp),%xmm1
  movapd %xmm0,%xmm5
  addsd  %xmm0,%xmm5
  addsd  %xmm1,%xmm1
  mulsd  %xmm0,%xmm5
  movapd %xmm1,%xmm0
  movsd  %xmm5,(%rsp)
  call   <sin@plt>
  test   %r14d,%r14d
  movsd  (%rsp),%xmm5
  movapd %xmm0,%xmm4
  js     L0
  xorpd  D34(%rip),%xmm4
  movsd  D14(%rip),%xmm0
L5:
  test   %r13,%r13
  js     L38
  movsd  D43(%rip),%xmm6
  lea    0x1(%r13),%rcx
  pxor   %xmm3,%xmm3
  xor    %eax,%eax
  lea    0x0(,%rbx,8),%rdx
  nop
L11:
  movsd  0x0(%rbp,%rax,8),%xmm7
  movsd  (%r12,%rax,8),%xmm2
  movsd  (%r12,%rdx,1),%xmm10
  movsd  0x0(%rbp,%rdx,1),%xmm1
  movapd %xmm7,%xmm8
  movapd %xmm2,%xmm9
  addsd  %xmm1,%xmm8
  addsd  %xmm10,%xmm9
  subsd  %xmm7,%xmm1
  subsd  %xmm10,%xmm2
  movapd %xmm9,%xmm7
  mulsd  %xmm3,%xmm9
  movapd %xmm1,%xmm10
  mulsd  %xmm0,%xmm1
  mulsd  %xmm0,%xmm7
  mulsd  %xmm3,%xmm10
  subsd  %xmm9,%xmm1
  addsd  %xmm10,%xmm7
  movapd %xmm1,%xmm9
  subsd  %xmm2,%xmm9
  addsd  %xmm1,%xmm2
  movapd %xmm8,%xmm1
  subsd  %xmm7,%xmm1
  addsd  %xmm8,%xmm7
  movsd  %xmm9,(%r12,%rdx,1)
  movsd  %xmm2,(%r12,%rax,8)
  movapd %xmm4,%xmm2
  movsd  %xmm1,0x0(%rbp,%rdx,1)
  movapd %xmm5,%xmm1
  sub    $0x8,%rdx
  mulsd  %xmm3,%xmm2
  movsd  %xmm7,0x0(%rbp,%rax,8)
  add    $0x1,%rax
  mulsd  %xmm0,%xmm1
  addsd  %xmm2,%xmm1
  movapd %xmm0,%xmm2
  mulsd  %xmm4,%xmm0
  subsd  %xmm1,%xmm2
  movapd %xmm5,%xmm1
  mulsd  %xmm3,%xmm1
  subsd  %xmm1,%xmm0
  addsd  %xmm3,%xmm0
  movapd %xmm2,%xmm3
  mulsd  %xmm2,%xmm3
  movapd %xmm0,%xmm1
  mulsd  %xmm0,%xmm1
  addsd  %xmm1,%xmm3
  movapd %xmm6,%xmm1
  divsd  %xmm3,%xmm1
  movapd %xmm0,%xmm3
  movapd %xmm2,%xmm0
  addsd  %xmm6,%xmm1
  mulsd  %xmm1,%xmm3
  mulsd  %xmm1,%xmm0
  cmp    %rcx,%rax
  jne    L11
L38:
  add    $0x10,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
  nopl   (%rax)
L0:
  movsd  0x0(%rbp),%xmm0
  movsd  %xmm0,0x0(%rbp,%rbx,8)
  movsd  (%r12),%xmm0
  movsd  %xmm0,(%r12,%rbx,8)
  movsd  D9(%rip),%xmm0
  jmp    L5
```

tool_calls: name: parse_data
arguments: label: D38
type_guessed: f64
name: parse_data
arguments: label: D34
type_guessed: byte[16]
name: parse_data
arguments: label: D14
type_guessed: byte[8]
name: parse_data
arguments: label: D43
type_guessed: f64
name: parse_data
arguments: label: D9
type_guessed: byte[8]

tool: D38: 1.5707963267948966

tool: {"D34": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00]}

tool: {"D14": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0xbf]}

tool: D43: 0.5

tool: {"D9": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0x3f]}

assistant: ```c
void realtr(double *a, double *b, long half_length, int isn) {
  long nh, j, k;
  double sd, cd, sn, cn, a1r, a1i, a2r, a2i, re, im;
  nh = half_length >> 1;
  sd = 1.57079632679489661923 / half_length;
  cd = 2.0 * sin(sd) * sin(sd);
  sd = sin(sd + sd);
  sn = 0;
  if (isn < 0) {
    cn = 1;
    a[half_length] = a[0];
    b[half_length] = b[0];
  } else {
    cn = -1;
    sd = -sd;
  }
  for (j = 0; j <= nh; j++) {
    k = half_length - j;
    a1r = a[j] + a[k];
    a2r = b[j] + b[k];
    a1i = b[j] - b[k];
    a2i = -a[j] + a[k];
    re = cn * a2r + sn * a2i;
    im = -sn * a2r + cn * a2i;
    b[k] = +im - a1i;
    b[j] = im + a1i;
    a[k] = a1r - re;
    a[j] = a1r + re;
    a1r = cn - (cd * cn + sd * sn);
    sn = (sd * cn - cd * sn) + sn;
    cn = 0.5 / (a1r * a1r + sn * sn) + 0.5;
    sn *= cn;
    cn *= a1r;
  }
}
```