FM / interval

/usr/include/boost/numeric/interval

Name	Size	Mode	Actions
compare/	-	0755	rm
detail/	-	0755	rm
ext/	-	0755	rm
arith.hpp	10353	0644	edit dl rm
arith2.hpp	10092	0644	edit dl rm
arith3.hpp	2150	0644	edit dl rm
checking.hpp	2824	0644	edit dl rm
compare.hpp	705	0644	edit dl rm
constants.hpp	3177	0644	edit dl rm
hw_rounding.hpp	2441	0644	edit dl rm
interval.hpp	12198	0644	edit dl rm
io.hpp	1174	0644	edit dl rm
limits.hpp	1747	0644	edit dl rm
policies.hpp	1771	0644	edit dl rm
rounded_arith.hpp	4783	0644	edit dl rm
rounded_transc.hpp	5269	0644	edit dl rm
rounding.hpp	2595	0644	edit dl rm
transc.hpp	7917	0644	edit dl rm
utility.hpp	10741	0644	edit dl rm
utility_fwd.hpp	4589	0644	edit dl rm

Edit: /usr/include/boost/numeric/interval/arith3.hpp (2150B)

#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

/* 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 #include namespace boost { namespace numeric { namespace interval_lib { template 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(x, y)) return I::empty(); typename Policies::rounding rnd; return I(rnd.add_down(x, y), rnd.add_up(x, y), true); } template 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(x, y)) return I::empty(); typename Policies::rounding rnd; return I(rnd.sub_down(x, y), rnd.sub_up(x, y), true); } template 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(x, y)) return I::empty(); typename Policies::rounding rnd; return I(rnd.mul_down(x, y), rnd.mul_up(x, y), true); } template 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(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