Lines Matching refs:vr
337 #define SUMDIFF_EIGTH(x, y, ix, iy, vr, vi, ur, ui) \
338 vr = (x[ 0 + ix]>>1) + (x[16 + ix]>>1); /* Re A + Re B */ \
342 y[ 0 + iy] = vr + vi; /* Re A' = ReA + ReB +ReC + ReD */ \
343 y[ 4 + iy] = vr - vi; /* Re C' = -(ReC+ReD) + (ReA+ReB) */ \
346 vr -= x[16 + ix]; /* Re A - Re B */ \
350 y[ 2 + iy] = ui + vr; /* Re B' = Im C - Im D + Re A - Re B */ \
351 y[ 6 + iy] = vr - ui; /* Re D' = -Im C + Im D + Re A - Re B */ \
360 FIXP_DBL vr, vi, ur, ui;
363 SUMDIFF_EIGTH(x, y, 0, 0, vr, vi, ur, ui);
364 SUMDIFF_EIGTH(x, y, 4, 8, vr, vi, ur, ui);
365 SUMDIFF_EIGTH(x, y, 2, 16, vr, vi, ur, ui);
366 SUMDIFF_EIGTH(x, y, 6, 24, vr, vi, ur, ui);
370 vr = y[ 8];
374 x[ 0] = ur + (vr>>1);
376 x[ 8] = ur - (vr>>1);
381 vr = y[13];
385 x[ 4] = ur + (vr>>1);
387 x[12] = ur - (vr>>1);
392 vr = y[24];
396 x[16] = ur + (vr>>1);
398 x[24] = ur - (vr>>1);
403 vr = y[29];
407 x[20] = ur + (vr>>1);
409 x[28] = ur - (vr>>1);
414 SUMDIFF_PIFOURTH(vi, vr, y[10], y[11])
415 //vr = fMultDiv2((y[11] + y[10]),W_PiFOURTH);
419 x[ 2] = (ur>>1) + vr;
421 x[10] = (ur>>1) - vr;
426 SUMDIFF_PIFOURTH(vr, vi, y[14], y[15])
429 x[ 6] = (ur>>1) + vr;
431 x[14] = (ur>>1) - vr;
436 SUMDIFF_PIFOURTH(vi, vr, y[26], y[27])
439 x[18] = (ur>>1) + vr;
441 x[26] = (ur>>1) - vr;
446 SUMDIFF_PIFOURTH(vr, vi, y[30], y[31])
449 x[22] = (ur>>1) + vr;
451 x[30] = (ur>>1) - vr;
456 vr = x[16];
460 x[ 0] = ur + (vr>>1);
462 x[16] = ur - (vr>>1);
468 vr = x[25];
471 x[ 8] = ur + (vr>>1);
473 x[24] = ur - (vr>>1);
478 cplxMultDiv2(&vi, &vr, x[19], x[18], fft16_w16[0]);
481 x[ 2] = (ur>>1) + vr;
483 x[18] = (ur>>1) - vr;
488 cplxMultDiv2(&vr, &vi, x[27], x[26], fft16_w16[0]);
491 x[10] = (ur>>1) + vr;
493 x[26] = (ur>>1) - vr;
498 SUMDIFF_PIFOURTH(vi, vr, x[20], x[21])
501 x[ 4] = (ur>>1) + vr;
503 x[20] = (ur>>1) - vr;
508 SUMDIFF_PIFOURTH(vr, vi, x[28], x[29])
511 x[12] = (ur>>1) + vr;
513 x[28] = (ur>>1) - vr;
518 cplxMultDiv2(&vi, &vr, x[23], x[22], fft16_w16[1]);
521 x[ 6] = (ur>>1) + vr;
523 x[22] = (ur>>1) - vr;
528 cplxMultDiv2(&vr, &vi, x[31], x[30], fft16_w16[1]);
531 x[14] = (ur>>1) + vr;
533 x[30] = (ur>>1) - vr;
550 FIXP_DBL vr,vi,ur,ui;
560 vr = (x[ 0] + x[32])>>1; /* Re A + Re B */
565 y[ 0] = vr + (vi>>1); /* Re A' = ReA + ReB +ReC + ReD */
566 y[ 4] = vr - (vi>>1); /* Re C' = -(ReC+ReD) + (ReA+ReB) */
570 vr -= x[32]; /* Re A - Re B */
575 y[ 2] = ui + vr; /* Re B' = Im C - Im D + Re A - Re B */
576 y[ 6] = vr - ui; /* Re D' = -Im C + Im D + Re A - Re B */
581 vr = (x[ 8] + x[40])>>1; /* Re A + Re B */
586 y[ 8] = vr + (vi>>1); /* Re A' = ReA + ReB +ReC + ReD */
587 y[12] = vr - (vi>>1); /* Re C' = -(ReC+ReD) + (ReA+ReB) */
591 vr -= x[40]; /* Re A - Re B */
596 y[10] = ui + vr; /* Re B' = Im C - Im D + Re A - Re B */
597 y[14] = vr - ui; /* Re D' = -Im C + Im D + Re A - Re B */
602 vr = (x[ 4] + x[36])>>1; /* Re A + Re B */
607 y[16] = vr + (vi>>1); /* Re A' = ReA + ReB +ReC + ReD */
608 y[20] = vr - (vi>>1); /* Re C' = -(ReC+ReD) + (ReA+ReB) */
612 vr -= x[36]; /* Re A - Re B */
617 y[18] = ui + vr; /* Re B' = Im C - Im D + Re A - Re B */
618 y[22] = vr - ui; /* Re D' = -Im C + Im D + Re A - Re B */
623 vr = (x[12] + x[44])>>1; /* Re A + Re B */
628 y[24] = vr + (vi>>1); /* Re A' = ReA + ReB +ReC + ReD */
629 y[28] = vr - (vi>>1); /* Re C' = -(ReC+ReD) + (ReA+ReB) */
633 vr -= x[44]; /* Re A - Re B */
638 y[26] = ui + vr; /* Re B' = Im C - Im D + Re A - Re B */
639 y[30] = vr - ui; /* Re D' = -Im C + Im D + Re A - Re B */
644 vr = (x[ 2] + x[34])>>1; /* Re A + Re B */
649 y[32] = vr + (vi>>1); /* Re A' = ReA + ReB +ReC + ReD */
650 y[36] = vr - (vi>>1); /* Re C' = -(ReC+ReD) + (ReA+ReB) */
654 vr -= x[34]; /* Re A - Re B */
659 y[34] = ui + vr; /* Re B' = Im C - Im D + Re A - Re B */
660 y[38] = vr - ui; /* Re D' = -Im C + Im D + Re A - Re B */
665 vr = (x[10] + x[42])>>1; /* Re A + Re B */
670 y[40] = vr + (vi>>1); /* Re A' = ReA + ReB +ReC + ReD */
671 y[44] = vr - (vi>>1); /* Re C' = -(ReC+ReD) + (ReA+ReB) */
675 vr -= x[42]; /* Re A - Re B */
680 y[42] = ui + vr; /* Re B' = Im C - Im D + Re A - Re B */
681 y[46] = vr - ui; /* Re D' = -Im C + Im D + Re A - Re B */
686 vr = (x[ 6] + x[38])>>1; /* Re A + Re B */
691 y[48] = vr + (vi>>1); /* Re A' = ReA + ReB +ReC + ReD */
692 y[52] = vr - (vi>>1); /* Re C' = -(ReC+ReD) + (ReA+ReB) */
696 vr -= x[38]; /* Re A - Re B */
701 y[50] = ui + vr; /* Re B' = Im C - Im D + Re A - Re B */
702 y[54] = vr - ui; /* Re D' = -Im C + Im D + Re A - Re B */
707 vr = (x[14] + x[46])>>1; /* Re A + Re B */
712 y[56] = vr + (vi>>1); /* Re A' = ReA + ReB +ReC + ReD */
713 y[60] = vr - (vi>>1); /* Re C' = -(ReC+ReD) + (ReA+ReB) */
717 vr -= x[46]; /* Re A - Re B */
722 y[58] = ui + vr; /* Re B' = Im C - Im D + Re A - Re B */
723 y[62] = vr - ui; /* Re D' = -Im C + Im D + Re A - Re B */
734 vr = yt[8];
738 xt[ 0] = ur + (vr>>1);
740 xt[ 8] = ur - (vr>>1);
743 vr = yt[13];
747 xt[ 4] = ur + (vr>>1);
749 xt[12] = ur - (vr>>1);
752 SUMDIFF_PIFOURTH(vi, vr, yt[10], yt[11])
755 xt[ 2] = (ur>>1) + vr;
757 xt[10] = (ur>>1) - vr;
760 SUMDIFF_PIFOURTH(vr, vi, yt[14], yt[15])
764 xt[ 6] = (ur>>1) + vr;
766 xt[14] = (ur>>1) - vr;
772 vr = x[16];
776 x[ 0] = ur + (vr>>1);
778 x[16] = ur - (vr>>1);
782 vr = x[25];
785 x[ 8] = ur + (vr>>1);
787 x[24] = ur - (vr>>1);
790 vr = x[48];
794 x[32] = ur + (vr>>1);
796 x[48] = ur - (vr>>1);
800 vr = x[57];
803 x[40] = ur + (vr>>1);
805 x[56] = ur - (vr>>1);
808 cplxMultDiv2(&vi, &vr, x[19], x[18], fft32_w32[0]);
811 x[ 2] = (ur>>1) + vr;
813 x[18] = (ur>>1) - vr;
816 cplxMultDiv2(&vr, &vi, x[27], x[26], fft32_w32[0]);
819 x[10] = (ur>>1) + vr;
821 x[26] = (ur>>1) - vr;
824 cplxMultDiv2(&vi, &vr, x[51], x[50], fft32_w32[0]);
827 x[34] = (ur>>1) + vr;
829 x[50] = (ur>>1) - vr;
832 cplxMultDiv2(&vr, &vi, x[59], x[58], fft32_w32[0]);
835 x[42] = (ur>>1) + vr;
837 x[58] = (ur>>1) - vr;
840 SUMDIFF_PIFOURTH(vi, vr, x[20], x[21])
843 x[ 4] = (ur>>1) + vr;
845 x[20] = (ur>>1) - vr;
848 SUMDIFF_PIFOURTH(vr, vi, x[28], x[29])
851 x[12] = (ur>>1) + vr;
853 x[28] = (ur>>1) - vr;
856 SUMDIFF_PIFOURTH(vi, vr, x[52], x[53])
859 x[36] = (ur>>1) + vr;
861 x[52] = (ur>>1) - vr;
864 SUMDIFF_PIFOURTH(vr, vi, x[60], x[61])
867 x[44] = (ur>>1) + vr;
869 x[60] = (ur>>1) - vr;
873 cplxMultDiv2(&vi, &vr, x[23], x[22], fft32_w32[1]);
876 x[ 6] = (ur>>1) + vr;
878 x[22] = (ur>>1) - vr;
881 cplxMultDiv2(&vr, &vi, x[31], x[30], fft32_w32[1]);
884 x[14] = (ur>>1) + vr;
886 x[30] = (ur>>1) - vr;
889 cplxMultDiv2(&vi, &vr, x[55], x[54], fft32_w32[1]);
892 x[38] = (ur>>1) + vr;
894 x[54] = (ur>>1) - vr;
897 cplxMultDiv2(&vr, &vi, x[63], x[62], fft32_w32[1]);
901 x[46] = (ur>>1) + vr;
903 x[62] = (ur>>1) - vr;
906 vr = x[32];
910 x[ 0] = ur + (vr>>1);
912 x[32] = ur - (vr>>1);
916 vr = x[49];
919 x[16] = ur + (vr>>1);
921 x[48] = ur - (vr>>1);
924 cplxMultDiv2(&vi, &vr, x[35], x[34], fft32_w32[2]);
927 x[ 2] = (ur>>1) + vr;
929 x[34] = (ur>>1) - vr;
932 cplxMultDiv2(&vr, &vi, x[51], x[50], fft32_w32[2]);
935 x[18] = (ur>>1) + vr;
937 x[50] = (ur>>1) - vr;
940 cplxMultDiv2(&vi, &vr, x[37], x[36], fft32_w32[0]);
943 x[ 4] = (ur>>1) + vr;
945 x[36] = (ur>>1) - vr;
948 cplxMultDiv2(&vr, &vi, x[53], x[52], fft32_w32[0]);
951 x[20] = (ur>>1) + vr;
953 x[52] = (ur>>1) - vr;
956 cplxMultDiv2(&vi, &vr, x[39], x[38], fft32_w32[3]);
959 x[ 6] = (ur>>1) + vr;
961 x[38] = (ur>>1) - vr;
964 cplxMultDiv2(&vr, &vi, x[55], x[54], fft32_w32[3]);
967 x[22] = (ur>>1) + vr;
969 x[54] = (ur>>1) - vr;
972 SUMDIFF_PIFOURTH(vi, vr, x[40], x[41])
975 x[ 8] = (ur>>1) + vr;
977 x[40] = (ur>>1) - vr;
980 SUMDIFF_PIFOURTH(vr, vi, x[56], x[57])
983 x[24] = (ur>>1) + vr;
985 x[56] = (ur>>1) - vr;
988 cplxMultDiv2(&vi, &vr, x[43], x[42], fft32_w32[4]);
992 x[10] = (ur>>1) + vr;
994 x[42] = (ur>>1) - vr;
997 cplxMultDiv2(&vr, &vi, x[59], x[58], fft32_w32[4]);
1000 x[26] = (ur>>1) + vr;
1002 x[58] = (ur>>1) - vr;
1005 cplxMultDiv2(&vi, &vr, x[45], x[44], fft32_w32[1]);
1008 x[12] = (ur>>1) + vr;
1010 x[44] = (ur>>1) - vr;
1013 cplxMultDiv2(&vr, &vi, x[61], x[60], fft32_w32[1]);
1016 x[28] = (ur>>1) + vr;
1018 x[60] = (ur>>1) - vr;
1021 cplxMultDiv2(&vi, &vr, x[47], x[46], fft32_w32[5]);
1024 x[14] = (ur>>1) + vr;
1026 x[46] = (ur>>1) - vr;
1029 cplxMultDiv2(&vr, &vi, x[63], x[62], fft32_w32[5]);
1032 x[30] = (ur>>1) + vr;
1034 x[62] = (ur>>1) - vr;
