Averager4.h

#ifndef _theplu_yat_statistics_averager4
#define _theplu_yat_statistics_averager4

// $Id: Averager4.h 2809 2012-08-06 13:44:47Z peter $

/*
  Copyright (C) 2012 Peter Johansson

  This file is part of the yat library, http://dev.thep.lu.se/yat

  The yat library 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; either version 3 of the
  License, or (at your option) any later version.

  The yat library 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 for more details.

  You should have received a copy of the GNU General Public License
  along with yat. If not, see <http://www.gnu.org/licenses/>.
*/

#include "averager_base.h"

namespace theplu {
namespace yat {
namespace statistics {

  class Averager4 : public averager_base4<Averager4>
  {
  public:
    Averager4(void);

    //const Averager& averager(void) const;

    template<class Derived>
    const Averager4& operator+=(const averager_base4<Derived>& rhs);
  private:
    friend class averager_base<Averager4>;
    // static virtual function in averager_base
    void add_impl(double, long int n);
    void rescale_impl(double);

    // using compiler generated copy and assignment
    // Averager4(const Averager4&)
    // Averager4& operator=(const Averager4& rhs);
  };

    template<class Derived>
    const Averager4& Averager4::operator+=(const averager_base4<Derived>& rhs)
    {
      this->add4(rhs.mean(), rhs.sum_xx_centered(),
                 rhs.sum_x3_centered(), rhs.sum_x4_centered(), rhs.n());
      return *this;
    }

}}} // of namespace statistics, yat, and theplu
#endif