Subversion Repositories Kolibri OS

Rev

Go to most recent revision | Blame | Last modification | View Log | Download | RSS feed

  1. /*
  2.         tabinit.c: initialize tables...
  3.  
  4.         copyright ?-2006 by the mpg123 project - free software under the terms of the LGPL 2.1
  5.         see COPYING and AUTHORS files in distribution or http://mpg123.de
  6.         initially written by Michael Hipp
  7. */
  8.  
  9. //#include <stdlib.h>
  10.  
  11. //#include "config.h"
  12. #include "mpg123.h"
  13. //#include "debug.h"
  14.  
  15. static unsigned char *conv16to8_buf = NULL;
  16. unsigned char *conv16to8;
  17.  
  18. #ifndef USE_MMX
  19. real decwin[512+32];
  20. #ifdef USE_ALTIVEC
  21. static real __attribute__ ((aligned (16))) cos64[16];
  22. static real __attribute__ ((aligned (16))) cos32[8];
  23. static real __attribute__ ((aligned (16))) cos16[4];
  24. static real __attribute__ ((aligned (16))) cos8[2];
  25. static real __attribute__ ((aligned (16))) cos4[1];
  26. #else
  27. static real cos64[16],cos32[8],cos16[4],cos8[2],cos4[1];
  28. #endif
  29.  
  30. real *pnts[] = { cos64,cos32,cos16,cos8,cos4 };
  31.  
  32.  
  33. static long intwinbase[] = {
  34.      0,    -1,    -1,    -1,    -1,    -1,    -1,    -2,    -2,    -2,
  35.     -2,    -3,    -3,    -4,    -4,    -5,    -5,    -6,    -7,    -7,
  36.     -8,    -9,   -10,   -11,   -13,   -14,   -16,   -17,   -19,   -21,
  37.    -24,   -26,   -29,   -31,   -35,   -38,   -41,   -45,   -49,   -53,
  38.    -58,   -63,   -68,   -73,   -79,   -85,   -91,   -97,  -104,  -111,
  39.   -117,  -125,  -132,  -139,  -147,  -154,  -161,  -169,  -176,  -183,
  40.   -190,  -196,  -202,  -208,  -213,  -218,  -222,  -225,  -227,  -228,
  41.   -228,  -227,  -224,  -221,  -215,  -208,  -200,  -189,  -177,  -163,
  42.   -146,  -127,  -106,   -83,   -57,   -29,     2,    36,    72,   111,
  43.    153,   197,   244,   294,   347,   401,   459,   519,   581,   645,
  44.    711,   779,   848,   919,   991,  1064,  1137,  1210,  1283,  1356,
  45.   1428,  1498,  1567,  1634,  1698,  1759,  1817,  1870,  1919,  1962,
  46.   2001,  2032,  2057,  2075,  2085,  2087,  2080,  2063,  2037,  2000,
  47.   1952,  1893,  1822,  1739,  1644,  1535,  1414,  1280,  1131,   970,
  48.    794,   605,   402,   185,   -45,  -288,  -545,  -814, -1095, -1388,
  49.  -1692, -2006, -2330, -2663, -3004, -3351, -3705, -4063, -4425, -4788,
  50.  -5153, -5517, -5879, -6237, -6589, -6935, -7271, -7597, -7910, -8209,
  51.  -8491, -8755, -8998, -9219, -9416, -9585, -9727, -9838, -9916, -9959,
  52.  -9966, -9935, -9863, -9750, -9592, -9389, -9139, -8840, -8492, -8092,
  53.  -7640, -7134, -6574, -5959, -5288, -4561, -3776, -2935, -2037, -1082,
  54.    -70,   998,  2122,  3300,  4533,  5818,  7154,  8540,  9975, 11455,
  55.  12980, 14548, 16155, 17799, 19478, 21189, 22929, 24694, 26482, 28289,
  56.  30112, 31947, 33791, 35640, 37489, 39336, 41176, 43006, 44821, 46617,
  57.  48390, 50137, 51853, 53534, 55178, 56778, 58333, 59838, 61289, 62684,
  58.  64019, 65290, 66494, 67629, 68692, 69679, 70590, 71420, 72169, 72835,
  59.  73415, 73908, 74313, 74630, 74856, 74992, 75038 };
  60.  
  61. #pragma warning(disable:4244)
  62. void make_decode_tables(long scaleval)
  63. {
  64.   int i,j,k,kr,divv;
  65.   real *costab;
  66.   int idx;
  67.  
  68.  
  69.   for(i=0;i<5;i++)
  70.   {
  71.     kr=0x10>>i; divv=0x40>>i;
  72.     costab = pnts[i];
  73.     for(k=0;k<kr;k++)
  74.       costab[k] = DOUBLE_TO_REAL(1.0 / (2.0 * cos(M_PI * ((double) k * 2.0 + 1.0) / (double) divv)));
  75.  
  76.   }
  77.  
  78.   idx = 0;
  79.   scaleval = -scaleval;
  80.   for(i=0,j=0;i<256;i++,j++,idx+=32)
  81.   {
  82.     if(idx < 512+16)
  83.       decwin[idx+16] = decwin[idx] = DOUBLE_TO_REAL((double) intwinbase[j] / 65536.0 * (double) scaleval);
  84.  
  85.     if(i % 32 == 31)
  86.       idx -= 1023;
  87.     if(i % 64 == 63)
  88.       scaleval = - scaleval;
  89.   }
  90.  
  91.   for( /* i=256 */ ;i<512;i++,j--,idx+=32)
  92.   {
  93.     if(idx < 512+16)
  94.       decwin[idx+16] = decwin[idx] = DOUBLE_TO_REAL((double) intwinbase[j] / 65536.0 * (double) scaleval);
  95.  
  96.     if(i % 32 == 31)
  97.       idx -= 1023;
  98.     if(i % 64 == 63)
  99.       scaleval = - scaleval;
  100.   }
  101. }
  102. #endif
  103.  
  104. #if 0
  105. int make_conv16to8_table(int mode)
  106. {
  107.   int i;
  108.  
  109.   /*
  110.    * ????: 8.0 is right but on SB cards '2.0' is a better value ???
  111.    */
  112.   const double mul = 8.0;
  113.  
  114.   if(!conv16to8_buf) {
  115.     conv16to8_buf = (unsigned char *) malloc(8192);
  116.     if(!conv16to8_buf) {
  117.       error("Can't allocate 16 to 8 converter table!");
  118.       return -1;
  119.     }
  120.     conv16to8 = conv16to8_buf + 4096;
  121.   }
  122.  
  123.   if(mode == AUDIO_FORMAT_ULAW_8) {
  124.     double m=127.0 / log(256.0);
  125.     int c1;
  126.  
  127.     for(i=-4096;i<4096;i++) {
  128. /* dunno whether this is a valid transformation rule ?!?!? */
  129.       if(i < 0)
  130.         c1 = 127 - (int) (log( 1.0 - 255.0 * (double) i*mul / 32768.0 ) * m);
  131.       else
  132.         c1 = 255 - (int) (log( 1.0 + 255.0 * (double) i*mul / 32768.0 ) * m);
  133.       if(c1 < 0 || c1 > 255)
  134.         fprintf(stderr,"Converror %d %d\n",i,c1);
  135.       if(c1 == 0)
  136.         c1 = 2;
  137.       conv16to8[i] = (unsigned char) c1;
  138.     }
  139.   }
  140.   else if(mode == AUDIO_FORMAT_SIGNED_8) {
  141.     for(i=-4096;i<4096;i++) {
  142.       conv16to8[i] = i>>5;
  143.     }
  144.   }
  145.   else if(mode == AUDIO_FORMAT_UNSIGNED_8) {
  146.     for(i=-4096;i<4096;i++) {
  147.       conv16to8[i] = (i>>5)+128;
  148.     }
  149.   }
  150.   else {
  151.     for(i=-4096;i<4096;i++) {
  152.       conv16to8[i] = 0;
  153.     }
  154.   }
  155.         return 0;
  156. }
  157.  
  158.  
  159. #endif