From: sgunderson@bigfoot.com <>
Date: Thu, 28 Dec 2006 13:59:38 +0000 (+0100)
Subject: Fix more bugs in estimate-skew.
X-Git-Url: https://git.sesse.net/?p=audiosync;a=commitdiff_plain;h=e9e02908106107a3265fabc92acd06eb9ee683c1;hp=eacd887c25d1ddde6f165513e9fd1f4a0c5b8afc

Fix more bugs in estimate-skew.
---

diff --git a/estimate-skew.c b/estimate-skew.c
index d2991fc..e8ba75e 100644
--- a/estimate-skew.c
+++ b/estimate-skew.c
@@ -11,7 +11,7 @@ double frac(double x)
 double filter(double x)
 {
 	static double l = 0.0;
-	l = 0.01 * x + 0.99 * l;
+	l = 0.001 * x + 0.999 * l;
 	return l;
 }
 
@@ -54,8 +54,8 @@ int main(int argc, char **argv)
 
 	short prev_sample, sample = 0;
 	int in_pos = -1;
-	double p = -1.0;
-	double speed = 1.0;
+	double p = 0.0;
+	double speed = 1.001;
 
 	while (!feof(in1) && !feof(in2)) {
 		short refs;
@@ -82,9 +82,9 @@ int main(int argc, char **argv)
 		// subtract the two samples and pull it through a filter. we assume
 		// the sines are normalized for now (and that there's no dc skew)
 		double offset = filter(intp_sample - refs);
-		speed -= 0.00000001 * offset;
+		speed += 1e-8 * offset;
 		fwrite(&speed, sizeof(double), 1, skew);
 
-		printf("%f (offset=%f)\n", speed, intp_sample - sample);
+		printf("%f (offset=%f / filt=%f)\n", speed, intp_sample - refs, offset);
 	}
 }