Centralizer.h

#ifndef _theplu_yat_normalizer_centralizer_
#define _theplu_yat_normalizer_centralizer_

// $Id: Centralizer.h 2415 2011-01-22 06:01:31Z peter $

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

#include "yat/statistics/Average.h"
#include "yat/utility/DataIterator.h"
#include "yat/utility/iterator_traits.h"
#include "yat/utility/WeightIterator.h"

#include <boost/concept_check.hpp>

#include <algorithm>
#include <functional>

namespace theplu {
namespace yat {
namespace normalizer {

  template<class UnaryFunction = statistics::Average>
  class Centralizer
  {
  public:
    Centralizer(void){}

    Centralizer(const UnaryFunction& uf)
      : func_(uf) {}

    template<class InputIterator, class ForwardIterator>
    void operator()(InputIterator first, InputIterator last,
                    ForwardIterator result) const
    {
      BOOST_CONCEPT_ASSERT((boost::InputIterator<InputIterator>));
      BOOST_CONCEPT_ASSERT((boost::Mutable_ForwardIterator<ForwardIterator>));
      // copy weight if result is weighted
      detail::copy_weight_if_weighted(first, last, result); 
      std::transform(utility::data_iterator(first), 
                     utility::data_iterator(last), 
                     utility::data_iterator(result),
                     std::bind2nd(std::minus<double>(),func_(first, last)));
    }

  private:
    UnaryFunction func_;
  };

}}} // end of namespace normalizer, yat and thep
#endif