dlls/winemp3.acm/tabinit.c

   1 /*
   2  * Copyright (c) Michael Hipp and other authors of the mpglib project.
   3  *
   4  * This library 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  * This library 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 this library; if not, write to the Free Software
  16  * Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA
  17  */
  18
  19 #include <stdlib.h>
  20
  21 #include "mpg123.h"
  22
  23 real decwin[512+32];
  24 static real cos64[16],cos32[8],cos16[4],cos8[2],cos4[1];
  25 real *pnts[] = { cos64,cos32,cos16,cos8,cos4 };
  26
  27 #if 0
  28 static unsigned char *conv16to8_buf = NULL;
  29 unsigned char *conv16to8;
  30 #endif
  31
  32 static const long intwinbase[] = {
  33      0,    -1,    -1,    -1,    -1,    -1,    -1,    -2,    -2,    -2,
  34     -2,    -3,    -3,    -4,    -4,    -5,    -5,    -6,    -7,    -7,
  35     -8,    -9,   -10,   -11,   -13,   -14,   -16,   -17,   -19,   -21,
  36    -24,   -26,   -29,   -31,   -35,   -38,   -41,   -45,   -49,   -53,
  37    -58,   -63,   -68,   -73,   -79,   -85,   -91,   -97,  -104,  -111,
  38   -117,  -125,  -132,  -139,  -147,  -154,  -161,  -169,  -176,  -183,
  39   -190,  -196,  -202,  -208,  -213,  -218,  -222,  -225,  -227,  -228,
  40   -228,  -227,  -224,  -221,  -215,  -208,  -200,  -189,  -177,  -163,
  41   -146,  -127,  -106,   -83,   -57,   -29,     2,    36,    72,   111,
  42    153,   197,   244,   294,   347,   401,   459,   519,   581,   645,
  43    711,   779,   848,   919,   991,  1064,  1137,  1210,  1283,  1356,
  44   1428,  1498,  1567,  1634,  1698,  1759,  1817,  1870,  1919,  1962,
  45   2001,  2032,  2057,  2075,  2085,  2087,  2080,  2063,  2037,  2000,
  46   1952,  1893,  1822,  1739,  1644,  1535,  1414,  1280,  1131,   970,
  47    794,   605,   402,   185,   -45,  -288,  -545,  -814, -1095, -1388,
  48  -1692, -2006, -2330, -2663, -3004, -3351, -3705, -4063, -4425, -4788,
  49  -5153, -5517, -5879, -6237, -6589, -6935, -7271, -7597, -7910, -8209,
  50  -8491, -8755, -8998, -9219, -9416, -9585, -9727, -9838, -9916, -9959,
  51  -9966, -9935, -9863, -9750, -9592, -9389, -9139, -8840, -8492, -8092,
  52  -7640, -7134, -6574, -5959, -5288, -4561, -3776, -2935, -2037, -1082,
  53    -70,   998,  2122,  3300,  4533,  5818,  7154,  8540,  9975, 11455,
  54  12980, 14548, 16155, 17799, 19478, 21189, 22929, 24694, 26482, 28289,
  55  30112, 31947, 33791, 35640, 37489, 39336, 41176, 43006, 44821, 46617,
  56  48390, 50137, 51853, 53534, 55178, 56778, 58333, 59838, 61289, 62684,
  57  64019, 65290, 66494, 67629, 68692, 69679, 70590, 71420, 72169, 72835,
  58  73415, 73908, 74313, 74630, 74856, 74992, 75038 };
  59
  60 void make_decode_tables(long scaleval)
  61 {
  62   int i,j,k,kr,divv;
  63   real *table,*costab;
  64
  65
  66   for(i=0;i<5;i++)
  67   {
  68     kr=0x10>>i; divv=0x40>>i;
  69     costab = pnts[i];
  70     for(k=0;k<kr;k++)
  71       costab[k] = 1.0 / (2.0 * cos(M_PI * ((double) k * 2.0 + 1.0) / (double) divv));
  72   }
  73
  74   table = decwin;
  75   scaleval = -scaleval;
  76   for(i=0,j=0;i<256;i++,j++,table+=32)
  77   {
  78     if(table < decwin+512+16)
  79       table[16] = table[0] = (double) intwinbase[j] / 65536.0 * (double) scaleval;
  80     if(i % 32 == 31)
  81       table -= 1023;
  82     if(i % 64 == 63)
  83       scaleval = - scaleval;
  84   }
  85
  86   for( /* i=256 */ ;i<512;i++,j--,table+=32)
  87   {
  88     if(table < decwin+512+16)
  89       table[16] = table[0] = (double) intwinbase[j] / 65536.0 * (double) scaleval;
  90     if(i % 32 == 31)
  91       table -= 1023;
  92     if(i % 64 == 63)
  93       scaleval = - scaleval;
  94   }
  95 }