utilities/src/randomgenerator.cpp

/********************************************************************************
 *  FARSA Genetic Algorithm Library                                             *
 *  Copyright (C) 2007-2008 Gianluca Massera <emmegian@yahoo.it>                *
 *                                                                              *
 *  This program 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.                                         *
 *                                                                              *
 *  This program 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., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA  *
 ********************************************************************************/

#include "randomgenerator.h"
#include <cmath>

#ifdef FARSA_USE_GSL
    #include "gsl/gsl_rng.h"
    #include "gsl/gsl_randist.h"
#endif

namespace farsa {

class FARSA_UTIL_INTERNAL RandomGeneratorPrivate {
public:
#ifdef FARSA_USE_GSL
    gsl_rng* rng;
    RandomGeneratorPrivate() {
        rng = gsl_rng_alloc( gsl_rng_taus2 );
    };
    ~RandomGeneratorPrivate() {
        gsl_rng_free( rng );
    };
#endif
};

RandomGenerator::RandomGenerator( int seed ) {
    prive = new RandomGeneratorPrivate();
#ifdef FARSA_USE_GSL
    gsl_rng_set( prive->rng, seed );
#endif
}

RandomGenerator::~RandomGenerator() {
    delete prive;
}

void RandomGenerator::setSeed( int seed ) {
    this->seedv = seed;
#ifdef FARSA_USE_GSL
    gsl_rng_set( prive->rng, seed );
#endif
}

int RandomGenerator::seed() {
    return seedv;
}

bool RandomGenerator::getBool( double trueProbability ) {
#ifdef FARSA_USE_GSL
    return ( gsl_rng_uniform( prive->rng ) < trueProbability );
#endif
}

int RandomGenerator::getInt( int min, int max ) {
#ifdef FARSA_USE_GSL
    return ( gsl_rng_uniform_int( prive->rng, qAbs( max-min )+1 ) + min );
#endif
}

double RandomGenerator::getDouble( double min, double max ) {
#ifdef FARSA_USE_GSL
    //--- FIXME this implementation never return max
    return gsl_ran_flat( prive->rng, min, max );
#endif
}

double RandomGenerator::getGaussian( double var, double mean ) {
#ifdef FARSA_USE_GSL
    return gsl_ran_gaussian( prive->rng, var ) + mean;
#endif
}

} // end namespace farsa