lib/window.c

   1 /********************************************************************
   2  *                                                                  *
   3  * THIS FILE IS PART OF THE OggVorbis SOFTWARE CODEC SOURCE CODE.   *
   4  * USE, DISTRIBUTION AND REPRODUCTION OF THIS SOURCE IS GOVERNED BY *
   5  * THE GNU LESSER/LIBRARY PUBLIC LICENSE, WHICH IS INCLUDED WITH    *
   6  * THIS SOURCE. PLEASE READ THESE TERMS BEFORE DISTRIBUTING.        *
   7  *                                                                  *
   8  * THE OggVorbis SOURCE CODE IS (C) COPYRIGHT 1994-2000             *
   9  * by Monty <monty@xiph.org> and the XIPHOPHORUS Company            *
  10  * http://www.xiph.org/                                             *
  11  *                                                                  *
  12  ********************************************************************
  13
  14  function: window functions
  15  last mod: $Id: window.c,v 1.10 2000/11/14 00:05:32 xiphmont Exp $
  16
  17  ********************************************************************/
  18
  19 #include <stdlib.h>
  20 #include <math.h>
  21 #include "os.h"
  22 #include "misc.h"
  23
  24 float *_vorbis_window(int type, int window,int left,int right){
  25   float *ret=_ogg_calloc(window,sizeof(float));
  26
  27   switch(type){
  28   case 0:
  29     /* The 'vorbis window' (window 0) is sin(sin(x)*sin(x)*2pi) */
  30     {
  31       int leftbegin=window/4-left/2;
  32       int rightbegin=window-window/4-right/2;
  33       int i;
  34
  35       for(i=0;i<left;i++){
  36         float x=(i+.5)/left*M_PI/2.;
  37         x=sin(x);
  38         x*=x;
  39         x*=M_PI/2.;
  40         x=sin(x);
  41         ret[i+leftbegin]=x;
  42       }
  43
  44       for(i=leftbegin+left;i<rightbegin;i++)
  45         ret[i]=1.;
  46
  47       for(i=0;i<right;i++){
  48         float x=(right-i-.5)/right*M_PI/2.;
  49         x=sin(x);
  50         x*=x;
  51         x*=M_PI/2.;
  52         x=sin(x);
  53         ret[i+rightbegin]=x;
  54       }
  55     }
  56     break;
  57   default:
  58     _ogg_free(ret);
  59     return(NULL);
  60   }
  61   return(ret);
  62 }
  63