From ee04c5c8742baa5deb4fe9944da1b6d4319b4b6a Mon Sep 17 00:00:00 2001
From: Gael Guennebaud <g.gael@free.fr>
Date: Tue, 17 Feb 2009 09:57:32 +0000
Subject: add tests showing bug in unitOrthogonal such that we don't forget it!

---
 test/geo_orthomethods.cpp | 10 +++++++++-
 1 file changed, 9 insertions(+), 1 deletion(-)

(limited to 'test/geo_orthomethods.cpp')

diff --git a/test/geo_orthomethods.cpp b/test/geo_orthomethods.cpp
index 08110e6d8..1485877d0 100644
--- a/test/geo_orthomethods.cpp
+++ b/test/geo_orthomethods.cpp
@@ -76,6 +76,15 @@ template<typename Scalar, int Size> void orthomethods(int size=Size)
   VERIFY_IS_MUCH_SMALLER_THAN(v0.unitOrthogonal().dot(v0), Scalar(1));
   VERIFY_IS_APPROX(v0.unitOrthogonal().norm(), Scalar(1));
 
+  if (size>3)
+  {
+    v0.template start<3>().setZero();
+    v0.end(size-3).setRandom();
+
+    VERIFY_IS_MUCH_SMALLER_THAN(v0.unitOrthogonal().dot(v0), Scalar(1));
+    VERIFY_IS_APPROX(v0.unitOrthogonal().norm(), Scalar(1));
+  }
+
   // colwise/rowwise cross product
   Vector3 vec3 = Vector3::Random();
   int i = ei_random<int>(0,size-1);
@@ -89,7 +98,6 @@ template<typename Scalar, int Size> void orthomethods(int size=Size)
   matN3.setRandom();
   mcrossN3 = matN3.rowwise().cross(vec3);
   VERIFY_IS_APPROX(mcrossN3.row(i), matN3.row(i).cross(vec3));
-
 }
 
 void test_geo_orthomethods()
-- 
cgit v1.2.3