source: trunk/lib/gslapi/utility.cc @ 298

// $Id: utility.cc 298 2005-04-29 12:16:21Z peter $

#include <c++_tools/gslapi/utility.h>

#include <c++_tools/gslapi/vector.h>
#include <c++_tools/utility/random_singleton.h>

#include <utility>


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