commit 0ce5bc0d14557891ecd39144007e8e11d7887489
author Gael Guennebaud <g.gael@free.fr> Wed Jan 27 23:23:59 2010 +0100
committer Gael Guennebaud <g.gael@free.fr> Wed Jan 27 23:23:59 2010 +0100
tree 79b60fe0eca2051b5d1a4c44f189f862860157ab
parent 7d92b7ad5f9910e194d13787c9fde9c9fd41154d

add support for global math function for array

test/array.cpp

diff --git a/test/array.cpp b/test/array.cpp
index 8a71b40..e51dbac 100644
--- a/test/array.cpp
+++ b/test/array.cpp
@@ -132,6 +132,30 @@
   VERIFY_IS_APPROX(((m1.abs()+1)>RealScalar(0.1)).rowwise().count(), ArrayXi::Constant(rows, cols));
 }
 
+template<typename MatrixType> void array_real(const MatrixType& m)
+{
+  typedef typename MatrixType::Scalar Scalar;
+  typedef typename NumTraits<Scalar>::Real RealScalar;
+
+  int rows = m.rows();
+  int cols = m.cols();
+
+  MatrixType m1 = MatrixType::Random(rows, cols),
+             m2 = MatrixType::Random(rows, cols),
+             m3(rows, cols);
+
+  VERIFY_IS_APPROX(m1.sin(), std::sin(m1));
+  VERIFY_IS_APPROX(m1.sin(), ei_sin(m1));
+  VERIFY_IS_APPROX(m1.cos(), ei_cos(m1));
+  VERIFY_IS_APPROX(m1.cos(), ei_cos(m1));
+  VERIFY_IS_APPROX(m1.abs().sqrt(), std::sqrt(std::abs(m1)));
+  VERIFY_IS_APPROX(m1.abs().sqrt(), ei_sqrt(ei_abs(m1)));
+  VERIFY_IS_APPROX(m1.abs().log(), std::log(std::abs(m1)));
+  VERIFY_IS_APPROX(m1.abs().log(), ei_log(ei_abs(m1)));
+  VERIFY_IS_APPROX(m1.exp(), std::exp(m1));
+  VERIFY_IS_APPROX(m1.exp(), ei_exp(m1));
+}
+
 void test_array()
 {
   for(int i = 0; i < g_repeat; i++) {
@@ -149,4 +173,10 @@
     CALL_SUBTEST_5( comparisons(ArrayXXf(8, 12)) );
     CALL_SUBTEST_6( comparisons(ArrayXXi(8, 12)) );
   }
+  for(int i = 0; i < g_repeat; i++) {
+    CALL_SUBTEST_1( array_real(Array<float, 1, 1>()) );
+    CALL_SUBTEST_2( array_real(Array22f()) );
+    CALL_SUBTEST_3( array_real(Array44d()) );
+    CALL_SUBTEST_5( array_real(ArrayXXf(8, 12)) );
+  }
 }