summaryrefslogtreecommitdiff log msg author committer range
blob: 18d06c00dff0e19cbfb165b4e4c905d518e791a8 (plain)
 ```1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 ``` ``````// Copyright (c) 2007 John Maddock // Use, modification and distribution are subject to 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_MATH_COS_PI_HPP #define BOOST_MATH_COS_PI_HPP #ifdef _MSC_VER #pragma once #endif #include #include #include #include #include #include namespace boost{ namespace math{ namespace detail{ template T cos_pi_imp(T x, const Policy& pol) { BOOST_MATH_STD_USING // ADL of std names // cos of pi*x: bool invert = false; if(fabs(x) < 0.5) return cos(constants::pi() * x); if(x < 1) { x = -x; } T rem = floor(x); if(itrunc(rem, pol) & 1) invert = !invert; rem = x - rem; if(rem > 0.5f) { rem = 1 - rem; invert = !invert; } if(rem == 0.5f) return 0; rem = cos(constants::pi() * rem); return invert ? T(-rem) : rem; } } // namespace detail template inline typename tools::promote_args::type cos_pi(T x, const Policy& pol) { typedef typename tools::promote_args::type result_type; return boost::math::detail::cos_pi_imp(x, pol); } template inline typename tools::promote_args::type cos_pi(T x) { return boost::math::cos_pi(x, policies::policy<>()); } } // namespace math } // namespace boost #endif ``````