Mercurial > SDL_sound_CoreAudio

comparison decoders/libmpg123/tabinit.c @ 562:7e08477b0fc1

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
MP3 decoder upgrade work. Ripped out SMPEG and mpglib support, replaced it with "mpg123.c" and libmpg123. libmpg123 is a much better version of mpglib, so it should solve all the problems about MP3's not seeking, or most modern MP3's not playing at all, etc. Since you no longer have to make a tradeoff with SMPEG for features, and SMPEG is basically rotting, I removed it from the project. There is still work to be done with libmpg123...there are MMX, 3DNow, SSE, Altivec, etc decoders which we don't have enabled at the moment, and the build system could use some work to make this compile more cleanly, etc. Still: huge win.
author Ryan C. Gordon <icculus@icculus.org>
date Fri, 30 Jan 2009 02:44:47 -0500
parents
children
comparison
equal deleted inserted replaced
561:f2985e08589c 562:7e08477b0fc1
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