Averager3.h

#ifndef _theplu_yat_statistics_averager3
#define _theplu_yat_statistics_averager3

// $Id: Averager3.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 Averager3 : public averager_base3<Averager3>
  {
  public:
    Averager3(void);

    template<class Derived>
    const Averager3& operator+=(const averager_base3<Derived>& rhs);
  private:
    friend class averager_base<Averager3>;
    // 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 Averager3& Averager3::operator+=(const averager_base3<Derived>& rhs)
  {
    this->add3(rhs.mean(), rhs.sum_xx_centered(), rhs.sum_x3_centered(),
               rhs.n());
    return *this;
  }

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