From 25544dbec3429848226c9a567ccd7e82973c04e7 Mon Sep 17 00:00:00 2001 From: Chen-Pang He Date: Wed, 10 Jul 2013 02:36:34 +0800 Subject: [PATCH] Add assertion against undefined matrix power. --- unsupported/Eigen/src/MatrixFunctions/MatrixPower.h | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/unsupported/Eigen/src/MatrixFunctions/MatrixPower.h b/unsupported/Eigen/src/MatrixFunctions/MatrixPower.h index e0a687978..e4f13c993 100644 --- a/unsupported/Eigen/src/MatrixFunctions/MatrixPower.h +++ b/unsupported/Eigen/src/MatrixFunctions/MatrixPower.h @@ -398,8 +398,11 @@ void MatrixPower::initialize() } m_nulls = rows() - m_rank; - if (m_nulls) + if (m_nulls) { + eigen_assert(m_T.bottomRightCorner(m_nulls, m_nulls).isZero() + && "Base of matrix power should be invertible or with a semisimple zero eigenvalue."); m_fT.bottomRows(m_nulls).fill(RealScalar(0)); + } } template