2012-04-10 03:54:38 +00:00
|
|
|
#ifndef SPROUT_NUMERIC_DFT_FIXED_DFT_HPP
|
|
|
|
#define SPROUT_NUMERIC_DFT_FIXED_DFT_HPP
|
2012-02-25 14:59:46 +00:00
|
|
|
|
|
|
|
#include <sprout/config.hpp>
|
|
|
|
#include <sprout/index_tuple.hpp>
|
2012-03-31 07:24:13 +00:00
|
|
|
#include <sprout/container/traits.hpp>
|
|
|
|
#include <sprout/container/functions.hpp>
|
2012-02-25 14:59:46 +00:00
|
|
|
#include <sprout/iterator/operation.hpp>
|
|
|
|
#include <sprout/algorithm/fixed/result_of.hpp>
|
2012-04-10 03:54:38 +00:00
|
|
|
#include <sprout/numeric/dft/dft_element.hpp>
|
2012-04-01 13:15:09 +00:00
|
|
|
#include HDR_ITERATOR_SSCRISK_CEL_OR_SPROUT
|
2012-02-25 14:59:46 +00:00
|
|
|
|
|
|
|
namespace sprout {
|
|
|
|
namespace fixed {
|
|
|
|
namespace detail {
|
2012-02-28 01:46:39 +00:00
|
|
|
template<typename InputIterator, typename Result, sprout::index_t... Indexes>
|
2012-04-04 13:23:41 +00:00
|
|
|
inline SPROUT_CONSTEXPR typename sprout::fixed::result_of::algorithm<Result>::type dft_impl(
|
2012-02-25 14:59:46 +00:00
|
|
|
InputIterator first,
|
|
|
|
InputIterator last,
|
|
|
|
Result const& result,
|
|
|
|
sprout::index_tuple<Indexes...>,
|
2012-03-31 07:24:13 +00:00
|
|
|
typename sprout::container_traits<Result>::difference_type offset,
|
|
|
|
typename sprout::container_traits<Result>::size_type size,
|
|
|
|
typename sprout::container_traits<Result>::size_type input_size
|
2012-02-25 14:59:46 +00:00
|
|
|
)
|
|
|
|
{
|
2012-03-31 07:24:13 +00:00
|
|
|
return sprout::remake<Result>(
|
2012-02-25 14:59:46 +00:00
|
|
|
result,
|
2012-04-11 14:28:29 +00:00
|
|
|
size,
|
2012-02-25 14:59:46 +00:00
|
|
|
(Indexes >= offset && Indexes < offset + size && Indexes < offset + input_size
|
|
|
|
? sprout::detail::dft_element_impl(first, last, Indexes - offset, input_size)
|
2012-03-31 07:24:13 +00:00
|
|
|
: *sprout::next(sprout::internal_begin(result), Indexes)
|
2012-02-25 14:59:46 +00:00
|
|
|
)...
|
|
|
|
);
|
|
|
|
}
|
|
|
|
template<typename InputIterator, typename Result>
|
2012-04-04 13:23:41 +00:00
|
|
|
inline SPROUT_CONSTEXPR typename sprout::fixed::result_of::algorithm<Result>::type dft(
|
2012-02-25 14:59:46 +00:00
|
|
|
InputIterator first,
|
|
|
|
InputIterator last,
|
|
|
|
Result const& result
|
|
|
|
)
|
|
|
|
{
|
|
|
|
return sprout::fixed::detail::dft_impl(
|
|
|
|
first,
|
|
|
|
last,
|
|
|
|
result,
|
2012-03-31 07:24:13 +00:00
|
|
|
typename sprout::index_range<0, sprout::container_traits<Result>::static_size>::type(),
|
|
|
|
sprout::internal_begin_offset(result),
|
2012-02-25 14:59:46 +00:00
|
|
|
sprout::size(result),
|
2012-04-01 13:15:09 +00:00
|
|
|
NS_SSCRISK_CEL_OR_SPROUT::distance(first, last)
|
2012-02-25 14:59:46 +00:00
|
|
|
);
|
|
|
|
}
|
|
|
|
} // namespace detail
|
|
|
|
//
|
|
|
|
// dft
|
|
|
|
//
|
|
|
|
template<typename InputIterator, typename Result>
|
2012-04-04 13:23:41 +00:00
|
|
|
inline SPROUT_CONSTEXPR typename sprout::fixed::result_of::algorithm<Result>::type dft(
|
2012-02-25 14:59:46 +00:00
|
|
|
InputIterator first,
|
|
|
|
InputIterator last,
|
|
|
|
Result const& result
|
|
|
|
)
|
|
|
|
{
|
|
|
|
return sprout::fixed::detail::dft(first, last, result);
|
|
|
|
}
|
|
|
|
} // namespace fixed
|
|
|
|
|
|
|
|
using sprout::fixed::dft;
|
|
|
|
} // namespace sprout
|
|
|
|
|
2012-04-10 03:54:38 +00:00
|
|
|
#endif // #ifndef SPROUT_NUMERIC_DFT_FIXED_DFT_HPP
|