asdf.hpp

Go to the documentation of this file.

#ifndef EDSP_ASDF_HPP
#define EDSP_ASDF_HPP

#include <iterator>

namespace edsp { namespace feature { inline namespace temporal {

        template <typename InputIt, typename OutputIt>
        constexpr void amdf(InputIt first, InputIt last, OutputIt d_first) {
            using value_type = typename std::iterator_traits<OutputIt>::value_type;
            const auto N      = std::distance(first, last);
            auto * array = &(*first);
            for (auto i = 0; i < N; ++i, ++d_first) {
                for (auto j = 0; j < (N - i); ++j) {
                    *d_first += std::abs(array[j] - array[j + i]);
                }
                *d_first /= static_cast<value_type>(N);
            }
        }
}}}
#endif //EDSP_ASDF_HPP