From a2feb6f3c7a7c5a88e187a21bbb85f3d6ae4f65b Mon Sep 17 00:00:00 2001
From: Jitse Niesen <jitse@maths.leeds.ac.uk>
Date: Sat, 3 Sep 2011 04:58:06 +0100
Subject: [PATCH] Add defensive assert to MatrixExponential,

---
 unsupported/Eigen/src/MatrixFunctions/MatrixExponential.h | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/unsupported/Eigen/src/MatrixFunctions/MatrixExponential.h b/unsupported/Eigen/src/MatrixFunctions/MatrixExponential.h
index ee8f042de..7fd36eca6 100644
--- a/unsupported/Eigen/src/MatrixFunctions/MatrixExponential.h
+++ b/unsupported/Eigen/src/MatrixFunctions/MatrixExponential.h
@@ -395,6 +395,9 @@ void MatrixExponential<MatrixType>::computeUV(long double)
     MatrixType A = m_M / pow(Scalar(2), m_squarings);
     pade17(A);
   }
+#else
+  // this case should be handled in compute()
+  eigen_assert(false && "Bug in MatrixExponential"); 
 #endif  // LDBL_MANT_DIG
 }