// File: j-faria[P]kima/eigen/Eigen/src/Core/products/TriangularSolverVector.h
// Function: Eigen::internal::triangular_solve_vector<double, double, long, 1, 2, false, 1>::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,RowMajor>, 0, OuterStride<> > LhsMap;
    const LhsMap lhs(_lhs,size,size,OuterStride<>(lhsStride));

    typedef const_blas_data_mapper<LhsScalar,Index,RowMajor> 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 r = IsLower ? pi : size - pi; // 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
        Index startRow = IsLower ? pi : pi-actualPanelWidth;
        Index startCol = IsLower ? 0 : pi;

        general_matrix_vector_product<Index,LhsScalar,LhsMapper,RowMajor,Conjugate,RhsScalar,RhsMapper,false>::run(
          actualPanelWidth, r,
          LhsMapper(&lhs.coeffRef(startRow,startCol), lhsStride),
          RhsMapper(rhs + startCol, 1),
          rhs + startRow, 1,
          RhsScalar(-1));
      }

      for(Index k=0; k<actualPanelWidth; ++k)
      {
        Index i = IsLower ? pi+k : pi-k-1;
        Index s = IsLower ? pi   : i+1;
        if (k>0)
          rhs[i] -= (cjLhs.row(i).segment(s,k).transpose().cwiseProduct(Map<const Matrix<RhsScalar,Dynamic,1> >(rhs+s,k))).sum();

        if((!(Mode & UnitDiag)) && numext::not_equal_strict(rhs[i],RhsScalar(0)))
          rhs[i] /= cjLhs(i,i);
      }
    }
  }

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x3b0, %rsp            # imm = 0x3B0
movq %rdi, -0x78(%rbp)
movq %rsi, -0x80(%rbp)
movq %rdx, -0x88(%rbp)
movq %rcx, -0x90(%rbp)
movq -0x80(%rbp), %rax
movq %rax, -0x388(%rbp)
movq -0x78(%rbp), %rax
movq %rax, -0x380(%rbp)
movq -0x78(%rbp), %rax
movq %rax, -0x378(%rbp)
movq -0x88(%rbp), %rsi
leaq -0xc8(%rbp), %rdi
callq 0x4e6e0
movq -0x388(%rbp), %rsi
movq -0x380(%rbp), %rdx
movq -0x378(%rbp), %rcx
leaq -0xb8(%rbp), %rdi
leaq -0xc8(%rbp), %r8
callq 0x53ed0
leaq -0xb8(%rbp), %rax
movq %rax, -0xd0(%rbp)
movq -0x78(%rbp), %rax
movq %rax, -0xd8(%rbp)
cmpq $0x0, -0xd8(%rbp)
jle 0x53e46
movq -0xd8(%rbp), %rax
movq %rax, -0xe8(%rbp)
leaq -0xe8(%rbp), %rdi
leaq 0x71e10(%rip), %rsi     # 0xc58b0
callq 0x7cf0
movq (%rax), %rax
movq %rax, -0xe0(%rbp)
movq -0x78(%rbp), %rax
subq -0xd8(%rbp), %rax
movq %rax, -0xf0(%rbp)
cmpq $0x0, -0xf0(%rbp)
jle 0x53c67
movq -0xd8(%rbp), %rax
subq -0xe0(%rbp), %rax
movq %rax, -0xf8(%rbp)
movq -0xd8(%rbp), %rax
movq %rax, -0x100(%rbp)
movq -0xe0(%rbp), %rax
movq %rax, -0x398(%rbp)
movq -0xf0(%rbp), %rax
movq %rax, -0x390(%rbp)
movq -0xf8(%rbp), %rsi
movq -0x100(%rbp), %rdx
leaq -0xb8(%rbp), %rdi
callq 0x55eb0
movq %rax, %rcx
movq -0x88(%rbp), %rax
leaq -0x110(%rbp), %rdx
movq %rdx, -0x48(%rbp)
movq %rcx, -0x50(%rbp)
movq %rax, -0x58(%rbp)
movq -0x48(%rbp), %rdx
movq -0x50(%rbp), %rcx
movq -0x58(%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 0x53b9e
leaq 0x670bd(%rip), %rdi     # 0xbac43
leaq 0x670be(%rip), %rsi     # 0xbac4b
movl $0xd1, %edx
leaq 0x6718f(%rip), %rcx     # 0xbad28
callq 0x6340
movq -0x90(%rbp), %rax
movq -0x100(%rbp), %rcx
shlq $0x3, %rcx
addq %rcx, %rax
leaq -0x120(%rbp), %rcx
movq %rcx, -0x60(%rbp)
movq %rax, -0x68(%rbp)
movq $0x1, -0x70(%rbp)
movq -0x60(%rbp), %rdx
movq -0x68(%rbp), %rcx
movq -0x70(%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 0x53c23
leaq 0x67038(%rip), %rdi     # 0xbac43
leaq 0x67039(%rip), %rsi     # 0xbac4b
movl $0xd1, %edx
leaq 0x67052(%rip), %rcx     # 0xbac70
callq 0x6340
movq -0x390(%rbp), %rsi
movq -0x398(%rbp), %rdi
movq -0x90(%rbp), %r8
movq -0xf8(%rbp), %rax
shlq $0x3, %rax
addq %rax, %r8
leaq -0x110(%rbp), %rdx
leaq -0x120(%rbp), %rcx
movl $0x1, %r9d
movsd 0x63a36(%rip), %xmm0    # 0xb7698
callq 0x53f30
movq $0x0, -0x128(%rbp)
movq -0x128(%rbp), %rax
cmpq -0xe0(%rbp), %rax
jge 0x53e2d
movq -0xd8(%rbp), %rax
subq -0x128(%rbp), %rax
subq $0x1, %rax
movq %rax, -0x130(%rbp)
movq -0x130(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x138(%rbp)
cmpq $0x0, -0x128(%rbp)
jle 0x53da9
movq -0xd0(%rbp), %rsi
movq -0x130(%rbp), %rdx
leaq -0x348(%rbp), %rdi
callq 0x55f30
movq -0x138(%rbp), %rdx
movq -0x128(%rbp), %rcx
leaq -0x2f0(%rbp), %rdi
leaq -0x348(%rbp), %rsi
callq 0x55f80
leaq -0x268(%rbp), %rdi
leaq -0x2f0(%rbp), %rsi
callq 0x55fe0
movq -0x90(%rbp), %rax
movq -0x138(%rbp), %rcx
shlq $0x3, %rcx
addq %rcx, %rax
movq %rax, -0x3a8(%rbp)
movq -0x128(%rbp), %rax
movq %rax, -0x3a0(%rbp)
leaq -0x362(%rbp), %rdi
callq 0x4d5a0
movq -0x3a8(%rbp), %rsi
movq -0x3a0(%rbp), %rdx
leaq -0x360(%rbp), %rdi
leaq -0x362(%rbp), %rcx
callq 0x56090
leaq -0x1e0(%rbp), %rdi
leaq -0x268(%rbp), %rsi
leaq -0x360(%rbp), %rdx
callq 0x56020
leaq -0x1e0(%rbp), %rdi
callq 0x560f0
movaps %xmm0, %xmm1
movq -0x90(%rbp), %rax
movq -0x130(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
subsd %xmm1, %xmm0
movsd %xmm0, (%rax,%rcx,8)
movq -0x90(%rbp), %rdi
movq -0x130(%rbp), %rax
shlq $0x3, %rax
addq %rax, %rdi
xorps %xmm0, %xmm0
movsd %xmm0, -0x370(%rbp)
leaq -0x370(%rbp), %rsi
callq 0x4e770
testb $0x1, %al
jne 0x53ddb
jmp 0x53e14
movq -0xd0(%rbp), %rdi
movq -0x130(%rbp), %rsi
movq -0x130(%rbp), %rdx
callq 0x56150
movaps %xmm0, %xmm1
movq -0x90(%rbp), %rax
movq -0x130(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
divsd %xmm1, %xmm0
movsd %xmm0, (%rax,%rcx,8)
jmp 0x53e16
movq -0x128(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x128(%rbp)
jmp 0x53c72
jmp 0x53e2f
movq -0xd8(%rbp), %rax
subq $0x8, %rax
movq %rax, -0xd8(%rbp)
jmp 0x53a76
addq $0x3b0, %rsp            # imm = 0x3B0
popq %rbp
retq
nop