1 //============================================================================== 2 // Copyright 2014 LRI UMR 8623 CNRS/Univ Paris Sud XI 3 // Copyright 2014 NumScale SAS 4 // 5 // Distributed under the Boost Software License, Version 1.0. 6 // See accompanying file LICENSE.txt or copy at 7 // http://www.boost.org/LICENSE_1_0.txt 8 //============================================================================== 9 #define BOOST_TEST_MODULE odeint_nt2_algebra_dispatcher 10 11 #include <boost/test/included/unit_test.hpp> 12 #include <boost/test/floating_point_comparison.hpp> 13 #include <boost/numeric/odeint/external/nt2/nt2_algebra_dispatcher.hpp> 14 #include <boost/numeric/odeint/algebra/default_operations.hpp> 15 #include <boost/mpl/list.hpp> 16 17 #include <boost/preprocessor/repetition.hpp> 18 #include <boost/preprocessor/arithmetic/mul.hpp> 19 20 #include <nt2/table.hpp> 21 #include <nt2/sdk/meta/as.hpp> 22 #include <nt2/include/functions/ones.hpp> 23 24 using namespace boost::unit_test; 25 using namespace boost::numeric::odeint; 26 27 typedef boost::mpl::list< float , double > fp_types; 28 29 #define TABLE(z,n,text) nt2::table<T> y ## n = \ 30 nt2::ones(1,2,nt2::meta::as_<T>() )*T(BOOST_PP_ADD(n,1)); 31 32 #define PARAMS(z,n,text) T(BOOST_PP_ADD(n,1)), 33 34 #define SUM(z,n,text) +BOOST_PP_MUL(BOOST_PP_ADD(n,3),BOOST_PP_ADD(n,2)) 35 36 #define TEST(z,n,text) BOOST_CHECK_SMALL( y0(BOOST_PP_ADD(n,1)) \ 37 -T( 2 BOOST_PP_REPEAT(text, SUM, text) ), T(1e-10) ); 38 39 #define TEST_CASE(z,n,text) BOOST_AUTO_TEST_CASE_TEMPLATE ( \ 40 BOOST_PP_CAT(odeint_foreach, n), T, fp_types ) \ 41 { \ 42 vector_space_algebra algebra; \ 43 BOOST_PP_REPEAT(BOOST_PP_ADD(n,2),TABLE,tt) \ 44 BOOST_PP_CAT(algebra.for_each,BOOST_PP_ADD(n,2))( \ 45 BOOST_PP_ENUM_PARAMS(BOOST_PP_ADD(n,2),y), default_operations:: \ 46 BOOST_PP_CAT(scale_sum,BOOST_PP_ADD(n,1)) <T>( \ 47 BOOST_PP_REPEAT(n, PARAMS, text ) T(BOOST_PP_ADD(n,1)))); \ 48 BOOST_PP_REPEAT(2,TEST,n) \ 49 } 50 51 BOOST_AUTO_TEST_SUITE( nt2_algebra ) 52 53 BOOST_PP_REPEAT(7,TEST_CASE,dummy) 54 55 BOOST_AUTO_TEST_SUITE_END() 56