// File: j-faria[P]kima/eigen/Eigen/src/Core/products/SelfadjointProduct.h
// Function: Eigen::selfadjoint_product_selector<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>, 1, false>::run(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> const&, double const&)
static void run(MatrixType& mat, const OtherType& other, const typename MatrixType::Scalar& alpha)
  {
    typedef typename MatrixType::Scalar Scalar;
    typedef internal::blas_traits<OtherType> OtherBlasTraits;
    typedef typename OtherBlasTraits::DirectLinearAccessType ActualOtherType;
    typedef typename internal::remove_all<ActualOtherType>::type _ActualOtherType;
    typename internal::add_const_on_value_type<ActualOtherType>::type actualOther = OtherBlasTraits::extract(other.derived());

    Scalar actualAlpha = alpha * OtherBlasTraits::extractScalarFactor(other.derived());

    enum {
      IsRowMajor = (internal::traits<MatrixType>::Flags&RowMajorBit) ? 1 : 0,
      OtherIsRowMajor = _ActualOtherType::Flags&RowMajorBit ? 1 : 0
    };

    Index size = mat.cols();
    Index depth = actualOther.cols();

    typedef internal::gemm_blocking_space<IsRowMajor ? RowMajor : ColMajor,Scalar,Scalar,
              MatrixType::MaxColsAtCompileTime, MatrixType::MaxColsAtCompileTime, _ActualOtherType::MaxColsAtCompileTime> BlockingType;

    BlockingType blocking(size, size, depth, 1, false);


    internal::general_matrix_matrix_triangular_product<Index,
      Scalar, OtherIsRowMajor ? RowMajor : ColMajor,   OtherBlasTraits::NeedToConjugate  && NumTraits<Scalar>::IsComplex,
      Scalar, OtherIsRowMajor ? ColMajor : RowMajor, (!OtherBlasTraits::NeedToConjugate) && NumTraits<Scalar>::IsComplex,
      IsRowMajor ? RowMajor : ColMajor, MatrixType::InnerStrideAtCompileTime, UpLo>
      ::run(size, depth,
            &actualOther.coeffRef(0,0), actualOther.outerStride(), &actualOther.coeffRef(0,0), actualOther.outerStride(),
            mat.data(), mat.innerStride(), mat.outerStride(), actualAlpha, blocking);
  }

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r14
pushq %rbx
subq $0x100, %rsp            # imm = 0x100
movq %rdi, -0x18(%rbp)
movq %rsi, -0x20(%rbp)
movq %rdx, -0x28(%rbp)
movq -0x20(%rbp), %rdi
callq 0x390f0
movq %rax, %rdi
callq 0x398c0
movq %rax, -0x30(%rbp)
movq -0x28(%rbp), %rax
movsd (%rax), %xmm0
movsd %xmm0, -0xb8(%rbp)
movq -0x20(%rbp), %rdi
callq 0x390f0
movq %rax, %rdi
callq 0x39920
movaps %xmm0, %xmm1
movsd -0xb8(%rbp), %xmm0
mulsd %xmm1, %xmm0
movsd %xmm0, -0x38(%rbp)
movq -0x18(%rbp), %rdi
callq 0x39210
movq %rax, -0x40(%rbp)
movq -0x30(%rbp), %rdi
callq 0x39210
movq %rax, -0x48(%rbp)
movq -0x40(%rbp), %rdx
movq -0x48(%rbp), %rcx
leaq -0x80(%rbp), %rdi
movl $0x1, %r8d
xorl %r9d, %r9d
movl %r9d, -0xac(%rbp)
movq %rdx, %rsi
callq 0x48980
movl -0xac(%rbp), %eax
movq -0x40(%rbp), %rcx
movq %rcx, -0xa8(%rbp)
movq -0x48(%rbp), %rcx
movq %rcx, -0xa0(%rbp)
movq -0x30(%rbp), %rdi
movl %eax, %edx
movq %rdx, %rsi
callq 0x48650
movq %rax, -0x98(%rbp)
jmp 0x48865
movq -0x30(%rbp), %rdi
callq 0x3b900
movq %rax, -0xc8(%rbp)
movq -0x30(%rbp), %rdi
xorl %eax, %eax
movl %eax, %edx
movq %rdx, %rsi
callq 0x48650
movq %rax, -0xc0(%rbp)
jmp 0x4888e
movq -0x30(%rbp), %rdi
callq 0x3b900
movq %rax, -0xd0(%rbp)
movq -0x18(%rbp), %rdi
callq 0x3d290
movq %rax, -0xd8(%rbp)
movq -0x18(%rbp), %rdi
callq 0x3d2b0
movq %rax, -0xe0(%rbp)
movq -0x18(%rbp), %rdi
callq 0x3b900
movq -0xe0(%rbp), %r11
movq -0xd8(%rbp), %r10
movq -0xa8(%rbp), %rdi
movq -0xa0(%rbp), %rsi
movq -0x98(%rbp), %rdx
movq -0xc8(%rbp), %rcx
movq -0xc0(%rbp), %r8
movq -0xd0(%rbp), %r9
movq %rax, %rbx
movq %rsp, %rax
leaq -0x80(%rbp), %r14
movq %r14, 0x20(%rax)
leaq -0x38(%rbp), %r14
movq %r14, 0x18(%rax)
movq %rbx, 0x10(%rax)
movq %r11, 0x8(%rax)
movq %r10, (%rax)
callq 0x48a40
jmp 0x48927
leaq -0x80(%rbp), %rdi
callq 0x495a0
addq $0x100, %rsp            # imm = 0x100
popq %rbx
popq %r14
popq %rbp
retq
movq %rax, %rcx
movl %edx, %eax
movq %rcx, -0x88(%rbp)
movl %eax, -0x8c(%rbp)
leaq -0x80(%rbp), %rdi
callq 0x495a0
movq -0x88(%rbp), %rdi
callq 0x6830
nopw %cs:(%rax,%rax)
nopl (%rax)