/* 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

