// File: /workspace/github_c_makefile_star2_O0/diaxen[P]fft-garden/ct_dit_dr_P.c
// Function: ct_fft_dr
static inline void
ct_fft_dr(const complex double * restrict in,
           complex double * restrict out,
           int s, int log2_N)
{
  unsigned l = 1 << (log2_N - s - 1);
  unsigned stride = 1 << s;

  switch (l)
    {
    case 1:                     /* size 2 base case  */
      ct_dit_bf2_0(1, stride, out, in);
      break;

    case 2:                     /* size 4 base case  */
      // ct_fft_dr(in, out, s + 1, log2_N);
      ct_dit_bf2_0(1, stride * 2, out, in);
      // ct_fft_dr(in + stride, out + l, s + 1, log2_N);
      ct_dit_bf2_0(1, stride * 2, out + 2, in + stride);

      ct_dit_bf2_0(2, 2, out, out);
      ct_dit_bf2_pi2(2, out + 1);
      break;

    default:
      ct_fft_dr(in, out, s + 1, log2_N);
      ct_fft_dr(in + stride, out + l, s + 1, log2_N);

      ct_dit_bf2_0(l, l, out, out);
      ct_dit_bf2_pi2(l, out + l / 2);
      ct_dit_bf2_pi4(l, out + l / 4);
      ct_dit_bf2_3pi4(l, out + 3 * l / 4);

      unsigned k;
      for (k = 1; k < l / 4; k++)
        ct_dit_bf2(l, out + k, get_twiddle(k << s, log2_N));
      for (k++; k < l / 2; k++)
        ct_dit_bf2(l, out + k, get_twiddle(k << s, log2_N));
      for (k++; k < 3 * l / 4; k++)
        ct_dit_bf2(l, out + k, get_twiddle(k << s, log2_N));
      for (k++; k < l; k++)
        ct_dit_bf2(l, out + k, get_twiddle(k << s, log2_N));
    }
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0xb0, %rsp
movq %rdi, -0x8(%rbp)
movq %rsi, -0x10(%rbp)
movl %edx, -0x14(%rbp)
movl %ecx, -0x18(%rbp)
movl -0x18(%rbp), %eax
movl -0x14(%rbp), %ecx
subl %ecx, %eax
movb %al, %cl
decb %cl
movl $0x1, %eax
movl %eax, %edx
shll %cl, %edx
movl %edx, %ecx
movl %ecx, -0x1c(%rbp)
movb -0x14(%rbp), %cl
shll %cl, %eax
movl %eax, -0x20(%rbp)
movl -0x1c(%rbp), %eax
movl %eax, -0x6c(%rbp)
subl $0x1, %eax
je 0x1ee5
jmp 0x1ed8
movl -0x6c(%rbp), %eax
subl $0x2, %eax
je 0x1eff
jmp 0x1f66
movl -0x20(%rbp), %esi
movq -0x10(%rbp), %rdx
movq -0x8(%rbp), %rcx
movl $0x1, %edi
callq 0x2210
jmp 0x21cd
movl -0x20(%rbp), %esi
shll %esi
movq -0x10(%rbp), %rdx
movq -0x8(%rbp), %rcx
movl $0x1, %edi
callq 0x2210
movl -0x20(%rbp), %esi
shll %esi
movq -0x10(%rbp), %rdx
addq $0x20, %rdx
movq -0x8(%rbp), %rcx
movl -0x20(%rbp), %eax
shlq $0x4, %rax
addq %rax, %rcx
movl $0x1, %edi
callq 0x2210
movq -0x10(%rbp), %rdx
movq -0x10(%rbp), %rcx
movl $0x2, %esi
movl %esi, %edi
callq 0x2210
movq -0x10(%rbp), %rsi
addq $0x10, %rsi
movl $0x2, %edi
callq 0x22c0
jmp 0x21cd
movq -0x8(%rbp), %rdi
movq -0x10(%rbp), %rsi
movl -0x14(%rbp), %edx
addl $0x1, %edx
movl -0x18(%rbp), %ecx
callq 0x1e90
movq -0x8(%rbp), %rdi
movl -0x20(%rbp), %eax
shlq $0x4, %rax
addq %rax, %rdi
movq -0x10(%rbp), %rsi
movl -0x1c(%rbp), %eax
shlq $0x4, %rax
addq %rax, %rsi
movl -0x14(%rbp), %edx
addl $0x1, %edx
movl -0x18(%rbp), %ecx
callq 0x1e90
movl -0x1c(%rbp), %edi
movl -0x1c(%rbp), %esi
movq -0x10(%rbp), %rdx
movq -0x10(%rbp), %rcx
callq 0x2210
movl -0x1c(%rbp), %edi
movq -0x10(%rbp), %rsi
movl -0x1c(%rbp), %eax
shrl %eax
movl %eax, %eax
shlq $0x4, %rax
addq %rax, %rsi
callq 0x22c0
movl -0x1c(%rbp), %edi
movq -0x10(%rbp), %rsi
movl -0x1c(%rbp), %eax
shrl $0x2, %eax
movl %eax, %eax
shlq $0x4, %rax
addq %rax, %rsi
callq 0x23a0
movl -0x1c(%rbp), %edi
movq -0x10(%rbp), %rsi
imull $0x3, -0x1c(%rbp), %eax
shrl $0x2, %eax
movl %eax, %eax
shlq $0x4, %rax
addq %rax, %rsi
callq 0x2490
movl $0x1, -0x24(%rbp)
movl -0x24(%rbp), %eax
movl -0x1c(%rbp), %ecx
shrl $0x2, %ecx
cmpl %ecx, %eax
jae 0x2071
movl -0x1c(%rbp), %eax
movl %eax, -0x7c(%rbp)
movq -0x10(%rbp), %rax
movl -0x24(%rbp), %ecx
shlq $0x4, %rcx
addq %rcx, %rax
movq %rax, -0x78(%rbp)
movl -0x24(%rbp), %edi
movl -0x14(%rbp), %ecx
shll %cl, %edi
movl -0x18(%rbp), %esi
callq 0x2710
movl -0x7c(%rbp), %edi
movq -0x78(%rbp), %rsi
movsd %xmm0, -0x38(%rbp)
movsd %xmm1, -0x30(%rbp)
movsd -0x38(%rbp), %xmm0
movsd -0x30(%rbp), %xmm1
callq 0x25b0
movl -0x24(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x24(%rbp)
jmp 0x2011
movl -0x24(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x24(%rbp)
movl -0x24(%rbp), %eax
movl -0x1c(%rbp), %ecx
shrl %ecx
cmpl %ecx, %eax
jae 0x20e5
movl -0x1c(%rbp), %eax
movl %eax, -0x8c(%rbp)
movq -0x10(%rbp), %rax
movl -0x24(%rbp), %ecx
shlq $0x4, %rcx
addq %rcx, %rax
movq %rax, -0x88(%rbp)
movl -0x24(%rbp), %edi
movl -0x14(%rbp), %ecx
shll %cl, %edi
movl -0x18(%rbp), %esi
callq 0x2710
movl -0x8c(%rbp), %edi
movq -0x88(%rbp), %rsi
movsd %xmm0, -0x48(%rbp)
movsd %xmm1, -0x40(%rbp)
movsd -0x48(%rbp), %xmm0
movsd -0x40(%rbp), %xmm1
callq 0x25b0
movl -0x24(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x24(%rbp)
jmp 0x207a
movl -0x24(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x24(%rbp)
movl -0x24(%rbp), %eax
imull $0x3, -0x1c(%rbp), %ecx
shrl $0x2, %ecx
cmpl %ecx, %eax
jae 0x215b
movl -0x1c(%rbp), %eax
movl %eax, -0x9c(%rbp)
movq -0x10(%rbp), %rax
movl -0x24(%rbp), %ecx
shlq $0x4, %rcx
addq %rcx, %rax
movq %rax, -0x98(%rbp)
movl -0x24(%rbp), %edi
movl -0x14(%rbp), %ecx
shll %cl, %edi
movl -0x18(%rbp), %esi
callq 0x2710
movl -0x9c(%rbp), %edi
movq -0x98(%rbp), %rsi
movsd %xmm0, -0x58(%rbp)
movsd %xmm1, -0x50(%rbp)
movsd -0x58(%rbp), %xmm0
movsd -0x50(%rbp), %xmm1
callq 0x25b0
movl -0x24(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x24(%rbp)
jmp 0x20ee
movl -0x24(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x24(%rbp)
movl -0x24(%rbp), %eax
cmpl -0x1c(%rbp), %eax
jae 0x21cb
movl -0x1c(%rbp), %eax
movl %eax, -0xac(%rbp)
movq -0x10(%rbp), %rax
movl -0x24(%rbp), %ecx
shlq $0x4, %rcx
addq %rcx, %rax
movq %rax, -0xa8(%rbp)
movl -0x24(%rbp), %edi
movl -0x14(%rbp), %ecx
shll %cl, %edi
movl -0x18(%rbp), %esi
callq 0x2710
movl -0xac(%rbp), %edi
movq -0xa8(%rbp), %rsi
movsd %xmm0, -0x68(%rbp)
movsd %xmm1, -0x60(%rbp)
movsd -0x68(%rbp), %xmm0
movsd -0x60(%rbp), %xmm1
callq 0x25b0
movl -0x24(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x24(%rbp)
jmp 0x2164
jmp 0x21cd
addq $0xb0, %rsp
popq %rbp
retq
nopw %cs:(%rax,%rax)