// File: Luthaf[P]zeopp/./Eigen/src/Core/products/TriangularSolverVector.h
// Function: Eigen::internal::triangular_solve_vector<double, double, long, 1, 2, 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));
    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(!(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 slighlty faster at runtime
        general_matrix_vector_product<Index,LhsScalar,ColMajor,Conjugate,RhsScalar,false>::run(
            r, actualPanelWidth,
            &lhs.coeffRef(endBlock,startBlock), lhsStride,
            rhs+startBlock, 1,
            rhs+endBlock, 1, RhsScalar(-1));
      }
    }
  }

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x240, %rsp            # imm = 0x240
movq %rdi, -0x8(%rbp)
movq %rsi, -0x10(%rbp)
movq %rdx, -0x18(%rbp)
movq %rcx, -0x20(%rbp)
movq -0x10(%rbp), %rax
movq %rax, -0x208(%rbp)
movq -0x8(%rbp), %rax
movq %rax, -0x200(%rbp)
movq -0x8(%rbp), %rax
movq %rax, -0x1f8(%rbp)
movq -0x18(%rbp), %rsi
leaq -0x58(%rbp), %rdi
callq 0x18cf50
movq -0x208(%rbp), %rsi
movq -0x200(%rbp), %rdx
movq -0x1f8(%rbp), %rcx
leaq -0x48(%rbp), %rdi
leaq -0x58(%rbp), %r8
callq 0x18cf80
leaq -0x48(%rbp), %rax
movq %rax, -0x60(%rbp)
movq -0x8(%rbp), %rax
movq %rax, -0x68(%rbp)
cmpq $0x0, -0x68(%rbp)
jle 0x18c7c5
movq -0x68(%rbp), %rax
movq %rax, -0x78(%rbp)
leaq -0x78(%rbp), %rdi
leaq 0x7c8ef(%rip), %rsi     # 0x208e58
callq 0x175fd0
movq (%rax), %rax
movq %rax, -0x70(%rbp)
movq -0x68(%rbp), %rax
subq -0x70(%rbp), %rax
movq %rax, -0x80(%rbp)
movq $0x0, -0x88(%rbp)
movq $0x0, -0x90(%rbp)
movq -0x90(%rbp), %rax
cmpq -0x70(%rbp), %rax
jge 0x18c716
movq -0x68(%rbp), %rax
subq -0x90(%rbp), %rax
subq $0x1, %rax
movq %rax, -0x98(%rbp)
movq -0x60(%rbp), %rdi
movq -0x98(%rbp), %rsi
movq -0x98(%rbp), %rdx
callq 0x18cfe0
movsd (%rax), %xmm1
movq -0x20(%rbp), %rax
movq -0x98(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
divsd %xmm1, %xmm0
movsd %xmm0, (%rax,%rcx,8)
movq -0x70(%rbp), %rax
subq -0x90(%rbp), %rax
subq $0x1, %rax
movq %rax, -0xa0(%rbp)
movq -0x98(%rbp), %rax
subq -0xa0(%rbp), %rax
movq %rax, -0xa8(%rbp)
cmpq $0x0, -0xa0(%rbp)
jle 0x18c6fd
movq -0x20(%rbp), %rax
movq -0x98(%rbp), %rcx
shlq $0x3, %rcx
addq %rcx, %rax
movq %rax, -0x220(%rbp)
movq -0x60(%rbp), %rsi
movq -0x98(%rbp), %rdx
leaq -0x1c8(%rbp), %rdi
callq 0x18d0a0
movq -0xa8(%rbp), %rdx
movq -0xa0(%rbp), %rcx
leaq -0x180(%rbp), %rdi
leaq -0x1c8(%rbp), %rsi
callq 0x18d0f0
movq -0x220(%rbp), %rsi
leaq -0x118(%rbp), %rdi
leaq -0x180(%rbp), %rdx
callq 0x18d060
movq -0x20(%rbp), %rax
movq -0xa8(%rbp), %rcx
shlq $0x3, %rcx
addq %rcx, %rax
movq %rax, -0x218(%rbp)
movq -0xa0(%rbp), %rax
movq %rax, -0x210(%rbp)
leaq -0x1e2(%rbp), %rdi
callq 0x17cec0
movq -0x218(%rbp), %rsi
movq -0x210(%rbp), %rdx
leaq -0x1e0(%rbp), %rdi
leaq -0x1e2(%rbp), %rcx
callq 0x18d140
leaq -0x1e0(%rbp), %rdi
leaq -0x118(%rbp), %rsi
callq 0x18d1a0
jmp 0x18c6ff
movq -0x90(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x90(%rbp)
jmp 0x18c597
movq -0x80(%rbp), %rax
movq %rax, -0x1f0(%rbp)
cmpq $0x0, -0x1f0(%rbp)
jle 0x18c7b2
movq -0x1f0(%rbp), %rax
movq %rax, -0x230(%rbp)
movq -0x70(%rbp), %rax
movq %rax, -0x228(%rbp)
movq -0x88(%rbp), %rsi
movq -0x80(%rbp), %rdx
leaq -0x48(%rbp), %rdi
callq 0x18e990
movq -0x230(%rbp), %rdi
movq -0x228(%rbp), %rsi
movq %rax, %rdx
movq -0x18(%rbp), %rcx
movq -0x20(%rbp), %r8
movq -0x80(%rbp), %rax
shlq $0x3, %rax
addq %rax, %r8
movq -0x20(%rbp), %rax
movq -0x88(%rbp), %r9
shlq $0x3, %r9
addq %r9, %rax
movl $0x1, %r9d
movsd 0x65fd8(%rip), %xmm0    # 0x1f2778
movq %rax, (%rsp)
movq $0x1, 0x8(%rsp)
callq 0x18d200
jmp 0x18c7b4
movq -0x68(%rbp), %rax
subq $0x8, %rax
movq %rax, -0x68(%rbp)
jmp 0x18c54b
addq $0x240, %rsp            # imm = 0x240
popq %rbp
retq
nop