Mercurial > libavcodec.hg
comparison ppc/fdct_altivec.c @ 2612:210cca8f5764 libavcodec
gcc 4 compilation fix
patch by Bret Hughes <bhughes at elevating dot com>
author | diego |
---|---|
date | Sun, 17 Apr 2005 21:41:16 +0000 |
parents | 6a4cfc5f9f96 |
children | 4f14a3062d0a |
comparison
equal
deleted
inserted
replaced
2611:00f608ae3791 | 2612:210cca8f5764 |
---|---|
212 POWERPC_PERF_START_COUNT(altivec_fdct, 1); | 212 POWERPC_PERF_START_COUNT(altivec_fdct, 1); |
213 | 213 |
214 | 214 |
215 /* setup constants {{{ */ | 215 /* setup constants {{{ */ |
216 /* mzero = -0.0 */ | 216 /* mzero = -0.0 */ |
217 vu32(mzero) = vec_splat_u32(-1); | 217 mzero = ((vector float)vec_splat_u32(-1)); |
218 vu32(mzero) = vec_sl(vu32(mzero), vu32(mzero)); | 218 mzero = ((vector float)vec_sl(vu32(mzero), vu32(mzero))); |
219 cp = fdctconsts; | 219 cp = fdctconsts; |
220 cnsts0 = vec_ld(0, cp); cp++; | 220 cnsts0 = vec_ld(0, cp); cp++; |
221 cnsts1 = vec_ld(0, cp); cp++; | 221 cnsts1 = vec_ld(0, cp); cp++; |
222 cnsts2 = vec_ld(0, cp); | 222 cnsts2 = vec_ld(0, cp); |
223 /* }}} */ | 223 /* }}} */ |
225 | 225 |
226 /* 8x8 matrix transpose (vector short[8]) {{{ */ | 226 /* 8x8 matrix transpose (vector short[8]) {{{ */ |
227 #define MERGE_S16(hl,a,b) vec_merge##hl(vs16(a), vs16(b)) | 227 #define MERGE_S16(hl,a,b) vec_merge##hl(vs16(a), vs16(b)) |
228 | 228 |
229 bp = (vector signed short*)block; | 229 bp = (vector signed short*)block; |
230 vs16(b00) = vec_ld(0, bp); | 230 b00 = ((vector float)vec_ld(0, bp)); |
231 vs16(b40) = vec_ld(16*4, bp); | 231 b40 = ((vector float)vec_ld(16*4, bp)); |
232 vs16(b01) = MERGE_S16(h, b00, b40); | 232 b01 = ((vector float)MERGE_S16(h, b00, b40)); |
233 vs16(b11) = MERGE_S16(l, b00, b40); | 233 b11 = ((vector float)MERGE_S16(l, b00, b40)); |
234 bp++; | 234 bp++; |
235 vs16(b10) = vec_ld(0, bp); | 235 b10 = ((vector float)vec_ld(0, bp)); |
236 vs16(b50) = vec_ld(16*4, bp); | 236 b50 = ((vector float)vec_ld(16*4, bp)); |
237 vs16(b21) = MERGE_S16(h, b10, b50); | 237 b21 = ((vector float)MERGE_S16(h, b10, b50)); |
238 vs16(b31) = MERGE_S16(l, b10, b50); | 238 b31 = ((vector float)MERGE_S16(l, b10, b50)); |
239 bp++; | 239 bp++; |
240 vs16(b20) = vec_ld(0, bp); | 240 b20 = ((vector float)vec_ld(0, bp)); |
241 vs16(b60) = vec_ld(16*4, bp); | 241 b60 = ((vector float)vec_ld(16*4, bp)); |
242 vs16(b41) = MERGE_S16(h, b20, b60); | 242 b41 = ((vector float)MERGE_S16(h, b20, b60)); |
243 vs16(b51) = MERGE_S16(l, b20, b60); | 243 b51 = ((vector float)MERGE_S16(l, b20, b60)); |
244 bp++; | 244 bp++; |
245 vs16(b30) = vec_ld(0, bp); | 245 b30 = ((vector float)vec_ld(0, bp)); |
246 vs16(b70) = vec_ld(16*4, bp); | 246 b70 = ((vector float)vec_ld(16*4, bp)); |
247 vs16(b61) = MERGE_S16(h, b30, b70); | 247 b61 = ((vector float)MERGE_S16(h, b30, b70)); |
248 vs16(b71) = MERGE_S16(l, b30, b70); | 248 b71 = ((vector float)MERGE_S16(l, b30, b70)); |
249 | 249 |
250 vs16(x0) = MERGE_S16(h, b01, b41); | 250 x0 = ((vector float)MERGE_S16(h, b01, b41)); |
251 vs16(x1) = MERGE_S16(l, b01, b41); | 251 x1 = ((vector float)MERGE_S16(l, b01, b41)); |
252 vs16(x2) = MERGE_S16(h, b11, b51); | 252 x2 = ((vector float)MERGE_S16(h, b11, b51)); |
253 vs16(x3) = MERGE_S16(l, b11, b51); | 253 x3 = ((vector float)MERGE_S16(l, b11, b51)); |
254 vs16(x4) = MERGE_S16(h, b21, b61); | 254 x4 = ((vector float)MERGE_S16(h, b21, b61)); |
255 vs16(x5) = MERGE_S16(l, b21, b61); | 255 x5 = ((vector float)MERGE_S16(l, b21, b61)); |
256 vs16(x6) = MERGE_S16(h, b31, b71); | 256 x6 = ((vector float)MERGE_S16(h, b31, b71)); |
257 vs16(x7) = MERGE_S16(l, b31, b71); | 257 x7 = ((vector float)MERGE_S16(l, b31, b71)); |
258 | 258 |
259 vs16(b00) = MERGE_S16(h, x0, x4); | 259 b00 = ((vector float)MERGE_S16(h, x0, x4)); |
260 vs16(b10) = MERGE_S16(l, x0, x4); | 260 b10 = ((vector float)MERGE_S16(l, x0, x4)); |
261 vs16(b20) = MERGE_S16(h, x1, x5); | 261 b20 = ((vector float)MERGE_S16(h, x1, x5)); |
262 vs16(b30) = MERGE_S16(l, x1, x5); | 262 b30 = ((vector float)MERGE_S16(l, x1, x5)); |
263 vs16(b40) = MERGE_S16(h, x2, x6); | 263 b40 = ((vector float)MERGE_S16(h, x2, x6)); |
264 vs16(b50) = MERGE_S16(l, x2, x6); | 264 b50 = ((vector float)MERGE_S16(l, x2, x6)); |
265 vs16(b60) = MERGE_S16(h, x3, x7); | 265 b60 = ((vector float)MERGE_S16(h, x3, x7)); |
266 vs16(b70) = MERGE_S16(l, x3, x7); | 266 b70 = ((vector float)MERGE_S16(l, x3, x7)); |
267 | 267 |
268 #undef MERGE_S16 | 268 #undef MERGE_S16 |
269 /* }}} */ | 269 /* }}} */ |
270 | 270 |
271 | 271 |
273 * conversion to vector float. The following code section takes advantage | 273 * conversion to vector float. The following code section takes advantage |
274 * of this. | 274 * of this. |
275 */ | 275 */ |
276 #if 1 | 276 #if 1 |
277 /* fdct rows {{{ */ | 277 /* fdct rows {{{ */ |
278 vs16(x0) = vec_add(vs16(b00), vs16(b70)); | 278 x0 = ((vector float)vec_add(vs16(b00), vs16(b70))); |
279 vs16(x7) = vec_sub(vs16(b00), vs16(b70)); | 279 x7 = ((vector float)vec_sub(vs16(b00), vs16(b70))); |
280 vs16(x1) = vec_add(vs16(b10), vs16(b60)); | 280 x1 = ((vector float)vec_add(vs16(b10), vs16(b60))); |
281 vs16(x6) = vec_sub(vs16(b10), vs16(b60)); | 281 x6 = ((vector float)vec_sub(vs16(b10), vs16(b60))); |
282 vs16(x2) = vec_add(vs16(b20), vs16(b50)); | 282 x2 = ((vector float)vec_add(vs16(b20), vs16(b50))); |
283 vs16(x5) = vec_sub(vs16(b20), vs16(b50)); | 283 x5 = ((vector float)vec_sub(vs16(b20), vs16(b50))); |
284 vs16(x3) = vec_add(vs16(b30), vs16(b40)); | 284 x3 = ((vector float)vec_add(vs16(b30), vs16(b40))); |
285 vs16(x4) = vec_sub(vs16(b30), vs16(b40)); | 285 x4 = ((vector float)vec_sub(vs16(b30), vs16(b40))); |
286 | 286 |
287 vs16(b70) = vec_add(vs16(x0), vs16(x3)); | 287 b70 = ((vector float)vec_add(vs16(x0), vs16(x3))); |
288 vs16(b10) = vec_add(vs16(x1), vs16(x2)); | 288 b10 = ((vector float)vec_add(vs16(x1), vs16(x2))); |
289 | 289 |
290 vs16(b00) = vec_add(vs16(b70), vs16(b10)); | 290 b00 = ((vector float)vec_add(vs16(b70), vs16(b10))); |
291 vs16(b40) = vec_sub(vs16(b70), vs16(b10)); | 291 b40 = ((vector float)vec_sub(vs16(b70), vs16(b10))); |
292 | 292 |
293 #define CTF0(n) \ | 293 #define CTF0(n) \ |
294 vs32(b##n##1) = vec_unpackl(vs16(b##n##0)); \ | 294 b##n##1 = ((vector float)vec_unpackl(vs16(b##n##0))); \ |
295 vs32(b##n##0) = vec_unpackh(vs16(b##n##0)); \ | 295 b##n##0 = ((vector float)vec_unpackh(vs16(b##n##0))); \ |
296 b##n##1 = vec_ctf(vs32(b##n##1), 0); \ | 296 b##n##1 = vec_ctf(vs32(b##n##1), 0); \ |
297 b##n##0 = vec_ctf(vs32(b##n##0), 0); | 297 b##n##0 = vec_ctf(vs32(b##n##0), 0); |
298 | 298 |
299 CTF0(0); | 299 CTF0(0); |
300 CTF0(4); | 300 CTF0(4); |
301 | 301 |
302 vs16(b20) = vec_sub(vs16(x0), vs16(x3)); | 302 b20 = ((vector float)vec_sub(vs16(x0), vs16(x3))); |
303 vs16(b60) = vec_sub(vs16(x1), vs16(x2)); | 303 b60 = ((vector float)vec_sub(vs16(x1), vs16(x2))); |
304 | 304 |
305 CTF0(2); | 305 CTF0(2); |
306 CTF0(6); | 306 CTF0(6); |
307 | 307 |
308 #undef CTF0 | 308 #undef CTF0 |
319 cnst = LD_W0; | 319 cnst = LD_W0; |
320 b60 = vec_madd(cnst, b60, x0); | 320 b60 = vec_madd(cnst, b60, x0); |
321 b61 = vec_madd(cnst, b61, x1); | 321 b61 = vec_madd(cnst, b61, x1); |
322 | 322 |
323 #define CTFX(x,b) \ | 323 #define CTFX(x,b) \ |
324 vs32(b##0) = vec_unpackh(vs16(x)); \ | 324 b##0 = ((vector float)vec_unpackh(vs16(x))); \ |
325 vs32(b##1) = vec_unpackl(vs16(x)); \ | 325 b##1 = ((vector float)vec_unpackl(vs16(x))); \ |
326 b##0 = vec_ctf(vs32(b##0), 0); \ | 326 b##0 = vec_ctf(vs32(b##0), 0); \ |
327 b##1 = vec_ctf(vs32(b##1), 0); \ | 327 b##1 = vec_ctf(vs32(b##1), 0); \ |
328 | 328 |
329 CTFX(x4, b7); | 329 CTFX(x4, b7); |
330 CTFX(x5, b5); | 330 CTFX(x5, b5); |
471 | 471 |
472 /* round, convert back to short {{{ */ | 472 /* round, convert back to short {{{ */ |
473 #define CTS(n) \ | 473 #define CTS(n) \ |
474 b##n##0 = vec_round(b##n##0); \ | 474 b##n##0 = vec_round(b##n##0); \ |
475 b##n##1 = vec_round(b##n##1); \ | 475 b##n##1 = vec_round(b##n##1); \ |
476 vs32(b##n##0) = vec_cts(b##n##0, 0); \ | 476 b##n##0 = ((vector float)vec_cts(b##n##0, 0)); \ |
477 vs32(b##n##1) = vec_cts(b##n##1, 0); \ | 477 b##n##1 = ((vector float)vec_cts(b##n##1, 0)); \ |
478 vs16(b##n##0) = vec_pack(vs32(b##n##0), vs32(b##n##1)); \ | 478 b##n##0 = ((vector float)vec_pack(vs32(b##n##0), vs32(b##n##1))); \ |
479 vec_st(vs16(b##n##0), 0, bp); | 479 vec_st(vs16(b##n##0), 0, bp); |
480 | 480 |
481 bp = (vector signed short*)block; | 481 bp = (vector signed short*)block; |
482 CTS(0); bp++; | 482 CTS(0); bp++; |
483 CTS(1); bp++; | 483 CTS(1); bp++; |