X-Git-Url: https://git.donarmstrong.com/?p=rsem.git;a=blobdiff_plain;f=sam%2Fkprobaln.c;fp=sam%2Fkprobaln.c;h=894a2ae625f9e0b2b9d38337b8651d8f02bcce4a;hp=5201c1af75b9c24636ee67481b0461ab7713e7a1;hb=946f9a6adb2a82048c8453d44693cd3838d32939;hpb=9712cd844059fd5cd25b52162ae10bc1d756d899 diff --git a/sam/kprobaln.c b/sam/kprobaln.c index 5201c1a..894a2ae 100644 --- a/sam/kprobaln.c +++ b/sam/kprobaln.c @@ -161,7 +161,7 @@ int kpa_glocal(const uint8_t *_ref, int l_ref, const uint8_t *_query, int l_quer double p = 1., Pr1 = 0.; for (i = 0; i <= l_query + 1; ++i) { p *= s[i]; - if (p < 1e-100) Pr += -4.343 * log(p), p = 1.; + if (p < 1e-100) Pr1 += -4.343 * log(p), p = 1.; } Pr1 += -4.343 * log(p * l_ref * l_query); Pr = (int)(Pr1 + .499);