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

```asm
<cfftf1>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x58,%rsp
  mov    %edi,0x14(%rsp)
  mov    %rsi,0x8(%rsp)
  mov    %r8,0x20(%rsp)
  mov    %r9d,0x1c(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x48(%rsp)
  mov    0x4(%r8),%eax
  mov    %eax,0x18(%rsp)
  test   %eax,%eax
  jle    L171
  mov    0x14(%rsp),%eax
  mov    0x8(%r8),%ebx
  mov    %rdx,%r10
  mov    %rcx,%r11
  lea    0x44(%rsp),%rcx
  mov    %r10,0x28(%rsp)
  xor    %r15d,%r15d
  xor    %r13d,%r13d
  cltd
  mov    %rcx,0x38(%rsp)
  mov    %ebx,%r14d
  mov    $0x2,%r12d
  idiv   %ebx
  mov    0x8(%rsp),%rdx
  mov    $0x1,%esi
  lea    (%rax,%rax,1),%ebp
  mov    %ebp,%edi
  jmp    L78
  cs nopw 0x0(%rax,%rax,1)
L39:
  cmp    $0x3,%ebx
  jne    L181
  sub    $0x8,%rsp
  lea    0x0(%rbp,%r13,1),%ecx
  mov    %ebp,%edi
  xor    $0x1,%r15d
  mov    0x24(%rsp),%eax
  movslq %ecx,%rcx
  lea    (%r11,%rcx,8),%r9
  mov    %r10,%rcx
  push   %rax
  mov    %r11,0x40(%rsp)
  call   <passf3@plt>
  pop    %rdx
  pop    %rcx
  mov    0x30(%rsp),%r11
L107:
  sub    $0x1,%ebx
  imul   %ebp,%ebx
  add    %ebx,%r13d
  cmp    %r12d,0x18(%rsp)
  jl     L176
L182:
  mov    0x20(%rsp),%rax
  mov    %r14d,%esi
  mov    0x4(%rax,%r12,4),%ebx
  mov    0x14(%rsp),%eax
  mov    %ebx,%ecx
  cltd
  imul   %r14d,%ecx
  idiv   %ecx
  mov    0x8(%rsp),%rdx
  lea    (%rax,%rax,1),%ebp
  mov    0x28(%rsp),%rax
  mov    %ebp,%edi
  imul   %r14d,%edi
  mov    %rax,%r10
  test   %r15d,%r15d
  mov    %ecx,%r14d
  cmovne %rdx,%r10
  cmovne %rax,%rdx
  add    $0x1,%r12
L78:
  movslq %r13d,%rcx
  lea    (%r11,%rcx,8),%r8
  cmp    $0x4,%ebx
  je     L145
  jg     L127
  cmp    $0x2,%ebx
  jne    L39
  sub    $0x1,%ebx
  mov    0x1c(%rsp),%r9d
  mov    %r10,%rcx
  mov    %ebp,%edi
  imul   %ebp,%ebx
  mov    %r11,0x30(%rsp)
  xor    $0x1,%r15d
  call   <passf2@plt>
  mov    0x30(%rsp),%r11
  add    %ebx,%r13d
  cmp    %r12d,0x18(%rsp)
  jge    L182
  xchg   %ax,%ax
L176:
  mov    0x28(%rsp),%r10
  test   %r15d,%r15d
  je     L171
  mov    0x14(%rsp),%edx
  add    %edx,%edx
  test   %edx,%edx
  jle    L171
  mov    0x8(%rsp),%rcx
  lea    0x8(%r10),%rax
  cmp    %rax,%rcx
  je     L72
  lea    -0x1(%rdx),%eax
  cmp    $0x2,%eax
  jbe    L72
  shr    %edx
  xor    %eax,%eax
  shl    $0x4,%rdx
  nopw   0x0(%rax,%rax,1)
L18:
  movupd (%r10,%rax,1),%xmm1
  movups %xmm1,(%rcx,%rax,1)
  add    $0x10,%rax
  cmp    %rax,%rdx
  jne    L18
L171:
  mov    0x48(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L43
  add    $0x58,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  cs nopw 0x0(%rax,%rax,1)
L127:
  cmp    $0x5,%ebx
  jne    L181
  lea    0x0(%rbp,%r13,1),%ecx
  sub    $0x8,%rsp
  xor    $0x1,%r15d
  lea    (%rcx,%rbp,1),%edi
  mov    0x24(%rsp),%eax
  movslq %ecx,%rcx
  lea    (%rdi,%rbp,1),%r9d
  movslq %edi,%rdi
  movslq %r9d,%r9
  lea    (%r11,%rdi,8),%rdi
  push   %rax
  lea    (%r11,%r9,8),%r9
  push   %r9
  lea    (%r11,%rcx,8),%r9
  mov    %r10,%rcx
  push   %rdi
  mov    %ebp,%edi
  mov    %r11,0x50(%rsp)
  call   <passf5@plt>
  add    $0x20,%rsp
  mov    0x30(%rsp),%r11
  jmp    L107
  nopw   0x0(%rax,%rax,1)
L145:
  lea    0x0(%rbp,%r13,1),%ecx
  mov    0x1c(%rsp),%eax
  xor    $0x1,%r15d
  movslq %ecx,%rdi
  add    %ebp,%ecx
  movslq %ecx,%rcx
  push   %rax
  lea    (%r11,%rdi,8),%r9
  mov    %ebp,%edi
  lea    (%r11,%rcx,8),%rcx
  push   %rcx
  mov    %r10,%rcx
  mov    %r11,0x40(%rsp)
  call   <passf4@plt>
  pop    %rsi
  pop    %rdi
  mov    0x30(%rsp),%r11
  jmp    L107
  nopw   0x0(%rax,%rax,1)
L181:
  mov    %r11,0x30(%rsp)
  sub    $0x8,%rsp
  mov    %rdx,%r9
  mov    %esi,%ecx
  mov    0x24(%rsp),%eax
  mov    %ebx,%edx
  mov    %ebp,%esi
  push   %rax
  push   %r8
  mov    %edi,%r8d
  push   %r10
  mov    0x58(%rsp),%rdi
  call   <passf@plt>
  mov    0x64(%rsp),%eax
  add    $0x20,%rsp
  mov    0x30(%rsp),%r11
  test   %eax,%eax
  je     L107
  xor    $0x1,%r15d
  jmp    L107
  cs nopw 0x0(%rax,%rax,1)
L72:
  mov    0x8(%rsp),%rcx
  xor    %eax,%eax
  nopw   0x0(%rax,%rax,1)
L121:
  movsd  (%r10,%rax,8),%xmm0
  movsd  %xmm0,(%rcx,%rax,8)
  add    $0x1,%rax
  cmp    %rdx,%rax
  jne    L121
  jmp    L171
L43:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
void cfftf1(int n, double c[], double ch[], const double wa[], const int ifac[13 + 2], int isign) {
  int idot, i;
  int k1, l1, l2;
  int na, nf, ip, iw, ix2, ix3, ix4, nac, ido, idl1;
  double *cinput, *coutput;
  nf = ifac[1];
  na = 0;
  l1 = 1;
  iw = 0;
  for (k1 = 2; k1 <= nf + 1; k1++) {
    ip = ifac[k1];
    l2 = ip * l1;
    ido = n / l2;
    idot = ido + ido;
    idl1 = idot * l1;
    if (na) {
      cinput = ch;
      coutput = c;
    } else {
      cinput = c;
      coutput = ch;
    }
    switch (ip) {
    case 4:
      ix2 = iw + idot;
      ix3 = ix2 + idot;
      passf4(idot, l1, cinput, coutput, &wa[iw], &wa[ix2], &wa[ix3], isign);
      na = !na;
      break;
    case 2:
      passf2(idot, l1, cinput, coutput, &wa[iw], isign);
      na = !na;
      break;
    case 3:
      ix2 = iw + idot;
      passf3(idot, l1, cinput, coutput, &wa[iw], &wa[ix2], isign);
      na = !na;
      break;
    case 5:
      ix2 = iw + idot;
      ix3 = ix2 + idot;
      ix4 = ix3 + idot;
      passf5(idot, l1, cinput, coutput, &wa[iw], &wa[ix2], &wa[ix3], &wa[ix4], isign);
      na = !na;
      break;
    default:
      passf(&nac, idot, ip, l1, idl1, cinput, coutput, &wa[iw], isign);
      if (nac != 0)
        na = !na;
    }
    l1 = l2;
    iw += (ip - 1) * idot;
  }
  if (na == 0)
    return;
  for (i = 0; i < 2 * n; i++)
    c[i] = ch[i];
}
```