Random.hpp

Go to the documentation of this file.

//
// This file is part of MorphoGraphX - http://www.MorphoGraphX.org
// Copyright (C) 2012-2015 Richard S. Smith and collaborators.
//
// If you use MorphoGraphX in your work, please cite:
//   http://dx.doi.org/10.7554/eLife.05864
//
// MorphoGraphX is free software, and is licensed under under the terms of the 
// GNU General (GPL) Public License version 2.0, http://www.gnu.org/licenses.
// 
#ifndef RANDOM_HPP
#define RANDOM_HPP
 
#include <Config.hpp>
 
#include <Vector.hpp>
 
#ifdef WIN32
static long int random() { return rand(); }
 
static void srandom(unsigned int seed) { return srand(seed); }
#endif
 
namespace mgx 
{
  mgx_EXPORT unsigned int sran_time();
  
  mgx_EXPORT void sran(unsigned int seed);
  
  mgx_EXPORT double ran(double M);
  mgx_EXPORT long double ran(long double M);
  mgx_EXPORT float ran(float M);
  
  template <typename T> double ran(T M) {
    return ran(double(M));
  }
  
  template <size_t dim, typename T> Vector<dim, T> ran(const Vector<dim, T>& V)
  {
    Vector<dim, T> result;
    for(size_t i = 0; i < dim; ++i)
      result[i] = ran(V[i]);
    return result;
  }
  
  mgx_EXPORT double gaussRan(double mean, double sigma);
  
  template <size_t dim> Vector<dim, double> gaussRan(const Vector<dim, double>& mean, const Vector<dim, double>& sigma)
  {
    Vector<dim, double> result;
    for(size_t i = 0; i < dim; ++i)
      result[i] = gaussRan(mean[i], sigma[i]);
    return result;
  }
  
  mgx_EXPORT long int ranInt();
  
  template <typename T> T ranInt()
  {
    long int i = ranInt();
    return (T)i;
  }
  
  template <size_t dim, typename T> Vector<dim, T> ranInt()
  {
    Vector<dim, long int> result;
    for(size_t i = 0; i < dim; ++i)
      result[i] = ranInt<T>();
    return result;
  }
  
  template <size_t dim> Vector<dim, long int> ranInt() {
    return ranInt<dim, long int>();
  }
  
  mgx_EXPORT long int ranInt(long int n);
  
  template <size_t dim, typename T> Vector<dim, T> ranInt(const Vector<dim, T>& n)
  {
    Vector<dim, T> result;
    for(size_t i = 0; i < dim; ++i)
      result[i] = ranInt<T>(n[i]);
    return result;
  }
}
#endif