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