From 41f80b26cbe6c67b0ca5da0cd91e284cfe9071f9 Mon Sep 17 00:00:00 2001
From: Gael Guennebaud <g.gael@free.fr>
Date: Wed, 4 Feb 2009 20:20:34 +0000
Subject: [PATCH] bugfix in LDLt for size==1

---
 Eigen/src/Cholesky/LDLT.h | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/Eigen/src/Cholesky/LDLT.h b/Eigen/src/Cholesky/LDLT.h
index d0fbddcb3..eeccf9725 100644
--- a/Eigen/src/Cholesky/LDLT.h
+++ b/Eigen/src/Cholesky/LDLT.h
@@ -127,6 +127,8 @@ void LDLT<MatrixType>::compute(const MatrixType& a)
   m_matrix = a;
 
   if (size <= 1) {
+    m_p.setZero();
+    m_transpositions.setZero();
     return;
   }