Big rework of the Traits system, unittests updated.

tvmet-1.7.1/testsuite/testtraits.cpp

@@ -0,0 +1,126 @@
+/* This file is part of Eigen, a C++ template library for linear algebra
+ * Copyright (C) 2007 Benoit Jacob <jacob@math.jussieu.fr>
+ *
+ * Based on Tvmet source code, http://tvmet.sourceforge.net,
+ * Copyright (C) 2001 - 2003 Olaf Petzold <opetzold@users.sourceforge.net>
+ *
+ * This library is free software; you can redistribute it and/or
+ * modify it under the terms of the GNU Lesser General Public
+ * License as published by the Free Software Foundation; either
+ * version 2.1 of the License, or (at your option) any later version.
+ *
+ * This library is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+ * Lesser General Public License for more details.
+ *
+ * You should have received a copy of the GNU Lesser General Public
+ * License along with this library; if not, write to the Free Software
+ * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
+ *
+ * $Id: SelfTest.cc,v 1.1 2004/04/24 11:55:15 opetzold Exp $
+ */
+
+#include "main.h"
+
+template<typename T> struct TestTraits
+{
+  void real()
+  {
+    T x = Traits<T>::random();
+    typedef typename Traits<T>::real_type real_type;
+    real_type r = Traits<T>::real(x);
+    TEST_APPROX(r, x);
+  }
+
+  void imag()
+  {
+    T x = Traits<T>::random();
+    typedef typename Traits<T>::real_type real_type;
+    real_type r = Traits<T>::imag(x);
+    TEST_ZERO(r);
+  }
+  
+  void conj()
+  {
+    T x = Traits<T>::random();
+    typedef typename Traits<T>::real_type conj_type;
+    conj_type r = Traits<T>::conj(x);
+    TEST_APPROX(r, x);
+  }
+
+  void abs()
+  {
+    T x = Traits<T>::random();
+    typedef typename Traits<T>::real_type value_type;
+    value_type r1 = Traits<T>::abs(x);
+    value_type r2 = Traits<T>::abs(-x);
+    TEST_APPROX(r1, r2);
+  }
+  
+  void sqrt()
+  {
+    T x = Traits<T>::random();
+    T a = Traits<T>::abs(x);
+    T b = Traits<T>::sqrt(a);
+    // T could be an integer type, so b*b=a is not necessarily true
+    TEST_LESSTHAN(b*b, a);
+    TEST_LESSTHAN(a, (b+1)*(b+1));
+  }
+  
+  void isApprox()
+  {
+    T x = Traits<T>::random();
+    TEST(Traits<T>::isApprox(x,x));
+    TEST(Traits<T>::isApprox(x,x+Traits<T>::epsilon()/10));
+    TEST(!Traits<T>::isApprox(x,x+1));
+  }
+  
+  void isNegligible()
+  {
+    T zero(0), one(1), x = Traits<T>::random();
+    TEST(Traits<T>::isNegligible(zero, zero));
+    TEST(Traits<T>::isNegligible(zero, one));
+    TEST(!Traits<T>::isNegligible(one, one));
+    if(x != zero) {
+      TEST(Traits<T>::isNegligible(zero, x));
+    }
+    TEST(!Traits<T>::isNegligible(one, x));
+  }
+  
+  void isZero()
+  {
+    T zero(0), one(1), x = Traits<T>::random(), y = Traits<T>::random();
+    TEST(Traits<T>::isZero(zero));
+    TEST(!Traits<T>::isZero(one));
+    TEST(Traits<T>::isZero((x+y)*(x+y)-x*x-y*y-x*y-y*x));
+  }
+  
+  void isLessThan()
+  {
+    T one(1), x = Traits<T>::random();
+    TEST(Traits<T>::isLessThan(x, x+one));
+    TEST(!Traits<T>::isLessThan(x+one, x));
+    TEST(Traits<T>::isLessThan(x, x+Traits<T>::epsilon()/10));
+  }
+  
+  TestTraits()
+  {
+    real();
+    imag();
+    conj();
+    abs();
+    sqrt();
+    isApprox();
+    isNegligible();
+    isZero();
+    isLessThan();
+  }
+};
+
+void TvmetTestSuite::testTraits()
+{
+  TestTraits<int>();
+  TestTraits<float>();
+  TestTraits<double>();
+}