X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=src%2Fdist_dna.c;h=b8797e955b0647d88ed547e68c22119c9a0112fd;hb=f295ab19440298e543db5a270e54f10a84382197;hp=a289f159e5bafad3c353845b7a364c3c99d3c6e2;hpb=82dd3702485179ba5408f1e3e57eb856d025e16c;p=ape.git diff --git a/src/dist_dna.c b/src/dist_dna.c index a289f15..b8797e9 100644 --- a/src/dist_dna.c +++ b/src/dist_dna.c @@ -1,4 +1,4 @@ -/* dist_dna.c 2010-07-14 */ +/* dist_dna.c 2010-09-16 */ /* Copyright 2005-2010 Emmanuel Paradis @@ -742,7 +742,7 @@ void distDNA_GG95_pairdel(unsigned char *x, int *n, int *s, double *d, #define COMPUTE_DIST_LogDet\ for (k = 0; k < 16; k++) Ftab[k] = ((double) Ntab[k]/L);\ - d[target] = (-log(detFourByFour(Ftab))/4 - LN4)/4;\ + d[target] = -log(detFourByFour(Ftab))/4 - LN4;\ if (*variance) {\ /* For the inversion, we first make U an identity matrix */\ for (k = 1; k < 15; k++) U[k] = 0;\