annotate timefilter.c @ 4656:c467710f2fe2 libavformat

Selftest code.
author michael
date Fri, 06 Mar 2009 00:26:51 +0000
parents d9f082019a79
children 6b24b7aa7351
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
4650
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
1 /*
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
2 * Delay Locked Loop based time filter
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
3 * Copyright (c) 2009 Samalyse
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
4 * Author: Olivier Guilyardi <olivier samalyse com>
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
5 *
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
6 * This file is part of FFmpeg.
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
7 *
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
8 * FFmpeg is free software; you can redistribute it and/or
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
9 * modify it under the terms of the GNU Lesser General Public
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
10 * License as published by the Free Software Foundation; either
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
11 * version 2.1 of the License, or (at your option) any later version.
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
12 *
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
13 * FFmpeg is distributed in the hope that it will be useful,
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
14 * but WITHOUT ANY WARRANTY; without even the implied warranty of
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
16 * Lesser General Public License for more details.
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
17 *
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
18 * You should have received a copy of the GNU Lesser General Public
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
19 * License along with FFmpeg; if not, write to the Free Software
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
20 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
21 */
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
22
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
23
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
24 #include "config.h"
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
25 #include "avformat.h"
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
26 #include "timefilter.h"
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
27
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
28 struct TimeFilter {
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
29 /// Delay Locked Loop data. These variables refer to mathematical
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
30 /// concepts described in: http://www.kokkinizita.net/papers/usingdll.pdf
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
31 double cycle_time;
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
32 double feedback2_factor;
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
33 double feedback3_factor;
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
34 double integrator2_state;
4654
d90df812a2a1 Adapt first factor.
michael
parents: 4653
diff changeset
35 int count;
4650
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
36 };
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
37
4652
28b80bbf6288 Allow timespans to vary.
michael
parents: 4651
diff changeset
38 TimeFilter * ff_timefilter_new(double feedback2_factor, double feedback3_factor)
4650
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
39 {
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
40 TimeFilter *self = av_mallocz(sizeof(TimeFilter));
4652
28b80bbf6288 Allow timespans to vary.
michael
parents: 4651
diff changeset
41 self->integrator2_state = 1.0;
4650
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
42 self->feedback2_factor = feedback2_factor;
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
43 self->feedback3_factor = feedback3_factor;
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
44 return self;
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
45 }
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
46
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
47 void ff_timefilter_destroy(TimeFilter *self)
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
48 {
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
49 av_freep(&self);
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
50 }
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
51
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
52 void ff_timefilter_reset(TimeFilter *self)
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
53 {
4654
d90df812a2a1 Adapt first factor.
michael
parents: 4653
diff changeset
54 self->count = 0;
4650
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
55 }
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
56
4653
34a843f958b7 Merge ff_timefilter_read() and ff_timefilter_update(), this simplifies API and
michael
parents: 4652
diff changeset
57 double ff_timefilter_update(TimeFilter *self, double system_time, double period)
4650
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
58 {
4654
d90df812a2a1 Adapt first factor.
michael
parents: 4653
diff changeset
59 self->count++;
4655
d9f082019a79 Fix bug with time==0 being special.
michael
parents: 4654
diff changeset
60 if (self->count==1) {
4650
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
61 /// init loop
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
62 self->cycle_time = system_time;
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
63 } else {
4651
2f6e99161450 Reorder operations, 20% improved performance, 1 variable less.
michael
parents: 4650
diff changeset
64 double loop_error;
4652
28b80bbf6288 Allow timespans to vary.
michael
parents: 4651
diff changeset
65 self->cycle_time+= self->integrator2_state * period;
4650
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
66 /// calculate loop error
4651
2f6e99161450 Reorder operations, 20% improved performance, 1 variable less.
michael
parents: 4650
diff changeset
67 loop_error = system_time - self->cycle_time;
4650
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
68
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
69 /// update loop
4654
d90df812a2a1 Adapt first factor.
michael
parents: 4653
diff changeset
70 self->cycle_time += FFMAX(self->feedback2_factor, 1.0/(self->count)) * loop_error;
4652
28b80bbf6288 Allow timespans to vary.
michael
parents: 4651
diff changeset
71 self->integrator2_state += self->feedback3_factor * loop_error / period;
4650
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
72 }
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
73 return self->cycle_time;
5090afc3cd19 Timefilter code by Olivier Guilyardi o list O samalyse o com
michael
parents:
diff changeset
74 }
4656
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
75
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
76 #ifdef TEST
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
77 main(){
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
78 double n0,n1;
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
79 #define SAMPLES 1000
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
80 double ideal[SAMPLES];
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
81 double samples[SAMPLES];
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
82 for(n0= 0; n0<40; n0=2*n0+1){
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
83 for(n1= 0; n1<10; n1=2*n1+1){
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
84 double best_error= 1000000000;
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
85 double bestpar0=1;
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
86 double bestpar1=0.001;
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
87 int better, i;
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
88
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
89 srandom(123);
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
90 for(i=0; i<SAMPLES; i++){
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
91 ideal[i] = 10 + i + n1*i/(1000);
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
92 samples[i]= ideal[i] + n0*(rand()-RAND_MAX/2)/(RAND_MAX*10LL);
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
93 }
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
94
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
95 do{
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
96 double par0, par1;
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
97 better=0;
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
98 for(par0= bestpar0*0.8; par0<=bestpar0*1.21; par0+=bestpar0*0.05){
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
99 for(par1= bestpar1*0.8; par1<=bestpar1*1.21; par1+=bestpar1*0.05){
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
100 double error=0;
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
101 TimeFilter *tf= ff_timefilter_new(par0, par1);
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
102 for(i=0; i<SAMPLES; i++){
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
103 double filtered;
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
104 filtered= ff_timefilter_update(tf, samples[i], 1);
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
105 error += (filtered - ideal[i]) * (filtered - ideal[i]);
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
106 }
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
107 ff_timefilter_destroy(tf);
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
108 if(error < best_error){
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
109 best_error= error;
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
110 bestpar0= par0;
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
111 bestpar1= par1;
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
112 better=1;
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
113 }
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
114 }
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
115 }
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
116 }while(better);
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
117 printf(" [%f %f %f]", bestpar0, bestpar1, best_error);
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
118 }
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
119 printf("\n");
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
120 }
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
121 }
c467710f2fe2 Selftest code.
michael
parents: 4655
diff changeset
122 #endif