// File: mouda[P]LTE_random_access-simulate/alglibinternal.cpp
// Function: alglib_impl::hermitianmatrixvectormultiply(alglib_impl::ae_matrix*, bool, long, long, alglib_impl::ae_vector*, alglib_impl::ae_complex, alglib_impl::ae_vector*, alglib_impl::ae_state*)
void hermitianmatrixvectormultiply(/* Complex */ ae_matrix* a,
     ae_bool isupper,
     ae_int_t i1,
     ae_int_t i2,
     /* Complex */ ae_vector* x,
     ae_complex alpha,
     /* Complex */ ae_vector* y,
     ae_state *_state)
{
    ae_int_t i;
    ae_int_t ba1;
    ae_int_t ba2;
    ae_int_t by1;
    ae_int_t by2;
    ae_int_t bx1;
    ae_int_t bx2;
    ae_int_t n;
    ae_complex v;


    n = i2-i1+1;
    if( n<=0 )
    {
        return;
    }
    
    /*
     * Let A = L + D + U, where
     *  L is strictly lower triangular (main diagonal is zero)
     *  D is diagonal
     *  U is strictly upper triangular (main diagonal is zero)
     *
     * A*x = L*x + D*x + U*x
     *
     * Calculate D*x first
     */
    for(i=i1; i<=i2; i++)
    {
        y->ptr.p_complex[i-i1+1] = ae_c_mul(a->ptr.pp_complex[i][i],x->ptr.p_complex[i-i1+1]);
    }
    
    /*
     * Add L*x + U*x
     */
    if( isupper )
    {
        for(i=i1; i<=i2-1; i++)
        {
            
            /*
             * Add L*x to the result
             */
            v = x->ptr.p_complex[i-i1+1];
            by1 = i-i1+2;
            by2 = n;
            ba1 = i+1;
            ba2 = i2;
            ae_v_caddc(&y->ptr.p_complex[by1], 1, &a->ptr.pp_complex[i][ba1], 1, "Conj", ae_v_len(by1,by2), v);
            
            /*
             * Add U*x to the result
             */
            bx1 = i-i1+2;
            bx2 = n;
            ba1 = i+1;
            ba2 = i2;
            v = ae_v_cdotproduct(&x->ptr.p_complex[bx1], 1, "N", &a->ptr.pp_complex[i][ba1], 1, "N", ae_v_len(bx1,bx2));
            y->ptr.p_complex[i-i1+1] = ae_c_add(y->ptr.p_complex[i-i1+1],v);
        }
    }
    else
    {
        for(i=i1+1; i<=i2; i++)
        {
            
            /*
             * Add L*x to the result
             */
            bx1 = 1;
            bx2 = i-i1;
            ba1 = i1;
            ba2 = i-1;
            v = ae_v_cdotproduct(&x->ptr.p_complex[bx1], 1, "N", &a->ptr.pp_complex[i][ba1], 1, "N", ae_v_len(bx1,bx2));
            y->ptr.p_complex[i-i1+1] = ae_c_add(y->ptr.p_complex[i-i1+1],v);
            
            /*
             * Add U*x to the result
             */
            v = x->ptr.p_complex[i-i1+1];
            by1 = 1;
            by2 = i-i1;
            ba1 = i1;
            ba2 = i-1;
            ae_v_caddc(&y->ptr.p_complex[by1], 1, &a->ptr.pp_complex[i][ba1], 1, "Conj", ae_v_len(by1,by2), v);
        }
    }
    ae_v_cmulc(&y->ptr.p_complex[1], 1, ae_v_len(1,n), alpha);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x1c0, %rsp            # imm = 0x1C0
movb %sil, %al
movq 0x10(%rbp), %rsi
movsd %xmm0, -0x10(%rbp)
movsd %xmm1, -0x8(%rbp)
movq %rdi, -0x18(%rbp)
andb $0x1, %al
movb %al, -0x19(%rbp)
movq %rdx, -0x28(%rbp)
movq %rcx, -0x30(%rbp)
movq %r8, -0x38(%rbp)
movq %r9, -0x40(%rbp)
movq -0x30(%rbp), %rax
subq -0x28(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x80(%rbp)
cmpq $0x0, -0x80(%rbp)
jg 0x24a51
jmp 0x250de
movq -0x28(%rbp), %rax
movq %rax, -0x48(%rbp)
movq -0x48(%rbp), %rax
cmpq -0x30(%rbp), %rax
jg 0x24b3d
movq -0x18(%rbp), %rax
movq 0x38(%rax), %rax
movq -0x48(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movq -0x48(%rbp), %rcx
shlq $0x4, %rcx
addq %rcx, %rax
movq (%rax), %rcx
movq %rcx, -0xb0(%rbp)
movq 0x8(%rax), %rax
movq %rax, -0xa8(%rbp)
movq -0x38(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x48(%rbp), %rcx
subq -0x28(%rbp), %rcx
addq $0x1, %rcx
shlq $0x4, %rcx
addq %rcx, %rax
movq (%rax), %rcx
movq %rcx, -0xc0(%rbp)
movq 0x8(%rax), %rax
movq %rax, -0xb8(%rbp)
movsd -0xb0(%rbp), %xmm0
movsd -0xa8(%rbp), %xmm1
movsd -0xc0(%rbp), %xmm2
movsd -0xb8(%rbp), %xmm3
callq 0x3f510
movsd %xmm0, -0xa0(%rbp)
movsd %xmm1, -0x98(%rbp)
movq -0x40(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x48(%rbp), %rcx
subq -0x28(%rbp), %rcx
addq $0x1, %rcx
shlq $0x4, %rcx
addq %rcx, %rax
movq -0xa0(%rbp), %rcx
movq %rcx, (%rax)
movq -0x98(%rbp), %rcx
movq %rcx, 0x8(%rax)
movq -0x48(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x48(%rbp)
jmp 0x24a59
testb $0x1, -0x19(%rbp)
je 0x24dec
movq -0x28(%rbp), %rax
movq %rax, -0x48(%rbp)
movq -0x48(%rbp), %rax
movq -0x30(%rbp), %rcx
subq $0x1, %rcx
cmpq %rcx, %rax
jg 0x24de7
movq -0x38(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x48(%rbp), %rcx
subq -0x28(%rbp), %rcx
addq $0x1, %rcx
shlq $0x4, %rcx
addq %rcx, %rax
movq (%rax), %rcx
movq %rcx, -0x90(%rbp)
movq 0x8(%rax), %rax
movq %rax, -0x88(%rbp)
movq -0x48(%rbp), %rax
subq -0x28(%rbp), %rax
addq $0x2, %rax
movq %rax, -0x60(%rbp)
movq -0x80(%rbp), %rax
movq %rax, -0x68(%rbp)
movq -0x48(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x50(%rbp)
movq -0x30(%rbp), %rax
movq %rax, -0x58(%rbp)
movq -0x40(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x60(%rbp), %rcx
shlq $0x4, %rcx
addq %rcx, %rax
movq %rax, -0x190(%rbp)
movq -0x18(%rbp), %rax
movq 0x38(%rax), %rax
movq -0x48(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movq -0x50(%rbp), %rcx
shlq $0x4, %rcx
addq %rcx, %rax
movq %rax, -0x188(%rbp)
movq -0x60(%rbp), %rdi
movq -0x68(%rbp), %rsi
callq 0x42370
movq -0x190(%rbp), %rdi
movq -0x188(%rbp), %rdx
movq %rax, %r9
movq -0x90(%rbp), %rax
movq %rax, -0xd0(%rbp)
movq -0x88(%rbp), %rax
movq %rax, -0xc8(%rbp)
movsd -0xd0(%rbp), %xmm0
movsd -0xc8(%rbp), %xmm1
movl $0x1, %ecx
leaq 0x330e0(%rip), %r8      # 0x57d32
movq %rcx, %rsi
callq 0x40fb0
movq -0x48(%rbp), %rax
subq -0x28(%rbp), %rax
addq $0x2, %rax
movq %rax, -0x70(%rbp)
movq -0x80(%rbp), %rax
movq %rax, -0x78(%rbp)
movq -0x48(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x50(%rbp)
movq -0x30(%rbp), %rax
movq %rax, -0x58(%rbp)
movq -0x38(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x70(%rbp), %rcx
shlq $0x4, %rcx
addq %rcx, %rax
movq %rax, -0x180(%rbp)
movq -0x18(%rbp), %rax
movq 0x38(%rax), %rax
movq -0x48(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movq -0x50(%rbp), %rcx
shlq $0x4, %rcx
addq %rcx, %rax
movq %rax, -0x178(%rbp)
movq -0x70(%rbp), %rdi
movq -0x78(%rbp), %rsi
callq 0x42370
movq -0x180(%rbp), %rdi
movq -0x178(%rbp), %rcx
movl $0x1, %r8d
leaq 0x3385e(%rip), %r9      # 0x58548
movq %r8, %rsi
movq %r9, %rdx
movq %rax, (%rsp)
callq 0x3faa0
movsd %xmm0, -0xe0(%rbp)
movsd %xmm1, -0xd8(%rbp)
movq -0xe0(%rbp), %rax
movq %rax, -0x90(%rbp)
movq -0xd8(%rbp), %rax
movq %rax, -0x88(%rbp)
movq -0x40(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x48(%rbp), %rcx
subq -0x28(%rbp), %rcx
addq $0x1, %rcx
shlq $0x4, %rcx
addq %rcx, %rax
movq (%rax), %rcx
movq %rcx, -0x100(%rbp)
movq 0x8(%rax), %rax
movq %rax, -0xf8(%rbp)
movq -0x90(%rbp), %rax
movq %rax, -0x110(%rbp)
movq -0x88(%rbp), %rax
movq %rax, -0x108(%rbp)
movsd -0x100(%rbp), %xmm0
movsd -0xf8(%rbp), %xmm1
movsd -0x110(%rbp), %xmm2
movsd -0x108(%rbp), %xmm3
callq 0x3f4c0
movsd %xmm0, -0xf0(%rbp)
movsd %xmm1, -0xe8(%rbp)
movq -0x40(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x48(%rbp), %rcx
subq -0x28(%rbp), %rcx
addq $0x1, %rcx
shlq $0x4, %rcx
addq %rcx, %rax
movq -0xf0(%rbp), %rcx
movq %rcx, (%rax)
movq -0xe8(%rbp), %rcx
movq %rcx, 0x8(%rax)
movq -0x48(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x48(%rbp)
jmp 0x24b4f
jmp 0x25083
movq -0x28(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x48(%rbp)
movq -0x48(%rbp), %rax
cmpq -0x30(%rbp), %rax
jg 0x25081
movq $0x1, -0x70(%rbp)
movq -0x48(%rbp), %rax
subq -0x28(%rbp), %rax
movq %rax, -0x78(%rbp)
movq -0x28(%rbp), %rax
movq %rax, -0x50(%rbp)
movq -0x48(%rbp), %rax
subq $0x1, %rax
movq %rax, -0x58(%rbp)
movq -0x38(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x70(%rbp), %rcx
shlq $0x4, %rcx
addq %rcx, %rax
movq %rax, -0x1b0(%rbp)
movq -0x18(%rbp), %rax
movq 0x38(%rax), %rax
movq -0x48(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movq -0x50(%rbp), %rcx
shlq $0x4, %rcx
addq %rcx, %rax
movq %rax, -0x1a8(%rbp)
movq -0x70(%rbp), %rdi
movq -0x78(%rbp), %rsi
callq 0x42370
movq -0x1b0(%rbp), %rdi
movq -0x1a8(%rbp), %rcx
movl $0x1, %r8d
leaq 0x336b6(%rip), %r9      # 0x58548
movq %r8, %rsi
movq %r9, %rdx
movq %rax, (%rsp)
callq 0x3faa0
movsd %xmm0, -0x120(%rbp)
movsd %xmm1, -0x118(%rbp)
movq -0x120(%rbp), %rax
movq %rax, -0x90(%rbp)
movq -0x118(%rbp), %rax
movq %rax, -0x88(%rbp)
movq -0x40(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x48(%rbp), %rcx
subq -0x28(%rbp), %rcx
addq $0x1, %rcx
shlq $0x4, %rcx
addq %rcx, %rax
movq (%rax), %rcx
movq %rcx, -0x140(%rbp)
movq 0x8(%rax), %rax
movq %rax, -0x138(%rbp)
movq -0x90(%rbp), %rax
movq %rax, -0x150(%rbp)
movq -0x88(%rbp), %rax
movq %rax, -0x148(%rbp)
movsd -0x140(%rbp), %xmm0
movsd -0x138(%rbp), %xmm1
movsd -0x150(%rbp), %xmm2
movsd -0x148(%rbp), %xmm3
callq 0x3f4c0
movsd %xmm0, -0x130(%rbp)
movsd %xmm1, -0x128(%rbp)
movq -0x40(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x48(%rbp), %rcx
subq -0x28(%rbp), %rcx
addq $0x1, %rcx
shlq $0x4, %rcx
addq %rcx, %rax
movq -0x130(%rbp), %rcx
movq %rcx, (%rax)
movq -0x128(%rbp), %rcx
movq %rcx, 0x8(%rax)
movq -0x38(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x48(%rbp), %rcx
subq -0x28(%rbp), %rcx
addq $0x1, %rcx
shlq $0x4, %rcx
addq %rcx, %rax
movq (%rax), %rcx
movq %rcx, -0x90(%rbp)
movq 0x8(%rax), %rax
movq %rax, -0x88(%rbp)
movq $0x1, -0x60(%rbp)
movq -0x48(%rbp), %rax
subq -0x28(%rbp), %rax
movq %rax, -0x68(%rbp)
movq -0x28(%rbp), %rax
movq %rax, -0x50(%rbp)
movq -0x48(%rbp), %rax
subq $0x1, %rax
movq %rax, -0x58(%rbp)
movq -0x40(%rbp), %rax
movq 0x28(%rax), %rax
movq -0x60(%rbp), %rcx
shlq $0x4, %rcx
addq %rcx, %rax
movq %rax, -0x1a0(%rbp)
movq -0x18(%rbp), %rax
movq 0x38(%rax), %rax
movq -0x48(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movq -0x50(%rbp), %rcx
shlq $0x4, %rcx
addq %rcx, %rax
movq %rax, -0x198(%rbp)
movq -0x60(%rbp), %rdi
movq -0x68(%rbp), %rsi
callq 0x42370
movq -0x1a0(%rbp), %rdi
movq -0x198(%rbp), %rdx
movq %rax, %r9
movq -0x90(%rbp), %rax
movq %rax, -0x160(%rbp)
movq -0x88(%rbp), %rax
movq %rax, -0x158(%rbp)
movsd -0x160(%rbp), %xmm0
movsd -0x158(%rbp), %xmm1
movl $0x1, %ecx
leaq 0x32cca(%rip), %r8      # 0x57d32
movq %rcx, %rsi
callq 0x40fb0
movq -0x48(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x48(%rbp)
jmp 0x24df8
jmp 0x25083
movq -0x40(%rbp), %rax
movq 0x28(%rax), %rax
addq $0x10, %rax
movq %rax, -0x1b8(%rbp)
movq -0x80(%rbp), %rsi
movl $0x1, %edi
callq 0x42370
movq -0x1b8(%rbp), %rdi
movq %rax, %rdx
movq -0x10(%rbp), %rax
movq %rax, -0x170(%rbp)
movq -0x8(%rbp), %rax
movq %rax, -0x168(%rbp)
movsd -0x170(%rbp), %xmm0
movsd -0x168(%rbp), %xmm1
movl $0x1, %esi
callq 0x41780
addq $0x1c0, %rsp            # imm = 0x1C0
popq %rbp
retq
nopw (%rax,%rax)