source: trunk/yat/utility/utility.cc @ 737

// $Id: utility.cc 737 2007-01-07 13:48:38Z peter $

/*
  Copyright (C) 2005 Jari Häkkinen, Markus Ringnér
  Copyright (C) 2006 Jari Häkkinen
  Copyright (C) 2007 Peter Johansson

  This file is part of the yat library, http://lev.thep.lu.se/trac/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 2 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 this program; if not, write to the Free Software
  Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA
  02111-1307, USA.
*/

#include "utility.h"

#include "matrix.h"
#include "stl_utility.h"
#include "vector.h"
#include "yat/random/random.h"

#include <sstream>
#include <string>

namespace theplu {
namespace yat {
namespace utility {

  bool is_double(const std::string& s)
  {
    std::stringstream ss(s);
    double a;
    ss>>a;
    if(ss.fail()) 
      return false;
    // Check that nothing is left on stream
    std::string b;
    ss >> b;
    return (b.size() ? false : true);
  }
  
  bool is_float(const std::string& s)
  {
    std::stringstream ss(s);
    float a;
    ss>>a;
    if(ss.fail()) 
      return false;
    // Check that nothing is left on stream
    std::string b;
    ss >> b;
    return (b.size() ? false : true);
  }

  bool is_int(const std::string& s)
  {
    std::stringstream ss(s);
    int a;
    ss >> a;
    if(ss.fail()) 
      return false;
    // Check that nothing is left on stream
    std::string b;
    ss >> b;
    return (b.size() ? false : true);
  }

  bool is_nan(const std::string& s)
  {
    std::stringstream ss(s);
    std::string s2;
    ss >> s2; // to trim surrounding whitespaces
    to_lower(s2);
    // Check that nothing is left on stream
    std::string s3;
    ss >> s3;
    if(s3.size())
      return false;
    std::string nan("nan");
    return (nan==s2);
  }


  std::pair<size_t, size_t> minmax_index(const vector& vec,
                                         const std::vector<size_t>& subset)
  {
    size_t min_index = subset[0];
    size_t max_index = subset[0];
    for (unsigned int i=1; i<subset.size(); i++)
      if (vec[subset[i]] < vec[min_index])
        min_index = subset[i];
      else if (vec[subset[i]] > vec[max_index])
        max_index = subset[i];
    return std::pair<size_t,size_t>(min_index, max_index);
  }


  void shuffle(vector& invec)
  {
    vector vec(invec);
    random::DiscreteUniform rnd;
    for (size_t i=0; i<vec.size(); i++){
      size_t index = rnd(vec.size()-i);
      invec[i]=vec(index);
      vec(index)=vec(vec.size()-i-1);
    }
  }


}}} // end of namespace utility, yat and thep