user: Demonstrate the C code that generates the assembly output below:

```asm
<drftf1>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x28,%rsp
  mov    %edi,0xc(%rsp)
  mov    0x4(%r8),%eax
  test   %eax,%eax
  jle    L12
  mov    %rdx,%r14
  mov    %rcx,%r15
  movslq %eax,%rdx
  lea    0x4(,%rdx,4),%rdx
  lea    (%r8,%rdx,1),%rbx
  lea    -0x4(%r8,%rdx,1),%rdx
  lea    -0x1(%rax),%eax
  shl    $0x2,%rax
  sub    %rax,%rdx
  mov    %rdx,0x18(%rsp)
  mov    %edi,%ebp
  mov    $0x1,%r8d
  mov    %rsi,0x10(%rsp)
  mov    %edi,%r13d
  jmp    L110
L62:
  lea    (%rdi,%rbp,1),%eax
  lea    (%rdi,%rax,1),%edx
  test   %r12d,%r12d
  je     L90
  cltq
  movslq %ebp,%rcx
  sub    $0x8,%rsp
  movslq %edx,%rdx
  lea    -0x8(%r15,%rdx,8),%rdx
  push   %rdx
  lea    -0x8(%r15,%rax,8),%r9
  lea    -0x8(%r15,%rcx,8),%r8
  mov    0x20(%rsp),%rcx
  mov    %r14,%rdx
  mov    %r10d,%esi
  call   <dradf4@plt>
  add    $0x10,%rsp
  mov    %r12d,%r8d
L4:
  sub    $0x4,%rbx
  cmp    0x18(%rsp),%rbx
  je     L29
L110:
  mov    (%rbx),%esi
  mov    %r13d,%ecx
  mov    %r13d,%eax
  cltd
  idiv   %esi
  mov    %eax,%r10d
  mov    %eax,%r13d
  mov    0xc(%rsp),%eax
  cltd
  idiv   %ecx
  mov    %eax,%edi
  lea    -0x1(%rsi),%eax
  imul   %edi,%eax
  sub    %eax,%ebp
  mov    $0x1,%r12d
  sub    %r8d,%r12d
  cmp    $0x4,%esi
  je     L62
  cmp    $0x2,%esi
  jne    L46
  test   %r12d,%r12d
  je     L24
  movslq %ebp,%rax
  lea    -0x8(%r15,%rax,8),%r8
  mov    0x10(%rsp),%rcx
  mov    %r14,%rdx
  mov    %r10d,%esi
  call   <dradf2@plt>
  mov    %r12d,%r8d
  jmp    L4
L90:
  cltq
  movslq %ebp,%rcx
  sub    $0x8,%rsp
  movslq %edx,%rdx
  lea    -0x8(%r15,%rdx,8),%rdx
  push   %rdx
  lea    -0x8(%r15,%rax,8),%r9
  lea    -0x8(%r15,%rcx,8),%r8
  mov    %r14,%rcx
  mov    0x20(%rsp),%rdx
  mov    %r10d,%esi
  call   <dradf4@plt>
  add    $0x10,%rsp
  mov    %r12d,%r8d
  jmp    L4
L46:
  cmp    $0x1,%edi
  cmove  %r8d,%r12d
  mov    %r10d,%ecx
  imul   %edi,%ecx
  test   %r12d,%r12d
  je     L52
  movslq %ebp,%rax
  lea    -0x8(%r15,%rax,8),%rax
  push   %rax
  mov    0x18(%rsp),%rax
  push   %rax
  push   %rax
  push   %r14
  mov    %r14,%r9
  mov    %r14,%r8
  mov    %r10d,%edx
  call   <dradfg@plt>
  add    $0x20,%rsp
  mov    $0x0,%r8d
  jmp    L4
L24:
  movslq %ebp,%rax
  lea    -0x8(%r15,%rax,8),%r8
  mov    %r14,%rcx
  mov    0x10(%rsp),%rdx
  mov    %r10d,%esi
  call   <dradf2@plt>
  mov    %r12d,%r8d
  jmp    L4
L52:
  movslq %ebp,%rax
  lea    -0x8(%r15,%rax,8),%rax
  push   %rax
  push   %r14
  push   %r14
  mov    0x28(%rsp),%r9
  push   %r9
  mov    %r9,%r8
  mov    %r10d,%edx
  call   <dradfg@plt>
  add    $0x20,%rsp
  mov    $0x1,%r8d
  jmp    L4
L29:
  mov    0x10(%rsp),%r13
  cmp    $0x1,%r8d
  je     L12
  mov    0xc(%rsp),%eax
  test   %eax,%eax
  jle    L12
  mov    %eax,%edx
  mov    $0x0,%eax
L94:
  movsd  (%r14,%rax,8),%xmm0
  movsd  %xmm0,0x0(%r13,%rax,8)
  add    $0x1,%rax
  cmp    %rax,%rdx
  jne    L94
L12:
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void drftf1(int n, double *c, double *ch, double *wa, int *ifac) {
  int i, k1, l1, l2;
  int na, kh, nf;
  int ip, iw, ido, idl1, ix2, ix3;
  nf = ifac[1];
  na = 1;
  l2 = n;
  iw = n;
  for (k1 = 0; k1 < nf; k1++) {
    kh = nf - k1;
    ip = ifac[kh + 1];
    l1 = l2 / ip;
    ido = n / l2;
    idl1 = ido * l1;
    iw -= (ip - 1) * ido;
    na = 1 - na;
    if (ip != 4)
      goto L102;
    ix2 = iw + ido;
    ix3 = ix2 + ido;
    if (na != 0)
      dradf4(ido, l1, ch, c, wa + iw - 1, wa + ix2 - 1, wa + ix3 - 1);
    else
      dradf4(ido, l1, c, ch, wa + iw - 1, wa + ix2 - 1, wa + ix3 - 1);
    goto L110;
  L102:
    if (ip != 2)
      goto L104;
    if (na != 0)
      goto L103;
    dradf2(ido, l1, c, ch, wa + iw - 1);
    goto L110;
  L103:
    dradf2(ido, l1, ch, c, wa + iw - 1);
    goto L110;
  L104:
    if (ido == 1)
      na = 1 - na;
    if (na != 0)
      goto L109;
    dradfg(ido, ip, l1, idl1, c, c, c, ch, ch, wa + iw - 1);
    na = 1;
    goto L110;
  L109:
    dradfg(ido, ip, l1, idl1, ch, ch, ch, c, c, wa + iw - 1);
    na = 0;
  L110:
    l2 = l1;
  }
  if (na == 1)
    return;
  for (i = 0; i < n; i++)
    c[i] = ch[i];
}
```