source: trunk/test/phred.cc @ 3325

// $Id: phred.cc 3325 2014-10-07 01:44:54Z peter $

/*
  Copyright (C) 2014 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 <config.h>

#include "Suite.h"

#include <yat/omic/phred.h>

#include <utility>
#include <vector>

using namespace theplu::yat;

int main(int argc, char* argv[])
{
  test::Suite suite(argc, argv);

  std::vector<std::pair<double, double> > data;
  data.push_back(std::make_pair(0.1, 10.0));
  data.push_back(std::make_pair(0.01, 20.0));
  data.push_back(std::make_pair(0.0001, 40.0));


  for (size_t i=0; i<data.size(); ++i) {
    double phred = omic::phred(data[i].first);
    if (!suite.equal(phred, data[i].second))
      suite.add(false);
    double x = omic::phred_inv(data[i].second);
    if (!suite.equal(x, data[i].first))
      suite.add(false);
  }

  return suite.return_value();
}