TFL2013sim/diploid_functions.hpp at master · ThorntonLab/TFL2013sim · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
#ifndef _DIPLOID_FUNCTIONS_HPP_
#define _DIPLOID_FUNCTIONS_HPP_

#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>

template< typename gamete_type,
	  typename vector_type_allocator,
	  template<typename,typename> class vector_type,
	  typename diploid_fitness_function>
double sample_diploid(gsl_rng * r, vector_type<gamete_type,vector_type_allocator > * gametes,
		      const unsigned & twoN, const diploid_fitness_function & ff,
		      const double & f = 0.);

template< typename gamete_type,
	  typename vector_type_allocator,
	  template<typename,typename> class vector_type,
	  typename diploid_fitness_function>
double sample_diploid_mc(gsl_rng * r, vector_type<gamete_type,vector_type_allocator > * gametes,
		      const unsigned & twoN, const diploid_fitness_function & ff,
		      const double & f = 0.);


template< typename gamete_type,
	  typename vector_type_allocator,
	  template<typename,typename> class vector_type,
	  typename diploid_fitness_function >
	  double sample_diploid(gsl_rng * r, vector_type<gamete_type,vector_type_allocator > * gametes,
				const unsigned & twoN_curr, const unsigned & twoN_next,
				const diploid_fitness_function & ff,
				const double & f = 0.);

template< typename gamete_type,
	  typename vector_type_allocator,
	  template<typename,typename> class vector_type,
	  typename diploid_fitness_function >
	  double sample_diploid_mc(gsl_rng * r, vector_type<gamete_type,vector_type_allocator > * gametes,
				const unsigned & twoN_curr, const unsigned & twoN_next,
				const diploid_fitness_function & ff,
				const double & f = 0.);


template< typename gamete_type,
	  typename vector_type_allocator,
 	  template<typename,typename> class vector_type,
	  typename map_function>
	  unsigned recombine(gsl_rng * r, vector_type<gamete_type,vector_type_allocator > * gametes,
			     const unsigned & twoN, const double & littler,const map_function & mf,
			     const bool single,
			     const double & f = 0);

#endif /* _DIPLOID_FUNCTIONS_HPP_ */
#include <diploid_functions.tcc>