From 5339db61645daa2c829898cbd2c3310b5f80f7a8 Mon Sep 17 00:00:00 2001
From: Benoit Jacob <jacob.benoit.1@gmail.com>
Date: Mon, 31 Aug 2009 22:08:43 -0400
Subject: [PATCH] add VERIFY_IS_UNITARY

---
 test/main.h | 8 ++++++++
 1 file changed, 8 insertions(+)

diff --git a/test/main.h b/test/main.h
index e3866c68b..8c84e525c 100644
--- a/test/main.h
+++ b/test/main.h
@@ -153,6 +153,8 @@ namespace Eigen
 #define VERIFY_IS_APPROX_OR_LESS_THAN(a, b) VERIFY(test_ei_isApproxOrLessThan(a, b))
 #define VERIFY_IS_NOT_APPROX_OR_LESS_THAN(a, b) VERIFY(!test_ei_isApproxOrLessThan(a, b))
 
+#define VERIFY_IS_UNITARY(a) VERIFY(test_isUnitary(a))
+
 #define CALL_SUBTEST(FUNC) do { \
     g_test_stack.push_back(EI_PP_MAKE_STRING(FUNC)); \
     FUNC; \
@@ -227,6 +229,12 @@ inline bool test_ei_isMuchSmallerThan(const MatrixBase<Derived>& m,
   return m.isMuchSmallerThan(s, test_precision<typename ei_traits<Derived>::Scalar>());
 }
 
+template<typename Derived>
+inline bool test_isUnitary(const MatrixBase<Derived>& m)
+{
+  return m.isUnitary(test_precision<typename ei_traits<Derived>::Scalar>());
+}
+
 template<typename MatrixType>
 void createRandomMatrixOfRank(int desired_rank, int rows, int cols, MatrixType& m)
 {