18 #include "MapFunction.h" 19 #include "MathConstant.h" 23 double DistanceToRecombination(
double distance)
25 return (1.0 - exp(-2.0 * distance)) * 0.5;
28 double RecombinationToDistance(
double recombination)
30 return (log(max(1.0 - 2 * recombination, 1e-7)) * -0.5);
33 double KosambiDistanceToRecombination(
double distance)
35 double e_to_4x = exp(4.0 * distance);
37 return (0.5 *(e_to_4x - 1.0) / (e_to_4x + 1.0));
40 double RecombinationToKosambiDistance(
double theta)
42 return 0.25 * log((1.0 + 2*theta) / max(1.0 - 2.0*theta, 1e-7));