diff options
author | Mans Rullgard <mans@mansr.com> | 2011-03-19 17:19:04 +0000 |
---|---|---|
committer | Mans Rullgard <mans@mansr.com> | 2011-03-19 19:49:27 +0000 |
commit | a45fbda994df54600101e845e29f0af60e456e95 (patch) | |
tree | 5b7967059edcf79ea6814d77db25adf6a3a387eb /libavcodec/mdct.c | |
parent | 26f548bb59177cfc8c45ff633dd37b60cfd23edf (diff) | |
download | ffmpeg-a45fbda994df54600101e845e29f0af60e456e95.tar.gz |
Move ff_kbd_window_init() to a separate file
This function is not tightly coupled to mdct, and it's in the way
of making a fixed-point mdct implementation.
Signed-off-by: Mans Rullgard <mans@mansr.com>
Diffstat (limited to 'libavcodec/mdct.c')
-rw-r--r-- | libavcodec/mdct.c | 25 |
1 files changed, 0 insertions, 25 deletions
diff --git a/libavcodec/mdct.c b/libavcodec/mdct.c index c99a6cfee2..f2f34f9eda 100644 --- a/libavcodec/mdct.c +++ b/libavcodec/mdct.c @@ -30,31 +30,6 @@ * MDCT/IMDCT transforms. */ -// Generate a Kaiser-Bessel Derived Window. -#define BESSEL_I0_ITER 50 // default: 50 iterations of Bessel I0 approximation -av_cold void ff_kbd_window_init(float *window, float alpha, int n) -{ - int i, j; - double sum = 0.0, bessel, tmp; - double local_window[FF_KBD_WINDOW_MAX]; - double alpha2 = (alpha * M_PI / n) * (alpha * M_PI / n); - - assert(n <= FF_KBD_WINDOW_MAX); - - for (i = 0; i < n; i++) { - tmp = i * (n - i) * alpha2; - bessel = 1.0; - for (j = BESSEL_I0_ITER; j > 0; j--) - bessel = bessel * tmp / (j * j) + 1; - sum += bessel; - local_window[i] = sum; - } - - sum++; - for (i = 0; i < n; i++) - window[i] = sqrt(local_window[i] / sum); -} - #include "mdct_tablegen.h" /** |