Chris@16: // Boost.Units - A C++ library for zero-overhead dimensional analysis and 
Chris@16: // unit/quantity manipulation and conversion
Chris@16: //
Chris@16: // Copyright (C) 2003-2008 Matthias Christian Schabel
Chris@16: // Copyright (C) 2008 Steven Watanabe
Chris@16: //
Chris@16: // Distributed under the Boost Software License, Version 1.0. (See
Chris@16: // accompanying file LICENSE_1_0.txt or copy at
Chris@16: // http://www.boost.org/LICENSE_1_0.txt)
Chris@16: 
Chris@16: #ifndef BOOST_UNITS_CMATH_IMPL_HPP 
Chris@16: #define BOOST_UNITS_CMATH_IMPL_HPP
Chris@16: 
Chris@16: #include <boost/config.hpp>
Chris@16: #include <boost/math/special_functions/fpclassify.hpp>
Chris@16: 
Chris@16: namespace boost {
Chris@16: namespace units {
Chris@16: namespace detail {
Chris@16: 
Chris@16: template<class Y> 
Chris@16: inline bool isgreater BOOST_PREVENT_MACRO_SUBSTITUTION(const Y& v1,const Y& v2)
Chris@16: {
Chris@16:     if((boost::math::isnan)(v1) || (boost::math::isnan)(v2)) return false;
Chris@16:     else return v1 > v2;
Chris@16: }
Chris@16: 
Chris@16: template<class Y> 
Chris@16: inline bool isgreaterequal BOOST_PREVENT_MACRO_SUBSTITUTION(const Y& v1,const Y& v2)
Chris@16: {
Chris@16:     if((boost::math::isnan)(v1) || (boost::math::isnan)(v2)) return false;
Chris@16:     else return v1 >= v2;
Chris@16: }
Chris@16: 
Chris@16: template<class Y> 
Chris@16: inline bool isless BOOST_PREVENT_MACRO_SUBSTITUTION(const Y& v1,const Y& v2)
Chris@16: {
Chris@16:     if((boost::math::isnan)(v1) || (boost::math::isnan)(v2)) return false;
Chris@16:     else return v1 < v2;
Chris@16: }
Chris@16: 
Chris@16: template<class Y> 
Chris@16: inline bool islessequal BOOST_PREVENT_MACRO_SUBSTITUTION(const Y& v1,const Y& v2)
Chris@16: {
Chris@16:     if((boost::math::isnan)(v1) || (boost::math::isnan)(v2)) return false;
Chris@16:     else return v1 <= v2;
Chris@16: }
Chris@16: 
Chris@16: template<class Y> 
Chris@16: inline bool islessgreater BOOST_PREVENT_MACRO_SUBSTITUTION(const Y& v1,const Y& v2)
Chris@16: {
Chris@16:     if((boost::math::isnan)(v1) || (boost::math::isnan)(v2)) return false;
Chris@16:     else return v1 < v2 || v1 > v2;
Chris@16: }
Chris@16: 
Chris@16: template<class Y> 
Chris@16: inline bool isunordered BOOST_PREVENT_MACRO_SUBSTITUTION(const Y& v1,const Y& v2)
Chris@16: {
Chris@16:     return (boost::math::isnan)(v1) || (boost::math::isnan)(v2);
Chris@16: }
Chris@16: 
Chris@16: template<class Y>
Chris@16: inline Y fdim BOOST_PREVENT_MACRO_SUBSTITUTION(const Y& v1,const Y& v2)
Chris@16: {
Chris@16:     if((boost::math::isnan)(v1)) return v1;
Chris@16:     else if((boost::math::isnan)(v2)) return v2;
Chris@16:     else if(v1 > v2) return(v1 - v2);
Chris@16:     else return(Y(0));
Chris@16: }
Chris@16: 
Chris@16: #if 0
Chris@16: 
Chris@16: template<class T>
Chris@16: struct fma_issue_warning {
Chris@16:     enum { value = false };
Chris@16: };
Chris@16: 
Chris@16: template<class Y>
Chris@16: inline Y fma(const Y& v1,const Y& v2,const Y& v3)
Chris@16: {
Chris@16:     //this implementation does *not* meet the
Chris@16:     //requirement of infinite intermediate precision
Chris@16:     BOOST_STATIC_WARNING((fma_issue_warning<Y>::value));
Chris@16: 
Chris@16:     return v1 * v2 + v3;
Chris@16: }
Chris@16: 
Chris@16: #endif
Chris@16: 
Chris@16: template<class Y>
Chris@16: inline Y fmax BOOST_PREVENT_MACRO_SUBSTITUTION(const Y& v1,const Y& v2)
Chris@16: {
Chris@16:     if((boost::math::isnan)(v1)) return(v2);
Chris@16:     else if((boost::math::isnan)(v2)) return(v1);
Chris@16:     else if(v1 > v2) return(v1);
Chris@16:     else return(v2);
Chris@16: }
Chris@16: 
Chris@16: template<class Y>
Chris@16: inline Y fmin BOOST_PREVENT_MACRO_SUBSTITUTION(const Y& v1,const Y& v2)
Chris@16: {
Chris@16:     if((boost::math::isnan)(v1)) return(v2);
Chris@16:     else if((boost::math::isnan)(v2)) return(v1);
Chris@16:     else if(v1 < v2) return(v1);
Chris@16:     else return(v2);
Chris@16: }
Chris@16: 
Chris@16: //template<class Y>
Chris@16: //inline long long llrint(const Y& val)
Chris@16: //{
Chris@16: //    return static_cast<long long>(rint(val));
Chris@16: //}
Chris@16: //
Chris@16: //template<class Y>
Chris@16: //inline long long llround(const Y& val)
Chris@16: //{
Chris@16: //    return static_cast<long long>(round(val));
Chris@16: //}
Chris@16: 
Chris@16: #if 0
Chris@16: 
Chris@16: template<class Y>
Chris@16: inline Y nearbyint(const Y& val)
Chris@16: {
Chris@16:     //this is not really correct.
Chris@16:     //the result should be according to the
Chris@16:     //current rounding mode.
Chris@16:     using boost::math::round;
Chris@16:     return round(val);
Chris@16: }
Chris@16: 
Chris@16: template<class Y>
Chris@16: inline Y rint(const Y& val)
Chris@16: {
Chris@16:     //I don't feel like trying to figure out
Chris@16:     //how to raise a floating pointer exception
Chris@16:     return nearbyint(val);
Chris@16: }
Chris@16: 
Chris@16: #endif
Chris@16: 
Chris@16: template<class Y>
Chris@16: inline Y trunc BOOST_PREVENT_MACRO_SUBSTITUTION(const Y& val)
Chris@16: {
Chris@16:     if(val > 0) return std::floor(val);
Chris@16:     else if(val < 0) return std::ceil(val);
Chris@16:     else return val;
Chris@16: }
Chris@16: 
Chris@16: }
Chris@16: }
Chris@16: }
Chris@16: 
Chris@16: #endif // BOOST_UNITS_CMATH_IMPL_HPP