Test refactoring: poly evaluation

Author: JohanMabille

test/CMakeLists.txt

@@ -153,6 +153,7 @@ set(XSIMD_TESTS
     test_hyperbolic.cpp
     test_interface.cpp
     test_memory.cpp
+    test_poly_evaluation.cpp
     test_utils.hpp
     #[[    xsimd_api_test.hpp
     xsimd_api_test.cpp

test/test_poly_evaluation.cpp

@@ -0,0 +1,61 @@
+/***************************************************************************
+* Copyright (c) Johan Mabille, Sylvain Corlay, Wolf Vollprecht and         *
+* Martin Renou                                                             *
+* Copyright (c) QuantStack                                                 *
+*                                                                          *
+* Distributed under the terms of the BSD 3-Clause License.                 *
+*                                                                          *
+* The full license is in the file LICENSE, distributed with this software. *
+****************************************************************************/
+
+#include "test_utils.hpp"
+
+template <class B>
+class poly_evaluation_test : public testing::Test
+{
+protected:
+
+    using batch_type = B;
+    using value_type = typename B::value_type;
+    static constexpr size_t size = B::size;
+    using vector_type = std::vector<value_type>;
+
+    size_t nb_input;
+    vector_type input;
+    vector_type horner_res;
+    vector_type estrin_res;
+
+    poly_evaluation_test()
+    {
+        nb_input = size * 10000;
+        input.resize(nb_input);
+        for (size_t i = 0; i < nb_input; ++i)
+        {
+            input[i] = value_type(i) / 4 + value_type(1.2) * std::sqrt(value_type(i + 0.25));
+        }
+        horner_res.resize(nb_input);
+        estrin_res.resize(nb_input);
+    }
+
+    void test_poly_evaluation()
+    {
+        batch_type in, out;
+        for (size_t i = 0; i < nb_input; i += size)
+        {
+            detail::load_batch(in, input, i);
+            out = xsimd::horner<batch_type, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16>(in);
+            detail::store_batch(out, horner_res, i);
+            out = xsimd::estrin<batch_type, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16>(in);
+            detail::store_batch(out, estrin_res, i);
+        }
+        size_t diff = detail::get_nb_diff(horner_res, estrin_res);
+        EXPECT_EQ(diff, 0) << print_function_name("estrin");
+    }
+};
+
+TYPED_TEST_SUITE(poly_evaluation_test, batch_float_types, simd_test_names);
+
+TYPED_TEST(poly_evaluation_test, poly_evaluation)
+{
+    this->test_poly_evaluation();
+}