// File: /workspace/github_cpp_makefile_star2_O2/mouda[P]LTE_random_access-simulate/alglibinternal.cpp
// Function: alglib_impl::ftbase_ftbasegenerateplanrec(long, long, alglib_impl::ftplan*, long*, long*, long*, long*, long*, long, alglib_impl::ae_state*)
static void ftbase_ftbasegenerateplanrec(ae_int_t n,
     ae_int_t tasktype,
     ftplan* plan,
     ae_int_t* plansize,
     ae_int_t* precomputedsize,
     ae_int_t* planarraysize,
     ae_int_t* tmpmemsize,
     ae_int_t* stackmemsize,
     ae_int_t stackptr,
     ae_state *_state)
{
    ae_int_t k;
    ae_int_t m;
    ae_int_t n1;
    ae_int_t n2;
    ae_int_t esize;
    ae_int_t entryoffset;


    
    /*
     * prepare
     */
    if( *plansize+ftbase_ftbaseplanentrysize>(*planarraysize) )
    {
        ftbase_fftarrayresize(&plan->plan, planarraysize, 8*(*planarraysize), _state);
    }
    entryoffset = *plansize;
    esize = ftbase_ftbaseplanentrysize;
    *plansize = *plansize+esize;
    
    /*
     * if N=1, generate empty plan and exit
     */
    if( n==1 )
    {
        plan->plan.ptr.p_int[entryoffset+0] = esize;
        plan->plan.ptr.p_int[entryoffset+1] = -1;
        plan->plan.ptr.p_int[entryoffset+2] = -1;
        plan->plan.ptr.p_int[entryoffset+3] = ftbase_fftemptyplan;
        plan->plan.ptr.p_int[entryoffset+4] = -1;
        plan->plan.ptr.p_int[entryoffset+5] = -1;
        plan->plan.ptr.p_int[entryoffset+6] = -1;
        plan->plan.ptr.p_int[entryoffset+7] = -1;
        return;
    }
    
    /*
     * generate plans
     */
    ftbasefactorize(n, tasktype, &n1, &n2, _state);
    if( tasktype==ftbase_ftbasecffttask||tasktype==ftbase_ftbaserffttask )
    {
        
        /*
         * complex FFT plans
         */
        if( n1!=1 )
        {
            
            /*
             * Cooley-Tukey plan (real or complex)
             *
             * Note that child plans are COMPLEX
             * (whether plan itself is complex or not).
             */
            *tmpmemsize = ae_maxint(*tmpmemsize, 2*n1*n2, _state);
            plan->plan.ptr.p_int[entryoffset+0] = esize;
            plan->plan.ptr.p_int[entryoffset+1] = n1;
            plan->plan.ptr.p_int[entryoffset+2] = n2;
            if( tasktype==ftbase_ftbasecffttask )
            {
                plan->plan.ptr.p_int[entryoffset+3] = ftbase_fftcooleytukeyplan;
            }
            else
            {
                plan->plan.ptr.p_int[entryoffset+3] = ftbase_fftrealcooleytukeyplan;
            }
            plan->plan.ptr.p_int[entryoffset+4] = 0;
            plan->plan.ptr.p_int[entryoffset+5] = *plansize;
            ftbase_ftbasegenerateplanrec(n1, ftbase_ftbasecffttask, plan, plansize, precomputedsize, planarraysize, tmpmemsize, stackmemsize, stackptr, _state);
            plan->plan.ptr.p_int[entryoffset+6] = *plansize;
            ftbase_ftbasegenerateplanrec(n2, ftbase_ftbasecffttask, plan, plansize, precomputedsize, planarraysize, tmpmemsize, stackmemsize, stackptr, _state);
            plan->plan.ptr.p_int[entryoffset+7] = -1;
            return;
        }
        else
        {
            if( ((n==2||n==3)||n==4)||n==5 )
            {
                
                /*
                 * hard-coded plan
                 */
                plan->plan.ptr.p_int[entryoffset+0] = esize;
                plan->plan.ptr.p_int[entryoffset+1] = n1;
                plan->plan.ptr.p_int[entryoffset+2] = n2;
                plan->plan.ptr.p_int[entryoffset+3] = ftbase_fftcodeletplan;
                plan->plan.ptr.p_int[entryoffset+4] = 0;
                plan->plan.ptr.p_int[entryoffset+5] = -1;
                plan->plan.ptr.p_int[entryoffset+6] = -1;
                plan->plan.ptr.p_int[entryoffset+7] = *precomputedsize;
                if( n==3 )
                {
                    *precomputedsize = *precomputedsize+2;
                }
                if( n==5 )
                {
                    *precomputedsize = *precomputedsize+5;
                }
                return;
            }
            else
            {
                
                /*
                 * Bluestein's plan
                 *
                 * Select such M that M>=2*N-1, M is composite, and M's
                 * factors are 2, 3, 5
                 */
                k = 2*n2-1;
                m = ftbasefindsmooth(k, _state);
                *tmpmemsize = ae_maxint(*tmpmemsize, 2*m, _state);
                plan->plan.ptr.p_int[entryoffset+0] = esize;
                plan->plan.ptr.p_int[entryoffset+1] = n2;
                plan->plan.ptr.p_int[entryoffset+2] = -1;
                plan->plan.ptr.p_int[entryoffset+3] = ftbase_fftbluesteinplan;
                plan->plan.ptr.p_int[entryoffset+4] = m;
                plan->plan.ptr.p_int[entryoffset+5] = *plansize;
                stackptr = stackptr+2*2*m;
                *stackmemsize = ae_maxint(*stackmemsize, stackptr, _state);
                ftbase_ftbasegenerateplanrec(m, ftbase_ftbasecffttask, plan, plansize, precomputedsize, planarraysize, tmpmemsize, stackmemsize, stackptr, _state);
                stackptr = stackptr-2*2*m;
                plan->plan.ptr.p_int[entryoffset+6] = -1;
                plan->plan.ptr.p_int[entryoffset+7] = *precomputedsize;
                *precomputedsize = *precomputedsize+2*m+2*n;
                return;
            }
        }
    }
    if( tasktype==ftbase_ftbaserfhttask )
    {
        
        /*
         * real FHT plans
         */
        if( n1!=1 )
        {
            
            /*
             * Cooley-Tukey plan
             *
             */
            *tmpmemsize = ae_maxint(*tmpmemsize, 2*n1*n2, _state);
            plan->plan.ptr.p_int[entryoffset+0] = esize;
            plan->plan.ptr.p_int[entryoffset+1] = n1;
            plan->plan.ptr.p_int[entryoffset+2] = n2;
            plan->plan.ptr.p_int[entryoffset+3] = ftbase_fhtcooleytukeyplan;
            plan->plan.ptr.p_int[entryoffset+4] = 0;
            plan->plan.ptr.p_int[entryoffset+5] = *plansize;
            ftbase_ftbasegenerateplanrec(n1, tasktype, plan, plansize, precomputedsize, planarraysize, tmpmemsize, stackmemsize, stackptr, _state);
            plan->plan.ptr.p_int[entryoffset+6] = *plansize;
            ftbase_ftbasegenerateplanrec(n2, tasktype, plan, plansize, precomputedsize, planarraysize, tmpmemsize, stackmemsize, stackptr, _state);
            plan->plan.ptr.p_int[entryoffset+7] = -1;
            return;
        }
        else
        {
            
            /*
             * N2 plan
             */
            plan->plan.ptr.p_int[entryoffset+0] = esize;
            plan->plan.ptr.p_int[entryoffset+1] = n1;
            plan->plan.ptr.p_int[entryoffset+2] = n2;
            plan->plan.ptr.p_int[entryoffset+3] = ftbase_fhtn2plan;
            plan->plan.ptr.p_int[entryoffset+4] = 0;
            plan->plan.ptr.p_int[entryoffset+5] = -1;
            plan->plan.ptr.p_int[entryoffset+6] = -1;
            plan->plan.ptr.p_int[entryoffset+7] = -1;
            if( ((n==2||n==3)||n==4)||n==5 )
            {
                
                /*
                 * hard-coded plan
                 */
                plan->plan.ptr.p_int[entryoffset+0] = esize;
                plan->plan.ptr.p_int[entryoffset+1] = n1;
                plan->plan.ptr.p_int[entryoffset+2] = n2;
                plan->plan.ptr.p_int[entryoffset+3] = ftbase_fhtcodeletplan;
                plan->plan.ptr.p_int[entryoffset+4] = 0;
                plan->plan.ptr.p_int[entryoffset+5] = -1;
                plan->plan.ptr.p_int[entryoffset+6] = -1;
                plan->plan.ptr.p_int[entryoffset+7] = *precomputedsize;
                if( n==3 )
                {
                    *precomputedsize = *precomputedsize+2;
                }
                if( n==5 )
                {
                    *precomputedsize = *precomputedsize+5;
                }
                return;
            }
            return;
        }
    }
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x68, %rsp
movq %r9, %r15
movq %r8, 0x8(%rsp)
movq %rcx, %r13
movq %rdx, %rbx
movq %rsi, 0x10(%rsp)
movq 0xb8(%rsp), %rbp
movq (%rcx), %r14
leaq 0x8(%r14), %rax
movq (%r9), %r12
cmpq %r12, %rax
movq %rdi, (%rsp)
jle 0x2085b
movq %r13, 0x18(%rsp)
leaq 0x20(%rsp), %rsi
movq %rbp, %rdi
callq 0x23815
movq %rbp, %r13
xorl %ebp, %ebp
leaq 0x38(%rsp), %r14
pushq $0x2
popq %rdx
pushq $0x1
popq %r8
movq %r14, %rdi
xorl %esi, %esi
movq %r13, %rcx
callq 0x239a5
movq (%r15), %rsi
movq %r14, %rdi
movq %r13, %r14
movq %r13, %rdx
callq 0x23add
cmpq (%r15), %rbp
jge 0x2080d
movq 0x28(%rbx), %rax
movq (%rax,%rbp,8), %rax
movq 0x60(%rsp), %rcx
movq %rax, (%rcx,%rbp,8)
incq %rbp
jmp 0x207f2
shlq $0x3, %r12
movq %rbx, %rdi
movq %r12, %rsi
movq %r14, %rbp
movq %r14, %rdx
callq 0x23add
xorl %eax, %eax
movq 0x18(%rsp), %r13
cmpq (%r15), %rax
jge 0x20844
movq 0x60(%rsp), %rcx
movq (%rcx,%rax,8), %rcx
movq 0x28(%rbx), %rdx
movq %rcx, (%rdx,%rax,8)
incq %rax
jmp 0x20829
movq %r12, (%r15)
movq %rbp, %rdi
callq 0x237b6
movq (%r13), %r14
leaq 0x8(%r14), %rax
movq (%rsp), %rdi
movq %rax, (%r13)
cmpq $0x1, %rdi
jne 0x208b5
movq 0x28(%rbx), %rax
movq $0x8, (%rax,%r14,8)
movq 0x28(%rbx), %rax
orq $-0x1, 0x8(%rax,%r14,8)
movq 0x28(%rbx), %rax
orq $-0x1, 0x10(%rax,%r14,8)
movq 0x28(%rbx), %rax
movq $0x6, 0x18(%rax,%r14,8)
movq 0x28(%rbx), %rax
orq $-0x1, 0x20(%rax,%r14,8)
movq 0x28(%rbx), %rax
orq $-0x1, 0x28(%rax,%r14,8)
movq 0x28(%rbx), %rax
orq $-0x1, 0x30(%rax,%r14,8)
jmp 0x20c14
leaq 0x38(%rsp), %rdx
leaq 0x20(%rsp), %rcx
callq 0x22ad5
testq $-0x3, 0x10(%rsp)
jne 0x2094c
movq 0x38(%rsp), %r12
cmpq $0x1, %r12
jne 0x20a44
movq (%rsp), %rsi
leaq -0x2(%rsi), %rax
cmpq $0x3, %rax
ja 0x20c2d
movq 0x28(%rbx), %rax
movq $0x8, (%rax,%r14,8)
movq 0x28(%rbx), %rax
movq $0x1, 0x8(%rax,%r14,8)
movq 0x20(%rsp), %rax
movq 0x28(%rbx), %rcx
movq %rax, 0x10(%rcx,%r14,8)
movq 0x28(%rbx), %rax
movq $0x2, 0x18(%rax,%r14,8)
movq 0x28(%rbx), %rax
andq $0x0, 0x20(%rax,%r14,8)
movq 0x28(%rbx), %rax
orq $-0x1, 0x28(%rax,%r14,8)
movq 0x28(%rbx), %rax
orq $-0x1, 0x30(%rax,%r14,8)
movq 0x8(%rsp), %rdx
jmp 0x20a1b
movq 0x38(%rsp), %r12
cmpq $0x1, %r12
jne 0x20b26
movq 0x28(%rbx), %rax
movq $0x8, (%rax,%r14,8)
movq 0x28(%rbx), %rax
movq $0x1, 0x8(%rax,%r14,8)
movq 0x20(%rsp), %rax
movq 0x28(%rbx), %rcx
movq %rax, 0x10(%rcx,%r14,8)
movq 0x28(%rbx), %rcx
movq $0x3e7, 0x18(%rcx,%r14,8) # imm = 0x3E7
movq 0x28(%rbx), %rcx
andq $0x0, 0x20(%rcx,%r14,8)
movq 0x28(%rbx), %rcx
orq $-0x1, 0x28(%rcx,%r14,8)
movq 0x28(%rbx), %rcx
orq $-0x1, 0x30(%rcx,%r14,8)
movq 0x28(%rbx), %rcx
orq $-0x1, 0x38(%rcx,%r14,8)
movq (%rsp), %rsi
leaq -0x2(%rsi), %rcx
cmpq $0x3, %rcx
movq 0x8(%rsp), %rdx
ja 0x20c1e
movq 0x28(%rbx), %rcx
movq $0x8, (%rcx,%r14,8)
movq 0x28(%rbx), %rcx
movq $0x1, 0x8(%rcx,%r14,8)
movq 0x28(%rbx), %rcx
movq %rax, 0x10(%rcx,%r14,8)
movq 0x28(%rbx), %rax
movq $0x4, 0x18(%rax,%r14,8)
movq 0x28(%rbx), %rax
andq $0x0, 0x20(%rax,%r14,8)
movq 0x28(%rbx), %rax
orq $-0x1, 0x28(%rax,%r14,8)
movq 0x28(%rbx), %rax
orq $-0x1, 0x30(%rax,%r14,8)
movq (%rdx), %rax
movq 0x28(%rbx), %rcx
movq %rax, 0x38(%rcx,%r14,8)
cmpq $0x5, %rsi
je 0x20d3a
cmpq $0x3, %rsi
jne 0x20c1e
addq $0x2, %rax
jmp 0x20d3e
movq 0xa0(%rsp), %rax
movq (%rax), %rdi
movq %rbp, %rdx
movq 0x20(%rsp), %rax
movq %r12, %rsi
imulq %rax, %rsi
movq %rax, %rbp
movq %rax, (%rsp)
addq %rsi, %rsi
callq 0x244a0
movq 0xa0(%rsp), %r10
movq %rax, (%r10)
movq 0x28(%rbx), %rax
movq $0x8, (%rax,%r14,8)
movq 0x28(%rbx), %rax
movq %r12, 0x8(%rax,%r14,8)
movq 0x28(%rbx), %rax
movq %rbp, 0x10(%rax,%r14,8)
xorl %eax, %eax
cmpq $0x0, 0x10(%rsp)
setne %al
movq 0x28(%rbx), %rcx
leaq (%rax,%rax,4), %rax
movq %rax, 0x18(%rcx,%r14,8)
movq 0x28(%rbx), %rax
andq $0x0, 0x20(%rax,%r14,8)
movq (%r13), %rax
movq 0x28(%rbx), %rcx
movq %rax, 0x28(%rcx,%r14,8)
movq %r12, %rdi
xorl %esi, %esi
movq %rbx, %rdx
movq %r13, %rcx
movq 0x8(%rsp), %r12
movq %r12, %r8
movq %r15, %r9
movq 0xb8(%rsp), %rbp
pushq %rbp
movq 0xb8(%rsp), %rax
pushq %rax
movq 0xb8(%rsp), %rax
pushq %rax
pushq %r10
callq 0x20771
addq $0x20, %rsp
movq (%r13), %rax
movq 0x28(%rbx), %rcx
movq %rax, 0x30(%rcx,%r14,8)
movq (%rsp), %rdi
xorl %esi, %esi
movq %rbx, %rdx
movq %r13, %rcx
movq %r12, %r8
movq %r15, %r9
pushq %rbp
jmp 0x20bf6
movq 0xa0(%rsp), %rax
movq (%rax), %rdi
movq %rbp, %rdx
movq 0x20(%rsp), %rbp
movq %r12, %rsi
imulq %rbp, %rsi
addq %rsi, %rsi
callq 0x244a0
movq 0xa0(%rsp), %r10
movq %rax, (%r10)
movq 0x28(%rbx), %rax
movq $0x8, (%rax,%r14,8)
movq 0x28(%rbx), %rax
movq %r12, 0x8(%rax,%r14,8)
movq 0x28(%rbx), %rax
movq %rbp, 0x10(%rax,%r14,8)
movq 0x28(%rbx), %rax
movq $0x3, 0x18(%rax,%r14,8)
movq 0x28(%rbx), %rax
andq $0x0, 0x20(%rax,%r14,8)
movq (%r13), %rax
movq 0x28(%rbx), %rcx
movq %rax, 0x28(%rcx,%r14,8)
pushq $0x1
popq %rsi
movq %r12, %rdi
movq %rbx, %rdx
movq %r13, %rcx
movq 0x8(%rsp), %r12
movq %r12, %r8
movq %r15, %r9
pushq 0xb8(%rsp)
movq 0xb8(%rsp), %rax
pushq %rax
movq 0xb8(%rsp), %rax
pushq %rax
pushq %r10
callq 0x20771
addq $0x20, %rsp
movq (%r13), %rax
movq 0x28(%rbx), %rcx
movq %rax, 0x30(%rcx,%r14,8)
movq %rbp, %rdi
pushq $0x1
popq %rsi
movq %rbx, %rdx
movq %r13, %rcx
movq %r12, %r8
movq %r15, %r9
pushq 0xb8(%rsp)
pushq 0xb8(%rsp)
pushq 0xb8(%rsp)
pushq 0xb8(%rsp)
callq 0x20771
addq $0x20, %rsp
movq 0x28(%rbx), %rax
orq $-0x1, 0x38(%rax,%r14,8)
addq $0x68, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movq %r13, 0x18(%rsp)
movq 0x20(%rsp), %r13
leaq -0x1(,%r13,2), %rdi
movq %rbp, %rsi
callq 0x22b9e
movq %rax, %r12
movq 0xa0(%rsp), %rax
movq (%rax), %rdi
leaq (%r12,%r12), %rsi
movq %rbp, %rdx
callq 0x244a0
movq 0xa0(%rsp), %rcx
movq %rax, (%rcx)
movq 0x28(%rbx), %rax
movq $0x8, (%rax,%r14,8)
movq 0x28(%rbx), %rax
movq %r13, 0x8(%rax,%r14,8)
movq 0x28(%rbx), %rax
orq $-0x1, 0x10(%rax,%r14,8)
movq 0x28(%rbx), %rax
movq $0x1, 0x18(%rax,%r14,8)
movq 0x28(%rbx), %rax
movq %r12, 0x20(%rax,%r14,8)
movq 0x18(%rsp), %r13
movq (%r13), %rax
movq 0x28(%rbx), %rcx
movq %rax, 0x28(%rcx,%r14,8)
movq 0xb0(%rsp), %rax
leaq (%rax,%r12,4), %rsi
movq %rsi, 0x10(%rsp)
movq 0xa8(%rsp), %rax
movq (%rax), %rdi
movq %rbp, %rdx
callq 0x244a0
movq 0xa8(%rsp), %r10
movq %rax, (%r10)
movq %r12, %rdi
xorl %esi, %esi
movq %rbx, %rdx
movq %r13, %rcx
movq 0x8(%rsp), %r13
movq %r13, %r8
movq %r15, %r9
pushq %rbp
pushq 0x18(%rsp)
pushq %r10
pushq 0xb8(%rsp)
callq 0x20771
addq $0x20, %rsp
movq 0x28(%rbx), %rax
orq $-0x1, 0x30(%rax,%r14,8)
movq (%r13), %rax
movq 0x28(%rbx), %rcx
movq %rax, 0x38(%rcx,%r14,8)
movq (%rsp), %rcx
addq %r12, %rcx
leaq (%rax,%rcx,2), %rax
movq %rax, (%r13)
jmp 0x20c1e
addq $0x5, %rax
movq %rax, (%rdx)
jmp 0x20c1e