← Back to branch summary

~ubuntu-branches/ubuntu/lucid/mpg123/lucid

« back to all changes in this revision

Viewing changes to tabinit.c

Tags: upstream-0.60
Import upstream version 0.60

Show diffs side-by-side

added added

removed removed

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