user: Converte this assembly snippet into its equivalent C code:

```asm
<cffti1>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  mov    %rsi,%r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rdx,%rbx
  sub    $0x88,%rsp
  mov    %rdi,0x40(%rsp)
  test   %rdi,%rdi
  js     L86
  pxor   %xmm0,%xmm0
  cvtsi2sd %rdi,%xmm0
L92:
  movsd  D20(%rip),%xmm1
  mov    0x40(%rsp),%rdi
  mov    %rbx,%rcx
  mov    $0x5,%edx
  lea    D25(%rip),%rsi
  divsd  %xmm0,%xmm1
  movsd  %xmm1,0x50(%rsp)
  call   <factorize@plt>
  mov    0x8(%rbx),%rax
  test   %rax,%rax
  je     L37
  lea    0x10(%rbx),%rdx
  lea    0x78(%rsp),%r14
  movq   $0x1,0x38(%rsp)
  lea    (%rdx,%rax,8),%rax
  mov    %rdx,0x60(%rsp)
  mov    %rax,0x68(%rsp)
  movq   $0x0,0x8(%rsp)
  nopl   0x0(%rax,%rax,1)
L63:
  mov    0x60(%rsp),%rax
  xor    %edx,%edx
  mov    (%rax),%rdi
  mov    0x38(%rsp),%rax
  mov    %rax,0x48(%rsp)
  mov    %rax,%rcx
  imul   %rdi,%rax
  mov    %rdi,0x20(%rsp)
  mov    %rax,%rsi
  mov    %rax,0x38(%rsp)
  mov    0x40(%rsp),%rax
  div    %rsi
  mov    %rax,%rbp
  cmp    $0x1,%rdi
  jbe    L75
  lea    0x1(%rbp),%rdx
  lea    -0x2(%rax,%rax,1),%rax
  mov    %rcx,0x18(%rsp)
  movq   $0x1,0x10(%rsp)
  cmp    $0x1,%rdx
  mov    $0x0,%edx
  cmovbe %rdx,%rax
  mov    %rax,0x58(%rsp)
  xchg   %ax,%ax
L16:
  mov    0x18(%rsp),%rax
  test   %rax,%rax
  js     L68
  pxor   %xmm0,%xmm0
  cvtsi2sd %rax,%xmm0
L24:
  mulsd  0x50(%rsp),%xmm0
  mov    0x8(%rsp),%rax
  mov    D29(%rip),%rsi
  shl    $0x3,%rax
  lea    0x8(%r13,%rax,1),%rdx
  lea    0x0(%r13,%rax,1),%rcx
  mov    %rsi,(%rcx)
  mov    %rdx,0x30(%rsp)
  movq   $0x0,(%rdx)
  mov    0x40(%rsp),%rdx
  mov    %rcx,0x28(%rsp)
  movsd  %xmm0,(%rsp)
  cmp    %rdx,0x38(%rsp)
  ja     L74
  lea    0x10(%r13,%rax,1),%r15
  mov    $0x1,%r12d
  lea    0x70(%rsp),%rbx
  nopl   0x0(%rax)
L2:
  pxor   %xmm0,%xmm0
  mov    %rbx,%rsi
  mov    %r14,%rdi
  add    $0x10,%r15
  cvtsi2sd %r12,%xmm0
  mulsd  (%rsp),%xmm0
  add    $0x1,%r12
  call   <sincos@plt>
  movsd  0x70(%rsp),%xmm0
  movsd  %xmm0,-0x10(%r15)
  movsd  0x78(%rsp),%xmm0
  movsd  %xmm0,-0x8(%r15)
  cmp    %r12,%rbp
  jae    L2
  mov    0x8(%rsp),%rax
  mov    0x58(%rsp),%rbx
  lea    0x2(%rax,%rbx,1),%rax
  mov    %rax,0x8(%rsp)
L74:
  cmpq   $0x6,0x20(%rsp)
  jbe    L44
  mov    0x8(%rsp),%rax
  mov    0x28(%rsp),%rdi
  movsd  0x0(%r13,%rax,8),%xmm0
  movsd  %xmm0,(%rdi)
  movsd  0x8(%r13,%rax,8),%xmm0
  mov    0x30(%rsp),%rax
  movsd  %xmm0,(%rax)
L44:
  addq   $0x1,0x10(%rsp)
  mov    0x48(%rsp),%rsi
  mov    0x10(%rsp),%rax
  add    %rsi,0x18(%rsp)
  cmp    %rax,0x20(%rsp)
  jne    L16
L75:
  addq   $0x8,0x60(%rsp)
  mov    0x60(%rsp),%rax
  cmp    %rax,0x68(%rsp)
  jne    L63
L37:
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L68:
  mov    0x18(%rsp),%rdi
  pxor   %xmm0,%xmm0
  mov    %rdi,%rax
  mov    %rdi,%rdx
  shr    %rax
  and    $0x1,%edx
  or     %rdx,%rax
  cvtsi2sd %rax,%xmm0
  addsd  %xmm0,%xmm0
  jmp    L24
L86:
  mov    0x40(%rsp),%rdx
  pxor   %xmm0,%xmm0
  mov    %rdx,%rax
  and    $0x1,%edx
  shr    %rax
  or     %rdx,%rax
  cvtsi2sd %rax,%xmm0
  addsd  %xmm0,%xmm0
  jmp    L92
```

tool_calls: name: extract_symbolic_data
arguments: tag: D20
type: f64
name: extract_symbolic_data
arguments: tag: D25
type: string
name: extract_symbolic_data
arguments: tag: D29
type: byte[8]

tool: D20: 6.283185307179586

tool: D25: 

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

assistant: ```c
void cffti1(size_t n, double wa[], size_t ifac[]) {
  static const size_t ntryh[5] = {4, 6, 3, 2, 5};
  static const double twopi = 6.28318530717958647692;
  size_t j, k, fi;
  double argh = twopi / n;
  size_t i = 0, l1 = 1;
  factorize(n, ntryh, 5, ifac);
  for (k = 1; k <= ifac[1]; k++) {
    size_t ip = ifac[k + 1];
    size_t ido = n / (l1 * ip);
    for (j = 1; j < ip; j++) {
      size_t is = i;
      double argld = j * l1 * argh;
      wa[i] = 1;
      wa[i + 1] = 0;
      for (fi = 1; fi <= ido; fi++) {
        double arg = fi * argld;
        i += 2;
        wa[i] = cos(arg);
        wa[i + 1] = sin(arg);
      }
      if (ip > 6) {
        wa[is] = wa[i];
        wa[is + 1] = wa[i + 1];
      }
    }
    l1 *= ip;
  }
}
```