70 lines
2.2 KiB
C++
70 lines
2.2 KiB
C++
|
/* Boost interval/arith3.hpp template implementation file
|
||
|
|
*
|
||
|
|
* This headers provides arithmetical functions
|
||
|
|
* which compute an interval given some base
|
||
|
|
* numbers. The resulting interval encloses the
|
||
|
|
* real result of the arithmetic operation.
|
||
|
|
*
|
||
|
|
* Copyright 2003 Guillaume Melquiond
|
||
|
|
*
|
||
|
|
* 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)
|
||
|
|
*/
|
||
|
|
|
||
|
|
#ifndef BOOST_NUMERIC_INTERVAL_ARITH3_HPP
|
||
|
|
#define BOOST_NUMERIC_INTERVAL_ARITH3_HPP
|
||
|
|
|
||
|
|
#include <boost/numeric/interval/detail/interval_prototype.hpp>
|
||
|
|
#include <boost/numeric/interval/detail/test_input.hpp>
|
||
|
|
|
||
|
|
namespace boost {
|
||
|
|
namespace numeric {
|
||
|
|
namespace interval_lib {
|
||
|
|
|
||
|
|
template<class I> inline
|
||
|
|
I add(const typename I::base_type& x, const typename I::base_type& y)
|
||
|
|
{
|
||
|
|
typedef typename I::traits_type Policies;
|
||
|
|
if (detail::test_input<typename I::base_type, Policies>(x, y))
|
||
|
|
return I::empty();
|
||
|
|
typename Policies::rounding rnd;
|
||
|
|
return I(rnd.add_down(x, y), rnd.add_up(x, y), true);
|
||
|
|
}
|
||
|
|
|
||
|
|
template<class I> inline
|
||
|
|
I sub(const typename I::base_type& x, const typename I::base_type& y)
|
||
|
|
{
|
||
|
|
typedef typename I::traits_type Policies;
|
||
|
|
if (detail::test_input<typename I::base_type, Policies>(x, y))
|
||
|
|
return I::empty();
|
||
|
|
typename Policies::rounding rnd;
|
||
|
|
return I(rnd.sub_down(x, y), rnd.sub_up(x, y), true);
|
||
|
|
}
|
||
|
|
|
||
|
|
template<class I> inline
|
||
|
|
I mul(const typename I::base_type& x, const typename I::base_type& y)
|
||
|
|
{
|
||
|
|
typedef typename I::traits_type Policies;
|
||
|
|
if (detail::test_input<typename I::base_type, Policies>(x, y))
|
||
|
|
return I::empty();
|
||
|
|
typename Policies::rounding rnd;
|
||
|
|
return I(rnd.mul_down(x, y), rnd.mul_up(x, y), true);
|
||
|
|
}
|
||
|
|
|
||
|
|
template<class I> inline
|
||
|
|
I div(const typename I::base_type& x, const typename I::base_type& y)
|
||
|
|
{
|
||
|
|
typedef typename I::traits_type Policies;
|
||
|
|
if (detail::test_input<typename I::base_type, Policies>(x, y) || user::is_zero(y))
|
||
|
|
return I::empty();
|
||
|
|
typename Policies::rounding rnd;
|
||
|
|
return I(rnd.div_down(x, y), rnd.div_up(x, y), true);
|
||
|
|
}
|
||
|
|
|
||
|
|
} // namespace interval_lib
|
||
|
|
} // namespace numeric
|
||
|
|
} // namespace boost
|
||
|
|
|
||
|
|
#endif // BOOST_NUMERIC_INTERVAL_ARITH3_HPP
|