// File: /workspace/github_cpp_makefile_star2_O2/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
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x48, %rsp
movq %rsi, %rbx
movq %rdi, %r14
movq (%rdi), %rdi
movq 0x8(%r14), %rax
movq 0x8(%rax), %r15
callq 0x25fba
cmpq 0x8(%rax), %r15
jne 0x25f9b
movq (%r14), %rdi
callq 0x25fe4
movq 0x8(%r14), %rcx
leaq 0x18(%rsp), %r15
movq %rax, (%r15)
movq %rcx, 0x8(%r15)
movq %rbx, %rdi
movq %r15, %rsi
callq 0x2614a
movq (%r14), %rdi
callq 0x26012
movq %rax, (%r15)
leaq 0x18(%rsp), %r13
movq %r13, %rdi
movq %rbx, %rsi
callq 0x2621c
movq (%r14), %rsi
leaq 0x38(%rsp), %r12
movq %r12, %rdi
callq 0x2603c
movq (%r12), %rax
movq %rax, (%r13)
leaq 0x18(%rsp), %rdi
callq 0x2606a
mulsd 0x1a248(%rip), %xmm0    # 0x40108
movsd 0x1a248(%rip), %xmm1    # 0x40110
maxsd %xmm0, %xmm1
movsd %xmm1, 0x8(%rsp)
xorl %r13d, %r13d
leaq 0x10(%rsp), %rbp
movq 0x38(%rsp), %rax
movq 0x8(%rax), %rcx
movq 0x10(%rax), %rax
cmpq %rcx, %rax
cmovlq %rax, %rcx
cmpq %rcx, %r13
jge 0x25f55
movq %r12, %rdi
movq %r13, %rsi
callq 0x26078
andpd 0x1a22a(%rip), %xmm0    # 0x40130
ucomisd 0x8(%rsp), %xmm0
jbe 0x25f3a
movq %r12, %rdi
movq %r13, %rsi
callq 0x26078
movsd %xmm0, 0x10(%rsp)
movq %r15, %rdi
movq %rbx, %rsi
movq %r13, %rdx
callq 0x26a06
movq %r15, %rdi
movq %rbp, %rsi
callq 0x260c6
jmp 0x25f50
movq %r15, %rdi
movq %rbx, %rsi
movq %r13, %rdx
callq 0x26a06
movq %r15, %rdi
callq 0x26112
incq %r13
jmp 0x25eda
movq (%r14), %rdi
callq 0x26122
leaq 0x18(%rsp), %r15
movq %rax, (%r15)
movq %r15, %rdi
movq %rbx, %rsi
callq 0x26c38
movq (%r14), %rdi
callq 0x25fe4
movq %rax, (%r15)
movq %rbx, 0x8(%r15)
leaq 0x18(%rsp), %rsi
movq %rbx, %rdi
callq 0x27ad4
addq $0x48, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
leaq 0x21676(%rip), %rdi     # 0x47618
leaq 0x1ef46(%rip), %rsi     # 0x44eef
leaq 0x21692(%rip), %rcx     # 0x47642
movl $0x1f9, %edx            # imm = 0x1F9
callq 0xe380