ffpopsim_generic.h

Go to the documentation of this file.

#ifndef FFPOPGEN_GENERIC_H_
#define FFPOPGEN_GENERIC_H_

#include <time.h>
#include <cmath>
#include <vector>
#include <bitset>
#include <string>
#include <sstream>
#include <iostream>
#include <fstream>
#include <iomanip>
#include <gsl/gsl_sf.h>
#include <gsl/gsl_randist.h>
#include <gsl/gsl_rng.h>
#include <gsl/gsl_histogram.h>
#include <gsl/gsl_histogram2d.h>
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_vector.h>
#include <gsl/gsl_eigen.h>
#include <boost/dynamic_bitset.hpp>

#define MIN(a,b) (a<b)?a:b
#define MAX(a,b) (a>b)?a:b
#define RNG gsl_rng_taus2               //choose the random number generator algorithm, see http://www.gnu.org/software/gsl/manual/html_node/Random-number-generator-algorithms.html

#define FREE_RECOMBINATION 1
#define CROSSOVERS 2
#define SINGLE_CROSSOVER 3

using namespace std;

struct index_value_pair_t {
        size_t index;
        double val;
        index_value_pair_t(int index_in=0, double val_in=0) : index(index_in), val(val_in) {};
};

struct genotype_value_pair_t {
        boost::dynamic_bitset<> genotype;
        double val;
        genotype_value_pair_t(boost::dynamic_bitset<> genotype_in=boost::dynamic_bitset<>(0), double val_in=0) : genotype(genotype_in), val(val_in) {};
};

struct stat_t {
        double mean;
        double variance;
        stat_t(double mean_in=0, double variance_in=0) : mean(mean_in), variance(variance_in) {};
};

#define SAMPLE_ERROR -12312154

class sample {
public:
        int number_of_values;
        double *values;
        double mean;
        double variance;

        int bins;
        bool mem_dis;
        bool mem_values;
        gsl_histogram *distribution;

        bool with_range;
        double range_min;
        double range_max;

        sample();
        virtual ~sample();
        int set_up(int n);
        int set_distribution(int bins=100);
        void set_range(double min, double max) {range_min=min; range_max=max; with_range=true;}
        int calc_mean();
        int calc_variance();
        int calc_distribution();
        int print_distribution(ostream &out);
};

#endif /* FFPOPGEN_GENERIC_H_ */