247 lines
7.8 KiB
C++
247 lines
7.8 KiB
C++
// 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()
|