Prevent int overflow with +100k genotypes by casting to double. Fixes s…

…amtools#894
caot · Oct 1, 2018 · cf72059 · cf72059
1 parent 984adb8
commit cf72059
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/plugins/fill-tags.c b/plugins/fill-tags.c
@@ -342,7 +342,7 @@ void calc_hwe(args_t *args, int nref, int nalt, int nhet, float *p_hwe, float *p
     double *probs = args->hwe_probs;
 
     // start at midpoint
-    int mid = nrare * (nref + nalt - nrare) / (nref + nalt);
+    int mid = (double)nrare * (nref + nalt - nrare) / (nref + nalt);
 
     // check to ensure that midpoint and rare alleles have same parity
     if ( (nrare & 1) ^ (mid & 1) ) mid++;