// File: Kingsford-Group[P]ghost2/include/Eigen/src/Core/products/SelfadjointMatrixVector.h
// Function: void Eigen::SelfadjointProductMatrix<Eigen::Block<Eigen::Matrix<double, -1, -1, 0, -1, -1>, -1, -1, false>, 17, false, Eigen::CwiseUnaryOp<Eigen::internal::scalar_multiple_op<double>, Eigen::Block<Eigen::Block<Eigen::Matrix<double, -1, -1, 0, -1, -1>, -1, 1, true>, -1, 1, false> const>, 0, true>::scaleAndAddTo<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>&, double const&) const
void scaleAndAddTo(Dest& dest, const Scalar& alpha) const
  {
    typedef typename Dest::Scalar ResScalar;
    typedef typename Base::RhsScalar RhsScalar;
    typedef Map<Matrix<ResScalar,Dynamic,1>, Aligned> MappedDest;
    
    eigen_assert(dest.rows()==m_lhs.rows() && dest.cols()==m_rhs.cols());

    typename internal::add_const_on_value_type<ActualLhsType>::type lhs = LhsBlasTraits::extract(m_lhs);
    typename internal::add_const_on_value_type<ActualRhsType>::type rhs = RhsBlasTraits::extract(m_rhs);

    Scalar actualAlpha = alpha * LhsBlasTraits::extractScalarFactor(m_lhs)
                               * RhsBlasTraits::extractScalarFactor(m_rhs);

    enum {
      EvalToDest = (Dest::InnerStrideAtCompileTime==1),
      UseRhs = (_ActualRhsType::InnerStrideAtCompileTime==1)
    };
    
    internal::gemv_static_vector_if<ResScalar,Dest::SizeAtCompileTime,Dest::MaxSizeAtCompileTime,!EvalToDest> static_dest;
    internal::gemv_static_vector_if<RhsScalar,_ActualRhsType::SizeAtCompileTime,_ActualRhsType::MaxSizeAtCompileTime,!UseRhs> static_rhs;

    ei_declare_aligned_stack_constructed_variable(ResScalar,actualDestPtr,dest.size(),
                                                  EvalToDest ? dest.data() : static_dest.data());
                                                  
    ei_declare_aligned_stack_constructed_variable(RhsScalar,actualRhsPtr,rhs.size(),
        UseRhs ? const_cast<RhsScalar*>(rhs.data()) : static_rhs.data());
    
    if(!EvalToDest)
    {
      #ifdef EIGEN_DENSE_STORAGE_CTOR_PLUGIN
      int size = dest.size();
      EIGEN_DENSE_STORAGE_CTOR_PLUGIN
      #endif
      MappedDest(actualDestPtr, dest.size()) = dest;
    }
      
    if(!UseRhs)
    {
      #ifdef EIGEN_DENSE_STORAGE_CTOR_PLUGIN
      int size = rhs.size();
      EIGEN_DENSE_STORAGE_CTOR_PLUGIN
      #endif
      Map<typename _ActualRhsType::PlainObject>(actualRhsPtr, rhs.size()) = rhs;
    }
      
      
    internal::selfadjoint_matrix_vector_product<Scalar, Index, (internal::traits<_ActualLhsType>::Flags&RowMajorBit) ? RowMajor : ColMajor, int(LhsUpLo), bool(LhsBlasTraits::NeedToConjugate), bool(RhsBlasTraits::NeedToConjugate)>::run
      (
        lhs.rows(),                             // size
        &lhs.coeffRef(0,0),  lhs.outerStride(), // lhs info
        actualRhsPtr, 1,                        // rhs info
        actualDestPtr,                          // result info
        actualAlpha                             // scale factor
      );
    
    if(!EvalToDest)
      dest = MappedDest(actualDestPtr, dest.size());
  }

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x170, %rsp            # imm = 0x170
movq %rdi, -0x18(%rbp)
movq %rsi, -0x20(%rbp)
movq %rdx, -0x28(%rbp)
movq -0x18(%rbp), %rax
movq %rax, -0xb0(%rbp)
movq -0x20(%rbp), %rdi
callq 0x52870
movq -0xb0(%rbp), %rdi
movq %rax, -0xa8(%rbp)
callq 0x53b40
movq -0xa8(%rbp), %rcx
movq %rax, %rdx
xorl %eax, %eax
cmpq %rdx, %rcx
movb %al, -0x99(%rbp)
jne 0x6e97b
movq -0x20(%rbp), %rdi
callq 0x528b0
movq -0xb0(%rbp), %rdi
movq %rax, -0xb8(%rbp)
addq $0x28, %rdi
callq 0x6edd0
movq %rax, %rcx
movq -0xb8(%rbp), %rax
cmpq %rcx, %rax
sete %al
movb %al, -0x99(%rbp)
movb -0x99(%rbp), %al
testb $0x1, %al
jne 0x6e987
jmp 0x6e989
jmp 0x6e9a8
leaq 0x5c21f(%rip), %rdi     # 0xcabaf
leaq 0x5c24f(%rip), %rsi     # 0xcabe6
movl $0xc4, %edx
leaq 0x5c27d(%rip), %rcx     # 0xcac20
callq 0x104a0
movq -0xb0(%rbp), %rdi
callq 0x540d0
movq -0xb0(%rbp), %rdi
movq %rax, -0x30(%rbp)
addq $0x28, %rdi
callq 0x6edf0
movq -0xb0(%rbp), %rdi
movq %rax, -0x38(%rbp)
movq -0x28(%rbp), %rax
movsd (%rax), %xmm0
movsd %xmm0, -0xc8(%rbp)
callq 0x54120
movq -0xb0(%rbp), %rdi
movaps %xmm0, %xmm1
movsd -0xc8(%rbp), %xmm0
mulsd %xmm1, %xmm0
movsd %xmm0, -0xc0(%rbp)
addq $0x28, %rdi
callq 0x6ee20
movaps %xmm0, %xmm1
movsd -0xc0(%rbp), %xmm0
mulsd %xmm1, %xmm0
movsd %xmm0, -0x40(%rbp)
movq -0x20(%rbp), %rdi
callq 0x52910
movq %rax, -0x8(%rbp)
movabsq $0x1fffffffffffffff, %rax # imm = 0x1FFFFFFFFFFFFFFF
cmpq %rax, -0x8(%rbp)
jbe 0x6ea45
callq 0x4a840
movq -0x20(%rbp), %rdi
callq 0x6ee70
cmpq $0x0, %rax
je 0x6ea66
movq -0x20(%rbp), %rdi
callq 0x6ee70
movq %rax, -0xd0(%rbp)
jmp 0x6eacc
movq -0x20(%rbp), %rdi
callq 0x52910
shlq $0x3, %rax
cmpq $0x4e20, %rax           # imm = 0x4E20
ja 0x6eaa2
movq -0x20(%rbp), %rdi
callq 0x52910
leaq 0xf(,%rax,8), %rcx
andq $-0x10, %rcx
movq %rsp, %rax
subq %rcx, %rax
movq %rax, %rsp
movq %rax, -0xd8(%rbp)
jmp 0x6eabe
movq -0x20(%rbp), %rdi
callq 0x52910
movq %rax, %rdi
shlq $0x3, %rdi
callq 0x4a8c0
movq %rax, -0xd8(%rbp)
movq -0xd8(%rbp), %rax
movq %rax, -0xd0(%rbp)
movq -0xd0(%rbp), %rax
movq %rax, -0x50(%rbp)
movq -0x20(%rbp), %rdi
callq 0x6ee70
cmpq $0x0, %rax
jne 0x6eaf3
movq -0x50(%rbp), %rax
movq %rax, -0xe0(%rbp)
jmp 0x6eafe
xorl %eax, %eax
movq %rax, -0xe0(%rbp)
jmp 0x6eafe
movq -0xe0(%rbp), %rax
movq %rax, -0xf8(%rbp)
movq -0x20(%rbp), %rdi
callq 0x52910
movq %rax, -0xf0(%rbp)
movq -0x20(%rbp), %rdi
callq 0x52910
movq -0xf8(%rbp), %rsi
movq -0xf0(%rbp), %rdx
shlq $0x3, %rax
subq $0x4e21, %rax           # imm = 0x4E21
setae %al
movzbl %al, %ecx
leaq -0x68(%rbp), %rdi
callq 0x541f0
movq -0x38(%rbp), %rdi
callq 0x6b670
movq %rax, -0xe8(%rbp)
jmp 0x6eb5e
movq -0xe8(%rbp), %rax
movq %rax, -0x10(%rbp)
movabsq $0x1fffffffffffffff, %rax # imm = 0x1FFFFFFFFFFFFFFF
cmpq %rax, -0x10(%rbp)
jbe 0x6eb82
callq 0x4a840
jmp 0x6eb80
jmp 0x6eb82
jmp 0x6eb84
movq -0x38(%rbp), %rdi
callq 0x6ee90
movq %rax, -0x100(%rbp)
jmp 0x6eb96
movq -0x100(%rbp), %rax
cmpq $0x0, %rax
je 0x6ebc8
movq -0x38(%rbp), %rdi
callq 0x6ee90
movq %rax, -0x108(%rbp)
jmp 0x6ebb5
movq -0x108(%rbp), %rax
movq %rax, -0x110(%rbp)
jmp 0x6ec6d
movq -0x38(%rbp), %rdi
callq 0x6b670
movq %rax, -0x118(%rbp)
jmp 0x6ebda
movq -0x118(%rbp), %rax
shlq $0x3, %rax
cmpq $0x4e20, %rax           # imm = 0x4E20
ja 0x6ec24
movq -0x38(%rbp), %rdi
callq 0x6b670
movq %rax, -0x120(%rbp)
jmp 0x6ebff
movq -0x120(%rbp), %rax
leaq 0xf(,%rax,8), %rcx
andq $-0x10, %rcx
movq %rsp, %rax
subq %rcx, %rax
movq %rax, %rsp
movq %rax, -0x128(%rbp)
jmp 0x6ec5f
movq -0x38(%rbp), %rdi
callq 0x6b670
movq %rax, -0x130(%rbp)
jmp 0x6ec36
movq -0x130(%rbp), %rdi
shlq $0x3, %rdi
callq 0x4a8c0
movq %rax, -0x138(%rbp)
jmp 0x6ec4f
movq -0x138(%rbp), %rax
movq %rax, -0x128(%rbp)
jmp 0x6ec5f
movq -0x128(%rbp), %rax
movq %rax, -0x110(%rbp)
movq -0x110(%rbp), %rax
movq %rax, -0x80(%rbp)
movq -0x38(%rbp), %rdi
callq 0x6ee90
movq %rax, -0x140(%rbp)
jmp 0x6ec8a
movq -0x140(%rbp), %rax
cmpq $0x0, %rax
jne 0x6eca4
movq -0x80(%rbp), %rax
movq %rax, -0x148(%rbp)
jmp 0x6ecaf
xorl %eax, %eax
movq %rax, -0x148(%rbp)
jmp 0x6ecaf
movq -0x148(%rbp), %rax
movq %rax, -0x158(%rbp)
movq -0x38(%rbp), %rdi
callq 0x6b670
movq %rax, -0x150(%rbp)
jmp 0x6eccf
movq -0x38(%rbp), %rdi
callq 0x6b670
movq %rax, -0x160(%rbp)
jmp 0x6ece1
movq -0x150(%rbp), %rdx
movq -0x158(%rbp), %rsi
movq -0x160(%rbp), %rax
shlq $0x3, %rax
subq $0x4e21, %rax           # imm = 0x4E21
setae %al
movzbl %al, %ecx
leaq -0x98(%rbp), %rdi
callq 0x541f0
jmp 0x6ed14
movq -0x30(%rbp), %rdi
callq 0x53b40
movq %rax, -0x170(%rbp)
movq -0x30(%rbp), %rdi
xorl %eax, %eax
movl %eax, %edx
movq %rdx, %rsi
callq 0x6fba0
movq %rax, -0x168(%rbp)
jmp 0x6ed3d
movq -0x30(%rbp), %rdi
callq 0x55af0
movq -0x170(%rbp), %rdi
movq -0x168(%rbp), %rsi
movq %rax, %rdx
movq -0x80(%rbp), %rcx
movq -0x50(%rbp), %r9
movsd -0x40(%rbp), %xmm0
movl $0x1, %r8d
callq 0x6eeb0
jmp 0x6ed71
leaq -0x98(%rbp), %rdi
callq 0x55c10
leaq -0x68(%rbp), %rdi
callq 0x55c10
movq %rbp, %rsp
popq %rbp
retq
movq %rax, %rcx
movl %edx, %eax
movq %rcx, -0x70(%rbp)
movl %eax, -0x74(%rbp)
jmp 0x6edb1
movq %rax, %rcx
movl %edx, %eax
movq %rcx, -0x70(%rbp)
movl %eax, -0x74(%rbp)
leaq -0x98(%rbp), %rdi
callq 0x55c10
leaq -0x68(%rbp), %rdi
callq 0x55c10
movq -0x70(%rbp), %rdi
callq 0x10cc0
nopw %cs:(%rax,%rax)
nopl (%rax)