// File: mouda[P]LTE_random_access-simulate/alglibinternal.cpp
// Function: alglib_impl::ftbase_ffttwcalc(alglib_impl::ae_vector*, long, long, long, alglib_impl::ae_state*)
static void ftbase_ffttwcalc(/* Real    */ ae_vector* a,
     ae_int_t aoffset,
     ae_int_t n1,
     ae_int_t n2,
     ae_state *_state)
{
    ae_int_t i;
    ae_int_t j;
    ae_int_t n;
    ae_int_t idx;
    ae_int_t offs;
    double x;
    double y;
    double twxm1;
    double twy;
    double twbasexm1;
    double twbasey;
    double twrowxm1;
    double twrowy;
    double tmpx;
    double tmpy;
    double v;


    n = n1*n2;
    v = -2*ae_pi/n;
    twbasexm1 = -2*ae_sqr(ae_sin(0.5*v, _state), _state);
    twbasey = ae_sin(v, _state);
    twrowxm1 = 0;
    twrowy = 0;
    for(i=0; i<=n2-1; i++)
    {
        twxm1 = 0;
        twy = 0;
        for(j=0; j<=n1-1; j++)
        {
            idx = i*n1+j;
            offs = aoffset+2*idx;
            x = a->ptr.p_double[offs+0];
            y = a->ptr.p_double[offs+1];
            tmpx = x*twxm1-y*twy;
            tmpy = x*twy+y*twxm1;
            a->ptr.p_double[offs+0] = x+tmpx;
            a->ptr.p_double[offs+1] = y+tmpy;
            
            /*
             * update Tw: Tw(new) = Tw(old)*TwRow
             */
            if( j<n1-1 )
            {
                if( j%ftbase_ftbaseupdatetw==0 )
                {
                    v = -2*ae_pi*i*(j+1)/n;
                    twxm1 = -2*ae_sqr(ae_sin(0.5*v, _state), _state);
                    twy = ae_sin(v, _state);
                }
                else
                {
                    tmpx = twrowxm1+twxm1*twrowxm1-twy*twrowy;
                    tmpy = twrowy+twxm1*twrowy+twy*twrowxm1;
                    twxm1 = twxm1+tmpx;
                    twy = twy+tmpy;
                }
            }
        }
        
        /*
         * update TwRow: TwRow(new) = TwRow(old)*TwBase
         */
        if( i<n2-1 )
        {
            if( j%ftbase_ftbaseupdatetw==0 )
            {
                v = -2*ae_pi*(i+1)/n;
                twrowxm1 = -2*ae_sqr(ae_sin(0.5*v, _state), _state);
                twrowy = ae_sin(v, _state);
            }
            else
            {
                tmpx = twbasexm1+twrowxm1*twbasexm1-twrowy*twbasey;
                tmpy = twbasey+twrowxm1*twbasey+twrowy*twbasexm1;
                twrowxm1 = twrowxm1+tmpx;
                twrowy = twrowy+tmpy;
            }
        }
    }
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0xb0, %rsp
movq %rdi, -0x8(%rbp)
movq %rsi, -0x10(%rbp)
movq %rdx, -0x18(%rbp)
movq %rcx, -0x20(%rbp)
movq %r8, -0x28(%rbp)
movq -0x18(%rbp), %rax
imulq -0x20(%rbp), %rax
movq %rax, -0x40(%rbp)
cvtsi2sdq -0x40(%rbp), %xmm1
movsd 0x1f56e(%rip), %xmm0    # 0x57898
divsd %xmm1, %xmm0
movsd %xmm0, -0xa8(%rbp)
movsd 0x1d112(%rip), %xmm0    # 0x55450
mulsd -0xa8(%rbp), %xmm0
movq -0x28(%rbp), %rdi
callq 0x3cb90
movq -0x28(%rbp), %rdi
callq 0x3c870
movaps %xmm0, %xmm1
movsd 0x1d2c5(%rip), %xmm0    # 0x55628
mulsd %xmm1, %xmm0
movsd %xmm0, -0x78(%rbp)
movsd -0xa8(%rbp), %xmm0
movq -0x28(%rbp), %rdi
callq 0x3cb90
movsd %xmm0, -0x80(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x88(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x90(%rbp)
movq $0x0, -0x30(%rbp)
movq -0x30(%rbp), %rax
movq -0x20(%rbp), %rcx
subq $0x1, %rcx
cmpq %rcx, %rax
jg 0x387bc
xorps %xmm0, %xmm0
movsd %xmm0, -0x68(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x70(%rbp)
movq $0x0, -0x38(%rbp)
movq -0x38(%rbp), %rax
movq -0x18(%rbp), %rcx
subq $0x1, %rcx
cmpq %rcx, %rax
jg 0x38644
movq -0x30(%rbp), %rax
imulq -0x18(%rbp), %rax
addq -0x38(%rbp), %rax
movq %rax, -0x48(%rbp)
movq -0x10(%rbp), %rax
movq -0x48(%rbp), %rcx
shlq %rcx
addq %rcx, %rax
movq %rax, -0x50(%rbp)
movq -0x8(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x50(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
movsd %xmm0, -0x58(%rbp)
movq -0x8(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x50(%rbp), %rcx
movsd 0x8(%rax,%rcx,8), %xmm0
movsd %xmm0, -0x60(%rbp)
movsd -0x58(%rbp), %xmm0
movsd -0x68(%rbp), %xmm2
movsd -0x60(%rbp), %xmm1
mulsd -0x70(%rbp), %xmm1
movq %xmm1, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x98(%rbp)
movsd -0x58(%rbp), %xmm0
movsd -0x70(%rbp), %xmm2
movsd -0x60(%rbp), %xmm1
mulsd -0x68(%rbp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0xa0(%rbp)
movsd -0x58(%rbp), %xmm0
addsd -0x98(%rbp), %xmm0
movq -0x8(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x50(%rbp), %rcx
movsd %xmm0, (%rax,%rcx,8)
movsd -0x60(%rbp), %xmm0
addsd -0xa0(%rbp), %xmm0
movq -0x8(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x50(%rbp), %rcx
movsd %xmm0, 0x8(%rax,%rcx,8)
movq -0x38(%rbp), %rax
movq -0x18(%rbp), %rcx
subq $0x1, %rcx
cmpq %rcx, %rax
jge 0x38631
movq -0x38(%rbp), %rax
cqto
idivq 0x2cea0(%rip)           # 0x65390
cmpq $0x0, %rdx
jne 0x38580
cvtsi2sdq -0x30(%rbp), %xmm1
movsd 0x1f390(%rip), %xmm0    # 0x57898
mulsd %xmm1, %xmm0
movq -0x38(%rbp), %rax
addq $0x1, %rax
cvtsi2sd %rax, %xmm1
mulsd %xmm1, %xmm0
cvtsi2sdq -0x40(%rbp), %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, -0xa8(%rbp)
movsd 0x1cf19(%rip), %xmm0    # 0x55450
mulsd -0xa8(%rbp), %xmm0
movq -0x28(%rbp), %rdi
callq 0x3cb90
movq -0x28(%rbp), %rdi
callq 0x3c870
movaps %xmm0, %xmm1
movsd 0x1d0cc(%rip), %xmm0    # 0x55628
mulsd %xmm1, %xmm0
movsd %xmm0, -0x68(%rbp)
movsd -0xa8(%rbp), %xmm0
movq -0x28(%rbp), %rdi
callq 0x3cb90
movsd %xmm0, -0x70(%rbp)
jmp 0x3862f
movsd -0x88(%rbp), %xmm0
movsd -0x68(%rbp), %xmm1
movsd -0x88(%rbp), %xmm2
mulsd %xmm2, %xmm1
addsd %xmm0, %xmm1
movsd -0x70(%rbp), %xmm0
movsd -0x90(%rbp), %xmm2
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x98(%rbp)
movsd -0x90(%rbp), %xmm0
movsd -0x68(%rbp), %xmm1
movsd -0x90(%rbp), %xmm2
mulsd %xmm2, %xmm1
addsd %xmm0, %xmm1
movsd -0x70(%rbp), %xmm0
movsd -0x88(%rbp), %xmm2
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0xa0(%rbp)
movsd -0x68(%rbp), %xmm0
addsd -0x98(%rbp), %xmm0
movsd %xmm0, -0x68(%rbp)
movsd -0x70(%rbp), %xmm0
addsd -0xa0(%rbp), %xmm0
movsd %xmm0, -0x70(%rbp)
jmp 0x38631
jmp 0x38633
movq -0x38(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x38(%rbp)
jmp 0x383cd
movq -0x30(%rbp), %rax
movq -0x20(%rbp), %rcx
subq $0x1, %rcx
cmpq %rcx, %rax
jge 0x387a9
movq -0x38(%rbp), %rax
cqto
idivq 0x2cd2a(%rip)           # 0x65390
cmpq $0x0, %rdx
jne 0x386f2
movq -0x30(%rbp), %rax
addq $0x1, %rax
cvtsi2sd %rax, %xmm1
movsd 0x1f213(%rip), %xmm0    # 0x57898
mulsd %xmm1, %xmm0
cvtsi2sdq -0x40(%rbp), %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, -0xa8(%rbp)
movsd 0x1cdad(%rip), %xmm0    # 0x55450
mulsd -0xa8(%rbp), %xmm0
movq -0x28(%rbp), %rdi
callq 0x3cb90
movq -0x28(%rbp), %rdi
callq 0x3c870
movaps %xmm0, %xmm1
movsd 0x1cf60(%rip), %xmm0    # 0x55628
mulsd %xmm1, %xmm0
movsd %xmm0, -0x88(%rbp)
movsd -0xa8(%rbp), %xmm0
movq -0x28(%rbp), %rdi
callq 0x3cb90
movsd %xmm0, -0x90(%rbp)
jmp 0x387a7
movsd -0x78(%rbp), %xmm0
movsd -0x88(%rbp), %xmm1
movsd -0x78(%rbp), %xmm2
mulsd %xmm2, %xmm1
addsd %xmm0, %xmm1
movsd -0x90(%rbp), %xmm0
movsd -0x80(%rbp), %xmm2
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x98(%rbp)
movsd -0x80(%rbp), %xmm0
movsd -0x88(%rbp), %xmm1
movsd -0x80(%rbp), %xmm2
mulsd %xmm2, %xmm1
addsd %xmm0, %xmm1
movsd -0x90(%rbp), %xmm0
movsd -0x78(%rbp), %xmm2
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0xa0(%rbp)
movsd -0x88(%rbp), %xmm0
addsd -0x98(%rbp), %xmm0
movsd %xmm0, -0x88(%rbp)
movsd -0x90(%rbp), %xmm0
addsd -0xa0(%rbp), %xmm0
movsd %xmm0, -0x90(%rbp)
jmp 0x387a9
jmp 0x387ab
movq -0x30(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x30(%rbp)
jmp 0x383a0
addq $0xb0, %rsp
popq %rbp
retq
nopw %cs:(%rax,%rax)