// File: /workspace/github_c_makefile_star2_O3/pavel-demin[P]ft8d-c/pffft.c
// Function: radb4_ps
static NEVER_INLINE(void) radb4_ps(int ido, int l1, const v4sf * RESTRICT cc, v4sf * RESTRICT ch,
                                   const float * RESTRICT wa1, const float * RESTRICT wa2, const float *RESTRICT wa3)
{
  static const float minus_sqrt2 = (float)-1.414213562373095;
  static const float two = 2.f;
  int i, k, l1ido = l1*ido;
  v4sf ci2, ci3, ci4, cr2, cr3, cr4, ti1, ti2, ti3, ti4, tr1, tr2, tr3, tr4;
  {
    const v4sf *RESTRICT cc_ = cc, * RESTRICT ch_end = ch + l1ido;
    v4sf *ch_ = ch;
    while (ch < ch_end) {
      v4sf a = cc[0], b = cc[4*ido-1];
      v4sf c = cc[2*ido], d = cc[2*ido-1];
      tr3 = SVMUL(two,d);
      tr2 = VADD(a,b);
      tr1 = VSUB(a,b);
      tr4 = SVMUL(two,c);
      ch[0*l1ido] = VADD(tr2, tr3);
      ch[2*l1ido] = VSUB(tr2, tr3);
      ch[1*l1ido] = VSUB(tr1, tr4);
      ch[3*l1ido] = VADD(tr1, tr4);

      cc += 4*ido; ch += ido;
    }
    cc = cc_; ch = ch_;
  }
  if (ido < 2) return;
  if (ido != 2) {
    for (k = 0; k < l1ido; k += ido) {
      const v4sf * RESTRICT pc = (v4sf*)(cc - 1 + 4*k);
      v4sf * RESTRICT ph = (v4sf*)(ch + k + 1);
      for (i = 2; i < ido; i += 2) {

        tr1 = VSUB(pc[i], pc[4*ido - i]);
        tr2 = VADD(pc[i], pc[4*ido - i]);
        ti4 = VSUB(pc[2*ido + i], pc[2*ido - i]);
        tr3 = VADD(pc[2*ido + i], pc[2*ido - i]);
        ph[0] = VADD(tr2, tr3);
        cr3 = VSUB(tr2, tr3);

        ti3 = VSUB(pc[2*ido + i + 1], pc[2*ido - i + 1]);
        tr4 = VADD(pc[2*ido + i + 1], pc[2*ido - i + 1]);
        cr2 = VSUB(tr1, tr4);
        cr4 = VADD(tr1, tr4);

        ti1 = VADD(pc[i + 1], pc[4*ido - i + 1]);
        ti2 = VSUB(pc[i + 1], pc[4*ido - i + 1]);

        ph[1] = VADD(ti2, ti3); ph += l1ido;
        ci3 = VSUB(ti2, ti3);
        ci2 = VADD(ti1, ti4);
        ci4 = VSUB(ti1, ti4);
        VCPLXMUL(cr2, ci2, LD_PS1(wa1[i-2]), LD_PS1(wa1[i-1]));
        ph[0] = cr2;
        ph[1] = ci2; ph += l1ido;
        VCPLXMUL(cr3, ci3, LD_PS1(wa2[i-2]), LD_PS1(wa2[i-1]));
        ph[0] = cr3;
        ph[1] = ci3; ph += l1ido;
        VCPLXMUL(cr4, ci4, LD_PS1(wa3[i-2]), LD_PS1(wa3[i-1]));
        ph[0] = cr4;
        ph[1] = ci4; ph = ph - 3*l1ido + 2;
      }
    }
    if (ido % 2 == 1) return;
  }
  for (k=0; k < l1ido; k+=ido) {
    int i0 = 4*k + ido;
    v4sf c = cc[i0-1], d = cc[i0 + 2*ido-1];
    v4sf a = cc[i0+0], b = cc[i0 + 2*ido+0];
    tr1 = VSUB(c,d);
    tr2 = VADD(c,d);
    ti1 = VADD(b,a);
    ti2 = VSUB(b,a);
    ch[ido-1 + k + 0*l1ido] = VADD(tr2,tr2);
    ch[ido-1 + k + 1*l1ido] = SVMUL(minus_sqrt2, VSUB(ti1, tr1));
    ch[ido-1 + k + 2*l1ido] = VADD(ti2, ti2);
    ch[ido-1 + k + 3*l1ido] = SVMUL(minus_sqrt2, VADD(ti1, tr1));
  }
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
movq %rdx, %rbp
movq %rdi, -0x50(%rsp)
imull %edi, %esi
movq %rsi, -0x58(%rsp)
movslq %esi, %rax
movq %rax, %r10
shlq $0x4, %r10
testl %eax, %eax
jle 0x5cdc
leaq (%rcx,%r10), %rdx
movq %rdx, -0x48(%rsp)
movq -0x50(%rsp), %rbx
leal (,%rbx,4), %r11d
leal (%rbx,%rbx), %edx
movslq %edx, %rdx
movq -0x58(%rsp), %rdi
leal (%rdi,%rdi), %esi
leal (%rdi,%rdi,2), %edi
movslq %r11d, %r11
movslq %ebx, %rbx
movq %rbp, %r12
leal -0x4(,%rbx,4), %ebp
movslq %ebp, %r14
shlq $0x4, %r14
shlq $0x4, %r11
shlq $0x4, %rdx
leal -0x2(,%rbx,2), %ebp
movslq %ebp, %r15
movq %r12, %rbp
shlq $0x4, %r15
shlq $0x4, %rbx
shlq $0x4, %rdi
shlq $0x4, %rsi
movq %rcx, %r12
movq %rbp, %r13
movaps (%r13), %xmm0
movaps 0x30(%r13,%r14), %xmm1
movaps (%r13,%rdx), %xmm2
movaps 0x10(%r13,%r15), %xmm3
addps %xmm3, %xmm3
movaps %xmm0, %xmm4
addps %xmm1, %xmm4
subps %xmm1, %xmm0
addps %xmm2, %xmm2
movaps %xmm4, %xmm1
addps %xmm3, %xmm1
movaps %xmm1, (%r12)
subps %xmm3, %xmm4
movaps %xmm4, (%r12,%rsi)
movaps %xmm0, %xmm1
subps %xmm2, %xmm1
movaps %xmm1, (%r12,%r10)
addps %xmm0, %xmm2
movaps %xmm2, (%r12,%rdi)
addq %r11, %r13
addq %rbx, %r12
cmpq -0x48(%rsp), %r12
jb 0x5c83
movq -0x50(%rsp), %rdx
cmpl $0x2, %edx
movq -0x58(%rsp), %rsi
jl 0x5fe0
jne 0x5cfe
testl %esi, %esi
jg 0x5f34
jmp 0x5fe0
testl %esi, %esi
jle 0x5fe0
movq 0x38(%rsp), %r11
leal (,%rdx,4), %edi
leal (%rdx,%rdx), %ebx
leal (%rsi,%rsi,2), %esi
movl %edx, %r15d
shlq $0x4, %rbx
addq %rbp, %rbx
movq %rbx, -0x20(%rsp)
leaq (,%r15,4), %rdx
movq %rdx, -0x28(%rsp)
shlq $0x4, %rdi
movq %rbp, -0x40(%rsp)
addq %rbp, %rdi
movq %rdi, -0x18(%rsp)
movq %r15, %rdx
shlq $0x4, %rdx
movq %rdx, -0x30(%rsp)
leaq (%r10,%r10,2), %rbp
shlq $0x4, %rsi
movq %rbp, %r12
subq %rsi, %r12
addq $0x20, %r12
shlq $0x5, %rax
xorl %ebx, %ebx
movq %rcx, -0x38(%rsp)
movq %rcx, %r13
xorl %ecx, %ecx
movq %rcx, -0x10(%rsp)
movq %rbx, -0x48(%rsp)
shrl $0x2, %ebx
shlq $0x6, %rbx
movq -0x20(%rsp), %rcx
leaq (%rcx,%rbx), %r14
movq -0x18(%rsp), %rcx
leaq (%rcx,%rbx), %rdi
addq -0x40(%rsp), %rbx
movl $0x2, %esi
movq $-0x8, %rdx
movq %r13, -0x8(%rsp)
movq %r13, %rcx
xorl %r13d, %r13d
movaps 0x10(%rbx,%r13,4), %xmm1
movaps -0x10(%rdi,%rdx,4), %xmm0
movaps %xmm1, %xmm2
subps %xmm0, %xmm2
addps %xmm0, %xmm1
movaps 0x10(%r14,%r13,4), %xmm0
movaps -0x10(%r14,%rdx,4), %xmm3
movaps %xmm0, %xmm5
subps %xmm3, %xmm5
addps %xmm3, %xmm0
movaps %xmm1, %xmm3
addps %xmm0, %xmm3
movaps %xmm3, 0x10(%rcx)
subps %xmm0, %xmm1
movaps 0x20(%r14,%r13,4), %xmm0
movaps (%r14,%rdx,4), %xmm3
movaps %xmm0, %xmm6
subps %xmm3, %xmm6
addps %xmm3, %xmm0
movaps %xmm2, %xmm4
subps %xmm0, %xmm4
addps %xmm2, %xmm0
movaps 0x20(%rbx,%r13,4), %xmm3
movaps (%rdi,%rdx,4), %xmm7
movaps %xmm3, %xmm2
addps %xmm7, %xmm2
subps %xmm7, %xmm3
movaps %xmm6, %xmm7
addps %xmm3, %xmm7
movaps %xmm7, 0x20(%rcx)
subps %xmm6, %xmm3
movaps %xmm5, %xmm6
addps %xmm2, %xmm6
subps %xmm5, %xmm2
movss (%r8,%r13), %xmm5
movss 0x4(%r8,%r13), %xmm7
shufps $0x0, %xmm7, %xmm7      # xmm7 = xmm7[0,0,0,0]
movaps %xmm4, %xmm8
mulps %xmm7, %xmm8
shufps $0x0, %xmm5, %xmm5      # xmm5 = xmm5[0,0,0,0]
mulps %xmm5, %xmm4
mulps %xmm6, %xmm7
subps %xmm7, %xmm4
mulps %xmm6, %xmm5
addps %xmm8, %xmm5
movaps %xmm4, 0x10(%rcx,%r10)
movaps %xmm5, 0x20(%rcx,%r10)
movss (%r9,%r13), %xmm4
movss 0x4(%r9,%r13), %xmm5
shufps $0x0, %xmm5, %xmm5      # xmm5 = xmm5[0,0,0,0]
movaps %xmm1, %xmm6
mulps %xmm5, %xmm6
shufps $0x0, %xmm4, %xmm4      # xmm4 = xmm4[0,0,0,0]
mulps %xmm4, %xmm1
mulps %xmm3, %xmm5
subps %xmm5, %xmm1
mulps %xmm3, %xmm4
addps %xmm6, %xmm4
movaps %xmm1, 0x10(%rcx,%rax)
movaps %xmm4, 0x20(%rcx,%rax)
movss (%r11,%r13), %xmm1
movss 0x4(%r11,%r13), %xmm3
shufps $0x0, %xmm3, %xmm3      # xmm3 = xmm3[0,0,0,0]
movaps %xmm0, %xmm4
mulps %xmm3, %xmm4
shufps $0x0, %xmm1, %xmm1      # xmm1 = xmm1[0,0,0,0]
mulps %xmm1, %xmm0
mulps %xmm2, %xmm3
subps %xmm3, %xmm0
mulps %xmm2, %xmm1
addps %xmm4, %xmm1
movaps %xmm0, 0x10(%rcx,%rbp)
movaps %xmm1, 0x20(%rcx,%rbp)
addq $0x2, %rsi
addq $-0x8, %rdx
addq $0x8, %r13
addq %r12, %rcx
cmpq %r15, %rsi
jb 0x5db2
movq -0x10(%rsp), %rcx
addq %r15, %rcx
movq -0x48(%rsp), %rbx
addq -0x28(%rsp), %rbx
movq -0x8(%rsp), %r13
addq -0x30(%rsp), %r13
cmpl %ecx, -0x58(%rsp)
jg 0x5d73
movq -0x58(%rsp), %rsi
testl %esi, %esi
setg %al
movq -0x50(%rsp), %rdx
testb $0x1, %dl
sete %cl
testb %al, %cl
movq -0x38(%rsp), %rcx
movq -0x40(%rsp), %rbp
je 0x5fe0
leal -0x1(%rdx), %edi
leal (%rsi,%rsi), %eax
leal (%rsi,%rsi,2), %r14d
movl %edx, %r15d
movl %esi, %r12d
shlq $0x4, %r14
shlq $0x4, %rdi
addq %rdi, %rcx
shlq $0x4, %r15
shlq $0x4, %rax
shlq $0x4, %r12
leal (,%rdx,4), %edi
leal (%rdx,%rdx,2), %r8d
xorl %r9d, %r9d
movaps 0x22a0(%rip), %xmm0     # 0x8210
movl %edx, %r10d
movl %r10d, %r11d
shlq $0x4, %r11
movl %r8d, %ebx
shlq $0x4, %rbx
movaps -0x10(%rbp,%r11), %xmm1
movaps (%rbp,%r11), %xmm2
movaps -0x10(%rbp,%rbx), %xmm3
movaps (%rbp,%rbx), %xmm4
movaps %xmm1, %xmm5
subps %xmm3, %xmm5
addps %xmm3, %xmm1
movaps %xmm2, %xmm3
addps %xmm4, %xmm3
subps %xmm2, %xmm4
addps %xmm1, %xmm1
movaps %xmm1, (%rcx)
movaps %xmm3, %xmm1
subps %xmm5, %xmm1
mulps %xmm0, %xmm1
movaps %xmm1, (%rcx,%r12)
addps %xmm4, %xmm4
movaps %xmm4, (%rcx,%rax)
addps %xmm5, %xmm3
mulps %xmm0, %xmm3
movaps %xmm3, (%rcx,%r14)
addq %r15, %rcx
addl %edi, %r10d
addl %edi, %r8d
addl %edx, %r9d
cmpl %r9d, %esi
jg 0x5f73
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq