Imported Upstream version 0.1.12a

[samtools.git] / bcftools / ld.c

diff --git a/bcftools/ld.c b/bcftools/ld.c
index aa7ec072b2c7fc2e60ea3d1fe8145711a8841d28..dc84d4b02d836f0a07a0736588dd1142485d10cf 100644 (file)
--- a/bcftools/ld.c
+++ b/bcftools/ld.c
@@ -80,7 +80,7 @@ double bcf_ld_freq(const bcf1_t *b0, const bcf1_t *b1, double f[4])
                memcpy(flast, f, 4 * sizeof(double));
                freq_iter(n_smpl, pdg, f);
                for (i = 0; i < 4; ++i) {
-                       double x = fabs(f[0] - flast[0]);
+                       double x = fabs(f[i] - flast[i]);
                        if (x > eps) eps = x;
                }
                if (eps < LD_ITER_EPS) break;