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

    typename add_const<ActualLhsType>::type actualLhs = LhsBlasTraits::extract(prod.lhs());
    typename add_const<ActualRhsType>::type 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...
      DirectlyUseRhs = _ActualRhsType::InnerStrideAtCompileTime==1
    };

    gemv_static_vector_if<RhsScalar,_ActualRhsType::SizeAtCompileTime,_ActualRhsType::MaxSizeAtCompileTime,!DirectlyUseRhs> static_rhs;

    ei_declare_aligned_stack_constructed_variable(RhsScalar,actualRhsPtr,actualRhs.size(),
        DirectlyUseRhs ? const_cast<RhsScalar*>(actualRhs.data()) : static_rhs.data());

    if(!DirectlyUseRhs)
    {
      #ifdef EIGEN_DENSE_STORAGE_CTOR_PLUGIN
      int size = actualRhs.size();
      EIGEN_DENSE_STORAGE_CTOR_PLUGIN
      #endif
      Map<typename _ActualRhsType::PlainObject>(actualRhsPtr, actualRhs.size()) = actualRhs;
    }

    general_matrix_vector_product
      <Index,LhsScalar,RowMajor,LhsBlasTraits::NeedToConjugate,RhsScalar,RhsBlasTraits::NeedToConjugate>::run(
        actualLhs.rows(), actualLhs.cols(),
        actualLhs.data(), actualLhs.outerStride(),
        actualRhsPtr, 1,
        dest.data(), dest.innerStride(),
        actualAlpha);
  }

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0xf0, %rsp
movq %rdi, -0x10(%rbp)
movq %rsi, -0x18(%rbp)
movq %rdx, -0x20(%rbp)
movq -0x10(%rbp), %rdi
callq 0x639d0
movq %rax, %rdi
callq 0x5ea30
movq %rax, -0x28(%rbp)
movq -0x10(%rbp), %rdi
callq 0x63a20
movq %rax, %rdi
callq 0x639e0
movq %rax, -0x30(%rbp)
movq -0x20(%rbp), %rax
movsd (%rax), %xmm0
movsd %xmm0, -0x80(%rbp)
movq -0x10(%rbp), %rdi
callq 0x639d0
movq %rax, %rdi
callq 0x63a40
movaps %xmm0, %xmm1
movsd -0x80(%rbp), %xmm0
mulsd %xmm1, %xmm0
movsd %xmm0, -0x78(%rbp)
movq -0x10(%rbp), %rdi
callq 0x63a20
movq %rax, %rdi
callq 0x63a70
movaps %xmm0, %xmm1
movsd -0x78(%rbp), %xmm0
mulsd %xmm1, %xmm0
movsd %xmm0, -0x38(%rbp)
leaq -0x30(%rbp), %rdi
callq 0x63aa0
movq %rax, -0x8(%rbp)
movabsq $0x1fffffffffffffff, %rax # imm = 0x1FFFFFFFFFFFFFFF
cmpq %rax, -0x8(%rbp)
jbe 0x636c4
callq 0x4a840
leaq -0x30(%rbp), %rdi
callq 0x63ae0
cmpq $0x0, %rax
je 0x636e5
leaq -0x30(%rbp), %rdi
callq 0x63ae0
movq %rax, -0x88(%rbp)
jmp 0x6374b
leaq -0x30(%rbp), %rdi
callq 0x63aa0
shlq $0x3, %rax
cmpq $0x4e20, %rax           # imm = 0x4E20
ja 0x63721
leaq -0x30(%rbp), %rdi
callq 0x63aa0
leaq 0xf(,%rax,8), %rcx
andq $-0x10, %rcx
movq %rsp, %rax
subq %rcx, %rax
movq %rax, %rsp
movq %rax, -0x90(%rbp)
jmp 0x6373d
leaq -0x30(%rbp), %rdi
callq 0x63aa0
movq %rax, %rdi
shlq $0x3, %rdi
callq 0x4a8c0
movq %rax, -0x90(%rbp)
movq -0x90(%rbp), %rax
movq %rax, -0x88(%rbp)
movq -0x88(%rbp), %rax
movq %rax, -0x48(%rbp)
leaq -0x30(%rbp), %rdi
callq 0x63ae0
cmpq $0x0, %rax
jne 0x63772
movq -0x48(%rbp), %rax
movq %rax, -0x98(%rbp)
jmp 0x6377d
xorl %eax, %eax
movq %rax, -0x98(%rbp)
jmp 0x6377d
movq -0x98(%rbp), %rax
movq %rax, -0xb0(%rbp)
leaq -0x30(%rbp), %rdi
movq %rdi, -0xb8(%rbp)
callq 0x63aa0
movq -0xb8(%rbp), %rdi
movq %rax, -0xa8(%rbp)
callq 0x63aa0
movq -0xb0(%rbp), %rsi
movq -0xa8(%rbp), %rdx
shlq $0x3, %rax
subq $0x4e21, %rax           # imm = 0x4E21
setae %al
movzbl %al, %ecx
leaq -0x60(%rbp), %rdi
callq 0x541f0
leaq -0x28(%rbp), %rdi
callq 0x5ea10
movq %rax, -0xa0(%rbp)
jmp 0x637e7
leaq -0x28(%rbp), %rdi
callq 0x5e9f0
movq %rax, -0xc0(%rbp)
jmp 0x637f9
leaq -0x28(%rbp), %rdi
callq 0x63b10
movq %rax, -0xc8(%rbp)
jmp 0x6380b
leaq -0x28(%rbp), %rdi
callq 0x5ed70
movq %rax, -0xd0(%rbp)
jmp 0x6381d
movq -0x48(%rbp), %rax
movq %rax, -0xe0(%rbp)
movq -0x18(%rbp), %rdi
callq 0x63b40
movq %rax, -0xd8(%rbp)
jmp 0x6383a
movq -0x18(%rbp), %rdi
callq 0x63b70
movq %rax, -0xe8(%rbp)
jmp 0x6384c
movq -0xe0(%rbp), %r8
movq -0xd0(%rbp), %rcx
movq -0xc8(%rbp), %rdx
movq -0xc0(%rbp), %rsi
movq -0xa0(%rbp), %rdi
movq -0xd8(%rbp), %r9
movq -0xe8(%rbp), %r10
movsd -0x38(%rbp), %xmm0
subq $0x10, %rsp
movq %rsp, %rax
movq %r10, 0x8(%rax)
movq %r9, (%rax)
movl $0x1, %r9d
callq 0x5ee20
addq $0x10, %rsp
jmp 0x638a1
leaq -0x60(%rbp), %rdi
callq 0x55c10
movq %rbp, %rsp
popq %rbp
retq
movq %rax, %rcx
movl %edx, %eax
movq %rcx, -0x68(%rbp)
movl %eax, -0x6c(%rbp)
leaq -0x60(%rbp), %rdi
callq 0x55c10
movq -0x68(%rbp), %rdi
callq 0x10cc0
nopl (%rax)