// File: j-faria[P]kima/eigen/Eigen/src/Core/products/TriangularSolverVector.h
// Function: Eigen::internal::triangular_solve_vector<double, double, long, 1, 1, false, 0>::run(long, double const*, long, double*)
static void run(Index size, const LhsScalar* _lhs, Index lhsStride, RhsScalar* rhs)
  {
    typedef Map<const Matrix<LhsScalar,Dynamic,Dynamic,ColMajor>, 0, OuterStride<> > LhsMap;
    const LhsMap lhs(_lhs,size,size,OuterStride<>(lhsStride));
    typedef const_blas_data_mapper<LhsScalar,Index,ColMajor> LhsMapper;
    typedef const_blas_data_mapper<RhsScalar,Index,ColMajor> RhsMapper;
    typename internal::conditional<Conjugate,
                                   const CwiseUnaryOp<typename internal::scalar_conjugate_op<LhsScalar>,LhsMap>,
                                   const LhsMap&
                                  >::type cjLhs(lhs);
    static const Index PanelWidth = EIGEN_TUNE_TRIANGULAR_PANEL_WIDTH;

    for(Index pi=IsLower ? 0 : size;
        IsLower ? pi<size : pi>0;
        IsLower ? pi+=PanelWidth : pi-=PanelWidth)
    {
      Index actualPanelWidth = (std::min)(IsLower ? size - pi : pi, PanelWidth);
      Index startBlock = IsLower ? pi : pi-actualPanelWidth;
      Index endBlock = IsLower ? pi + actualPanelWidth : 0;

      for(Index k=0; k<actualPanelWidth; ++k)
      {
        Index i = IsLower ? pi+k : pi-k-1;
        if(numext::not_equal_strict(rhs[i],RhsScalar(0)))
        {
          if(!(Mode & UnitDiag))
            rhs[i] /= cjLhs.coeff(i,i);

          Index r = actualPanelWidth - k - 1; // remaining size
          Index s = IsLower ? i+1 : i-r;
          if (r>0)
            Map<Matrix<RhsScalar,Dynamic,1> >(rhs+s,r) -= rhs[i] * cjLhs.col(i).segment(s,r);
        }
      }
      Index r = IsLower ? size - endBlock : startBlock; // remaining size
      if (r > 0)
      {
        // let's directly call the low level product function because:
        // 1 - it is faster to compile
        // 2 - it is slightly faster at runtime
        general_matrix_vector_product<Index,LhsScalar,LhsMapper,ColMajor,Conjugate,RhsScalar,RhsMapper,false>::run(
            r, actualPanelWidth,
            LhsMapper(&lhs.coeffRef(endBlock,startBlock), lhsStride),
            RhsMapper(rhs+startBlock, 1),
            rhs+endBlock, 1, RhsScalar(-1));
      }
    }
  }

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x340, %rsp            # imm = 0x340
movq %rdi, -0x78(%rbp)
movq %rsi, -0x80(%rbp)
movq %rdx, -0x88(%rbp)
movq %rcx, -0x90(%rbp)
movq -0x80(%rbp), %rax
movq %rax, -0x310(%rbp)
movq -0x78(%rbp), %rax
movq %rax, -0x308(%rbp)
movq -0x78(%rbp), %rax
movq %rax, -0x300(%rbp)
movq -0x88(%rbp), %rsi
leaq -0xc8(%rbp), %rdi
callq 0x4e6e0
movq -0x310(%rbp), %rsi
movq -0x308(%rbp), %rdx
movq -0x300(%rbp), %rcx
leaq -0xb8(%rbp), %rdi
leaq -0xc8(%rbp), %r8
callq 0x4e710
leaq -0xb8(%rbp), %rax
movq %rax, -0xd0(%rbp)
movq $0x0, -0xd8(%rbp)
movq -0xd8(%rbp), %rax
cmpq -0x78(%rbp), %rax
jge 0x4daf9
movq -0x78(%rbp), %rax
subq -0xd8(%rbp), %rax
movq %rax, -0xe8(%rbp)
leaq -0xe8(%rbp), %rdi
leaq 0x78161(%rip), %rsi     # 0xc58a8
callq 0x7cf0
movq (%rax), %rax
movq %rax, -0xe0(%rbp)
movq -0xd8(%rbp), %rax
movq %rax, -0xf0(%rbp)
movq -0xd8(%rbp), %rax
addq -0xe0(%rbp), %rax
movq %rax, -0xf8(%rbp)
movq $0x0, -0x100(%rbp)
movq -0x100(%rbp), %rax
cmpq -0xe0(%rbp), %rax
jge 0x4d94b
movq -0xd8(%rbp), %rax
addq -0x100(%rbp), %rax
movq %rax, -0x108(%rbp)
movq -0x90(%rbp), %rdi
movq -0x108(%rbp), %rax
shlq $0x3, %rax
addq %rax, %rdi
xorps %xmm0, %xmm0
movsd %xmm0, -0x110(%rbp)
leaq -0x110(%rbp), %rsi
callq 0x4e770
testb $0x1, %al
jne 0x4d7e2
jmp 0x4d932
movq -0xd0(%rbp), %rdi
movq -0x108(%rbp), %rsi
movq -0x108(%rbp), %rdx
callq 0x4e7a0
movsd (%rax), %xmm1
movq -0x90(%rbp), %rax
movq -0x108(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
divsd %xmm1, %xmm0
movsd %xmm0, (%rax,%rcx,8)
movq -0xe0(%rbp), %rax
subq -0x100(%rbp), %rax
subq $0x1, %rax
movq %rax, -0x118(%rbp)
movq -0x108(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x120(%rbp)
cmpq $0x0, -0x118(%rbp)
jle 0x4d930
movq -0x90(%rbp), %rax
movq -0x108(%rbp), %rcx
shlq $0x3, %rcx
addq %rcx, %rax
movq %rax, -0x328(%rbp)
movq -0xd0(%rbp), %rsi
movq -0x108(%rbp), %rdx
leaq -0x2b0(%rbp), %rdi
callq 0x4e8c0
movq -0x120(%rbp), %rdx
movq -0x118(%rbp), %rcx
leaq -0x258(%rbp), %rdi
leaq -0x2b0(%rbp), %rsi
callq 0x4e910
movq -0x328(%rbp), %rsi
leaq -0x1d0(%rbp), %rdi
leaq -0x258(%rbp), %rdx
callq 0x4e820
movq -0x90(%rbp), %rax
movq -0x120(%rbp), %rcx
shlq $0x3, %rcx
addq %rcx, %rax
movq %rax, -0x320(%rbp)
movq -0x118(%rbp), %rax
movq %rax, -0x318(%rbp)
leaq -0x2ca(%rbp), %rdi
callq 0x4d5a0
movq -0x320(%rbp), %rsi
movq -0x318(%rbp), %rdx
leaq -0x2c8(%rbp), %rdi
leaq -0x2ca(%rbp), %rcx
callq 0x4e970
leaq -0x2c8(%rbp), %rdi
leaq -0x1d0(%rbp), %rsi
callq 0x4e9d0
jmp 0x4d932
jmp 0x4d934
movq -0x100(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x100(%rbp)
jmp 0x4d784
movq -0x78(%rbp), %rax
subq -0xf8(%rbp), %rax
movq %rax, -0x2d8(%rbp)
cmpq $0x0, -0x2d8(%rbp)
jle 0x4dae0
movq -0x2d8(%rbp), %rax
movq %rax, -0x338(%rbp)
movq -0xe0(%rbp), %rax
movq %rax, -0x330(%rbp)
movq -0xf8(%rbp), %rsi
movq -0xf0(%rbp), %rdx
leaq -0xb8(%rbp), %rdi
callq 0x50980
movq %rax, %rcx
movq -0x88(%rbp), %rax
leaq -0x2e8(%rbp), %rdx
movq %rdx, -0x60(%rbp)
movq %rcx, -0x68(%rbp)
movq %rax, -0x70(%rbp)
movq -0x60(%rbp), %rdx
movq -0x68(%rbp), %rcx
movq -0x70(%rbp), %rax
movq %rdx, -0x8(%rbp)
movq %rcx, -0x10(%rbp)
movq %rax, -0x18(%rbp)
movq $0x1, -0x20(%rbp)
movq -0x8(%rbp), %rax
movq -0x10(%rbp), %rcx
movq %rcx, (%rax)
movq -0x18(%rbp), %rcx
movq %rcx, 0x8(%rax)
cmpq $0x1, -0x20(%rbp)
je 0x4da17
leaq 0x6d244(%rip), %rdi     # 0xbac43
leaq 0x6d245(%rip), %rsi     # 0xbac4b
movl $0xd1, %edx
leaq 0x6d25e(%rip), %rcx     # 0xbac70
callq 0x6340
movq -0x90(%rbp), %rax
movq -0xf0(%rbp), %rcx
shlq $0x3, %rcx
addq %rcx, %rax
leaq -0x2f8(%rbp), %rcx
movq %rcx, -0x48(%rbp)
movq %rax, -0x50(%rbp)
movq $0x1, -0x58(%rbp)
movq -0x48(%rbp), %rdx
movq -0x50(%rbp), %rcx
movq -0x58(%rbp), %rax
movq %rdx, -0x28(%rbp)
movq %rcx, -0x30(%rbp)
movq %rax, -0x38(%rbp)
movq $0x1, -0x40(%rbp)
movq -0x28(%rbp), %rax
movq -0x30(%rbp), %rcx
movq %rcx, (%rax)
movq -0x38(%rbp), %rcx
movq %rcx, 0x8(%rax)
cmpq $0x1, -0x40(%rbp)
je 0x4da9c
leaq 0x6d1bf(%rip), %rdi     # 0xbac43
leaq 0x6d1c0(%rip), %rsi     # 0xbac4b
movl $0xd1, %edx
leaq 0x6d1d9(%rip), %rcx     # 0xbac70
callq 0x6340
movq -0x330(%rbp), %rsi
movq -0x338(%rbp), %rdi
movq -0x90(%rbp), %r8
movq -0xf8(%rbp), %rax
shlq $0x3, %rax
addq %rax, %r8
leaq -0x2e8(%rbp), %rdx
leaq -0x2f8(%rbp), %rcx
movl $0x1, %r9d
movsd 0x69bbd(%rip), %xmm0    # 0xb7698
callq 0x4ea30
jmp 0x4dae2
movq -0xd8(%rbp), %rax
addq $0x8, %rax
movq %rax, -0xd8(%rbp)
jmp 0x4d716
addq $0x340, %rsp            # imm = 0x340
popq %rbp
retq
nopw %cs:(%rax,%rax)
nopl (%rax)