annotate resample2.c @ 7059:54e5978dd41a libavcodec

Change iquant tables to int16.
author michael
date Tue, 17 Jun 2008 14:03:00 +0000
parents 01647ac078a7
children 3f819263176e
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
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
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2835
diff changeset
21
2082
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
22 /**
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
23 * @file resample2.c
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
a8794f26ea98 AUDIOPHILE_KIDDY_MODE
michael
parents: 4701
diff changeset
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
a8794f26ea98 AUDIOPHILE_KIDDY_MODE
michael
parents: 4701
diff changeset
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
a8794f26ea98 AUDIOPHILE_KIDDY_MODE
michael
parents: 4701
diff changeset
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
a8794f26ea98 AUDIOPHILE_KIDDY_MODE
michael
parents: 4701
diff changeset
49 #else
a8794f26ea98 AUDIOPHILE_KIDDY_MODE
michael
parents: 4701
diff changeset
50 #define FILTER_SHIFT 0
a8794f26ea98 AUDIOPHILE_KIDDY_MODE
michael
parents: 4701
diff changeset
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
a8794f26ea98 AUDIOPHILE_KIDDY_MODE
michael
parents: 4701
diff changeset
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
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2835
diff changeset
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
7ea9c6a179e8 Remove unused variable.
diego
parents: 4709
diff changeset
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
a8794f26ea98 AUDIOPHILE_KIDDY_MODE
michael
parents: 4701
diff changeset
133 #ifdef CONFIG_RESAMPLE_AUDIOPHILE_KIDDY_MODE
a8794f26ea98 AUDIOPHILE_KIDDY_MODE
michael
parents: 4701
diff changeset
134 filter[ph * tap_count + i] = tab[i] / norm;
a8794f26ea98 AUDIOPHILE_KIDDY_MODE
michael
parents: 4701
diff changeset
135 #else
a8794f26ea98 AUDIOPHILE_KIDDY_MODE
michael
parents: 4701
diff changeset
136 filter[ph * tap_count + i] = av_clip(lrintf(tab[i] * scale / norm), FELEM_MIN, FELEM_MAX);
a8794f26ea98 AUDIOPHILE_KIDDY_MODE
michael
parents: 4701
diff changeset
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
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
140 #if 0
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
141 {
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
142 #define LEN 1024
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
143 int j,k;
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
144 double sine[LEN + tap_count];
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
145 double filtered[LEN];
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
146 double maxff=-2, minff=2, maxsf=-2, minsf=2;
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
147 for(i=0; i<LEN; i++){
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
148 double ss=0, sf=0, ff=0;
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
149 for(j=0; j<LEN+tap_count; j++)
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
150 sine[j]= cos(i*j*M_PI/LEN);
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
151 for(j=0; j<LEN; j++){
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
152 double sum=0;
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
153 ph=0;
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
154 for(k=0; k<tap_count; k++)
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
155 sum += filter[ph * tap_count + k] * sine[k+j];
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
156 filtered[j]= sum / (1<<FILTER_SHIFT);
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
157 ss+= sine[j + center] * sine[j + center];
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
158 ff+= filtered[j] * filtered[j];
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
159 sf+= sine[j + center] * filtered[j];
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
160 }
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
161 ss= sqrt(2*ss/LEN);
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
162 ff= sqrt(2*ff/LEN);
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
163 sf= 2*sf/LEN;
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
164 maxff= FFMAX(maxff, ff);
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
165 minff= FFMIN(minff, ff);
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
166 maxsf= FFMAX(maxsf, sf);
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
167 minsf= FFMIN(minsf, sf);
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
168 if(i%11==0){
4702
a8794f26ea98 AUDIOPHILE_KIDDY_MODE
michael
parents: 4701
diff changeset
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
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
170 minff=minsf= 2;
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
171 maxff=maxsf= -2;
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
172 }
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
173 }
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
174 }
8155dd60b50b code for testing the accuracy of the resampler
michael
parents: 4694
diff changeset
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
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
178 /**
5129
0244bba24b43 misc typo fixes
diego
parents: 4962
diff changeset
179 * Initializes an audio resampler.
0244bba24b43 misc typo fixes
diego
parents: 4962
diff changeset
180 * Note, if either rate is not an integer then simply scale both rates up so they are.
2082
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
181 */
2308
de38526a1f3f user selectble cutoff frequency
michael
parents: 2307
diff changeset
182 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
183 AVResampleContext *c= av_mallocz(sizeof(AVResampleContext));
2308
de38526a1f3f user selectble cutoff frequency
michael
parents: 2307
diff changeset
184 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
185 int phase_count= 1<<phase_shift;
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2835
diff changeset
186
2306
072fc321fbe6 make most resample filter parameters selectable at runtime
michael
parents: 2305
diff changeset
187 c->phase_shift= phase_shift;
072fc321fbe6 make most resample filter parameters selectable at runtime
michael
parents: 2305
diff changeset
188 c->phase_mask= phase_count-1;
072fc321fbe6 make most resample filter parameters selectable at runtime
michael
parents: 2305
diff changeset
189 c->linear= linear;
2082
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
190
2835
d4c4b84e0fac minor fixes for invalid audio data patch by (Wolfram Gloger: wmglo, dent med uni-muenchen de)
michael
parents: 2426
diff changeset
191 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
192 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
193 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
194 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
195 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
196
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
197 c->src_incr= out_rate;
2306
072fc321fbe6 make most resample filter parameters selectable at runtime
michael
parents: 2305
diff changeset
198 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
199 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
200
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
201 return c;
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
202 }
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 void av_resample_close(AVResampleContext *c){
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
205 av_freep(&c->filter_bank);
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
206 av_freep(&c);
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
207 }
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
208
2346
424e6b29de74 av_resample_compensate() doxy
michael
parents: 2308
diff changeset
209 /**
424e6b29de74 av_resample_compensate() doxy
michael
parents: 2308
diff changeset
210 * Compensates samplerate/timestamp drift. The compensation is done by changing
6524
01647ac078a7 typo: occured --> occurred
diego
parents: 5523
diff changeset
211 * the resampler parameters, so no audible clicks or similar distortions occur
2346
424e6b29de74 av_resample_compensate() doxy
michael
parents: 2308
diff changeset
212 * @param compensation_distance distance in output samples over which the compensation should be performed
424e6b29de74 av_resample_compensate() doxy
michael
parents: 2308
diff changeset
213 * @param sample_delta number of output samples which should be output less
424e6b29de74 av_resample_compensate() doxy
michael
parents: 2308
diff changeset
214 *
424e6b29de74 av_resample_compensate() doxy
michael
parents: 2308
diff changeset
215 * example: av_resample_compensate(c, 10, 500)
424e6b29de74 av_resample_compensate() doxy
michael
parents: 2308
diff changeset
216 * here instead of 510 samples only 500 samples would be output
424e6b29de74 av_resample_compensate() doxy
michael
parents: 2308
diff changeset
217 *
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2835
diff changeset
218 * note, due to rounding the actual compensation might be slightly different,
2346
424e6b29de74 av_resample_compensate() doxy
michael
parents: 2308
diff changeset
219 * especially if the compensation_distance is large and the in_rate used during init is small
424e6b29de74 av_resample_compensate() doxy
michael
parents: 2308
diff changeset
220 */
2082
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
221 void av_resample_compensate(AVResampleContext *c, int sample_delta, int compensation_distance){
2083
76cdbe832239 avoid useless normalization and 10l fix
michael
parents: 2082
diff changeset
222 // 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
223 c->compensation_distance= compensation_distance;
2083
76cdbe832239 avoid useless normalization and 10l fix
michael
parents: 2082
diff changeset
224 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
225 }
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
226
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
227 /**
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
228 * resamples.
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
229 * @param src an array of unconsumed samples
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
230 * @param consumed the number of samples of src which have been consumed are returned here
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
231 * @param src_size the number of unconsumed samples available
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
232 * @param dst_size the amount of space in samples available in dst
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
233 * @param update_ctx if this is 0 then the context wont be modified, that way several channels can be resampled with the same context
6524
01647ac078a7 typo: occured --> occurred
diego
parents: 5523
diff changeset
234 * @return the number of samples written in dst or -1 if an error occurred
2082
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
235 */
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
236 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
237 int dst_index, i;
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
238 int index= c->index;
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
239 int frac= c->frac;
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
240 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
241 int dst_incr= c->dst_incr / c->src_incr;
2277
b88e05b9b445 dont just resample half of the data sometimes
michael
parents: 2171
diff changeset
242 int compensation_distance= c->compensation_distance;
2405
21c5ef8a4ff7 faster and slightly less accurate nearest neighbor resampler
michael
parents: 2403
diff changeset
243
2403
6d53608cf2cf special case for filter_length==1
michael
parents: 2346
diff changeset
244 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
245 int64_t index2= ((int64_t)index)<<32;
21c5ef8a4ff7 faster and slightly less accurate nearest neighbor resampler
michael
parents: 2403
diff changeset
246 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
247 dst_size= FFMIN(dst_size, (src_size-1-index) * (int64_t)c->src_incr / c->dst_incr);
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2835
diff changeset
248
2403
6d53608cf2cf special case for filter_length==1
michael
parents: 2346
diff changeset
249 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
250 dst[dst_index] = src[index2>>32];
21c5ef8a4ff7 faster and slightly less accurate nearest neighbor resampler
michael
parents: 2403
diff changeset
251 index2 += incr;
2403
6d53608cf2cf special case for filter_length==1
michael
parents: 2346
diff changeset
252 }
2405
21c5ef8a4ff7 faster and slightly less accurate nearest neighbor resampler
michael
parents: 2403
diff changeset
253 frac += dst_index * dst_incr_frac;
21c5ef8a4ff7 faster and slightly less accurate nearest neighbor resampler
michael
parents: 2403
diff changeset
254 index += dst_index * dst_incr;
21c5ef8a4ff7 faster and slightly less accurate nearest neighbor resampler
michael
parents: 2403
diff changeset
255 index += frac / c->src_incr;
21c5ef8a4ff7 faster and slightly less accurate nearest neighbor resampler
michael
parents: 2403
diff changeset
256 frac %= c->src_incr;
2403
6d53608cf2cf special case for filter_length==1
michael
parents: 2346
diff changeset
257 }else{
2082
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
258 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
259 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
260 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
261 FELEM2 val=0;
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2835
diff changeset
262
2082
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
263 if(sample_index < 0){
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
264 for(i=0; i<c->filter_length; i++)
4001
34fdffe98bd0 Rename ABS macro to FFABS.
diego
parents: 3947
diff changeset
265 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
266 }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
267 break;
2306
072fc321fbe6 make most resample filter parameters selectable at runtime
michael
parents: 2305
diff changeset
268 }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
269 FELEM2 v2=0;
2082
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
270 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
271 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
272 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
273 }
4707
441cded6920c optimize linear filter coeff interpolation code, this also makes the code less prone to overflows
michael
parents: 4703
diff changeset
274 val+=(v2-val)*(FELEML)frac / c->src_incr;
2306
072fc321fbe6 make most resample filter parameters selectable at runtime
michael
parents: 2305
diff changeset
275 }else{
2082
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
276 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
277 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
278 }
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
279 }
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
280
4702
a8794f26ea98 AUDIOPHILE_KIDDY_MODE
michael
parents: 4701
diff changeset
281 #ifdef CONFIG_RESAMPLE_AUDIOPHILE_KIDDY_MODE
5523
c2ab2ac31edb use av_clip_int16() where it makes sense
aurel
parents: 5215
diff changeset
282 dst[dst_index] = av_clip_int16(lrintf(val));
4702
a8794f26ea98 AUDIOPHILE_KIDDY_MODE
michael
parents: 4701
diff changeset
283 #else
2082
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
284 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
285 dst[dst_index] = (unsigned)(val + 32768) > 65535 ? (val>>31) ^ 32767 : val;
4702
a8794f26ea98 AUDIOPHILE_KIDDY_MODE
michael
parents: 4701
diff changeset
286 #endif
2082
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
287
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
288 frac += dst_incr_frac;
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
289 index += dst_incr;
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
290 if(frac >= c->src_incr){
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
291 frac -= c->src_incr;
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
292 index++;
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
293 }
2277
b88e05b9b445 dont just resample half of the data sometimes
michael
parents: 2171
diff changeset
294
b88e05b9b445 dont just resample half of the data sometimes
michael
parents: 2171
diff changeset
295 if(dst_index + 1 == compensation_distance){
b88e05b9b445 dont just resample half of the data sometimes
michael
parents: 2171
diff changeset
296 compensation_distance= 0;
b88e05b9b445 dont just resample half of the data sometimes
michael
parents: 2171
diff changeset
297 dst_incr_frac= c->ideal_dst_incr % c->src_incr;
b88e05b9b445 dont just resample half of the data sometimes
michael
parents: 2171
diff changeset
298 dst_incr= c->ideal_dst_incr / c->src_incr;
b88e05b9b445 dont just resample half of the data sometimes
michael
parents: 2171
diff changeset
299 }
2082
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
300 }
2403
6d53608cf2cf special case for filter_length==1
michael
parents: 2346
diff changeset
301 }
2306
072fc321fbe6 make most resample filter parameters selectable at runtime
michael
parents: 2305
diff changeset
302 *consumed= FFMAX(index, 0) >> c->phase_shift;
2307
df3e4a8e6aac 100l (dont reset the sample position at the block end)
michael
parents: 2306
diff changeset
303 if(index>=0) index &= c->phase_mask;
2109
84637c6ca74e various resampling fixes
michael
parents: 2083
diff changeset
304
2277
b88e05b9b445 dont just resample half of the data sometimes
michael
parents: 2171
diff changeset
305 if(compensation_distance){
b88e05b9b445 dont just resample half of the data sometimes
michael
parents: 2171
diff changeset
306 compensation_distance -= dst_index;
b88e05b9b445 dont just resample half of the data sometimes
michael
parents: 2171
diff changeset
307 assert(compensation_distance > 0);
b88e05b9b445 dont just resample half of the data sometimes
michael
parents: 2171
diff changeset
308 }
2082
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
309 if(update_ctx){
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
310 c->frac= frac;
2109
84637c6ca74e various resampling fixes
michael
parents: 2083
diff changeset
311 c->index= index;
2277
b88e05b9b445 dont just resample half of the data sometimes
michael
parents: 2171
diff changeset
312 c->dst_incr= dst_incr_frac + c->src_incr*dst_incr;
b88e05b9b445 dont just resample half of the data sometimes
michael
parents: 2171
diff changeset
313 c->compensation_distance= compensation_distance;
2082
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
314 }
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2835
diff changeset
315 #if 0
2083
76cdbe832239 avoid useless normalization and 10l fix
michael
parents: 2082
diff changeset
316 if(update_ctx && !c->compensation_distance){
76cdbe832239 avoid useless normalization and 10l fix
michael
parents: 2082
diff changeset
317 #undef rand
76cdbe832239 avoid useless normalization and 10l fix
michael
parents: 2082
diff changeset
318 av_resample_compensate(c, rand() % (8000*2) - 8000, 8000*2);
76cdbe832239 avoid useless normalization and 10l fix
michael
parents: 2082
diff changeset
319 av_log(NULL, AV_LOG_DEBUG, "%d %d %d\n", c->dst_incr, c->ideal_dst_incr, c->compensation_distance);
76cdbe832239 avoid useless normalization and 10l fix
michael
parents: 2082
diff changeset
320 }
76cdbe832239 avoid useless normalization and 10l fix
michael
parents: 2082
diff changeset
321 #endif
2967
ef2149182f1c COSMETICS: Remove all trailing whitespace.
diego
parents: 2835
diff changeset
322
2082
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
323 return dst_index;
3dc9bbe1b152 polyphase kaiser windowed sinc and blackman nuttall windowed sinc audio resample filters
michael
parents:
diff changeset
324 }