///////////////////////////////////////////////////////////////////////////////////////
// SDRdaemon - send I/Q samples read from a SDR device over the network via UDP      //
//             with FEC protection. GNUradio interface.                              //
//                                                                                   //
// http://stackoverflow.com/questions/10990618/calculate-rolling-moving-average-in-c //
//                                                                                   //
// Copyright (C) 2016 Edouard Griffiths, F4EXB                                       //
//                                                                                   //
// This program is free software; you can redistribute it and/or modify              //
// it under the terms of the GNU General Public License as published by              //
// the Free Software Foundation as version 3 of the License, or                      //
//                                                                                   //
// This program is distributed in the hope that it will be useful,                   //
// but WITHOUT ANY WARRANTY; without even the implied warranty of                    //
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the                      //
// GNU General Public License V3 for more details.                                   //
//                                                                                   //
// You should have received a copy of the GNU General Public License                 //
// along with this program. If not, see <http://www.gnu.org/licenses/>.              //
///////////////////////////////////////////////////////////////////////////////////////

#ifndef _UTIL_MOVINGAVERAGE_H_
#define _UTIL_MOVINGAVERAGE_H_

#include <algorithm>

template <typename T, typename Total, int N>
class MovingAverage
{
  public:
    MovingAverage()
      : m_num_samples(0), m_total(0)
    { }

    void operator()(T sample)
    {
        if (m_num_samples < N)
        {
            m_samples[m_num_samples++] = sample;
            m_total += sample;
        }
        else
        {
            T& oldest = m_samples[m_num_samples++ % N];
            m_total += sample - oldest;
            oldest = sample;
        }
    }

    operator double() const { return m_num_samples > 0 ? m_total / std::min(m_num_samples, N) : 0.0d; }
    operator float() const { return m_num_samples > 0 ? m_total / std::min(m_num_samples, N) : 0.0f; }

  private:
    T m_samples[N];
    int m_num_samples;
    Total m_total;
};

#endif /* GR_SDRDAEMONFEC_LIB_MOVINGAVERAGE_H_ */