X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=bam_mcns.c;h=d337e5a290aa76a1f1faa2c8a687db4639962a16;hb=d742c10178d8981ce062cb1f993149cfe9876613;hp=d98d48030560a42fd65c00bda4f1633bc90ac169;hpb=e23eb12b5f2c9b50b7db68b9aa972e50155540eb;p=samtools.git

diff --git a/bam_mcns.c b/bam_mcns.c
index d98d480..d337e5a 100644
--- a/bam_mcns.c
+++ b/bam_mcns.c
@@ -264,8 +264,8 @@ static void mc_add_afs(mc_aux_t *ma, double PD, double *f_map, double *p_map)
 			if (ma->afs1[k] > max) max = ma->afs1[k], max_k = k;
 			e += k * ma->afs1[k];
 		}
-		*f_map = .5 * max_k / ma->n; *p_map = max;
-		printf(" * %.3lg:%.3lg:%.3lg:%.3lg * ", sum, 1.-.5*max_k/ma->n, max, 1.-.5*e/ma->n);
+		*f_map = .5 * max_k / ma->n; *p_map = max; // e should equal mc_rst_t::f_exp
+//		printf(" * %.3lg:%.3lg:%.3lg:%.3lg * ", sum, 1.-.5*max_k/ma->n, max, 1.-.5*e/ma->n);
 	}
 }
 
@@ -297,3 +297,13 @@ int mc_cal(int ref, int *n, const bam_pileup1_t **plp, mc_aux_t *ma, mc_rst_t *r
 		mc_add_afs(ma, rst->PD, &rst->f_map, &rst->p_map);
 	return tot;
 }
+
+void mc_dump_afs(mc_aux_t *ma)
+{
+	int k;
+	fprintf(stderr, "[afs]");
+	for (k = 0; k <= ma->M; ++k)
+		fprintf(stderr, " %d:%.3lf", k, ma->afs[ma->M - k]);
+	fprintf(stderr, "\n");
+	memset(ma->afs, 0, sizeof(double) * (ma->M + 1));
+}