/////////////////////////////////////////////////////////////////////////////// | |
/// \file complex.hpp | |
/// | |
// Copyright 2005 Eric Niebler. 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_FUNCTIONAL_COMPLEX_HPP_EAN_01_17_2006 | |
#define BOOST_NUMERIC_FUNCTIONAL_COMPLEX_HPP_EAN_01_17_2006 | |
#ifdef BOOST_NUMERIC_FUNCTIONAL_HPP_INCLUDED | |
# error Include this file before boost/accumulators/numeric/functional.hpp | |
#endif | |
#include <complex> | |
#include <boost/mpl/or.hpp> | |
#include <boost/type_traits/is_same.hpp> | |
#include <boost/utility/enable_if.hpp> | |
#include <boost/typeof/std/complex.hpp> | |
#include <boost/accumulators/numeric/functional_fwd.hpp> | |
namespace boost { namespace numeric { namespace operators | |
{ | |
// So that the stats compile when Sample type is std::complex | |
template<typename T, typename U> | |
typename | |
disable_if< | |
mpl::or_<is_same<T, U>, is_same<std::complex<T>, U> > | |
, std::complex<T> | |
>::type | |
operator *(std::complex<T> ri, U const &u) | |
{ | |
// BUGBUG promote result to typeof(T()*u) ? | |
return ri *= static_cast<T>(u); | |
} | |
template<typename T, typename U> | |
typename | |
disable_if< | |
mpl::or_<is_same<T, U>, is_same<std::complex<T>, U> > | |
, std::complex<T> | |
>::type | |
operator /(std::complex<T> ri, U const &u) | |
{ | |
// BUGBUG promote result to typeof(T()*u) ? | |
return ri /= static_cast<T>(u); | |
} | |
}}} // namespace boost::numeric::operators | |
namespace boost { namespace numeric | |
{ | |
namespace detail | |
{ | |
template<typename T> | |
struct one_complex | |
{ | |
static std::complex<T> const value; | |
}; | |
template<typename T> | |
std::complex<T> const one_complex<T>::value | |
= std::complex<T>(numeric::one<T>::value, numeric::one<T>::value); | |
} | |
/// INTERNAL ONLY | |
/// | |
template<typename T> | |
struct one<std::complex<T> > | |
: detail::one_complex<T> | |
{ | |
typedef one type; | |
typedef std::complex<T> value_type; | |
operator value_type const & () const | |
{ | |
return detail::one_complex<T>::value; | |
} | |
}; | |
}} // namespace boost::numeric | |
#endif |