Rev 1905 | Details | Compare with Previous | Last modification | View Log | RSS feed
| Rev | Author | Line No. | Line |
|---|---|---|---|
| 1905 | serge | 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 | /* That altivec alignment part here should not hurt generic code, I hope */ |
||
| 13 | #ifdef OPT_ALTIVEC |
||
| 14 | static ALIGNED(16) real cos64[16]; |
||
| 15 | static ALIGNED(16) real cos32[8]; |
||
| 16 | static ALIGNED(16) real cos16[4]; |
||
| 17 | static ALIGNED(16) real cos8[2]; |
||
| 18 | static ALIGNED(16) real cos4[1]; |
||
| 19 | #elif defined(REAL_IS_FIXED) && defined(PRECALC_TABLES) |
||
| 20 | static real cos64[16] = |
||
| 21 | { |
||
| 22 | 8398725,8480395,8647771,8909416,9279544,9780026,10443886,11321405, |
||
| 23 | 12491246,14081950,16316987,19619946,24900150,34523836,57170182,170959967 |
||
| 24 | }; |
||
| 25 | static real cos32[8] = |
||
| 26 | { |
||
| 27 | 8429197,8766072,9511743,10851869,13223040,17795219,28897867,85583072 |
||
| 28 | }; |
||
| 29 | static real cos16[4] = |
||
| 30 | { |
||
| 31 | 8552951,10088893,15099095,42998586 |
||
| 32 | }; |
||
| 33 | static real cos8[2] = |
||
| 34 | { |
||
| 35 | 9079764,21920489 |
||
| 36 | }; |
||
| 37 | static real cos4[1] = |
||
| 38 | { |
||
| 39 | 11863283 |
||
| 40 | }; |
||
| 41 | #else |
||
| 42 | static real cos64[16],cos32[8],cos16[4],cos8[2],cos4[1]; |
||
| 43 | #endif |
||
| 44 | |||
| 45 | real *pnts[] = { cos64,cos32,cos16,cos8,cos4 }; |
||
| 46 | |||
| 47 | |||
| 48 | static long intwinbase[] = { |
||
| 49 | 0, -1, -1, -1, -1, -1, -1, -2, -2, -2, |
||
| 50 | -2, -3, -3, -4, -4, -5, -5, -6, -7, -7, |
||
| 51 | -8, -9, -10, -11, -13, -14, -16, -17, -19, -21, |
||
| 52 | -24, -26, -29, -31, -35, -38, -41, -45, -49, -53, |
||
| 53 | -58, -63, -68, -73, -79, -85, -91, -97, -104, -111, |
||
| 54 | -117, -125, -132, -139, -147, -154, -161, -169, -176, -183, |
||
| 55 | -190, -196, -202, -208, -213, -218, -222, -225, -227, -228, |
||
| 56 | -228, -227, -224, -221, -215, -208, -200, -189, -177, -163, |
||
| 57 | -146, -127, -106, -83, -57, -29, 2, 36, 72, 111, |
||
| 58 | 153, 197, 244, 294, 347, 401, 459, 519, 581, 645, |
||
| 59 | 711, 779, 848, 919, 991, 1064, 1137, 1210, 1283, 1356, |
||
| 60 | 1428, 1498, 1567, 1634, 1698, 1759, 1817, 1870, 1919, 1962, |
||
| 61 | 2001, 2032, 2057, 2075, 2085, 2087, 2080, 2063, 2037, 2000, |
||
| 62 | 1952, 1893, 1822, 1739, 1644, 1535, 1414, 1280, 1131, 970, |
||
| 63 | 794, 605, 402, 185, -45, -288, -545, -814, -1095, -1388, |
||
| 64 | -1692, -2006, -2330, -2663, -3004, -3351, -3705, -4063, -4425, -4788, |
||
| 65 | -5153, -5517, -5879, -6237, -6589, -6935, -7271, -7597, -7910, -8209, |
||
| 66 | -8491, -8755, -8998, -9219, -9416, -9585, -9727, -9838, -9916, -9959, |
||
| 67 | -9966, -9935, -9863, -9750, -9592, -9389, -9139, -8840, -8492, -8092, |
||
| 68 | -7640, -7134, -6574, -5959, -5288, -4561, -3776, -2935, -2037, -1082, |
||
| 69 | -70, 998, 2122, 3300, 4533, 5818, 7154, 8540, 9975, 11455, |
||
| 70 | 12980, 14548, 16155, 17799, 19478, 21189, 22929, 24694, 26482, 28289, |
||
| 71 | 30112, 31947, 33791, 35640, 37489, 39336, 41176, 43006, 44821, 46617, |
||
| 72 | 48390, 50137, 51853, 53534, 55178, 56778, 58333, 59838, 61289, 62684, |
||
| 73 | 64019, 65290, 66494, 67629, 68692, 69679, 70590, 71420, 72169, 72835, |
||
| 74 | 73415, 73908, 74313, 74630, 74856, 74992, 75038 }; |
||
| 75 | |||
| 76 | void prepare_decode_tables() |
||
| 77 | { |
||
| 78 | #if !defined(REAL_IS_FIXED) || !defined(PRECALC_TABLES) |
||
| 79 | int i,k,kr,divv; |
||
| 80 | real *costab; |
||
| 81 | |||
| 82 | for(i=0;i<5;i++) |
||
| 83 | { |
||
| 84 | kr=0x10>>i; divv=0x40>>i; |
||
| 85 | costab = pnts[i]; |
||
| 86 | for(k=0;k |
||
| 87 | costab[k] = DOUBLE_TO_REAL(1.0 / (2.0 * cos(M_PI * ((double) k * 2.0 + 1.0) / (double) divv))); |
||
| 88 | } |
||
| 89 | #endif |
||
| 90 | } |
||
| 91 | |||
| 92 | #ifdef OPT_MMXORSSE |
||
| 3960 | Serge | 93 | #if !defined(OPT_X86_64) && !defined(OPT_NEON) |
| 1905 | serge | 94 | void make_decode_tables_mmx_asm(long scaleval, float* decwin_mmx, float *decwins); |
| 95 | void make_decode_tables_mmx(mpg123_handle *fr) |
||
| 96 | { |
||
| 97 | debug("MMX decode tables"); |
||
| 98 | /* Take care: The scale should be like before, when we didn't have float output all around. */ |
||
| 99 | make_decode_tables_mmx_asm((long)((fr->lastscale < 0 ? fr->p.outscale : fr->lastscale)*SHORT_SCALE), fr->decwin_mmx, fr->decwins); |
||
| 100 | debug("MMX decode tables done"); |
||
| 101 | } |
||
| 102 | #else |
||
| 103 | |||
| 104 | /* This mimics round() as defined in C99. We stay C89. */ |
||
| 105 | static int rounded(double f) |
||
| 106 | { |
||
| 107 | return (int)(f>0 ? floor(f+0.5) : ceil(f-0.5)); |
||
| 108 | } |
||
| 109 | |||
| 110 | /* x86-64 doesn't use asm version */ |
||
| 111 | void make_decode_tables_mmx(mpg123_handle *fr) |
||
| 112 | { |
||
| 113 | int i,j,val; |
||
| 114 | int idx = 0; |
||
| 115 | short *ptr = (short *)fr->decwins; |
||
| 116 | /* Scale is always based on 1.0 . */ |
||
| 117 | double scaleval = -0.5*(fr->lastscale < 0 ? fr->p.outscale : fr->lastscale); |
||
| 118 | debug1("MMX decode tables with scaleval %g", scaleval); |
||
| 119 | for(i=0,j=0;i<256;i++,j++,idx+=32) |
||
| 120 | { |
||
| 121 | if(idx < 512+16) |
||
| 122 | fr->decwin_mmx[idx+16] = fr->decwin_mmx[idx] = DOUBLE_TO_REAL((double) intwinbase[j] * scaleval); |
||
| 123 | |||
| 124 | if(i % 32 == 31) |
||
| 125 | idx -= 1023; |
||
| 126 | if(i % 64 == 63) |
||
| 127 | scaleval = - scaleval; |
||
| 128 | } |
||
| 129 | |||
| 130 | for( /* i=256 */ ;i<512;i++,j--,idx+=32) |
||
| 131 | { |
||
| 132 | if(idx < 512+16) |
||
| 133 | fr->decwin_mmx[idx+16] = fr->decwin_mmx[idx] = DOUBLE_TO_REAL((double) intwinbase[j] * scaleval); |
||
| 134 | |||
| 135 | if(i % 32 == 31) |
||
| 136 | idx -= 1023; |
||
| 137 | if(i % 64 == 63) |
||
| 138 | scaleval = - scaleval; |
||
| 139 | } |
||
| 140 | |||
| 141 | for(i=0; i<512; i++) { |
||
| 142 | if(i&1) val = rounded(fr->decwin_mmx[i]*0.5); |
||
| 143 | else val = rounded(fr->decwin_mmx[i]*-0.5); |
||
| 144 | if(val > 32767) val = 32767; |
||
| 145 | else if(val < -32768) val = -32768; |
||
| 146 | ptr[i] = val; |
||
| 147 | } |
||
| 148 | for(i=512; i<512+32; i++) { |
||
| 149 | if(i&1) val = rounded(fr->decwin_mmx[i]*0.5); |
||
| 150 | else val = 0; |
||
| 151 | if(val > 32767) val = 32767; |
||
| 152 | else if(val < -32768) val = -32768; |
||
| 153 | ptr[i] = val; |
||
| 154 | } |
||
| 155 | for(i=0; i<512; i++) { |
||
| 156 | val = rounded(fr->decwin_mmx[511-i]*-0.5); |
||
| 157 | if(val > 32767) val = 32767; |
||
| 158 | else if(val < -32768) val = -32768; |
||
| 159 | ptr[512+32+i] = val; |
||
| 160 | } |
||
| 161 | debug("decode tables done"); |
||
| 162 | } |
||
| 163 | #endif |
||
| 164 | #endif |
||
| 165 | |||
| 166 | void make_decode_tables(mpg123_handle *fr) |
||
| 167 | { |
||
| 168 | int i,j; |
||
| 169 | int idx = 0; |
||
| 170 | /* Scale is always based on 1.0 . */ |
||
| 171 | double scaleval = -0.5*(fr->lastscale < 0 ? fr->p.outscale : fr->lastscale); |
||
| 172 | debug1("decode tables with scaleval %g", scaleval); |
||
| 173 | #ifdef REAL_IS_FIXED |
||
| 174 | long scaleval_long = DOUBLE_TO_REAL_15(scaleval); |
||
| 175 | #endif |
||
| 176 | for(i=0,j=0;i<256;i++,j++,idx+=32) |
||
| 177 | { |
||
| 178 | if(idx < 512+16) |
||
| 179 | #ifdef REAL_IS_FIXED |
||
| 180 | fr->decwin[idx+16] = fr->decwin[idx] = REAL_SCALE_WINDOW(intwinbase[j] * scaleval_long); |
||
| 181 | #else |
||
| 182 | fr->decwin[idx+16] = fr->decwin[idx] = DOUBLE_TO_REAL((double) intwinbase[j] * scaleval); |
||
| 183 | #endif |
||
| 184 | |||
| 185 | if(i % 32 == 31) |
||
| 186 | idx -= 1023; |
||
| 187 | if(i % 64 == 63) |
||
| 188 | #ifdef REAL_IS_FIXED |
||
| 189 | scaleval_long = - scaleval_long; |
||
| 190 | #else |
||
| 191 | scaleval = - scaleval; |
||
| 192 | #endif |
||
| 193 | } |
||
| 194 | |||
| 195 | for( /* i=256 */ ;i<512;i++,j--,idx+=32) |
||
| 196 | { |
||
| 197 | if(idx < 512+16) |
||
| 198 | #ifdef REAL_IS_FIXED |
||
| 199 | fr->decwin[idx+16] = fr->decwin[idx] = REAL_SCALE_WINDOW(intwinbase[j] * scaleval_long); |
||
| 200 | #else |
||
| 201 | fr->decwin[idx+16] = fr->decwin[idx] = DOUBLE_TO_REAL((double) intwinbase[j] * scaleval); |
||
| 202 | #endif |
||
| 203 | |||
| 204 | if(i % 32 == 31) |
||
| 205 | idx -= 1023; |
||
| 206 | if(i % 64 == 63) |
||
| 207 | #ifdef REAL_IS_FIXED |
||
| 208 | scaleval_long = - scaleval_long; |
||
| 209 | #else |
||
| 210 | scaleval = - scaleval; |
||
| 211 | #endif |
||
| 212 | } |
||
| 3960 | Serge | 213 | #if defined(OPT_X86_64) || defined(OPT_ALTIVEC) || defined(OPT_SSE) || defined(OPT_ARM) || defined(OPT_NEON) |
| 214 | if(fr->cpu_opts.type == x86_64 || fr->cpu_opts.type == altivec || fr->cpu_opts.type == sse || fr->cpu_opts.type == arm || fr->cpu_opts.type == neon) |
||
| 1905 | serge | 215 | { /* for float SSE / AltiVec / ARM decoder */ |
| 216 | for(i=512; i<512+32; i++) |
||
| 217 | { |
||
| 218 | fr->decwin[i] = (i&1) ? fr->decwin[i] : 0; |
||
| 219 | } |
||
| 220 | for(i=0; i<512; i++) |
||
| 221 | { |
||
| 222 | fr->decwin[512+32+i] = -fr->decwin[511-i]; |
||
| 223 | } |
||
| 3960 | Serge | 224 | #ifdef OPT_NEON |
| 225 | if(fr->cpu_opts.type == neon) |
||
| 226 | { |
||
| 227 | for(i=0; i<512; i+=2) |
||
| 228 | { |
||
| 229 | fr->decwin[i] = -fr->decwin[i]; |
||
| 230 | } |
||
| 231 | } |
||
| 232 | #endif |
||
| 1905 | serge | 233 | } |
| 234 | #endif |
||
| 235 | debug("decode tables done"); |
||
| 236 | } |
||
| 237 | |||
| 238 | #ifndef NO_8BIT |
||
| 239 | int make_conv16to8_table(mpg123_handle *fr) |
||
| 240 | { |
||
| 241 | int i; |
||
| 242 | int mode = fr->af.encoding; |
||
| 243 | |||
| 244 | /* |
||
| 245 | * ????: 8.0 is right but on SB cards '2.0' is a better value ??? |
||
| 246 | */ |
||
| 247 | const double mul = 8.0; |
||
| 248 | |||
| 249 | if(!fr->conv16to8_buf){ |
||
| 250 | fr->conv16to8_buf = (unsigned char *) malloc(8192); |
||
| 251 | if(!fr->conv16to8_buf) { |
||
| 252 | fr->err = MPG123_ERR_16TO8TABLE; |
||
| 253 | if(NOQUIET) error("Can't allocate 16 to 8 converter table!"); |
||
| 254 | return -1; |
||
| 255 | } |
||
| 256 | fr->conv16to8 = fr->conv16to8_buf + 4096; |
||
| 257 | } |
||
| 258 | |||
| 259 | if(fr->af.encoding == MPG123_ENC_ULAW_8){ |
||
| 260 | double m=127.0 / log(256.0); |
||
| 261 | int c1; |
||
| 262 | |||
| 263 | for(i=-4096;i<4096;i++) { |
||
| 264 | /* dunno whether this is a valid transformation rule ?!?!? */ |
||
| 265 | if(i < 0) |
||
| 266 | c1 = 127 - (int) (log( 1.0 - 255.0 * (double) i*mul / 32768.0 ) * m); |
||
| 267 | else |
||
| 268 | c1 = 255 - (int) (log( 1.0 + 255.0 * (double) i*mul / 32768.0 ) * m); |
||
| 269 | if((c1 < 0 || c1 > 255) && NOQUIET) error2("Converror %d %d",i,c1); |
||
| 270 | |||
| 271 | if(c1 == 0) |
||
| 272 | c1 = 2; |
||
| 273 | fr->conv16to8[i] = (unsigned char) c1; |
||
| 274 | } |
||
| 275 | } |
||
| 276 | else if(mode == MPG123_ENC_SIGNED_8) { |
||
| 277 | for(i=-4096;i<4096;i++) { |
||
| 278 | fr->conv16to8[i] = i>>5; |
||
| 279 | } |
||
| 280 | } |
||
| 281 | else if(mode == MPG123_ENC_UNSIGNED_8) { |
||
| 282 | for(i=-4096;i<4096;i++) { |
||
| 283 | fr->conv16to8[i] = (i>>5)+128; |
||
| 284 | } |
||
| 285 | } |
||
| 286 | else { |
||
| 287 | for(i=-4096;i<4096;i++) { |
||
| 288 | fr->conv16to8[i] = 0; |
||
| 289 | } |
||
| 290 | } |
||
| 291 | return 0; |
||
| 292 | } |
||
| 293 | #endif |
||
| 294 |