// File: mouda[P]LTE_random_access-simulate/ap.cpp
// Function: alglib_impl::_ialglib_cmv(long, long, double const*, double const*, alglib_impl::ae_complex*, double*, long, alglib_impl::ae_complex, alglib_impl::ae_complex)
void _ialglib_cmv(ae_int_t m, ae_int_t n, const double *a, const double *x, ae_complex *cy, double *dy, ae_int_t stride, ae_complex alpha, ae_complex beta)
{
    ae_int_t i, j;
    const double *pa, *parow, *pb;

    parow = a;
    for(i=0; i<m; i++)
    {
        double v0 = 0, v1 = 0;
        pa = parow;
        pb = x;
        for(j=0; j<n; j++)
        {
            v0 += pa[0]*pb[0];
            v1 += pa[0]*pb[1];
            v0 -= pa[1]*pb[1];
            v1 += pa[1]*pb[0];

            pa  += 2;
            pb  += 2;
        }
        if( cy!=NULL )
        {
            double tx = (beta.x*cy->x-beta.y*cy->y)+(alpha.x*v0-alpha.y*v1);
            double ty = (beta.x*cy->y+beta.y*cy->x)+(alpha.x*v1+alpha.y*v0);
            cy->x = tx;
            cy->y = ty;
            cy+=stride;
        }
        else
        {
            double tx = (beta.x*dy[0]-beta.y*dy[1])+(alpha.x*v0-alpha.y*v1);
            double ty = (beta.x*dy[1]+beta.y*dy[0])+(alpha.x*v1+alpha.y*v0);
            dy[0] = tx;
            dy[1] = ty;
            dy += 2*stride;
        }
        parow += 2*alglib_c_block;
    }
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x28, %rsp
movq 0x10(%rbp), %rax
movsd %xmm0, -0x10(%rbp)
movsd %xmm1, -0x8(%rbp)
movsd %xmm2, -0x20(%rbp)
movsd %xmm3, -0x18(%rbp)
movq %rdi, -0x28(%rbp)
movq %rsi, -0x30(%rbp)
movq %rdx, -0x38(%rbp)
movq %rcx, -0x40(%rbp)
movq %r8, -0x48(%rbp)
movq %r9, -0x50(%rbp)
movq -0x38(%rbp), %rax
movq %rax, -0x70(%rbp)
movq $0x0, -0x58(%rbp)
movq -0x58(%rbp), %rax
cmpq -0x28(%rbp), %rax
jge 0x4db0f
xorps %xmm0, %xmm0
movsd %xmm0, -0x80(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x88(%rbp)
movq -0x70(%rbp), %rax
movq %rax, -0x68(%rbp)
movq -0x40(%rbp), %rax
movq %rax, -0x78(%rbp)
movq $0x0, -0x60(%rbp)
movq -0x60(%rbp), %rax
cmpq -0x30(%rbp), %rax
jge 0x4d8e7
movq -0x68(%rbp), %rax
movsd (%rax), %xmm0
movq -0x78(%rbp), %rax
movsd (%rax), %xmm2
movsd -0x80(%rbp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x80(%rbp)
movq -0x68(%rbp), %rax
movsd (%rax), %xmm0
movq -0x78(%rbp), %rax
movsd 0x8(%rax), %xmm2
movsd -0x88(%rbp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x88(%rbp)
movq -0x68(%rbp), %rax
movsd 0x8(%rax), %xmm0
movq -0x78(%rbp), %rax
movsd 0x8(%rax), %xmm2
movsd -0x80(%rbp), %xmm1
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x80(%rbp)
movq -0x68(%rbp), %rax
movsd 0x8(%rax), %xmm0
movq -0x78(%rbp), %rax
movsd (%rax), %xmm2
movsd -0x88(%rbp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x88(%rbp)
movq -0x68(%rbp), %rax
addq $0x10, %rax
movq %rax, -0x68(%rbp)
movq -0x78(%rbp), %rax
addq $0x10, %rax
movq %rax, -0x78(%rbp)
movq -0x60(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x60(%rbp)
jmp 0x4d801
cmpq $0x0, -0x48(%rbp)
je 0x4d9f2
movsd -0x20(%rbp), %xmm0
movq -0x48(%rbp), %rax
movsd (%rax), %xmm2
movsd -0x18(%rbp), %xmm1
movq -0x48(%rbp), %rax
mulsd 0x8(%rax), %xmm1
movq %xmm1, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd -0x10(%rbp), %xmm1
movsd -0x80(%rbp), %xmm3
movsd -0x8(%rbp), %xmm2
mulsd -0x88(%rbp), %xmm2
movq %xmm2, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
addsd %xmm1, %xmm0
movsd %xmm0, -0x90(%rbp)
movsd -0x20(%rbp), %xmm0
movq -0x48(%rbp), %rax
movsd 0x8(%rax), %xmm2
movsd -0x18(%rbp), %xmm1
movq -0x48(%rbp), %rax
mulsd (%rax), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd -0x10(%rbp), %xmm1
movsd -0x88(%rbp), %xmm3
movsd -0x8(%rbp), %xmm2
mulsd -0x80(%rbp), %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
addsd %xmm1, %xmm0
movsd %xmm0, -0x98(%rbp)
movsd -0x90(%rbp), %xmm0
movq -0x48(%rbp), %rax
movsd %xmm0, (%rax)
movsd -0x98(%rbp), %xmm0
movq -0x48(%rbp), %rax
movsd %xmm0, 0x8(%rax)
movq 0x10(%rbp), %rax
shlq $0x4, %rax
addq -0x48(%rbp), %rax
movq %rax, -0x48(%rbp)
jmp 0x4daf0
movsd -0x20(%rbp), %xmm0
movq -0x50(%rbp), %rax
movsd (%rax), %xmm2
movsd -0x18(%rbp), %xmm1
movq -0x50(%rbp), %rax
mulsd 0x8(%rax), %xmm1
movq %xmm1, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd -0x10(%rbp), %xmm1
movsd -0x80(%rbp), %xmm3
movsd -0x8(%rbp), %xmm2
mulsd -0x88(%rbp), %xmm2
movq %xmm2, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
addsd %xmm1, %xmm0
movsd %xmm0, -0xa0(%rbp)
movsd -0x20(%rbp), %xmm0
movq -0x50(%rbp), %rax
movsd 0x8(%rax), %xmm2
movsd -0x18(%rbp), %xmm1
movq -0x50(%rbp), %rax
mulsd (%rax), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd -0x10(%rbp), %xmm1
movsd -0x88(%rbp), %xmm3
movsd -0x8(%rbp), %xmm2
mulsd -0x80(%rbp), %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
addsd %xmm1, %xmm0
movsd %xmm0, -0xa8(%rbp)
movsd -0xa0(%rbp), %xmm0
movq -0x50(%rbp), %rax
movsd %xmm0, (%rax)
movsd -0xa8(%rbp), %xmm0
movq -0x50(%rbp), %rax
movsd %xmm0, 0x8(%rax)
movq 0x10(%rbp), %rax
shlq %rax
shlq $0x3, %rax
addq -0x50(%rbp), %rax
movq %rax, -0x50(%rbp)
movq -0x70(%rbp), %rax
addq $0x180, %rax            # imm = 0x180
movq %rax, -0x70(%rbp)
movq -0x58(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x58(%rbp)
jmp 0x4d7c8
addq $0x28, %rsp
popq %rbp
retq
nopw %cs:(%rax,%rax)