# HG changeset patch
# User michael
# Date 1236372511 0
# Node ID 48d9e631f630bc7f0220d105ac42c4b93b5729a2
# Parent  5910b77d9502f7ffca68f799ff0a8b6b30bf14d2
Add code to the self test to print jitter & drift.

diff -r 5910b77d9502 -r 48d9e631f630 timefilter.c
--- a/timefilter.c	Fri Mar 06 03:12:33 2009 +0000
+++ b/timefilter.c	Fri Mar 06 20:48:31 2009 +0000
@@ -81,8 +81,14 @@
 #define SAMPLES 1000
     double ideal[SAMPLES];
     double samples[SAMPLES];
+#if 1
     for(n0= 0; n0<40; n0=2*n0+1){
         for(n1= 0; n1<10; n1=2*n1+1){
+#else
+    {{
+        n0=7;
+        n1=1;
+#endif
             double best_error= 1000000000;
             double bestpar0=1;
             double bestpar1=0.001;
@@ -116,7 +122,19 @@
                     }
                 }
             }while(better);
+#if 0
+            double lastfil=9;
+            TimeFilter *tf= ff_timefilter_new(1, bestpar0, bestpar1);
+            for(i=0; i<SAMPLES; i++){
+                double filtered;
+                filtered=  ff_timefilter_update(tf, samples[i], 1);
+                printf("%f %f %f %f\n", i - samples[i] + 10, filtered - samples[i], samples[FFMAX(i, 1)] - samples[FFMAX(i-1, 0)], filtered - lastfil);
+                lastfil= filtered;
+            }
+            ff_timefilter_destroy(tf);
+#else
             printf(" [%f %f %f]", bestpar0, bestpar1, best_error);
+#endif
         }
         printf("\n");
     }