// File: mouda[P]LTE_random_access-simulate/alglibinternal.cpp
// Function: alglib_impl::matrixvectormultiply(alglib_impl::ae_matrix*, long, long, long, long, bool, alglib_impl::ae_vector*, long, long, double, alglib_impl::ae_vector*, long, long, double, alglib_impl::ae_state*)
void matrixvectormultiply(/* Real    */ ae_matrix* a,
     ae_int_t i1,
     ae_int_t i2,
     ae_int_t j1,
     ae_int_t j2,
     ae_bool trans,
     /* Real    */ ae_vector* x,
     ae_int_t ix1,
     ae_int_t ix2,
     double alpha,
     /* Real    */ ae_vector* y,
     ae_int_t iy1,
     ae_int_t iy2,
     double beta,
     ae_state *_state)
{
    ae_int_t i;
    double v;


    if( !trans )
    {
        
        /*
         * y := alpha*A*x + beta*y;
         */
        if( i1>i2||j1>j2 )
        {
            return;
        }
        ae_assert(j2-j1==ix2-ix1, "MatrixVectorMultiply: A and X dont match!", _state);
        ae_assert(i2-i1==iy2-iy1, "MatrixVectorMultiply: A and Y dont match!", _state);
        
        /*
         * beta*y
         */
        if( ae_fp_eq(beta,0) )
        {
            for(i=iy1; i<=iy2; i++)
            {
                y->ptr.p_double[i] = 0;
            }
        }
        else
        {
            ae_v_muld(&y->ptr.p_double[iy1], 1, ae_v_len(iy1,iy2), beta);
        }
        
        /*
         * alpha*A*x
         */
        for(i=i1; i<=i2; i++)
        {
            v = ae_v_dotproduct(&a->ptr.pp_double[i][j1], 1, &x->ptr.p_double[ix1], 1, ae_v_len(j1,j2));
            y->ptr.p_double[iy1+i-i1] = y->ptr.p_double[iy1+i-i1]+alpha*v;
        }
    }
    else
    {
        
        /*
         * y := alpha*A'*x + beta*y;
         */
        if( i1>i2||j1>j2 )
        {
            return;
        }
        ae_assert(i2-i1==ix2-ix1, "MatrixVectorMultiply: A and X dont match!", _state);
        ae_assert(j2-j1==iy2-iy1, "MatrixVectorMultiply: A and Y dont match!", _state);
        
        /*
         * beta*y
         */
        if( ae_fp_eq(beta,0) )
        {
            for(i=iy1; i<=iy2; i++)
            {
                y->ptr.p_double[i] = 0;
            }
        }
        else
        {
            ae_v_muld(&y->ptr.p_double[iy1], 1, ae_v_len(iy1,iy2), beta);
        }
        
        /*
         * alpha*A'*x
         */
        for(i=i1; i<=i2; i++)
        {
            v = alpha*x->ptr.p_double[ix1+i-i1];
            ae_v_addd(&y->ptr.p_double[iy1], 1, &a->ptr.pp_double[i][j1], 1, ae_v_len(iy1,iy2), v);
        }
    }
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x80, %rsp
movb %r9b, %al
movq 0x40(%rbp), %r9
movq 0x38(%rbp), %r9
movq 0x30(%rbp), %r9
movq 0x28(%rbp), %r9
movq 0x20(%rbp), %r9
movq 0x18(%rbp), %r9
movq 0x10(%rbp), %r9
movq %rdi, -0x8(%rbp)
movq %rsi, -0x10(%rbp)
movq %rdx, -0x18(%rbp)
movq %rcx, -0x20(%rbp)
movq %r8, -0x28(%rbp)
andb $0x1, %al
movb %al, -0x29(%rbp)
movsd %xmm0, -0x38(%rbp)
movsd %xmm1, -0x40(%rbp)
testb $0x1, -0x29(%rbp)
jne 0x23b4b
movq -0x10(%rbp), %rax
cmpq -0x18(%rbp), %rax
jg 0x2399b
movq -0x20(%rbp), %rax
cmpq -0x28(%rbp), %rax
jle 0x239a0
jmp 0x23ceb
movq -0x28(%rbp), %rax
subq -0x20(%rbp), %rax
movq 0x20(%rbp), %rcx
subq 0x18(%rbp), %rcx
cmpq %rcx, %rax
sete %al
movq 0x40(%rbp), %rdx
movzbl %al, %edi
andl $0x1, %edi
leaq 0x342e9(%rip), %rsi     # 0x57cb0
callq 0x3b240
movq -0x18(%rbp), %rax
subq -0x10(%rbp), %rax
movq 0x38(%rbp), %rcx
subq 0x30(%rbp), %rcx
cmpq %rcx, %rax
sete %al
movq 0x40(%rbp), %rdx
movzbl %al, %edi
andl $0x1, %edi
leaq 0x342e7(%rip), %rsi     # 0x57cda
callq 0x3b240
movsd -0x40(%rbp), %xmm0
xorps %xmm1, %xmm1
callq 0x3c3d0
testb $0x1, %al
jne 0x23a0b
jmp 0x23a41
movq 0x30(%rbp), %rax
movq %rax, -0x48(%rbp)
movq -0x48(%rbp), %rax
cmpq 0x38(%rbp), %rax
jg 0x23a3f
movq 0x28(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x48(%rbp), %rcx
xorps %xmm0, %xmm0
movsd %xmm0, (%rax,%rcx,8)
movq -0x48(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x48(%rbp)
jmp 0x23a13
jmp 0x23a7b
movq 0x28(%rbp), %rax
movq 0x28(%rax), %rax
movq 0x30(%rbp), %rcx
shlq $0x3, %rcx
addq %rcx, %rax
movq %rax, -0x58(%rbp)
movq 0x30(%rbp), %rdi
movq 0x38(%rbp), %rsi
callq 0x42370
movq -0x58(%rbp), %rdi
movq %rax, %rdx
movsd -0x40(%rbp), %xmm0
movl $0x1, %esi
callq 0x422d0
movq -0x10(%rbp), %rax
movq %rax, -0x48(%rbp)
movq -0x48(%rbp), %rax
cmpq -0x18(%rbp), %rax
jg 0x23b46
movq -0x8(%rbp), %rax
movq 0x38(%rax), %rax
movq -0x48(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movq -0x20(%rbp), %rcx
shlq $0x3, %rcx
addq %rcx, %rax
movq %rax, -0x68(%rbp)
movq 0x10(%rbp), %rax
movq 0x28(%rax), %rax
movq 0x18(%rbp), %rcx
shlq $0x3, %rcx
addq %rcx, %rax
movq %rax, -0x60(%rbp)
movq -0x20(%rbp), %rdi
movq -0x28(%rbp), %rsi
callq 0x42370
movq -0x68(%rbp), %rdi
movq -0x60(%rbp), %rdx
movq %rax, %r8
movl $0x1, %ecx
movq %rcx, %rsi
callq 0x41940
movsd %xmm0, -0x50(%rbp)
movq 0x28(%rbp), %rax
movq 0x28(%rax), %rax
movq 0x30(%rbp), %rcx
addq -0x48(%rbp), %rcx
subq -0x10(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm1
movsd -0x38(%rbp), %xmm0
movsd -0x50(%rbp), %xmm2
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0x28(%rbp), %rax
movq 0x28(%rax), %rax
movq 0x30(%rbp), %rcx
addq -0x48(%rbp), %rcx
subq -0x10(%rbp), %rcx
movsd %xmm0, (%rax,%rcx,8)
movq -0x48(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x48(%rbp)
jmp 0x23a83
jmp 0x23ceb
movq -0x10(%rbp), %rax
cmpq -0x18(%rbp), %rax
jg 0x23b5f
movq -0x20(%rbp), %rax
cmpq -0x28(%rbp), %rax
jle 0x23b64
jmp 0x23ceb
movq -0x18(%rbp), %rax
subq -0x10(%rbp), %rax
movq 0x20(%rbp), %rcx
subq 0x18(%rbp), %rcx
cmpq %rcx, %rax
sete %al
movq 0x40(%rbp), %rdx
movzbl %al, %edi
andl $0x1, %edi
leaq 0x34125(%rip), %rsi     # 0x57cb0
callq 0x3b240
movq -0x28(%rbp), %rax
subq -0x20(%rbp), %rax
movq 0x38(%rbp), %rcx
subq 0x30(%rbp), %rcx
cmpq %rcx, %rax
sete %al
movq 0x40(%rbp), %rdx
movzbl %al, %edi
andl $0x1, %edi
leaq 0x34123(%rip), %rsi     # 0x57cda
callq 0x3b240
movsd -0x40(%rbp), %xmm0
xorps %xmm1, %xmm1
callq 0x3c3d0
testb $0x1, %al
jne 0x23bcf
jmp 0x23c05
movq 0x30(%rbp), %rax
movq %rax, -0x48(%rbp)
movq -0x48(%rbp), %rax
cmpq 0x38(%rbp), %rax
jg 0x23c03
movq 0x28(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x48(%rbp), %rcx
xorps %xmm0, %xmm0
movsd %xmm0, (%rax,%rcx,8)
movq -0x48(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x48(%rbp)
jmp 0x23bd7
jmp 0x23c3f
movq 0x28(%rbp), %rax
movq 0x28(%rax), %rax
movq 0x30(%rbp), %rcx
shlq $0x3, %rcx
addq %rcx, %rax
movq %rax, -0x70(%rbp)
movq 0x30(%rbp), %rdi
movq 0x38(%rbp), %rsi
callq 0x42370
movq -0x70(%rbp), %rdi
movq %rax, %rdx
movsd -0x40(%rbp), %xmm0
movl $0x1, %esi
callq 0x422d0
movq -0x10(%rbp), %rax
movq %rax, -0x48(%rbp)
movq -0x48(%rbp), %rax
cmpq -0x18(%rbp), %rax
jg 0x23ce9
movsd -0x38(%rbp), %xmm0
movq 0x10(%rbp), %rax
movq 0x28(%rax), %rax
movq 0x18(%rbp), %rcx
addq -0x48(%rbp), %rcx
subq -0x10(%rbp), %rcx
mulsd (%rax,%rcx,8), %xmm0
movsd %xmm0, -0x50(%rbp)
movq 0x28(%rbp), %rax
movq 0x28(%rax), %rax
movq 0x30(%rbp), %rcx
shlq $0x3, %rcx
addq %rcx, %rax
movq %rax, -0x80(%rbp)
movq -0x8(%rbp), %rax
movq 0x38(%rax), %rax
movq -0x48(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movq -0x20(%rbp), %rcx
shlq $0x3, %rcx
addq %rcx, %rax
movq %rax, -0x78(%rbp)
movq 0x30(%rbp), %rdi
movq 0x38(%rbp), %rsi
callq 0x42370
movq -0x80(%rbp), %rdi
movq -0x78(%rbp), %rdx
movq %rax, %r8
movsd -0x50(%rbp), %xmm0
movl $0x1, %ecx
movq %rcx, %rsi
callq 0x41fe0
movq -0x48(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x48(%rbp)
jmp 0x23c47
jmp 0x23ceb
addq $0x80, %rsp
popq %rbp
retq
nopw %cs:(%rax,%rax)