// File: Kingsford-Group[P]ghost2/include/Eigen/src/Cholesky/LDLT.h
// Function: void Eigen::internal::solve_retval<Eigen::LDLT<Eigen::Matrix<double, -1, -1, 0, -1, -1>, 1>, Eigen::Matrix<double, -1, 1, 0, -1, 1>>::evalTo<Eigen::Matrix<double, -1, 1, 0, -1, 1>>(Eigen::Matrix<double, -1, 1, 0, -1, 1>&) const
void evalTo(Dest& dst) const
  {
    eigen_assert(rhs().rows() == dec().matrixLDLT().rows());
    // dst = P b
    dst = dec().transpositionsP() * rhs();

    // dst = L^-1 (P b)
    dec().matrixL().solveInPlace(dst);

    // dst = D^-1 (L^-1 P b)
    // more precisely, use pseudo-inverse of D (see bug 241)
    using std::abs;
    using std::max;
    typedef typename LDLTType::MatrixType MatrixType;
    typedef typename LDLTType::Scalar Scalar;
    typedef typename LDLTType::RealScalar RealScalar;
    const Diagonal<const MatrixType> vectorD = dec().vectorD();
    RealScalar tolerance = (max)(vectorD.array().abs().maxCoeff() * NumTraits<Scalar>::epsilon(),
				 RealScalar(1) / NumTraits<RealScalar>::highest()); // motivated by LAPACK's xGELSS
    for (Index i = 0; i < vectorD.size(); ++i) {
      if(abs(vectorD(i)) > tolerance)
	dst.row(i) /= vectorD(i);
      else
	dst.row(i).setZero();
    }

    // dst = L^-T (D^-1 L^-1 P b)
    dec().matrixU().solveInPlace(dst);

    // dst = P^-1 (L^-T D^-1 L^-1 P b) = A^-1 b
    dst = dec().transpositionsP().transpose() * dst;
  }

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x110, %rsp            # imm = 0x110
movq %rdi, -0x8(%rbp)
movq %rsi, -0x10(%rbp)
movq -0x8(%rbp), %rdi
movq %rdi, -0xf8(%rbp)
callq 0x58a60
movq %rax, %rdi
callq 0x253b0
movq -0xf8(%rbp), %rdi
movq %rax, -0xf0(%rbp)
callq 0x58a80
movq %rax, %rdi
callq 0x58aa0
movq %rax, %rdi
callq 0x24690
movq %rax, %rcx
movq -0xf0(%rbp), %rax
cmpq %rcx, %rax
jne 0x587ff
jmp 0x5881e
leaq 0x6bbc2(%rip), %rdi     # 0xc43c8
leaq 0x689ef(%rip), %rsi     # 0xc11fc
movl $0x1f9, %edx            # imm = 0x1F9
leaq 0x6bbd9(%rip), %rcx     # 0xc43f2
callq 0x104a0
movq -0xf8(%rbp), %rdi
callq 0x58a80
movq %rax, %rdi
callq 0x58b50
movq -0xf8(%rbp), %rdi
movq %rax, -0x108(%rbp)
callq 0x58a60
movq -0x108(%rbp), %rdi
movq %rax, %rsi
callq 0x58b00
movq %rax, -0x20(%rbp)
movq %rdx, -0x18(%rbp)
movq -0x10(%rbp), %rdi
leaq -0x20(%rbp), %rsi
callq 0x58bb0
movq -0xf8(%rbp), %rdi
callq 0x58a80
movq %rax, %rdi
callq 0x58be0
movq %rax, -0x28(%rbp)
movq -0x10(%rbp), %rsi
leaq -0x28(%rbp), %rdi
callq 0x58c50
movq -0xf8(%rbp), %rdi
callq 0x58a80
movq %rax, %rsi
leaq -0x38(%rbp), %rdi
callq 0x58c80
leaq -0x70(%rbp), %rdi
leaq -0x38(%rbp), %rsi
callq 0x58d40
leaq -0x60(%rbp), %rdi
leaq -0x70(%rbp), %rsi
callq 0x58d80
leaq -0x60(%rbp), %rdi
callq 0x58dd0
movsd %xmm0, -0x100(%rbp)
callq 0x4dd20
movaps %xmm0, %xmm1
movsd -0x100(%rbp), %xmm0
mulsd %xmm1, %xmm0
movsd %xmm0, -0x48(%rbp)
callq 0x58e00
movaps %xmm0, %xmm1
movsd 0x63776(%rip), %xmm0    # 0xbc070
divsd %xmm1, %xmm0
movsd %xmm0, -0x78(%rbp)
leaq -0x48(%rbp), %rdi
leaq -0x78(%rbp), %rsi
callq 0x58d00
movsd (%rax), %xmm0
movsd %xmm0, -0x40(%rbp)
movq $0x0, -0x80(%rbp)
movq -0x80(%rbp), %rax
movq %rax, -0x110(%rbp)
leaq -0x38(%rbp), %rdi
callq 0x58e10
movq %rax, %rcx
movq -0x110(%rbp), %rax
cmpq %rcx, %rax
jge 0x589d2
movq -0x80(%rbp), %rsi
leaq -0x38(%rbp), %rdi
callq 0x58e50
callq 0x4dd00
ucomisd -0x40(%rbp), %xmm0
jbe 0x5899f
movq -0x80(%rbp), %rsi
leaq -0x38(%rbp), %rdi
callq 0x58e50
movsd %xmm0, -0x88(%rbp)
movq -0x10(%rbp), %rsi
movq -0x80(%rbp), %rdx
leaq -0xa8(%rbp), %rdi
callq 0x58ee0
leaq -0xa8(%rbp), %rdi
leaq -0x88(%rbp), %rsi
callq 0x58f30
jmp 0x589bf
movq -0x10(%rbp), %rsi
movq -0x80(%rbp), %rdx
leaq -0xc8(%rbp), %rdi
callq 0x58ee0
leaq -0xc8(%rbp), %rdi
callq 0x58fd0
jmp 0x589c1
movq -0x80(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x80(%rbp)
jmp 0x58921
movq -0xf8(%rbp), %rdi
callq 0x58a80
movq %rax, %rdi
callq 0x59000
movq %rax, -0xd0(%rbp)
movq -0x10(%rbp), %rsi
leaq -0xd0(%rbp), %rdi
callq 0x59070
movq -0xf8(%rbp), %rdi
callq 0x58a80
movq %rax, %rdi
callq 0x58b50
movq %rax, %rdi
callq 0x590a0
movq %rax, -0xe8(%rbp)
movq -0x10(%rbp), %rsi
leaq -0xe8(%rbp), %rdi
callq 0x590d0
movq %rax, -0xe0(%rbp)
movq %rdx, -0xd8(%rbp)
movq -0x10(%rbp), %rdi
leaq -0xe0(%rbp), %rsi
callq 0x59120
addq $0x110, %rsp            # imm = 0x110
popq %rbp
retq
nopw (%rax,%rax)