From 6c02cbbb0f1786f5f22285342850250c28845272 Mon Sep 17 00:00:00 2001
From: Gael Guennebaud <g.gael@free.fr>
Date: Tue, 1 Dec 2015 09:45:56 +0100
Subject: Fix matrix to quaternion (and angleaxis) conversion for matrix
 expression.

---
 test/geo_quaternion.cpp | 6 ++++++
 1 file changed, 6 insertions(+)

(limited to 'test/geo_quaternion.cpp')

diff --git a/test/geo_quaternion.cpp b/test/geo_quaternion.cpp
index 17229be4a..761bb52b4 100644
--- a/test/geo_quaternion.cpp
+++ b/test/geo_quaternion.cpp
@@ -49,6 +49,7 @@ template<typename Scalar, int Options> void quaternion(void)
   */
   using std::abs;
   typedef Matrix<Scalar,3,1> Vector3;
+  typedef Matrix<Scalar,3,3> Matrix3;
   typedef Matrix<Scalar,4,1> Vector4;
   typedef Quaternion<Scalar,Options> Quaternionx;
   typedef AngleAxis<Scalar> AngleAxisx;
@@ -101,6 +102,11 @@ template<typename Scalar, int Options> void quaternion(void)
   q2 = q1.toRotationMatrix();
   VERIFY_IS_APPROX(q1*v1,q2*v1);
 
+  Matrix3 rot1(q1);
+  VERIFY_IS_APPROX(q1*v1,rot1*v1);
+  Quaternionx q3(rot1.transpose()*rot1);
+  VERIFY_IS_APPROX(q3*v1,v1);
+
 
   // angle-axis conversion
   AngleAxisx aa = AngleAxisx(q1);
-- 
cgit v1.2.3