Mercurial > libavcodec.hg
annotate resample2.c @ 8991:ca768cb2bfb6 libavcodec
Use last decoded SPS as current SPS in order to parse picture timing SEI
correctly. This works around an apparent H.264 standard deficiency.
Patch by Ivan Schreter, schreter gmx net
author | cehoyos |
---|---|
date | Fri, 20 Feb 2009 16:20:01 +0000 |
parents | bdb11e2330d1 |
children | fff66291d84d |
rev | line source |
---|---|
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
1 /* |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
2 * audio resampling |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
3 * Copyright (c) 2004 Michael Niedermayer <michaelni@gmx.at> |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
4 * |
3947
c8c591fe26f8
Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents:
3347
diff
changeset
|
5 * This file is part of FFmpeg. |
c8c591fe26f8
Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents:
3347
diff
changeset
|
6 * |
c8c591fe26f8
Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents:
3347
diff
changeset
|
7 * FFmpeg is free software; you can redistribute it and/or |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
8 * modify it under the terms of the GNU Lesser General Public |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
9 * License as published by the Free Software Foundation; either |
3947
c8c591fe26f8
Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents:
3347
diff
changeset
|
10 * version 2.1 of the License, or (at your option) any later version. |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
11 * |
3947
c8c591fe26f8
Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents:
3347
diff
changeset
|
12 * FFmpeg is distributed in the hope that it will be useful, |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
15 * Lesser General Public License for more details. |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
16 * |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
17 * You should have received a copy of the GNU Lesser General Public |
3947
c8c591fe26f8
Change license headers to say 'FFmpeg' instead of 'this program/this library'
diego
parents:
3347
diff
changeset
|
18 * License along with FFmpeg; if not, write to the Free Software |
3036
0b546eab515d
Update licensing information: The FSF changed postal address.
diego
parents:
2967
diff
changeset
|
19 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
20 */ |
2967 | 21 |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
22 /** |
8718
e9d9d946f213
Use full internal pathname in doxygen @file directives.
diego
parents:
7377
diff
changeset
|
23 * @file libavcodec/resample2.c |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
24 * audio resampling |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
25 * @author Michael Niedermayer <michaelni@gmx.at> |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
26 */ |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
27 |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
28 #include "avcodec.h" |
2114
7f88c429ad65
ugly missing lrintf workaround by ("Steven M. Schultz" <sms at 2bsd dot com>)
michael
parents:
2109
diff
changeset
|
29 #include "dsputil.h" |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
30 |
4697
8e460c6a85a7
make high precision mode accessible at compile time
michael
parents:
4696
diff
changeset
|
31 #ifndef CONFIG_RESAMPLE_HP |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
32 #define FILTER_SHIFT 15 |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
33 |
2304
09b2ef0ac97d
make filter size, int32/int16 and a few other things selectable at compiletime
michael
parents:
2303
diff
changeset
|
34 #define FELEM int16_t |
09b2ef0ac97d
make filter size, int32/int16 and a few other things selectable at compiletime
michael
parents:
2303
diff
changeset
|
35 #define FELEM2 int32_t |
4702 | 36 #define FELEML int64_t |
2304
09b2ef0ac97d
make filter size, int32/int16 and a few other things selectable at compiletime
michael
parents:
2303
diff
changeset
|
37 #define FELEM_MAX INT16_MAX |
09b2ef0ac97d
make filter size, int32/int16 and a few other things selectable at compiletime
michael
parents:
2303
diff
changeset
|
38 #define FELEM_MIN INT16_MIN |
4700
e210bbc7dd19
select more sensible default windows (= attenuation beyond the dynamic range of your input is silly if it negatively affects other parameters)
michael
parents:
4699
diff
changeset
|
39 #define WINDOW_TYPE 9 |
4702 | 40 #elif !defined(CONFIG_RESAMPLE_AUDIOPHILE_KIDDY_MODE) |
4708
0ebf322d148e
change high precision mode to 30 bits again after making the code less overflow sensitive
michael
parents:
4707
diff
changeset
|
41 #define FILTER_SHIFT 30 |
2304
09b2ef0ac97d
make filter size, int32/int16 and a few other things selectable at compiletime
michael
parents:
2303
diff
changeset
|
42 |
09b2ef0ac97d
make filter size, int32/int16 and a few other things selectable at compiletime
michael
parents:
2303
diff
changeset
|
43 #define FELEM int32_t |
09b2ef0ac97d
make filter size, int32/int16 and a few other things selectable at compiletime
michael
parents:
2303
diff
changeset
|
44 #define FELEM2 int64_t |
4702 | 45 #define FELEML int64_t |
2304
09b2ef0ac97d
make filter size, int32/int16 and a few other things selectable at compiletime
michael
parents:
2303
diff
changeset
|
46 #define FELEM_MAX INT32_MAX |
09b2ef0ac97d
make filter size, int32/int16 and a few other things selectable at compiletime
michael
parents:
2303
diff
changeset
|
47 #define FELEM_MIN INT32_MIN |
4700
e210bbc7dd19
select more sensible default windows (= attenuation beyond the dynamic range of your input is silly if it negatively affects other parameters)
michael
parents:
4699
diff
changeset
|
48 #define WINDOW_TYPE 12 |
4702 | 49 #else |
50 #define FILTER_SHIFT 0 | |
51 | |
4709
6768647cc0fd
replace long double in the audiophile kiddy mode with doubles, its faster (and more portable) and the audiophile placebo 120db stopband attenuation isnt affected
michael
parents:
4708
diff
changeset
|
52 #define FELEM double |
6768647cc0fd
replace long double in the audiophile kiddy mode with doubles, its faster (and more portable) and the audiophile placebo 120db stopband attenuation isnt affected
michael
parents:
4708
diff
changeset
|
53 #define FELEM2 double |
6768647cc0fd
replace long double in the audiophile kiddy mode with doubles, its faster (and more portable) and the audiophile placebo 120db stopband attenuation isnt affected
michael
parents:
4708
diff
changeset
|
54 #define FELEML double |
4702 | 55 #define WINDOW_TYPE 24 |
2304
09b2ef0ac97d
make filter size, int32/int16 and a few other things selectable at compiletime
michael
parents:
2303
diff
changeset
|
56 #endif |
09b2ef0ac97d
make filter size, int32/int16 and a few other things selectable at compiletime
michael
parents:
2303
diff
changeset
|
57 |
09b2ef0ac97d
make filter size, int32/int16 and a few other things selectable at compiletime
michael
parents:
2303
diff
changeset
|
58 |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
59 typedef struct AVResampleContext{ |
2304
09b2ef0ac97d
make filter size, int32/int16 and a few other things selectable at compiletime
michael
parents:
2303
diff
changeset
|
60 FELEM *filter_bank; |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
61 int filter_length; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
62 int ideal_dst_incr; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
63 int dst_incr; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
64 int index; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
65 int frac; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
66 int src_incr; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
67 int compensation_distance; |
2306
072fc321fbe6
make most resample filter parameters selectable at runtime
michael
parents:
2305
diff
changeset
|
68 int phase_shift; |
072fc321fbe6
make most resample filter parameters selectable at runtime
michael
parents:
2305
diff
changeset
|
69 int phase_mask; |
072fc321fbe6
make most resample filter parameters selectable at runtime
michael
parents:
2305
diff
changeset
|
70 int linear; |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
71 }AVResampleContext; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
72 |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
73 /** |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
74 * 0th order modified bessel function of the first kind. |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
75 */ |
3347
82277c821113
Add const to (mostly) char* and make some functions static, which aren't used
diego
parents:
3036
diff
changeset
|
76 static double bessel(double x){ |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
77 double v=1; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
78 double t=1; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
79 int i; |
2967 | 80 |
4701
5fa551fb7640
optimize bessel function instead of trusting gcc to do trivial optimizations (as gcc doesnt ...)
michael
parents:
4700
diff
changeset
|
81 x= x*x/4; |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
82 for(i=1; i<50; i++){ |
4701
5fa551fb7640
optimize bessel function instead of trusting gcc to do trivial optimizations (as gcc doesnt ...)
michael
parents:
4700
diff
changeset
|
83 t *= x/(i*i); |
5fa551fb7640
optimize bessel function instead of trusting gcc to do trivial optimizations (as gcc doesnt ...)
michael
parents:
4700
diff
changeset
|
84 v += t; |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
85 } |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
86 return v; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
87 } |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
88 |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
89 /** |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
90 * builds a polyphase filterbank. |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
91 * @param factor resampling factor |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
92 * @param scale wanted sum of coefficients for each filter |
4699
4040ca52670a
make kaiser windows with other beta than 16 available
michael
parents:
4697
diff
changeset
|
93 * @param type 0->cubic, 1->blackman nuttall windowed sinc, 2..16->kaiser windowed sinc beta=2..16 |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
94 */ |
2304
09b2ef0ac97d
make filter size, int32/int16 and a few other things selectable at compiletime
michael
parents:
2303
diff
changeset
|
95 void av_build_filter(FELEM *filter, double factor, int tap_count, int phase_count, int scale, int type){ |
4822 | 96 int ph, i; |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
97 double x, y, w, tab[tap_count]; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
98 const int center= (tap_count-1)/2; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
99 |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
100 /* if upsampling, only need to interpolate, no filter */ |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
101 if (factor > 1.0) |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
102 factor = 1.0; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
103 |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
104 for(ph=0;ph<phase_count;ph++) { |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
105 double norm = 0; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
106 for(i=0;i<tap_count;i++) { |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
107 x = M_PI * ((double)(i - center) - (double)ph / phase_count) * factor; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
108 if (x == 0) y = 1.0; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
109 else y = sin(x) / x; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
110 switch(type){ |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
111 case 0:{ |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
112 const float d= -0.5; //first order derivative = -0.5 |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
113 x = fabs(((double)(i - center) - (double)ph / phase_count) * factor); |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
114 if(x<1.0) y= 1 - 3*x*x + 2*x*x*x + d*( -x*x + x*x*x); |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
115 else y= d*(-4 + 8*x - 5*x*x + x*x*x); |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
116 break;} |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
117 case 1: |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
118 w = 2.0*x / (factor*tap_count) + M_PI; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
119 y *= 0.3635819 - 0.4891775 * cos(w) + 0.1365995 * cos(2*w) - 0.0106411 * cos(3*w); |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
120 break; |
4699
4040ca52670a
make kaiser windows with other beta than 16 available
michael
parents:
4697
diff
changeset
|
121 default: |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
122 w = 2.0*x / (factor*tap_count*M_PI); |
4699
4040ca52670a
make kaiser windows with other beta than 16 available
michael
parents:
4697
diff
changeset
|
123 y *= bessel(type*sqrt(FFMAX(1-w*w, 0))); |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
124 break; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
125 } |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
126 |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
127 tab[i] = y; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
128 norm += y; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
129 } |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
130 |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
131 /* normalize so that an uniform color remains the same */ |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
132 for(i=0;i<tap_count;i++) { |
4702 | 133 #ifdef CONFIG_RESAMPLE_AUDIOPHILE_KIDDY_MODE |
134 filter[ph * tap_count + i] = tab[i] / norm; | |
135 #else | |
136 filter[ph * tap_count + i] = av_clip(lrintf(tab[i] * scale / norm), FELEM_MIN, FELEM_MAX); | |
137 #endif | |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
138 } |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
139 } |
4695 | 140 #if 0 |
141 { | |
142 #define LEN 1024 | |
143 int j,k; | |
144 double sine[LEN + tap_count]; | |
145 double filtered[LEN]; | |
146 double maxff=-2, minff=2, maxsf=-2, minsf=2; | |
147 for(i=0; i<LEN; i++){ | |
148 double ss=0, sf=0, ff=0; | |
149 for(j=0; j<LEN+tap_count; j++) | |
150 sine[j]= cos(i*j*M_PI/LEN); | |
151 for(j=0; j<LEN; j++){ | |
152 double sum=0; | |
153 ph=0; | |
154 for(k=0; k<tap_count; k++) | |
155 sum += filter[ph * tap_count + k] * sine[k+j]; | |
156 filtered[j]= sum / (1<<FILTER_SHIFT); | |
157 ss+= sine[j + center] * sine[j + center]; | |
158 ff+= filtered[j] * filtered[j]; | |
159 sf+= sine[j + center] * filtered[j]; | |
160 } | |
161 ss= sqrt(2*ss/LEN); | |
162 ff= sqrt(2*ff/LEN); | |
163 sf= 2*sf/LEN; | |
164 maxff= FFMAX(maxff, ff); | |
165 minff= FFMIN(minff, ff); | |
166 maxsf= FFMAX(maxsf, sf); | |
167 minsf= FFMIN(minsf, sf); | |
168 if(i%11==0){ | |
4702 | 169 av_log(NULL, AV_LOG_ERROR, "i:%4d ss:%f ff:%13.6e-%13.6e sf:%13.6e-%13.6e\n", i, ss, maxff, minff, maxsf, minsf); |
4695 | 170 minff=minsf= 2; |
171 maxff=maxsf= -2; | |
172 } | |
173 } | |
174 } | |
175 #endif | |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
176 } |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
177 |
2308 | 178 AVResampleContext *av_resample_init(int out_rate, int in_rate, int filter_size, int phase_shift, int linear, double cutoff){ |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
179 AVResampleContext *c= av_mallocz(sizeof(AVResampleContext)); |
2308 | 180 double factor= FFMIN(out_rate * cutoff / in_rate, 1.0); |
2306
072fc321fbe6
make most resample filter parameters selectable at runtime
michael
parents:
2305
diff
changeset
|
181 int phase_count= 1<<phase_shift; |
2967 | 182 |
2306
072fc321fbe6
make most resample filter parameters selectable at runtime
michael
parents:
2305
diff
changeset
|
183 c->phase_shift= phase_shift; |
072fc321fbe6
make most resample filter parameters selectable at runtime
michael
parents:
2305
diff
changeset
|
184 c->phase_mask= phase_count-1; |
072fc321fbe6
make most resample filter parameters selectable at runtime
michael
parents:
2305
diff
changeset
|
185 c->linear= linear; |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
186 |
2835
d4c4b84e0fac
minor fixes for invalid audio data patch by (Wolfram Gloger: wmglo, dent med uni-muenchen de)
michael
parents:
2426
diff
changeset
|
187 c->filter_length= FFMAX((int)ceil(filter_size/factor), 1); |
2306
072fc321fbe6
make most resample filter parameters selectable at runtime
michael
parents:
2305
diff
changeset
|
188 c->filter_bank= av_mallocz(c->filter_length*(phase_count+1)*sizeof(FELEM)); |
4700
e210bbc7dd19
select more sensible default windows (= attenuation beyond the dynamic range of your input is silly if it negatively affects other parameters)
michael
parents:
4699
diff
changeset
|
189 av_build_filter(c->filter_bank, factor, c->filter_length, phase_count, 1<<FILTER_SHIFT, WINDOW_TYPE); |
2306
072fc321fbe6
make most resample filter parameters selectable at runtime
michael
parents:
2305
diff
changeset
|
190 memcpy(&c->filter_bank[c->filter_length*phase_count+1], c->filter_bank, (c->filter_length-1)*sizeof(FELEM)); |
072fc321fbe6
make most resample filter parameters selectable at runtime
michael
parents:
2305
diff
changeset
|
191 c->filter_bank[c->filter_length*phase_count]= c->filter_bank[c->filter_length - 1]; |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
192 |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
193 c->src_incr= out_rate; |
2306
072fc321fbe6
make most resample filter parameters selectable at runtime
michael
parents:
2305
diff
changeset
|
194 c->ideal_dst_incr= c->dst_incr= in_rate * phase_count; |
072fc321fbe6
make most resample filter parameters selectable at runtime
michael
parents:
2305
diff
changeset
|
195 c->index= -phase_count*((c->filter_length-1)/2); |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
196 |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
197 return c; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
198 } |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
199 |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
200 void av_resample_close(AVResampleContext *c){ |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
201 av_freep(&c->filter_bank); |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
202 av_freep(&c); |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
203 } |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
204 |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
205 void av_resample_compensate(AVResampleContext *c, int sample_delta, int compensation_distance){ |
2083 | 206 // sample_delta += (c->ideal_dst_incr - c->dst_incr)*(int64_t)c->compensation_distance / c->ideal_dst_incr; |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
207 c->compensation_distance= compensation_distance; |
2083 | 208 c->dst_incr = c->ideal_dst_incr - c->ideal_dst_incr * (int64_t)sample_delta / compensation_distance; |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
209 } |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
210 |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
211 int av_resample(AVResampleContext *c, short *dst, short *src, int *consumed, int src_size, int dst_size, int update_ctx){ |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
212 int dst_index, i; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
213 int index= c->index; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
214 int frac= c->frac; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
215 int dst_incr_frac= c->dst_incr % c->src_incr; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
216 int dst_incr= c->dst_incr / c->src_incr; |
2277 | 217 int compensation_distance= c->compensation_distance; |
2405
21c5ef8a4ff7
faster and slightly less accurate nearest neighbor resampler
michael
parents:
2403
diff
changeset
|
218 |
2403 | 219 if(compensation_distance == 0 && c->filter_length == 1 && c->phase_shift==0){ |
2405
21c5ef8a4ff7
faster and slightly less accurate nearest neighbor resampler
michael
parents:
2403
diff
changeset
|
220 int64_t index2= ((int64_t)index)<<32; |
21c5ef8a4ff7
faster and slightly less accurate nearest neighbor resampler
michael
parents:
2403
diff
changeset
|
221 int64_t incr= (1LL<<32) * c->dst_incr / c->src_incr; |
21c5ef8a4ff7
faster and slightly less accurate nearest neighbor resampler
michael
parents:
2403
diff
changeset
|
222 dst_size= FFMIN(dst_size, (src_size-1-index) * (int64_t)c->src_incr / c->dst_incr); |
2967 | 223 |
2403 | 224 for(dst_index=0; dst_index < dst_size; dst_index++){ |
2405
21c5ef8a4ff7
faster and slightly less accurate nearest neighbor resampler
michael
parents:
2403
diff
changeset
|
225 dst[dst_index] = src[index2>>32]; |
21c5ef8a4ff7
faster and slightly less accurate nearest neighbor resampler
michael
parents:
2403
diff
changeset
|
226 index2 += incr; |
2403 | 227 } |
2405
21c5ef8a4ff7
faster and slightly less accurate nearest neighbor resampler
michael
parents:
2403
diff
changeset
|
228 frac += dst_index * dst_incr_frac; |
21c5ef8a4ff7
faster and slightly less accurate nearest neighbor resampler
michael
parents:
2403
diff
changeset
|
229 index += dst_index * dst_incr; |
21c5ef8a4ff7
faster and slightly less accurate nearest neighbor resampler
michael
parents:
2403
diff
changeset
|
230 index += frac / c->src_incr; |
21c5ef8a4ff7
faster and slightly less accurate nearest neighbor resampler
michael
parents:
2403
diff
changeset
|
231 frac %= c->src_incr; |
2403 | 232 }else{ |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
233 for(dst_index=0; dst_index < dst_size; dst_index++){ |
2306
072fc321fbe6
make most resample filter parameters selectable at runtime
michael
parents:
2305
diff
changeset
|
234 FELEM *filter= c->filter_bank + c->filter_length*(index & c->phase_mask); |
072fc321fbe6
make most resample filter parameters selectable at runtime
michael
parents:
2305
diff
changeset
|
235 int sample_index= index >> c->phase_shift; |
2304
09b2ef0ac97d
make filter size, int32/int16 and a few other things selectable at compiletime
michael
parents:
2303
diff
changeset
|
236 FELEM2 val=0; |
2967 | 237 |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
238 if(sample_index < 0){ |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
239 for(i=0; i<c->filter_length; i++) |
4001 | 240 val += src[FFABS(sample_index + i) % src_size] * filter[i]; |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
241 }else if(sample_index + c->filter_length > src_size){ |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
242 break; |
2306
072fc321fbe6
make most resample filter parameters selectable at runtime
michael
parents:
2305
diff
changeset
|
243 }else if(c->linear){ |
4707
441cded6920c
optimize linear filter coeff interpolation code, this also makes the code less prone to overflows
michael
parents:
4703
diff
changeset
|
244 FELEM2 v2=0; |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
245 for(i=0; i<c->filter_length; i++){ |
4707
441cded6920c
optimize linear filter coeff interpolation code, this also makes the code less prone to overflows
michael
parents:
4703
diff
changeset
|
246 val += src[sample_index + i] * (FELEM2)filter[i]; |
441cded6920c
optimize linear filter coeff interpolation code, this also makes the code less prone to overflows
michael
parents:
4703
diff
changeset
|
247 v2 += src[sample_index + i] * (FELEM2)filter[i + c->filter_length]; |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
248 } |
4707
441cded6920c
optimize linear filter coeff interpolation code, this also makes the code less prone to overflows
michael
parents:
4703
diff
changeset
|
249 val+=(v2-val)*(FELEML)frac / c->src_incr; |
2306
072fc321fbe6
make most resample filter parameters selectable at runtime
michael
parents:
2305
diff
changeset
|
250 }else{ |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
251 for(i=0; i<c->filter_length; i++){ |
2304
09b2ef0ac97d
make filter size, int32/int16 and a few other things selectable at compiletime
michael
parents:
2303
diff
changeset
|
252 val += src[sample_index + i] * (FELEM2)filter[i]; |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
253 } |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
254 } |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
255 |
4702 | 256 #ifdef CONFIG_RESAMPLE_AUDIOPHILE_KIDDY_MODE |
5523 | 257 dst[dst_index] = av_clip_int16(lrintf(val)); |
4702 | 258 #else |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
259 val = (val + (1<<(FILTER_SHIFT-1)))>>FILTER_SHIFT; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
260 dst[dst_index] = (unsigned)(val + 32768) > 65535 ? (val>>31) ^ 32767 : val; |
4702 | 261 #endif |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
262 |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
263 frac += dst_incr_frac; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
264 index += dst_incr; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
265 if(frac >= c->src_incr){ |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
266 frac -= c->src_incr; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
267 index++; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
268 } |
2277 | 269 |
270 if(dst_index + 1 == compensation_distance){ | |
271 compensation_distance= 0; | |
272 dst_incr_frac= c->ideal_dst_incr % c->src_incr; | |
273 dst_incr= c->ideal_dst_incr / c->src_incr; | |
274 } | |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
275 } |
2403 | 276 } |
2306
072fc321fbe6
make most resample filter parameters selectable at runtime
michael
parents:
2305
diff
changeset
|
277 *consumed= FFMAX(index, 0) >> c->phase_shift; |
2307
df3e4a8e6aac
100l (dont reset the sample position at the block end)
michael
parents:
2306
diff
changeset
|
278 if(index>=0) index &= c->phase_mask; |
2109 | 279 |
2277 | 280 if(compensation_distance){ |
281 compensation_distance -= dst_index; | |
282 assert(compensation_distance > 0); | |
283 } | |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
284 if(update_ctx){ |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
285 c->frac= frac; |
2109 | 286 c->index= index; |
2277 | 287 c->dst_incr= dst_incr_frac + c->src_incr*dst_incr; |
288 c->compensation_distance= compensation_distance; | |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
289 } |
2967 | 290 #if 0 |
2083 | 291 if(update_ctx && !c->compensation_distance){ |
292 #undef rand | |
293 av_resample_compensate(c, rand() % (8000*2) - 8000, 8000*2); | |
294 av_log(NULL, AV_LOG_DEBUG, "%d %d %d\n", c->dst_incr, c->ideal_dst_incr, c->compensation_distance); | |
295 } | |
296 #endif | |
2967 | 297 |
2082
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
298 return dst_index; |
3dc9bbe1b152
polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff
changeset
|
299 } |