fix bug of m2 & k2 in infer_allele

hxj5 · hxj5 · commit 06a8bf8b0898 · 2020-11-17T12:58:07.000+08:00
diff --git a/cellsnp_util.h b/cellsnp_util.h
@@ -852,7 +852,7 @@ static inline void csp_infer_allele(size_t *bc, int8_t *ref_idx, int8_t *alt_idx
     if (bc[0] < bc[1]) { m1 = bc[1]; m2 = bc[0]; k1 = 1; k2 = 0; }
     else { m1 = bc[0]; m2 = bc[1]; k1 = 0; k2 = 1; }
     for (i = 2; i < 5; i++) {
-        if (bc[i] > m1) { m1 = bc[i]; k1 = i; }
+        if (bc[i] > m1) { m2 = m1; k2 = k1; m1 = bc[i]; k1 = i; }
         else if (bc[i] > m2) { m2 = bc[i]; k2 = i; }
     }
     *ref_idx = k1; *alt_idx = k2;

Original file line number	Diff line number	Diff line change
`@@ -852,7 +852,7 @@ static inline void csp_infer_allele(size_t bc, int8_t ref_idx, int8_t *alt_idx`
`852`	`852`	`if (bc[0] < bc[1]) { m1 = bc[1]; m2 = bc[0]; k1 = 1; k2 = 0; }`
`853`	`853`	`else { m1 = bc[0]; m2 = bc[1]; k1 = 0; k2 = 1; }`
`854`	`854`	`for (i = 2; i < 5; i++) {`
`855`		`- if (bc[i] > m1) { m1 = bc[i]; k1 = i; }`
	`855`	`+ if (bc[i] > m1) { m2 = m1; k2 = k1; m1 = bc[i]; k1 = i; }`
`856`	`856`	`else if (bc[i] > m2) { m2 = bc[i]; k2 = i; }`
`857`	`857`	`}`
`858`	`858`	`ref_idx = k1; alt_idx = k2;`