Sprout/sprout/math/equal_to.hpp

75 lines
2.4 KiB
C++
Raw Normal View History

2013-08-08 09:54:33 +00:00
/*=============================================================================
2015-01-10 10:13:57 +00:00
Copyright (c) 2011-2015 Bolero MURAKAMI
2013-08-08 09:54:33 +00:00
https://github.com/bolero-MURAKAMI/Sprout
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)
=============================================================================*/
2012-07-06 14:10:49 +00:00
#ifndef SPROUT_MATH_EQUAL_TO_HPP
#define SPROUT_MATH_EQUAL_TO_HPP
#include <type_traits>
#include <sprout/config.hpp>
2013-08-06 15:15:09 +00:00
#include <sprout/limits.hpp>
#include <sprout/math/fabs.hpp>
2012-07-07 05:58:46 +00:00
#include <sprout/type_traits/arithmetic_promote.hpp>
2012-07-06 14:10:49 +00:00
#include <sprout/type_traits/enabler_if.hpp>
2013-11-10 13:13:25 +00:00
#include <sprout/algorithm/detail/max.hpp>
2012-07-06 14:10:49 +00:00
namespace sprout {
namespace math {
namespace detail {
template<typename T>
inline SPROUT_CONSTEXPR T
max3(T x, T y, T z) {
return sprout::max(sprout::max(x, y), z);
}
2012-07-07 05:58:46 +00:00
template<typename FloatType>
2012-07-06 14:10:49 +00:00
inline SPROUT_CONSTEXPR bool
2012-07-07 05:58:46 +00:00
equal_to_impl(FloatType x, FloatType y) {
2012-07-06 14:10:49 +00:00
return x == y
|| sprout::fabs(x - y)
2013-08-06 15:15:09 +00:00
<= sprout::numeric_limits<FloatType>::epsilon() * sprout::math::detail::max3(sprout::fabs(x), sprout::fabs(y), FloatType(1))
2012-07-06 14:10:49 +00:00
;
}
template<
2012-07-07 05:58:46 +00:00
typename FloatType1, typename FloatType2,
typename sprout::enabler_if<
std::is_floating_point<typename sprout::arithmetic_promote<FloatType1, FloatType2>::type>::value
>::type = sprout::enabler
>
inline SPROUT_CONSTEXPR bool
equal_to(FloatType1 x, FloatType2 y) {
typedef typename sprout::arithmetic_promote<FloatType1, FloatType2>::type type;
return sprout::math::detail::equal_to_impl<type>(x, y);
}
template<
typename IntType1, typename IntType2,
typename sprout::enabler_if<
std::is_integral<typename sprout::arithmetic_promote<IntType1, IntType2>::type>::value
>::type = sprout::enabler
2012-07-06 14:10:49 +00:00
>
inline SPROUT_CONSTEXPR bool
2012-07-07 05:58:46 +00:00
equal_to(IntType1 x, IntType2 y) {
2012-07-06 14:10:49 +00:00
return x == y;
}
} // namespace detail
//
// equal_to
//
template<
2012-07-07 05:58:46 +00:00
typename T, typename U,
typename sprout::enabler_if<std::is_arithmetic<T>::value && std::is_arithmetic<U>::value>::type = sprout::enabler
2012-07-06 14:10:49 +00:00
>
inline SPROUT_CONSTEXPR bool
2012-07-07 05:58:46 +00:00
equal_to(T x, U y) {
return sprout::math::detail::equal_to(x, y);
2012-07-06 14:10:49 +00:00
}
} // namespace math
} // namespace sprout
#endif // #ifndef SPROUT_MATH_EQUAL_TO_HPP