source: trunk/yat/statistics/VectorFunction.cc @ 1703

// $Id: VectorFunction.cc 1487 2008-09-10 08:41:36Z jari $

/*
  Copyright (C) 2007 Jari Häkkinen, Peter Johansson
  Copyright (C) 2008 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 "VectorFunction.h"
#include "utility.h"

#include <algorithm>
#include <cassert>
#include <numeric>
#include <vector>

namespace theplu {
namespace yat {
namespace statistics {  

  VectorFunction::~VectorFunction(){}


  double Max::operator()(const std::vector<double>& vec) const
  {
    assert(vec.size());
    return *(std::max_element(vec.begin(), vec.end()));
  }
  

  double Mean::operator()(const std::vector<double>& vec) const
  {
    assert(vec.size());
    return std::accumulate(vec.begin(), vec.end(), 0.0)/vec.size();
  }
  

  double Median::operator()(const std::vector<double>& vec) const
  {
    assert(vec.size());
    return median(vec.begin(), vec.end());
  }
  

  double Min::operator()(const std::vector<double>& vec) const
  {
    assert(vec.size());
    return *(std::min_element(vec.begin(), vec.end()));
  }
  

  Nth_Element::Nth_Element(size_t N)
    : N_(N) {}


  double Nth_Element::operator()(const std::vector<double>& vec) const
  {
    assert(N_-1<vec.size());
    std::vector<double> vec_copy(vec);
    std::vector<double>::iterator nth = vec_copy.begin()+(N_-1);
    std::nth_element(vec_copy.begin(), nth, vec_copy.end());
    return *(nth);
  }
  



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