From cf89d9371ac848b42ef0e960f278995e3b3a5213 Mon Sep 17 00:00:00 2001
From: Gael Guennebaud <g.gael@free.fr>
Date: Tue, 27 Jan 2009 23:01:53 +0000
Subject: [PATCH] LLT: makes the non positive definite test less strict, but we
 still need something better.

---
 Eigen/src/Cholesky/LDLT.h | 2 +-
 Eigen/src/Cholesky/LLT.h  | 2 +-
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/Eigen/src/Cholesky/LDLT.h b/Eigen/src/Cholesky/LDLT.h
index 4729c3621..4177000fd 100644
--- a/Eigen/src/Cholesky/LDLT.h
+++ b/Eigen/src/Cholesky/LDLT.h
@@ -97,7 +97,7 @@ void LDLT<MatrixType>::compute(const MatrixType& a)
   const int size = a.rows();
   m_matrix.resize(size, size);
   m_isPositiveDefinite = true;
-  const RealScalar eps = ei_sqrt(precision<Scalar>());
+  const RealScalar eps = precision<Scalar>();
 
   if (size<=1)
   {
diff --git a/Eigen/src/Cholesky/LLT.h b/Eigen/src/Cholesky/LLT.h
index c4e7fdffd..e02b82551 100644
--- a/Eigen/src/Cholesky/LLT.h
+++ b/Eigen/src/Cholesky/LLT.h
@@ -97,7 +97,7 @@ void LLT<MatrixType>::compute(const MatrixType& a)
   assert(a.rows()==a.cols());
   const int size = a.rows();
   m_matrix.resize(size, size);
-  const RealScalar eps = ei_sqrt(precision<Scalar>());
+  const RealScalar eps = precision<Scalar>();
 
   RealScalar x;
   x = ei_real(a.coeff(0,0));