summaryrefslogtreecommitdiff log msg author committer range
blob: 669a2c87ae327448318af59f6be7923e7f5b49c7 (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 70 71 72 73 74 75 76 77 78 79 80 81 82 ``` ``````// 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.25) return cos(constants::pi() * x); if(x < 0) { 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; if(rem > 0.25f) { rem = 0.5f - rem; rem = sin(constants::pi() * rem); } else 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&) { typedef typename tools::promote_args::type result_type; typedef typename policies::evaluation::type value_type; typedef typename policies::normalise< Policy, policies::promote_float, policies::promote_double, policies::discrete_quantile<>, policies::assert_undefined<> >::type forwarding_policy; return policies::checked_narrowing_cast(boost::math::detail::cos_pi_imp(x, forwarding_policy()), "cos_pi"); } template inline typename tools::promote_args::type cos_pi(T x) { return boost::math::cos_pi(x, policies::policy<>()); } } // namespace math } // namespace boost #endif ``````