user: 请从给定的汇编指令推导出C源代码：

```asm
<cffti1>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %rbx
  sub    $0x88,%rsp
  mov    %rdi,-0x78(%rbp)
  mov    %rsi,-0x80(%rbp)
  mov    %rdx,-0x88(%rbp)
  movsd  D27(%rip),%xmm1
  mov    -0x78(%rbp),%rax
  test   %rax,%rax
  js     L123
  pxor   %xmm0,%xmm0
  cvtsi2sd %rax,%xmm0
  jmp    L12
L123:
  mov    %rax,%rdx
  shr    %rdx
  and    $0x1,%eax
  or     %rax,%rdx
  pxor   %xmm0,%xmm0
  cvtsi2sd %rdx,%xmm0
  addsd  %xmm0,%xmm0
L12:
  divsd  %xmm0,%xmm1
  movsd  %xmm1,-0x40(%rbp)
  movq   $0x0,-0x50(%rbp)
  movq   $0x1,-0x48(%rbp)
  mov    -0x88(%rbp),%rdx
  mov    -0x78(%rbp),%rax
  mov    %rdx,%rcx
  mov    $0x5,%edx
  lea    D2(%rip),%rsi
  mov    %rax,%rdi
  call   <factorize@plt>
  movq   $0x1,-0x60(%rbp)
  jmp    L48
L51:
  mov    -0x60(%rbp),%rax
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  mov    %rax,-0x38(%rbp)
  mov    -0x48(%rbp),%rax
  imul   -0x38(%rbp),%rax
  mov    %rax,%rdi
  mov    -0x78(%rbp),%rax
  mov    $0x0,%edx
  div    %rdi
  mov    %rax,-0x30(%rbp)
  movq   $0x1,-0x68(%rbp)
  jmp    L45
L30:
  mov    -0x50(%rbp),%rax
  mov    %rax,-0x28(%rbp)
  mov    -0x68(%rbp),%rax
  imul   -0x48(%rbp),%rax
  test   %rax,%rax
  js     L71
  pxor   %xmm0,%xmm0
  cvtsi2sd %rax,%xmm0
  jmp    L98
L71:
  mov    %rax,%rdx
  shr    %rdx
  and    $0x1,%eax
  or     %rax,%rdx
  pxor   %xmm0,%xmm0
  cvtsi2sd %rdx,%xmm0
  addsd  %xmm0,%xmm0
L98:
  movsd  -0x40(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0x50(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  D14(%rip),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x50(%rbp),%rax
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rax)
  movq   $0x1,-0x58(%rbp)
  jmp    L79
L49:
  mov    -0x58(%rbp),%rax
  test   %rax,%rax
  js     L95
  pxor   %xmm0,%xmm0
  cvtsi2sd %rax,%xmm0
  jmp    L93
L95:
  mov    %rax,%rdx
  shr    %rdx
  and    $0x1,%eax
  or     %rax,%rdx
  pxor   %xmm0,%xmm0
  cvtsi2sd %rdx,%xmm0
  addsd  %xmm0,%xmm0
L93:
  movsd  -0x20(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  addq   $0x2,-0x50(%rbp)
  mov    -0x50(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  lea    (%rdx,%rax,1),%rbx
  mov    -0x18(%rbp),%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  movq   %xmm0,%rax
  mov    %rax,(%rbx)
  mov    -0x50(%rbp),%rax
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  lea    (%rdx,%rax,1),%rbx
  mov    -0x18(%rbp),%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movq   %xmm0,%rax
  mov    %rax,(%rbx)
  addq   $0x1,-0x58(%rbp)
L79:
  mov    -0x58(%rbp),%rax
  cmp    -0x30(%rbp),%rax
  jbe    L49
  cmpq   $0x6,-0x38(%rbp)
  jbe    L72
  mov    -0x50(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rax,%rdx
  mov    -0x28(%rbp),%rax
  lea    0x0(,%rax,8),%rcx
  mov    -0x80(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x50(%rbp),%rax
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rax,%rdx
  mov    -0x28(%rbp),%rax
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rcx
  mov    -0x80(%rbp),%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
L72:
  addq   $0x1,-0x68(%rbp)
L45:
  mov    -0x68(%rbp),%rax
  cmp    -0x38(%rbp),%rax
  jb     L30
  mov    -0x48(%rbp),%rax
  imul   -0x38(%rbp),%rax
  mov    %rax,-0x48(%rbp)
  addq   $0x1,-0x60(%rbp)
L48:
  mov    -0x88(%rbp),%rax
  add    $0x8,%rax
  mov    (%rax),%rax
  cmp    %rax,-0x60(%rbp)
  jbe    L51
  nop
  nop
  mov    -0x8(%rbp),%rbx
  leave
  ret
```

tool_calls: name: unpack_data
arguments: tag: D27
data_type: f64
name: unpack_data
arguments: tag: D2
data_type: string
name: unpack_data
arguments: tag: D14
data_type: byte[8]

tool: D27: 6.283185307179586

tool: D2: 

tool: {"D14": [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;
  }
}
```