// File: Kingsford-Group[P]ghost2/include/Eigen/src/Core/GeneralProduct.h
// Function: void Eigen::internal::gemv_selector<2, 0, true>::run<Eigen::GeneralProduct<Eigen::Block<Eigen::Matrix<double, -1, -1, 0, -1, -1>, -1, -1, false>, Eigen::Block<Eigen::Matrix<double, -1, 1, 0, -1, 1>, -1, 1, false>, 4>, Eigen::Block<Eigen::Matrix<double, -1, -1, 0, -1, -1>, -1, 1, false>>(Eigen::GeneralProduct<Eigen::Block<Eigen::Matrix<double, -1, -1, 0, -1, -1>, -1, -1, false>, Eigen::Block<Eigen::Matrix<double, -1, 1, 0, -1, 1>, -1, 1, false>, 4> const&, Eigen::Block<Eigen::Matrix<double, -1, -1, 0, -1, -1>, -1, 1, false>&, Eigen::GeneralProduct<Eigen::Block<Eigen::Matrix<double, -1, -1, 0, -1, -1>, -1, -1, false>, Eigen::Block<Eigen::Matrix<double, -1, 1, 0, -1, 1>, -1, 1, false>, 4>::Scalar const&)
static inline void run(const ProductType& prod, Dest& dest, const typename ProductType::Scalar& alpha)
  {
    typedef typename ProductType::Index Index;
    typedef typename ProductType::LhsScalar   LhsScalar;
    typedef typename ProductType::RhsScalar   RhsScalar;
    typedef typename ProductType::Scalar      ResScalar;
    typedef typename ProductType::RealScalar  RealScalar;
    typedef typename ProductType::ActualLhsType ActualLhsType;
    typedef typename ProductType::ActualRhsType ActualRhsType;
    typedef typename ProductType::LhsBlasTraits LhsBlasTraits;
    typedef typename ProductType::RhsBlasTraits RhsBlasTraits;
    typedef Map<Matrix<ResScalar,Dynamic,1>, Aligned> MappedDest;

    ActualLhsType actualLhs = LhsBlasTraits::extract(prod.lhs());
    ActualRhsType actualRhs = RhsBlasTraits::extract(prod.rhs());

    ResScalar actualAlpha = alpha * LhsBlasTraits::extractScalarFactor(prod.lhs())
                                  * RhsBlasTraits::extractScalarFactor(prod.rhs());

    enum {
      // FIXME find a way to allow an inner stride on the result if packet_traits<Scalar>::size==1
      // on, the other hand it is good for the cache to pack the vector anyways...
      EvalToDestAtCompileTime = Dest::InnerStrideAtCompileTime==1,
      ComplexByReal = (NumTraits<LhsScalar>::IsComplex) && (!NumTraits<RhsScalar>::IsComplex),
      MightCannotUseDest = (Dest::InnerStrideAtCompileTime!=1) || ComplexByReal
    };

    gemv_static_vector_if<ResScalar,Dest::SizeAtCompileTime,Dest::MaxSizeAtCompileTime,MightCannotUseDest> static_dest;

    bool alphaIsCompatible = (!ComplexByReal) || (numext::imag(actualAlpha)==RealScalar(0));
    bool evalToDest = EvalToDestAtCompileTime && alphaIsCompatible;
    
    RhsScalar compatibleAlpha = get_factor<ResScalar,RhsScalar>::run(actualAlpha);

    ei_declare_aligned_stack_constructed_variable(ResScalar,actualDestPtr,dest.size(),
                                                  evalToDest ? dest.data() : static_dest.data());
    
    if(!evalToDest)
    {
      #ifdef EIGEN_DENSE_STORAGE_CTOR_PLUGIN
      int size = dest.size();
      EIGEN_DENSE_STORAGE_CTOR_PLUGIN
      #endif
      if(!alphaIsCompatible)
      {
        MappedDest(actualDestPtr, dest.size()).setZero();
        compatibleAlpha = RhsScalar(1);
      }
      else
        MappedDest(actualDestPtr, dest.size()) = dest;
    }

    general_matrix_vector_product
      <Index,LhsScalar,ColMajor,LhsBlasTraits::NeedToConjugate,RhsScalar,RhsBlasTraits::NeedToConjugate>::run(
        actualLhs.rows(), actualLhs.cols(),
        actualLhs.data(), actualLhs.outerStride(),
        actualRhs.data(), actualRhs.innerStride(),
        actualDestPtr, 1,
        compatibleAlpha);

    if (!evalToDest)
    {
      if(!alphaIsCompatible)
        dest += actualAlpha * MappedDest(actualDestPtr, dest.size());
      else
        dest = MappedDest(actualDestPtr, dest.size());
    }
  }

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x1e0, %rsp            # imm = 0x1E0
movq %rdi, -0x10(%rbp)
movq %rsi, -0x18(%rbp)
movq %rdx, -0x20(%rbp)
movq -0x10(%rbp), %rdi
callq 0x540e0
movq %rax, %rdi
callq 0x540d0
movq %rax, -0x28(%rbp)
movq -0x10(%rbp), %rdi
callq 0x54100
movq %rax, %rdi
callq 0x540f0
movq %rax, -0x30(%rbp)
movq -0x20(%rbp), %rax
movsd (%rax), %xmm0
movsd %xmm0, -0x120(%rbp)
movq -0x10(%rbp), %rdi
callq 0x540e0
movq %rax, %rdi
callq 0x54120
movaps %xmm0, %xmm1
movsd -0x120(%rbp), %xmm0
mulsd %xmm1, %xmm0
movsd %xmm0, -0x118(%rbp)
movq -0x10(%rbp), %rdi
callq 0x54100
movq %rax, %rdi
callq 0x54140
movaps %xmm0, %xmm1
movsd -0x118(%rbp), %xmm0
mulsd %xmm1, %xmm0
movsd %xmm0, -0x38(%rbp)
movb $0x1, -0x3a(%rbp)
movb -0x3a(%rbp), %al
andb $0x1, %al
movb %al, -0x3b(%rbp)
leaq -0x38(%rbp), %rdi
callq 0x54160
movsd %xmm0, -0x48(%rbp)
movq -0x18(%rbp), %rdi
callq 0x54180
movq %rax, -0x8(%rbp)
movabsq $0x1fffffffffffffff, %rax # imm = 0x1FFFFFFFFFFFFFFF
cmpq %rax, -0x8(%rbp)
jbe 0x53c7a
callq 0x4a840
testb $0x1, -0x3b(%rbp)
je 0x53c92
movq -0x18(%rbp), %rdi
callq 0x541c0
movq %rax, -0x128(%rbp)
jmp 0x53ca2
leaq -0x39(%rbp), %rdi
callq 0x541e0
movq %rax, -0x128(%rbp)
movq -0x128(%rbp), %rax
cmpq $0x0, %rax
je 0x53ce7
testb $0x1, -0x3b(%rbp)
je 0x53cc7
movq -0x18(%rbp), %rdi
callq 0x541c0
movq %rax, -0x130(%rbp)
jmp 0x53cd7
leaq -0x39(%rbp), %rdi
callq 0x541e0
movq %rax, -0x130(%rbp)
movq -0x130(%rbp), %rax
movq %rax, -0x138(%rbp)
jmp 0x53d4d
movq -0x18(%rbp), %rdi
callq 0x54180
shlq $0x3, %rax
cmpq $0x4e20, %rax           # imm = 0x4E20
ja 0x53d23
movq -0x18(%rbp), %rdi
callq 0x54180
leaq 0xf(,%rax,8), %rcx
andq $-0x10, %rcx
movq %rsp, %rax
subq %rcx, %rax
movq %rax, %rsp
movq %rax, -0x140(%rbp)
jmp 0x53d3f
movq -0x18(%rbp), %rdi
callq 0x54180
movq %rax, %rdi
shlq $0x3, %rdi
callq 0x4a8c0
movq %rax, -0x140(%rbp)
movq -0x140(%rbp), %rax
movq %rax, -0x138(%rbp)
movq -0x138(%rbp), %rax
movq %rax, -0x50(%rbp)
testb $0x1, -0x3b(%rbp)
je 0x53d70
movq -0x18(%rbp), %rdi
callq 0x541c0
movq %rax, -0x148(%rbp)
jmp 0x53d80
leaq -0x39(%rbp), %rdi
callq 0x541e0
movq %rax, -0x148(%rbp)
movq -0x148(%rbp), %rax
cmpq $0x0, %rax
jne 0x53d9a
movq -0x50(%rbp), %rax
movq %rax, -0x150(%rbp)
jmp 0x53da5
xorl %eax, %eax
movq %rax, -0x150(%rbp)
jmp 0x53da5
movq -0x150(%rbp), %rax
movq %rax, -0x160(%rbp)
movq -0x18(%rbp), %rdi
callq 0x54180
movq %rax, -0x158(%rbp)
movq -0x18(%rbp), %rdi
callq 0x54180
movq -0x160(%rbp), %rsi
movq -0x158(%rbp), %rdx
shlq $0x3, %rax
cmpq $0x4e20, %rax           # imm = 0x4E20
seta %al
movzbl %al, %ecx
andl $0x1, %ecx
leaq -0x68(%rbp), %rdi
callq 0x541f0
testb $0x1, -0x3b(%rbp)
jne 0x53f00
testb $0x1, -0x3a(%rbp)
jne 0x53e92
movq -0x50(%rbp), %rax
movq %rax, -0x170(%rbp)
movq -0x18(%rbp), %rdi
callq 0x54180
movq %rax, -0x168(%rbp)
jmp 0x53e27
leaq -0x8e(%rbp), %rdi
callq 0x54230
jmp 0x53e35
movq -0x168(%rbp), %rdx
movq -0x170(%rbp), %rsi
leaq -0x80(%rbp), %rdi
leaq -0x8e(%rbp), %rcx
callq 0x54270
jmp 0x53e55
leaq -0x80(%rbp), %rdi
callq 0x542d0
jmp 0x53e60
movsd 0x68208(%rip), %xmm0    # 0xbc070
movsd %xmm0, -0x48(%rbp)
jmp 0x53efe
movq %rax, %rcx
movl %edx, %eax
movq %rcx, -0x88(%rbp)
movl %eax, -0x8c(%rbp)
leaq -0x68(%rbp), %rdi
callq 0x55c10
jmp 0x540b5
movq -0x18(%rbp), %rdi
movq %rdi, -0x188(%rbp)
movq -0x50(%rbp), %rax
movq %rax, -0x180(%rbp)
callq 0x54180
movq %rax, -0x178(%rbp)
jmp 0x53eb6
leaq -0xaa(%rbp), %rdi
callq 0x54230
jmp 0x53ec4
movq -0x178(%rbp), %rdx
movq -0x180(%rbp), %rsi
leaq -0xa8(%rbp), %rdi
leaq -0xaa(%rbp), %rcx
callq 0x54270
jmp 0x53ee7
movq -0x188(%rbp), %rsi
leaq -0xa8(%rbp), %rdi
callq 0x54300
jmp 0x53efc
jmp 0x53efe
jmp 0x53f00
movq -0x28(%rbp), %rdi
callq 0x53b40
movq %rax, -0x1a0(%rbp)
movq -0x28(%rbp), %rdi
callq 0x53950
movq %rax, -0x198(%rbp)
movq -0x28(%rbp), %rdi
callq 0x55ad0
movq %rax, -0x190(%rbp)
jmp 0x53f32
movq -0x28(%rbp), %rdi
callq 0x55af0
movq %rax, -0x1a8(%rbp)
jmp 0x53f44
movq -0x30(%rbp), %rdi
callq 0x55b10
movq %rax, -0x1b0(%rbp)
jmp 0x53f56
movq -0x30(%rbp), %rdi
callq 0x52be0
movq %rax, -0x1b8(%rbp)
jmp 0x53f68
movq -0x1b8(%rbp), %r9
movq -0x1b0(%rbp), %r8
movq -0x1a8(%rbp), %rcx
movq -0x190(%rbp), %rdx
movq -0x198(%rbp), %rsi
movq -0x1a0(%rbp), %rdi
movq -0x50(%rbp), %r10
movsd -0x48(%rbp), %xmm0
subq $0x10, %rsp
movq %rsp, %rax
movq %r10, (%rax)
movq $0x1, 0x8(%rax)
callq 0x54340
addq $0x10, %rsp
jmp 0x53fb8
testb $0x1, -0x3b(%rbp)
jne 0x540a7
testb $0x1, -0x3a(%rbp)
jne 0x54043
movq -0x50(%rbp), %rax
movq %rax, -0x1c8(%rbp)
movq -0x18(%rbp), %rdi
callq 0x54180
movq %rax, -0x1c0(%rbp)
jmp 0x53fe5
leaq -0xea(%rbp), %rdi
callq 0x54230
jmp 0x53ff3
movq -0x1c0(%rbp), %rdx
movq -0x1c8(%rbp), %rsi
leaq -0xe8(%rbp), %rdi
leaq -0xea(%rbp), %rcx
callq 0x54270
jmp 0x54016
leaq -0xd0(%rbp), %rdi
leaq -0x38(%rbp), %rsi
leaq -0xe8(%rbp), %rdx
callq 0x55b30
jmp 0x5402f
movq -0x18(%rbp), %rdi
leaq -0xd0(%rbp), %rsi
callq 0x55b70
jmp 0x54041
jmp 0x540a5
movq -0x50(%rbp), %rax
movq %rax, -0x1d8(%rbp)
movq -0x18(%rbp), %rdi
callq 0x54180
movq %rax, -0x1d0(%rbp)
jmp 0x54060
leaq -0x10a(%rbp), %rdi
callq 0x54230
jmp 0x5406e
movq -0x1d0(%rbp), %rdx
movq -0x1d8(%rbp), %rsi
leaq -0x108(%rbp), %rdi
leaq -0x10a(%rbp), %rcx
callq 0x54270
jmp 0x54091
movq -0x18(%rbp), %rdi
leaq -0x108(%rbp), %rsi
callq 0x55bd0
jmp 0x540a3
jmp 0x540a5
jmp 0x540a7
leaq -0x68(%rbp), %rdi
callq 0x55c10
movq %rbp, %rsp
popq %rbp
retq
movq -0x88(%rbp), %rdi
callq 0x10cc0
nopw %cs:(%rax,%rax)
nopl (%rax,%rax)