# HG changeset patch # User superdump # Date 1214147547 0 # Node ID b0820b8bd4dd225290e03875c81f2c418b3e980e # Parent 544da38cb2c9521853fac8ff50465a3493672f46 Add generic ff_sine_window_init function and implement in codecs appropriately diff -r 544da38cb2c9 -r b0820b8bd4dd cook.c --- a/cook.c Sun Jun 22 13:20:41 2008 +0000 +++ b/cook.c Sun Jun 22 15:12:27 2008 +0000 @@ -239,9 +239,9 @@ return -1; /* Initialize the MLT window: simple sine window. */ - alpha = M_PI / (2.0 * (float)mlt_size); + ff_sine_window_init(q->mlt_window, mlt_size); for(j=0 ; jmlt_window[j] = sin((j + 0.5) * alpha) * sqrt(2.0 / q->samples_per_channel); + q->mlt_window[j] *= sqrt(2.0 / q->samples_per_channel); /* Initialize the MDCT. */ if (ff_mdct_init(&q->mdct_ctx, av_log2(mlt_size)+1, 1)) { diff -r 544da38cb2c9 -r b0820b8bd4dd dsputil.h --- a/dsputil.h Sun Jun 22 13:20:41 2008 +0000 +++ b/dsputil.h Sun Jun 22 15:12:27 2008 +0000 @@ -649,6 +649,13 @@ */ void ff_kbd_window_init(float *window, float alpha, int n); +/** + * Generate a sine window. + * @param window pointer to half window + * @param n size of half window + */ +void ff_sine_window_init(float *window, int n); + int ff_mdct_init(MDCTContext *s, int nbits, int inverse); void ff_imdct_calc(MDCTContext *s, FFTSample *output, const FFTSample *input, FFTSample *tmp); diff -r 544da38cb2c9 -r b0820b8bd4dd imc.c --- a/imc.c Sun Jun 22 13:20:41 2008 +0000 +++ b/imc.c Sun Jun 22 15:12:27 2008 +0000 @@ -102,8 +102,9 @@ q->old_floor[i] = 1.0; /* Build mdct window, a simple sine window normalized with sqrt(2) */ + ff_sine_window_init(q->mdct_sine_window, COEFFS); for(i = 0; i < COEFFS; i++) - q->mdct_sine_window[i] = sin((i + 0.5) / 512.0 * M_PI) * sqrt(2.0); + q->mdct_sine_window[i] *= sqrt(2.0); for(i = 0; i < COEFFS/2; i++){ q->post_cos[i] = cos(i / 256.0 * M_PI); q->post_sin[i] = sin(i / 256.0 * M_PI); diff -r 544da38cb2c9 -r b0820b8bd4dd mdct.c --- a/mdct.c Sun Jun 22 13:20:41 2008 +0000 +++ b/mdct.c Sun Jun 22 15:12:27 2008 +0000 @@ -48,6 +48,13 @@ window[i] = sqrt(local_window[i] / sum); } +// Generate a sine window. +void ff_sine_window_init(float *window, int n) { + int i; + for(i = 0; i < n; i++) + window[i] = sin((i + 0.5) / (2 * n) * M_PI); +} + /** * init MDCT or IMDCT computation. */ diff -r 544da38cb2c9 -r b0820b8bd4dd nellymoserdec.c --- a/nellymoserdec.c Sun Jun 22 13:20:41 2008 +0000 +++ b/nellymoserdec.c Sun Jun 22 15:12:27 2008 +0000 @@ -148,9 +148,7 @@ /* Generate overlap window */ if (!sine_window[0]) - for (i=0 ; i<128; i++) { - sine_window[i] = sin((i + 0.5) / 256.0 * M_PI); - } + ff_sine_window_init(sine_window, 128); return 0; }