summaryrefslogtreecommitdiff log msg author committer range
blob: ae6b3e74423345333c07074c2dd2b4cc4b225f7a (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 ``` ``````// 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_SIN_PI_HPP #define BOOST_MATH_SIN_PI_HPP #ifdef _MSC_VER #pragma once #endif #include #include #include #include #include #include namespace boost{ namespace math{ namespace detail{ template T sin_pi_imp(T x, const Policy& pol) { BOOST_MATH_STD_USING // ADL of std names if(x < 0) return -sin_pi(-x); // sin of pi*x: bool invert; if(x < 0.5) return sin(constants::pi() * x); if(x < 1) { invert = true; x = -x; } else invert = false; T rem = floor(x); if(itrunc(rem, pol) & 1) invert = !invert; rem = x - rem; if(rem > 0.5f) rem = 1 - rem; if(rem == 0.5f) return static_cast(invert ? -1 : 1); rem = sin(constants::pi() * rem); return invert ? T(-rem) : rem; } } // namespace detail template inline typename tools::promote_args::type sin_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::sin_pi_imp(x, forwarding_policy()), "cos_pi"); } template inline typename tools::promote_args::type sin_pi(T x) { return boost::math::sin_pi(x, policies::policy<>()); } } // namespace math } // namespace boost #endif ``````