user: 下面的汇编指令对应的C源代码是什么:

```asm
<rfftf1>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x90,%rsp
  mov    %rdi,-0x68(%rbp)
  mov    %rsi,-0x70(%rbp)
  mov    %rdx,-0x78(%rbp)
  mov    %rcx,-0x80(%rbp)
  mov    %r8,-0x88(%rbp)
  mov    -0x68(%rbp),%rax
  mov    %rax,-0x50(%rbp)
  mov    -0x88(%rbp),%rax
  add    $0x8,%rax
  mov    (%rax),%rax
  mov    %rax,-0x30(%rbp)
  mov    -0x68(%rbp),%rax
  sub    $0x1,%rax
  mov    %rax,-0x48(%rbp)
  mov    -0x78(%rbp),%rax
  mov    %rax,-0x40(%rbp)
  mov    -0x70(%rbp),%rax
  mov    %rax,-0x38(%rbp)
  movq   $0x1,-0x58(%rbp)
  jmp    L116
L48:
  mov    -0x30(%rbp),%rax
  sub    -0x58(%rbp),%rax
  add    $0x2,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  mov    %rax,-0x28(%rbp)
  mov    -0x68(%rbp),%rax
  mov    $0x0,%edx
  divq   -0x50(%rbp)
  mov    %rax,-0x20(%rbp)
  mov    -0x50(%rbp),%rax
  mov    $0x0,%edx
  divq   -0x28(%rbp)
  mov    %rax,-0x50(%rbp)
  mov    -0x28(%rbp),%rax
  sub    $0x1,%rax
  imul   -0x20(%rbp),%rax
  sub    %rax,-0x48(%rbp)
  mov    -0x40(%rbp),%rax
  mov    %rax,-0x18(%rbp)
  mov    -0x38(%rbp),%rax
  mov    %rax,-0x40(%rbp)
  mov    -0x18(%rbp),%rax
  mov    %rax,-0x38(%rbp)
  cmpq   $0x4,-0x28(%rbp)
  jne    L118
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  lea    (%rdx,%rax,1),%rdi
  mov    -0x38(%rbp),%rcx
  mov    -0x40(%rbp),%rdx
  mov    -0x50(%rbp),%rsi
  mov    -0x20(%rbp),%rax
  mov    %rdi,%r8
  mov    %rax,%rdi
  call   <radf4@plt>
  jmp    L100
L118:
  cmpq   $0x2,-0x28(%rbp)
  jne    L117
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  lea    (%rdx,%rax,1),%rdi
  mov    -0x38(%rbp),%rcx
  mov    -0x40(%rbp),%rdx
  mov    -0x50(%rbp),%rsi
  mov    -0x20(%rbp),%rax
  mov    %rdi,%r8
  mov    %rax,%rdi
  call   <radf2@plt>
  jmp    L100
L117:
  cmpq   $0x3,-0x28(%rbp)
  jne    L62
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  lea    (%rdx,%rax,1),%rdi
  mov    -0x38(%rbp),%rcx
  mov    -0x40(%rbp),%rdx
  mov    -0x50(%rbp),%rsi
  mov    -0x20(%rbp),%rax
  mov    %rdi,%r8
  mov    %rax,%rdi
  call   <radf3@plt>
  jmp    L100
L62:
  cmpq   $0x5,-0x28(%rbp)
  jne    L61
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  lea    (%rdx,%rax,1),%rdi
  mov    -0x38(%rbp),%rcx
  mov    -0x40(%rbp),%rdx
  mov    -0x50(%rbp),%rsi
  mov    -0x20(%rbp),%rax
  mov    %rdi,%r8
  mov    %rax,%rdi
  call   <radf5@plt>
  jmp    L100
L61:
  cmpq   $0x1,-0x20(%rbp)
  jne    L21
  mov    -0x40(%rbp),%rax
  mov    %rax,-0x10(%rbp)
  mov    -0x38(%rbp),%rax
  mov    %rax,-0x40(%rbp)
  mov    -0x10(%rbp),%rax
  mov    %rax,-0x38(%rbp)
L21:
  mov    -0x48(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  lea    (%rdx,%rax,1),%rdi
  mov    -0x20(%rbp),%rax
  imul   -0x50(%rbp),%rax
  mov    %rax,%rcx
  mov    -0x38(%rbp),%r9
  mov    -0x40(%rbp),%r8
  mov    -0x50(%rbp),%rdx
  mov    -0x28(%rbp),%rsi
  mov    -0x20(%rbp),%rax
  sub    $0x8,%rsp
  push   %rdi
  mov    %rax,%rdi
  call   <radfg@plt>
  add    $0x10,%rsp
  mov    -0x40(%rbp),%rax
  mov    %rax,-0x8(%rbp)
  mov    -0x38(%rbp),%rax
  mov    %rax,-0x40(%rbp)
  mov    -0x8(%rbp),%rax
  mov    %rax,-0x38(%rbp)
L100:
  addq   $0x1,-0x58(%rbp)
L116:
  mov    -0x58(%rbp),%rax
  cmp    -0x30(%rbp),%rax
  jbe    L48
  mov    -0x40(%rbp),%rax
  cmp    -0x70(%rbp),%rax
  jne    L71
  mov    -0x68(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x78(%rbp),%rcx
  mov    -0x70(%rbp),%rax
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <memcpy@plt>
L71:
  nop
  leave
  ret
```

assistant: ```c
void rfftf1(size_t n, double c[], double ch[], const double wa[], const size_t ifac[]) {
  size_t k1, l1 = n, nf = ifac[1], iw = n - 1;
  double *p1 = ch, *p2 = c;
  for (k1 = 1; k1 <= nf; ++k1) {
    size_t ip = ifac[nf - k1 + 2];
    size_t ido = n / l1;
    l1 /= ip;
    iw -= (ip - 1) * ido;
    do {
      double *tmp_ = (p1);
      (p1) = (p2);
      (p2) = tmp_;
    } while (0);
    if (ip == 4)
      radf4(ido, l1, p1, p2, wa + iw);
    else if (ip == 2)
      radf2(ido, l1, p1, p2, wa + iw);
    else if (ip == 3)
      radf3(ido, l1, p1, p2, wa + iw);
    else if (ip == 5)
      radf5(ido, l1, p1, p2, wa + iw);
    else {
      if (ido == 1)
        do {
          double *tmp_ = (p1);
          (p1) = (p2);
          (p2) = tmp_;
        } while (0);
      radfg(ido, ip, l1, ido * l1, p1, p2, wa + iw);
      do {
        double *tmp_ = (p1);
        (p1) = (p2);
        (p2) = tmp_;
      } while (0);
    }
  }
  if (p1 == c)
    memcpy(c, ch, n * sizeof(double));
}
```