git.sesse.net Git - ffmpeg/blob - libavcodec/kbdwin.c

   1 /*
   2  * This file is part of FFmpeg.
   3  *
   4  * FFmpeg is free software; you can redistribute it and/or
   5  * modify it under the terms of the GNU Lesser General Public
   6  * License as published by the Free Software Foundation; either
   7  * version 2.1 of the License, or (at your option) any later version.
   8  *
   9  * FFmpeg is distributed in the hope that it will be useful,
  10  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  12  * Lesser General Public License for more details.
  13  *
  14  * You should have received a copy of the GNU Lesser General Public
  15  * License along with FFmpeg; if not, write to the Free Software
  16  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  17  */
  18
  19
  20 #include <libavutil/avassert.h>
  21 #include <libavutil/mathematics.h>
  22 #include "libavutil/attributes.h"
  23 #include "kbdwin.h"
  24
  25 #define BESSEL_I0_ITER 50 // default: 50 iterations of Bessel I0 approximation
  26
  27 av_cold void ff_kbd_window_init(float *window, float alpha, int n)
  28 {
  29    int i, j;
  30    double sum = 0.0, bessel, tmp;
  31    double local_window[FF_KBD_WINDOW_MAX];
  32    double alpha2 = (alpha * M_PI / n) * (alpha * M_PI / n);
  33
  34    av_assert0(n <= FF_KBD_WINDOW_MAX);
  35
  36    for (i = 0; i < n; i++) {
  37        tmp = i * (n - i) * alpha2;
  38        bessel = 1.0;
  39        for (j = BESSEL_I0_ITER; j > 0; j--)
  40            bessel = bessel * tmp / (j * j) + 1;
  41        sum += bessel;
  42        local_window[i] = sum;
  43    }
  44
  45    sum++;
  46    for (i = 0; i < n; i++)
  47        window[i] = sqrt(local_window[i] / sum);
  48 }