Mercurial > libavcodec.hg
comparison jfdctint.c @ 1567:e08df4d22d27 libavcodec
* introducing dct248 into the DSP context.
* simple/accurate implementation of dct248
* DV encoding now supports 2-4-8 DCT
* DV encoding gets a bit faster (but still miles away
from what I think it could do)
* misc. DV codec cleanups
author | romansh |
---|---|
date | Fri, 24 Oct 2003 18:28:01 +0000 |
parents | 1e39f273ecd6 |
children | eb26d190cf5a |
comparison
equal
deleted
inserted
replaced
1566:396e8526e82c | 1567:e08df4d22d27 |
---|---|
293 CONST_BITS+PASS1_BITS); | 293 CONST_BITS+PASS1_BITS); |
294 | 294 |
295 dataptr++; /* advance pointer to next column */ | 295 dataptr++; /* advance pointer to next column */ |
296 } | 296 } |
297 } | 297 } |
298 | |
299 /* | |
300 * The secret of DCT2-4-8 is really simple -- you do the usual 1-DCT | |
301 * on the rows and then, instead of doing even and odd, part on the colums | |
302 * you do even part two times. | |
303 */ | |
304 GLOBAL(void) | |
305 ff_fdct248_islow (DCTELEM * data) | |
306 { | |
307 int32_t tmp0, tmp1, tmp2, tmp3, tmp4, tmp5, tmp6, tmp7; | |
308 int32_t tmp10, tmp11, tmp12, tmp13; | |
309 int32_t z1, z2, z3, z4, z5; | |
310 DCTELEM *dataptr; | |
311 int ctr; | |
312 SHIFT_TEMPS | |
313 | |
314 /* Pass 1: process rows. */ | |
315 /* Note results are scaled up by sqrt(8) compared to a true DCT; */ | |
316 /* furthermore, we scale the results by 2**PASS1_BITS. */ | |
317 | |
318 dataptr = data; | |
319 for (ctr = DCTSIZE-1; ctr >= 0; ctr--) { | |
320 tmp0 = dataptr[0] + dataptr[7]; | |
321 tmp7 = dataptr[0] - dataptr[7]; | |
322 tmp1 = dataptr[1] + dataptr[6]; | |
323 tmp6 = dataptr[1] - dataptr[6]; | |
324 tmp2 = dataptr[2] + dataptr[5]; | |
325 tmp5 = dataptr[2] - dataptr[5]; | |
326 tmp3 = dataptr[3] + dataptr[4]; | |
327 tmp4 = dataptr[3] - dataptr[4]; | |
328 | |
329 /* Even part per LL&M figure 1 --- note that published figure is faulty; | |
330 * rotator "sqrt(2)*c1" should be "sqrt(2)*c6". | |
331 */ | |
332 | |
333 tmp10 = tmp0 + tmp3; | |
334 tmp13 = tmp0 - tmp3; | |
335 tmp11 = tmp1 + tmp2; | |
336 tmp12 = tmp1 - tmp2; | |
337 | |
338 dataptr[0] = (DCTELEM) ((tmp10 + tmp11) << PASS1_BITS); | |
339 dataptr[4] = (DCTELEM) ((tmp10 - tmp11) << PASS1_BITS); | |
340 | |
341 z1 = MULTIPLY(tmp12 + tmp13, FIX_0_541196100); | |
342 dataptr[2] = (DCTELEM) DESCALE(z1 + MULTIPLY(tmp13, FIX_0_765366865), | |
343 CONST_BITS-PASS1_BITS); | |
344 dataptr[6] = (DCTELEM) DESCALE(z1 + MULTIPLY(tmp12, - FIX_1_847759065), | |
345 CONST_BITS-PASS1_BITS); | |
346 | |
347 /* Odd part per figure 8 --- note paper omits factor of sqrt(2). | |
348 * cK represents cos(K*pi/16). | |
349 * i0..i3 in the paper are tmp4..tmp7 here. | |
350 */ | |
351 | |
352 z1 = tmp4 + tmp7; | |
353 z2 = tmp5 + tmp6; | |
354 z3 = tmp4 + tmp6; | |
355 z4 = tmp5 + tmp7; | |
356 z5 = MULTIPLY(z3 + z4, FIX_1_175875602); /* sqrt(2) * c3 */ | |
357 | |
358 tmp4 = MULTIPLY(tmp4, FIX_0_298631336); /* sqrt(2) * (-c1+c3+c5-c7) */ | |
359 tmp5 = MULTIPLY(tmp5, FIX_2_053119869); /* sqrt(2) * ( c1+c3-c5+c7) */ | |
360 tmp6 = MULTIPLY(tmp6, FIX_3_072711026); /* sqrt(2) * ( c1+c3+c5-c7) */ | |
361 tmp7 = MULTIPLY(tmp7, FIX_1_501321110); /* sqrt(2) * ( c1+c3-c5-c7) */ | |
362 z1 = MULTIPLY(z1, - FIX_0_899976223); /* sqrt(2) * (c7-c3) */ | |
363 z2 = MULTIPLY(z2, - FIX_2_562915447); /* sqrt(2) * (-c1-c3) */ | |
364 z3 = MULTIPLY(z3, - FIX_1_961570560); /* sqrt(2) * (-c3-c5) */ | |
365 z4 = MULTIPLY(z4, - FIX_0_390180644); /* sqrt(2) * (c5-c3) */ | |
366 | |
367 z3 += z5; | |
368 z4 += z5; | |
369 | |
370 dataptr[7] = (DCTELEM) DESCALE(tmp4 + z1 + z3, CONST_BITS-PASS1_BITS); | |
371 dataptr[5] = (DCTELEM) DESCALE(tmp5 + z2 + z4, CONST_BITS-PASS1_BITS); | |
372 dataptr[3] = (DCTELEM) DESCALE(tmp6 + z2 + z3, CONST_BITS-PASS1_BITS); | |
373 dataptr[1] = (DCTELEM) DESCALE(tmp7 + z1 + z4, CONST_BITS-PASS1_BITS); | |
374 | |
375 dataptr += DCTSIZE; /* advance pointer to next row */ | |
376 } | |
377 | |
378 /* Pass 2: process columns. | |
379 * We remove the PASS1_BITS scaling, but leave the results scaled up | |
380 * by an overall factor of 8. | |
381 */ | |
382 | |
383 dataptr = data; | |
384 for (ctr = DCTSIZE-1; ctr >= 0; ctr--) { | |
385 tmp0 = dataptr[DCTSIZE*0] + dataptr[DCTSIZE*1]; | |
386 tmp1 = dataptr[DCTSIZE*2] + dataptr[DCTSIZE*3]; | |
387 tmp2 = dataptr[DCTSIZE*4] + dataptr[DCTSIZE*5]; | |
388 tmp3 = dataptr[DCTSIZE*6] + dataptr[DCTSIZE*7]; | |
389 tmp4 = dataptr[DCTSIZE*0] - dataptr[DCTSIZE*1]; | |
390 tmp5 = dataptr[DCTSIZE*2] - dataptr[DCTSIZE*3]; | |
391 tmp6 = dataptr[DCTSIZE*4] - dataptr[DCTSIZE*5]; | |
392 tmp7 = dataptr[DCTSIZE*6] - dataptr[DCTSIZE*7]; | |
393 | |
394 tmp10 = tmp0 + tmp3; | |
395 tmp11 = tmp1 + tmp2; | |
396 tmp12 = tmp1 - tmp2; | |
397 tmp13 = tmp0 - tmp3; | |
398 | |
399 dataptr[DCTSIZE*0] = (DCTELEM) DESCALE(tmp10 + tmp11, PASS1_BITS); | |
400 dataptr[DCTSIZE*4] = (DCTELEM) DESCALE(tmp10 - tmp11, PASS1_BITS); | |
401 | |
402 z1 = MULTIPLY(tmp12 + tmp13, FIX_0_541196100); | |
403 dataptr[DCTSIZE*2] = (DCTELEM) DESCALE(z1 + MULTIPLY(tmp13, FIX_0_765366865), | |
404 CONST_BITS+PASS1_BITS); | |
405 dataptr[DCTSIZE*6] = (DCTELEM) DESCALE(z1 + MULTIPLY(tmp12, - FIX_1_847759065), | |
406 CONST_BITS+PASS1_BITS); | |
407 | |
408 tmp10 = tmp4 + tmp7; | |
409 tmp11 = tmp5 + tmp6; | |
410 tmp12 = tmp5 - tmp6; | |
411 tmp13 = tmp4 - tmp7; | |
412 | |
413 dataptr[DCTSIZE*1] = (DCTELEM) DESCALE(tmp10 + tmp11, PASS1_BITS); | |
414 dataptr[DCTSIZE*5] = (DCTELEM) DESCALE(tmp10 - tmp11, PASS1_BITS); | |
415 | |
416 z1 = MULTIPLY(tmp12 + tmp13, FIX_0_541196100); | |
417 dataptr[DCTSIZE*3] = (DCTELEM) DESCALE(z1 + MULTIPLY(tmp13, FIX_0_765366865), | |
418 CONST_BITS+PASS1_BITS); | |
419 dataptr[DCTSIZE*7] = (DCTELEM) DESCALE(z1 + MULTIPLY(tmp12, - FIX_1_847759065), | |
420 CONST_BITS+PASS1_BITS); | |
421 | |
422 dataptr++; /* advance pointer to next column */ | |
423 } | |
424 } |