123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246 |
- // Copyright (c) 2018-2019 Cem Bassoy
- //
- // Distributed under the Boost Software License, Version 1.0. (See
- // accompanying file LICENSE_1_0.txt or copy at
- // http://www.boost.org/LICENSE_1_0.txt)
- //
- // The authors gratefully acknowledge the support of
- // Fraunhofer and Google in producing this work
- // which started as a Google Summer of Code project.
- //
- #include <boost/numeric/ublas/tensor/operators_comparison.hpp>
- #include <boost/numeric/ublas/tensor/operators_arithmetic.hpp>
- #include <boost/numeric/ublas/tensor/tensor.hpp>
- #include <boost/test/unit_test.hpp>
- #include <boost/multiprecision/cpp_bin_float.hpp>
- #include "utility.hpp"
- using double_extended = boost::multiprecision::cpp_bin_float_double_extended;
- using test_types = zip<int,long,float,double,double_extended>::with_t<boost::numeric::ublas::first_order, boost::numeric::ublas::last_order>;
- struct fixture {
- using extents_type = boost::numeric::ublas::basic_extents<std::size_t>;
- fixture()
- : extents{
- extents_type{}, // 0
- extents_type{1,1}, // 1
- extents_type{1,2}, // 2
- extents_type{2,1}, // 3
- extents_type{2,3}, // 4
- extents_type{2,3,1}, // 5
- extents_type{4,1,3}, // 6
- extents_type{1,2,3}, // 7
- extents_type{4,2,3}, // 8
- extents_type{4,2,3,5}} // 9
- {
- }
- std::vector<extents_type> extents;
- };
- BOOST_AUTO_TEST_SUITE(test_tensor_comparison, * boost::unit_test::depends_on("test_tensor"))
- BOOST_FIXTURE_TEST_CASE_TEMPLATE( test_tensor_comparison, value, test_types, fixture)
- {
- using namespace boost::numeric;
- using value_type = typename value::first_type;
- using layout_type = typename value::second_type;
- using tensor_type = ublas::tensor<value_type, layout_type>;
- auto check = [](auto const& e)
- {
- auto t = tensor_type (e);
- auto t2 = tensor_type (e);
- auto v = value_type {};
- std::iota(t.begin(), t.end(), v);
- std::iota(t2.begin(), t2.end(), v+2);
- BOOST_CHECK( t == t );
- BOOST_CHECK( t != t2 );
- if(t.empty())
- return;
- BOOST_CHECK(!(t < t));
- BOOST_CHECK(!(t > t));
- BOOST_CHECK( t < t2 );
- BOOST_CHECK( t2 > t );
- BOOST_CHECK( t <= t );
- BOOST_CHECK( t >= t );
- BOOST_CHECK( t <= t2 );
- BOOST_CHECK( t2 >= t );
- BOOST_CHECK( t2 >= t2 );
- BOOST_CHECK( t2 >= t );
- };
- for(auto const& e : extents)
- check(e);
- auto e0 = extents.at(0);
- auto e1 = extents.at(1);
- auto e2 = extents.at(2);
- auto b = false;
- BOOST_CHECK_NO_THROW ( b = (tensor_type(e0) == tensor_type(e0)));
- BOOST_CHECK_NO_THROW ( b = (tensor_type(e1) == tensor_type(e2)));
- BOOST_CHECK_NO_THROW ( b = (tensor_type(e0) == tensor_type(e2)));
- BOOST_CHECK_NO_THROW ( b = (tensor_type(e1) != tensor_type(e2)));
- BOOST_CHECK_THROW ( b = (tensor_type(e1) >= tensor_type(e2)), std::runtime_error );
- BOOST_CHECK_THROW ( b = (tensor_type(e1) <= tensor_type(e2)), std::runtime_error );
- BOOST_CHECK_THROW ( b = (tensor_type(e1) < tensor_type(e2)), std::runtime_error );
- BOOST_CHECK_THROW ( b = (tensor_type(e1) > tensor_type(e2)), std::runtime_error );
- }
- BOOST_FIXTURE_TEST_CASE_TEMPLATE( test_tensor_comparison_with_tensor_expressions, value, test_types, fixture)
- {
- using namespace boost::numeric;
- using value_type = typename value::first_type;
- using layout_type = typename value::second_type;
- using tensor_type = ublas::tensor<value_type, layout_type>;
- auto check = [](auto const& e)
- {
- auto t = tensor_type (e);
- auto t2 = tensor_type (e);
- auto v = value_type {};
- std::iota(t.begin(), t.end(), v);
- std::iota(t2.begin(), t2.end(), v+2);
- BOOST_CHECK( t == t );
- BOOST_CHECK( t != t2 );
- if(t.empty())
- return;
- BOOST_CHECK( !(t < t) );
- BOOST_CHECK( !(t > t) );
- BOOST_CHECK( t < (t2+t) );
- BOOST_CHECK( (t2+t) > t );
- BOOST_CHECK( t <= (t+t) );
- BOOST_CHECK( (t+t2) >= t );
- BOOST_CHECK( (t2+t2+2) >= t);
- BOOST_CHECK( 2*t2 > t );
- BOOST_CHECK( t < 2*t2 );
- BOOST_CHECK( 2*t2 > t);
- BOOST_CHECK( 2*t2 >= t2 );
- BOOST_CHECK( t2 <= 2*t2);
- BOOST_CHECK( 3*t2 >= t );
- };
- for(auto const& e : extents)
- check(e);
- auto e0 = extents.at(0);
- auto e1 = extents.at(1);
- auto e2 = extents.at(2);
- auto b = false;
- BOOST_CHECK_NO_THROW (b = tensor_type(e0) == (tensor_type(e0) + tensor_type(e0)) );
- BOOST_CHECK_NO_THROW (b = tensor_type(e1) == (tensor_type(e2) + tensor_type(e2)) );
- BOOST_CHECK_NO_THROW (b = tensor_type(e0) == (tensor_type(e2) + 2) );
- BOOST_CHECK_NO_THROW (b = tensor_type(e1) != (2 + tensor_type(e2)) );
- BOOST_CHECK_NO_THROW (b = (tensor_type(e0) + tensor_type(e0)) == tensor_type(e0) );
- BOOST_CHECK_NO_THROW (b = (tensor_type(e2) + tensor_type(e2)) == tensor_type(e1) );
- BOOST_CHECK_NO_THROW (b = (tensor_type(e2) + 2) == tensor_type(e0) );
- BOOST_CHECK_NO_THROW (b = (2 + tensor_type(e2)) != tensor_type(e1) );
- BOOST_CHECK_THROW (b = tensor_type(e1) >= (tensor_type(e2) + tensor_type(e2)), std::runtime_error );
- BOOST_CHECK_THROW (b = tensor_type(e1) <= (tensor_type(e2) + tensor_type(e2)), std::runtime_error );
- BOOST_CHECK_THROW (b = tensor_type(e1) < (tensor_type(e2) + tensor_type(e2)), std::runtime_error );
- BOOST_CHECK_THROW (b = tensor_type(e1) > (tensor_type(e2) + tensor_type(e2)), std::runtime_error );
- BOOST_CHECK_THROW (b = tensor_type(e1) >= (tensor_type(e2) + 2), std::runtime_error );
- BOOST_CHECK_THROW (b = tensor_type(e1) <= (2 + tensor_type(e2)), std::runtime_error );
- BOOST_CHECK_THROW (b = tensor_type(e1) < (tensor_type(e2) + 3), std::runtime_error );
- BOOST_CHECK_THROW (b = tensor_type(e1) > (4 + tensor_type(e2)), std::runtime_error );
- }
- BOOST_FIXTURE_TEST_CASE_TEMPLATE( test_tensor_comparison_with_scalar, value, test_types, fixture)
- {
- using namespace boost::numeric;
- using value_type = typename value::first_type;
- using layout_type = typename value::second_type;
- using tensor_type = ublas::tensor<value_type, layout_type>;
- auto check = [](auto const& e)
- {
- BOOST_CHECK( tensor_type(e,value_type{2}) == tensor_type(e,value_type{2}) );
- BOOST_CHECK( tensor_type(e,value_type{2}) != tensor_type(e,value_type{1}) );
- if(e.empty())
- return;
- BOOST_CHECK( !(tensor_type(e,2) < 2) );
- BOOST_CHECK( !(tensor_type(e,2) > 2) );
- BOOST_CHECK( (tensor_type(e,2) >= 2) );
- BOOST_CHECK( (tensor_type(e,2) <= 2) );
- BOOST_CHECK( (tensor_type(e,2) == 2) );
- BOOST_CHECK( (tensor_type(e,2) != 3) );
- BOOST_CHECK( !(2 > tensor_type(e,2)) );
- BOOST_CHECK( !(2 < tensor_type(e,2)) );
- BOOST_CHECK( (2 <= tensor_type(e,2)) );
- BOOST_CHECK( (2 >= tensor_type(e,2)) );
- BOOST_CHECK( (2 == tensor_type(e,2)) );
- BOOST_CHECK( (3 != tensor_type(e,2)) );
- BOOST_CHECK( !( tensor_type(e,2)+3 < 5) );
- BOOST_CHECK( !( tensor_type(e,2)+3 > 5) );
- BOOST_CHECK( ( tensor_type(e,2)+3 >= 5) );
- BOOST_CHECK( ( tensor_type(e,2)+3 <= 5) );
- BOOST_CHECK( ( tensor_type(e,2)+3 == 5) );
- BOOST_CHECK( ( tensor_type(e,2)+3 != 6) );
- BOOST_CHECK( !( 5 > tensor_type(e,2)+3) );
- BOOST_CHECK( !( 5 < tensor_type(e,2)+3) );
- BOOST_CHECK( ( 5 >= tensor_type(e,2)+3) );
- BOOST_CHECK( ( 5 <= tensor_type(e,2)+3) );
- BOOST_CHECK( ( 5 == tensor_type(e,2)+3) );
- BOOST_CHECK( ( 6 != tensor_type(e,2)+3) );
- BOOST_CHECK( !( tensor_type(e,2)+tensor_type(e,3) < 5) );
- BOOST_CHECK( !( tensor_type(e,2)+tensor_type(e,3) > 5) );
- BOOST_CHECK( ( tensor_type(e,2)+tensor_type(e,3) >= 5) );
- BOOST_CHECK( ( tensor_type(e,2)+tensor_type(e,3) <= 5) );
- BOOST_CHECK( ( tensor_type(e,2)+tensor_type(e,3) == 5) );
- BOOST_CHECK( ( tensor_type(e,2)+tensor_type(e,3) != 6) );
- BOOST_CHECK( !( 5 > tensor_type(e,2)+tensor_type(e,3)) );
- BOOST_CHECK( !( 5 < tensor_type(e,2)+tensor_type(e,3)) );
- BOOST_CHECK( ( 5 >= tensor_type(e,2)+tensor_type(e,3)) );
- BOOST_CHECK( ( 5 <= tensor_type(e,2)+tensor_type(e,3)) );
- BOOST_CHECK( ( 5 == tensor_type(e,2)+tensor_type(e,3)) );
- BOOST_CHECK( ( 6 != tensor_type(e,2)+tensor_type(e,3)) );
- };
- for(auto const& e : extents)
- check(e);
- }
- BOOST_AUTO_TEST_SUITE_END()
|