From 22e4ba8319d8af4acc909a6f36cbc9a47d4b9309 Mon Sep 17 00:00:00 2001
From: Heng Li <lh3@live.co.uk>
Date: Tue, 23 Nov 2010 14:05:08 +0000
Subject: [PATCH] fixed a typo/bug in r^2 computation

---
 bcftools/ld.c | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/bcftools/ld.c b/bcftools/ld.c
index aa7ec07..dc84d4b 100644
--- 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;
-- 
2.39.2