Symaptic: os/ossrv/ossrv_pub/boost_apis/boost/numeric/interval/arith3.hpp@260cb5ec6c19

     1 /* Boost interval/arith3.hpp template implementation file

     2  *

     3  * This headers provides arithmetical functions

     4  * which compute an interval given some base

     5  * numbers. The resulting interval encloses the

     6  * real result of the arithmetic operation.

     7  *

     8  * Copyright 2003 Guillaume Melquiond

     9  *

    10  * Distributed under the Boost Software License, Version 1.0.

    11  * (See accompanying file LICENSE_1_0.txt or

    12  * copy at http://www.boost.org/LICENSE_1_0.txt)

    13  */

    15 #ifndef BOOST_NUMERIC_INTERVAL_ARITH3_HPP

    16 #define BOOST_NUMERIC_INTERVAL_ARITH3_HPP

    18 #include <boost/numeric/interval/detail/interval_prototype.hpp>

    19 #include <boost/numeric/interval/detail/test_input.hpp>

    21 namespace boost {

    22 namespace numeric {

    23 namespace interval_lib {

    25 template<class I> inline

    26 I add(const typename I::base_type& x, const typename I::base_type& y)

    27 {

    28   typedef typename I::traits_type Policies;

    29   if (detail::test_input<typename I::base_type, Policies>(x, y))

    30     return I::empty();

    31   typename Policies::rounding rnd;

    32   return I(rnd.add_down(x, y), rnd.add_up(x, y), true);

    33 }

    35 template<class I> inline

    36 I sub(const typename I::base_type& x, const typename I::base_type& y)

    37 {

    38   typedef typename I::traits_type Policies;

    39   if (detail::test_input<typename I::base_type, Policies>(x, y))

    40     return I::empty();

    41   typename Policies::rounding rnd;

    42   return I(rnd.sub_down(x, y), rnd.sub_up(x, y), true);

    43 }

    45 template<class I> inline

    46 I mul(const typename I::base_type& x, const typename I::base_type& y)

    47 {

    48   typedef typename I::traits_type Policies;

    49   if (detail::test_input<typename I::base_type, Policies>(x, y))

    50     return I::empty();

    51   typename Policies::rounding rnd;

    52   return I(rnd.mul_down(x, y), rnd.mul_up(x, y), true);

    53 }

    55 template<class I> inline

    56 I div(const typename I::base_type& x, const typename I::base_type& y)

    57 {

    58   typedef typename I::traits_type Policies;

    59   if (detail::test_input<typename I::base_type, Policies>(x, y) || user::is_zero(y))

    60     return I::empty();

    61   typename Policies::rounding rnd;

    62   return I(rnd.div_down(x, y), rnd.div_up(x, y), true);

    63 }

    65 } // namespace interval_lib

    66 } // namespace numeric

    67 } // namespace boost

    69 #endif // BOOST_NUMERIC_INTERVAL_ARITH3_HPP

author	sl
	Tue, 10 Jun 2014 14:32:02 +0200
changeset 1	260cb5ec6c19
permissions	-rw-r--r--