Skip to content

Commit bd0deed

Browse files
Do not weight recombination by n in LS HMM
1 parent 8bd218d commit bd0deed

File tree

2 files changed

+3
-2
lines changed

2 files changed

+3
-2
lines changed

Makefile

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -7,7 +7,7 @@ CFLAGS=-std=c99 -g -O3 -march=native -funroll-loops -ffast-math \
77
all: _tsinfer.cpython-34m.so
88

99
_tsinfer.cpython-34m.so: _tsinfermodule.c
10-
CC="${CC}" CFLAGS="${CFLAGS}" python setup.py build_ext --inplace
10+
CC="${CC}" CFLAGS="${CFLAGS}" python3 setup.py build_ext --inplace
1111

1212
ctags:
1313
ctags lib/*.c lib/*.h tsinfer/*.py

lib/ancestor_matcher.c

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -297,7 +297,8 @@ ancestor_matcher_update_site_likelihood_values(ancestor_matcher_t *self,
297297
double max_L, p_last, p_no_recomb, p_recomb, p_t, p_e;
298298
const double rho = self->recombination_rate[site];
299299
const double mu = self->mismatch_rate[site];
300-
const double n = (double) self->tree_sequence_builder->num_match_nodes;
300+
/* const double n = (double) self->tree_sequence_builder->num_match_nodes; */
301+
const double n = 1;
301302
const double num_alleles
302303
= (double) self->tree_sequence_builder->sites.num_alleles[site];
303304

0 commit comments

Comments
 (0)