Smoother.h

#ifndef _theplu_yat_statistics_smoother_
#define _theplu_yat_statistics_smoother_

// $Id: Smoother.h 3550 2017-01-03 05:41:02Z peter $

/*
  Copyright (C) 2008 Jari Häkkinen, Peter Johansson
  Copyright (C) 2009, 2010, 2016 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 "yat/utility/concept_check.h"
#include "yat/utility/iterator_traits.h"

#include <boost/concept_check.hpp>
#include <boost/iterator/iterator_concepts.hpp>

#include <vector>

namespace theplu {
namespace yat {
  namespace regression {
    class Kernel;
  }
namespace statistics {

  class Smoother
  {
  public:
    Smoother(const regression::Kernel&, double width,
             const std::vector<double>& values);

    Smoother(const regression::Kernel& kernel, double width,
             double xmin, double xmax, size_t n);

    void add(double x, double weight=1.0);

    const std::vector<double>& density(void) const;

    void rescale(double factor);

    void reset(void);

    const std::vector<double>& value(void) const;

  private:
    std::vector<double> density_;
    const regression::Kernel& kernel_;
    double width_;
    std::vector<double> x_;
  };

  template<typename Iterator>
  void add(Smoother& h, Iterator first, Iterator last)
  {
    BOOST_CONCEPT_ASSERT((utility::DataIteratorConcept<Iterator>));
    BOOST_CONCEPT_ASSERT((boost_concepts::SinglePassIterator<Iterator>));
    utility::iterator_traits<Iterator> traits;
    while (first!=last) {
      h.add(traits.data(first), traits.weight(first));
      ++first;
    }
  }

  std::ostream& operator<<(std::ostream& s,const Smoother&);

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

#endif