/*============================================================================= Copyright (c) 2011-2013 Bolero MURAKAMI https://github.com/bolero-MURAKAMI/Sprout 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 SPROUT_NUMERIC_DFT_FIXED_SQUARE_HPP #define SPROUT_NUMERIC_DFT_FIXED_SQUARE_HPP #include <sprout/config.hpp> #include <sprout/index_tuple/metafunction.hpp> #include <sprout/container/traits.hpp> #include <sprout/container/functions.hpp> #include <sprout/container/indexes.hpp> #include <sprout/algorithm/fixed/result_of.hpp> #include <sprout/pit/pit.hpp> #include <sprout/math/comparison.hpp> #include <sprout/numeric/dft/fixed/sawtooth.hpp> #include <sprout/detail/container_complate.hpp> namespace sprout { namespace fixed { namespace detail { template<typename T> inline SPROUT_CONSTEXPR T square_value(T const& t, T const& d) { return sprout::fixed::detail::sawtooth_value(t) - sprout::fixed::detail::sawtooth_value(t - d); } template<typename Container, sprout::index_t... Indexes> inline SPROUT_CONSTEXPR typename sprout::fixed::result_of::algorithm<Container>::type square_impl( Container const& cont, typename sprout::container_traits<Container>::value_type const& frequency, typename sprout::container_traits<Container>::value_type const& amplitude, typename sprout::container_traits<Container>::value_type const& phase, typename sprout::container_traits<Container>::value_type const& duty, sprout::index_tuple<Indexes...>, typename sprout::container_traits<Container>::difference_type offset, typename sprout::container_traits<Container>::size_type size ) { typedef typename sprout::container_traits<Container>::value_type value_type; return sprout::remake<Container>( cont, size, (Indexes >= offset && Indexes < offset + size ? amplitude * sprout::fixed::detail::square_value(frequency * value_type(Indexes) + phase, duty) : *sprout::next(sprout::internal_begin(cont), Indexes) )... ); } template<typename Container> inline SPROUT_CONSTEXPR typename sprout::fixed::result_of::algorithm<Container>::type square( Container const& cont, typename sprout::container_traits<Container>::value_type const& frequency, typename sprout::container_traits<Container>::value_type const& amplitude, typename sprout::container_traits<Container>::value_type const& phase, typename sprout::container_traits<Container>::value_type const& duty ) { return sprout::fixed::detail::square_impl( cont, frequency, amplitude, phase, duty, sprout::container_indexes<Container>::make(), sprout::internal_begin_offset(cont), sprout::size(cont) ); } } // namespace detail // // square // template<typename Container> inline SPROUT_CONSTEXPR typename sprout::fixed::result_of::algorithm<Container>::type square( Container const& cont, typename sprout::container_traits<Container>::value_type const& frequency = 1, typename sprout::container_traits<Container>::value_type const& amplitude = 1, typename sprout::container_traits<Container>::value_type const& phase = 0, typename sprout::container_traits<Container>::value_type const& duty = 0.5 ) { return sprout::fixed::detail::square(cont, frequency, amplitude, phase, duty); } template<typename Container> inline SPROUT_CONSTEXPR typename sprout::fixed::result_of::algorithm<Container>::type square( typename sprout::container_traits<Container>::value_type const& frequency = 1, typename sprout::container_traits<Container>::value_type const& amplitude = 1, typename sprout::container_traits<Container>::value_type const& phase = 0, typename sprout::container_traits<Container>::value_type const& duty = 0.5 ) { return sprout::fixed::square(sprout::pit<Container>(), frequency, amplitude, phase, duty); } } // namespace fixed using sprout::fixed::square; } // namespace sprout #endif // #ifndef SPROUT_NUMERIC_DFT_FIXED_SQUARE_HPP