Gauss.h

#ifndef _theplu_yat_normalizer_gauss_
#define _theplu_yat_normalizer_gauss_

// $Id: Gauss.h 2263 2010-05-26 22:13:07Z peter $

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

#include <boost/concept_check.hpp>

#include <gsl/gsl_cdf.h>

#include <iterator>

namespace theplu {
namespace yat {
namespace normalizer {

  class Gauss
  {
  public:
    template<typename RandomAccessIter1, typename RandomAccessIter2>
    void operator()(RandomAccessIter1 first, RandomAccessIter1 last,
                    RandomAccessIter2 result) const
    {
      BOOST_CONCEPT_ASSERT((utility::DataIteratorConcept<RandomAccessIter1>));
      BOOST_CONCEPT_ASSERT((utility::DataIteratorConcept<RandomAccessIter2>));
      BOOST_CONCEPT_ASSERT((boost::RandomAccessIterator<RandomAccessIter1>));
      BOOST_CONCEPT_ASSERT((boost::Mutable_RandomAccessIterator<RandomAccessIter2>));
      Spearman spearman;
      spearman(first, last, result);
      RandomAccessIter2 end = result + std::distance(first, last);
      utility::iterator_traits<RandomAccessIter2> trait;
      while (result != end) {
        trait.data(result) = gsl_cdf_ugaussian_Pinv(trait.data(result));
        ++result;
      }
    }


  private:
  };

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