source: trunk/src/gslapi_utility.cc @ 259

// $Id: gslapi_utility.cc 259 2005-03-04 01:09:55Z jari $

#include "gslapi_utility.h"
#include <utility>

#include "random_singleton.h"
#include "vector.h"


namespace theplu {
namespace gslapi {


  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);
    theplu::cpptools::random_singleton* 
      rnd=theplu::cpptools::random_singleton::get_instance();
    for (size_t i=0; i<vec.size(); i++){
      size_t index = rnd->get_uniform_int(vec.size()-i);
      invec[i]=vec(index);
      vec(index)=vec(vec.size()-i-1);
    }
  }

}} // of namespace gslapi and namespace thep