Fix LDLT with semi-definite complex matrices: owing to round-off errors, the diagonal was not real. Also exploit the fact that the diagonal is real in the rest of LDLT

1 files changed, 1 insertions, 1 deletions

diff --git a/test/cholesky.cpp b/test/cholesky.cpp
index 9aa7da7bb..1d147bd7a 100644
--- a/test/cholesky.cpp
+++ b/test/cholesky.cpp
@@ -209,7 +209,7 @@ template<typename MatrixType> void cholesky(const MatrixType& m)
       VERIFY_IS_APPROX(A * vecX, vecB);
     }
     
-    // check matrices with wide spectrum
+    // check matrices with a wide spectrum
     if(rows>=3)
     {
       RealScalar s = (std::min)(16,std::numeric_limits<RealScalar>::max_exponent10/8);