fixed cholesky with 0 sized matrix (cf. #785)
diff --git a/Eigen/src/Cholesky/LLT_LAPACKE.h b/Eigen/src/Cholesky/LLT_LAPACKE.h
index d7d75bd..bde9bcd 100644
--- a/Eigen/src/Cholesky/LLT_LAPACKE.h
+++ b/Eigen/src/Cholesky/LLT_LAPACKE.h
@@ -35,7 +35,7 @@
#include "./InternalHeaderCheck.h"
-namespace Eigen {
+namespace Eigen {
namespace internal {
@@ -106,7 +106,11 @@
template<typename MatrixType>
static Index blocked(MatrixType& m)
{
- eigen_assert(m.rows()==m.cols());
+ eigen_assert(m.rows() == m.cols());
+ if(m.rows() == 0) {
+ return -1;
+ }
+
/* Set up parameters for ?potrf */
lapack_int size = convert_index<lapack_int>(m.rows());
lapack_int StorageOrder = MatrixType::Flags&RowMajorBit?RowMajor:ColMajor;