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-2001             *
   9  * by the XIPHOPHORUS Company http://www.xiph.org/                  *
  10  *                                                                  *
  11  ********************************************************************
  12
  13  function: window functions
  14  last mod: $Id: window.c,v 1.12 2001/02/02 03:51:58 xiphmont Exp $
  15
  16  ********************************************************************/
  17
  18 #include <stdlib.h>
  19 #include <math.h>
  20 #include "os.h"
  21 #include "misc.h"
  22
  23 float *_vorbis_window(int type, int window,int left,int right){
  24   float *ret=_ogg_calloc(window,sizeof(float));
  25
  26   switch(type){
  27   case 0:
  28     /* The 'vorbis window' (window 0) is sin(sin(x)*sin(x)*2pi) */
  29     {
  30       int leftbegin=window/4-left/2;
  31       int rightbegin=window-window/4-right/2;
  32       int i;
  33
  34       for(i=0;i<left;i++){
  35         float x=(i+.5f)/left*M_PI/2.;
  36         x=sin(x);
  37         x*=x;
  38         x*=M_PI/2.f;
  39         x=sin(x);
  40         ret[i+leftbegin]=x;
  41       }
  42
  43       for(i=leftbegin+left;i<rightbegin;i++)
  44         ret[i]=1.f;
  45
  46       for(i=0;i<right;i++){
  47         float x=(right-i-.5f)/right*M_PI/2.;
  48         x=sin(x);
  49         x*=x;
  50         x*=M_PI/2.f;
  51         x=sin(x);
  52         ret[i+rightbegin]=x;
  53       }
  54     }
  55     break;
  56   default:
  57     _ogg_free(ret);
  58     return(NULL);
  59   }
  60   return(ret);
  61 }
  62