Lines Matching refs:vi
337 #define SUMDIFF_EIGTH(x, y, ix, iy, vr, vi, ur, ui) \
339 vi = (x[ 8 + ix]>>1) + (x[24 + ix]>>1); /* Re C + Re D */ \
342 y[ 0 + iy] = vr + vi; /* Re A' = ReA + ReB +ReC + ReD */ \
343 y[ 4 + iy] = vr - vi; /* Re C' = -(ReC+ReD) + (ReA+ReB) */ \
347 vi = vi - x[24 + ix]; /* Re C - Re D */ \
352 y[ 3 + iy] = ur - vi; /* Im B'= -Re C + Re D + Im A - Im B */ \
353 y[ 7 + iy] = vi + ur; /* Im D'= Re C - Re D + Im A - Im 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);
371 vi = y[ 9];
375 x[ 1] = ui + (vi>>1);
377 x[ 9] = ui - (vi>>1);
382 vi = y[12];
386 x[ 5] = ui - (vi>>1);
388 x[13] = ui + (vi>>1);
393 vi = y[25];
397 x[17] = ui + (vi>>1);
399 x[25] = ui - (vi>>1);
404 vi = y[28];
408 x[21] = ui - (vi>>1);
410 x[29] = ui + (vi>>1);
414 SUMDIFF_PIFOURTH(vi, vr, y[10], y[11])
416 //vi = fMultDiv2((y[11] - y[10]),W_PiFOURTH);
420 x[ 3] = (ui>>1) + vi;
422 x[11] = (ui>>1) - vi;
426 SUMDIFF_PIFOURTH(vr, vi, y[14], y[15])
430 x[ 7] = (ui>>1) - vi;
432 x[15] = (ui>>1) + vi;
436 SUMDIFF_PIFOURTH(vi, vr, y[26], y[27])
440 x[19] = (ui>>1) + vi;
442 x[27] = (ui>>1) - vi;
446 SUMDIFF_PIFOURTH(vr, vi, y[30], y[31])
450 x[23] = (ui>>1) - vi;
452 x[31] = (ui>>1) + vi;
457 vi = x[17];
461 x[ 1] = ui + (vi>>1);
463 x[17] = ui - (vi>>1);
467 vi = x[24];
472 x[ 9] = ui - (vi>>1);
474 x[25] = ui + (vi>>1);
478 cplxMultDiv2(&vi, &vr, x[19], x[18], fft16_w16[0]);
482 x[ 3] = (ui>>1) + vi;
484 x[19] = (ui>>1) - vi;
488 cplxMultDiv2(&vr, &vi, x[27], x[26], fft16_w16[0]);
492 x[11] = (ui>>1) - vi;
494 x[27] = (ui>>1) + vi;
498 SUMDIFF_PIFOURTH(vi, vr, x[20], x[21])
502 x[ 5] = (ui>>1) + vi;
504 x[21] = (ui>>1) - vi;
508 SUMDIFF_PIFOURTH(vr, vi, x[28], x[29])
512 x[13] = (ui>>1) - vi;
514 x[29] = (ui>>1) + vi;
518 cplxMultDiv2(&vi, &vr, x[23], x[22], fft16_w16[1]);
522 x[ 7] = (ui>>1) + vi;
524 x[23] = (ui>>1) - vi;
528 cplxMultDiv2(&vr, &vi, x[31], x[30], fft16_w16[1]);
532 x[15] = (ui>>1) - vi;
534 x[31] = (ui>>1) + vi;
550 FIXP_DBL vr,vi,ur,ui;
561 vi = (x[16] + x[48]); /* Re C + Re D */
565 y[ 0] = vr + (vi>>1); /* Re A' = ReA + ReB +ReC + ReD */
566 y[ 4] = vr - (vi>>1); /* Re C' = -(ReC+ReD) + (ReA+ReB) */
571 vi = (vi>>1) - x[48]; /* Re C - Re D */
577 y[ 3] = ur - vi; /* Im B'= -Re C + Re D + Im A - Im B */
578 y[ 7] = vi + ur; /* Im D'= Re C - Re D + Im A - Im B */
582 vi = (x[24] + x[56]); /* Re C + Re D */
586 y[ 8] = vr + (vi>>1); /* Re A' = ReA + ReB +ReC + ReD */
587 y[12] = vr - (vi>>1); /* Re C' = -(ReC+ReD) + (ReA+ReB) */
592 vi = (vi>>1) - x[56]; /* Re C - Re D */
598 y[11] = ur - vi; /* Im B'= -Re C + Re D + Im A - Im B */
599 y[15] = vi + ur; /* Im D'= Re C - Re D + Im A - Im B */
603 vi = (x[20] + x[52]); /* Re C + Re D */
607 y[16] = vr + (vi>>1); /* Re A' = ReA + ReB +ReC + ReD */
608 y[20] = vr - (vi>>1); /* Re C' = -(ReC+ReD) + (ReA+ReB) */
613 vi = (vi>>1) - x[52]; /* Re C - Re D */
619 y[19] = ur - vi; /* Im B'= -Re C + Re D + Im A - Im B */
620 y[23] = vi + ur; /* Im D'= Re C - Re D + Im A - Im B */
624 vi = (x[28] + x[60]); /* Re C + Re D */
628 y[24] = vr + (vi>>1); /* Re A' = ReA + ReB +ReC + ReD */
629 y[28] = vr - (vi>>1); /* Re C' = -(ReC+ReD) + (ReA+ReB) */
634 vi = (vi>>1) - x[60]; /* Re C - Re D */
640 y[27] = ur - vi; /* Im B'= -Re C + Re D + Im A - Im B */
641 y[31] = vi + ur; /* Im D'= Re C - Re D + Im A - Im B */
645 vi = (x[18] + x[50]); /* Re C + Re D */
649 y[32] = vr + (vi>>1); /* Re A' = ReA + ReB +ReC + ReD */
650 y[36] = vr - (vi>>1); /* Re C' = -(ReC+ReD) + (ReA+ReB) */
655 vi = (vi>>1) - x[50]; /* Re C - Re D */
661 y[35] = ur - vi; /* Im B'= -Re C + Re D + Im A - Im B */
662 y[39] = vi + ur; /* Im D'= Re C - Re D + Im A - Im B */
666 vi = (x[26] + x[58]); /* Re C + Re D */
670 y[40] = vr + (vi>>1); /* Re A' = ReA + ReB +ReC + ReD */
671 y[44] = vr - (vi>>1); /* Re C' = -(ReC+ReD) + (ReA+ReB) */
676 vi = (vi>>1) - x[58]; /* Re C - Re D */
682 y[43] = ur - vi; /* Im B'= -Re C + Re D + Im A - Im B */
683 y[47] = vi + ur; /* Im D'= Re C - Re D + Im A - Im B */
687 vi = (x[22] + x[54]); /* Re C + Re D */
691 y[48] = vr + (vi>>1); /* Re A' = ReA + ReB +ReC + ReD */
692 y[52] = vr - (vi>>1); /* Re C' = -(ReC+ReD) + (ReA+ReB) */
697 vi = (vi>>1) - x[54]; /* Re C - Re D */
703 y[51] = ur - vi; /* Im B'= -Re C + Re D + Im A - Im B */
704 y[55] = vi + ur; /* Im D'= Re C - Re D + Im A - Im B */
708 vi = (x[30] + x[62]); /* Re C + Re D */
712 y[56] = vr + (vi>>1); /* Re A' = ReA + ReB +ReC + ReD */
713 y[60] = vr - (vi>>1); /* Re C' = -(ReC+ReD) + (ReA+ReB) */
718 vi = (vi>>1) - x[62]; /* Re C - Re D */
724 y[59] = ur - vi; /* Im B'= -Re C + Re D + Im A - Im B */
725 y[63] = vi + ur; /* Im D'= Re C - Re D + Im A - Im B */
735 vi = yt[9];
739 xt[ 1] = ui + (vi>>1);
741 xt[ 9] = ui - (vi>>1);
744 vi = yt[12];
748 xt[ 5] = ui - (vi>>1);
750 xt[13] = ui + (vi>>1);
752 SUMDIFF_PIFOURTH(vi, vr, yt[10], yt[11])
756 xt[ 3] = (ui>>1) + vi;
758 xt[11] = (ui>>1) - vi;
760 SUMDIFF_PIFOURTH(vr, vi, yt[14], yt[15])
765 xt[ 7] = (ui>>1) - vi;
767 xt[15] = (ui>>1) + vi;
773 vi = x[17];
777 x[ 1] = ui + (vi>>1);
779 x[17] = ui - (vi>>1);
781 vi = x[24];
786 x[ 9] = ui - (vi>>1);
788 x[25] = ui + (vi>>1);
791 vi = x[49];
795 x[33] = ui + (vi>>1);
797 x[49] = ui - (vi>>1);
799 vi = x[56];
804 x[41] = ui - (vi>>1);
806 x[57] = ui + (vi>>1);
808 cplxMultDiv2(&vi, &vr, x[19], x[18], fft32_w32[0]);
812 x[ 3] = (ui>>1) + vi;
814 x[19] = (ui>>1) - vi;
816 cplxMultDiv2(&vr, &vi, x[27], x[26], fft32_w32[0]);
820 x[11] = (ui>>1) - vi;
822 x[27] = (ui>>1) + vi;
824 cplxMultDiv2(&vi, &vr, x[51], x[50], fft32_w32[0]);
828 x[35] = (ui>>1) + vi;
830 x[51] = (ui>>1) - vi;
832 cplxMultDiv2(&vr, &vi, x[59], x[58], fft32_w32[0]);
836 x[43] = (ui>>1) - vi;
838 x[59] = (ui>>1) + vi;
840 SUMDIFF_PIFOURTH(vi, vr, x[20], x[21])
844 x[ 5] = (ui>>1) + vi;
846 x[21] = (ui>>1) - vi;
848 SUMDIFF_PIFOURTH(vr, vi, x[28], x[29])
852 x[13] = (ui>>1) - vi;
854 x[29] = (ui>>1) + vi;
856 SUMDIFF_PIFOURTH(vi, vr, x[52], x[53])
860 x[37] = (ui>>1) + vi;
862 x[53] = (ui>>1) - vi;
864 SUMDIFF_PIFOURTH(vr, vi, x[60], x[61])
868 x[45] = (ui>>1) - vi;
870 x[61] = (ui>>1) + vi;
873 cplxMultDiv2(&vi, &vr, x[23], x[22], fft32_w32[1]);
877 x[ 7] = (ui>>1) + vi;
879 x[23] = (ui>>1) - vi;
881 cplxMultDiv2(&vr, &vi, x[31], x[30], fft32_w32[1]);
885 x[15] = (ui>>1) - vi;
887 x[31] = (ui>>1) + vi;
889 cplxMultDiv2(&vi, &vr, x[55], x[54], fft32_w32[1]);
893 x[39] = (ui>>1) + vi;
895 x[55] = (ui>>1) - vi;
897 cplxMultDiv2(&vr, &vi, x[63], x[62], fft32_w32[1]);
902 x[47] = (ui>>1) - vi;
904 x[63] = (ui>>1) + vi;
907 vi = x[33];
911 x[ 1] = ui + (vi>>1);
913 x[33] = ui - (vi>>1);
915 vi = x[48];
920 x[17] = ui - (vi>>1);
922 x[49] = ui + (vi>>1);
924 cplxMultDiv2(&vi, &vr, x[35], x[34], fft32_w32[2]);
928 x[ 3] = (ui>>1) + vi;
930 x[35] = (ui>>1) - vi;
932 cplxMultDiv2(&vr, &vi, x[51], x[50], fft32_w32[2]);
936 x[19] = (ui>>1) - vi;
938 x[51] = (ui>>1) + vi;
940 cplxMultDiv2(&vi, &vr, x[37], x[36], fft32_w32[0]);
944 x[ 5] = (ui>>1) + vi;
946 x[37] = (ui>>1) - vi;
948 cplxMultDiv2(&vr, &vi, x[53], x[52], fft32_w32[0]);
952 x[21] = (ui>>1) - vi;
954 x[53] = (ui>>1) + vi;
956 cplxMultDiv2(&vi, &vr, x[39], x[38], fft32_w32[3]);
960 x[ 7] = (ui>>1) + vi;
962 x[39] = (ui>>1) - vi;
964 cplxMultDiv2(&vr, &vi, x[55], x[54], fft32_w32[3]);
968 x[23] = (ui>>1) - vi;
970 x[55] = (ui>>1) + vi;
972 SUMDIFF_PIFOURTH(vi, vr, x[40], x[41])
976 x[ 9] = (ui>>1) + vi;
978 x[41] = (ui>>1) - vi;
980 SUMDIFF_PIFOURTH(vr, vi, x[56], x[57])
984 x[25] = (ui>>1) - vi;
986 x[57] = (ui>>1) + vi;
988 cplxMultDiv2(&vi, &vr, x[43], x[42], fft32_w32[4]);
993 x[11] = (ui>>1) + vi;
995 x[43] = (ui>>1) - vi;
997 cplxMultDiv2(&vr, &vi, x[59], x[58], fft32_w32[4]);
1001 x[27] = (ui>>1) - vi;
1003 x[59] = (ui>>1) + vi;
1005 cplxMultDiv2(&vi, &vr, x[45], x[44], fft32_w32[1]);
1009 x[13] = (ui>>1) + vi;
1011 x[45] = (ui>>1) - vi;
1013 cplxMultDiv2(&vr, &vi, x[61], x[60], fft32_w32[1]);
1017 x[29] = (ui>>1) - vi;
1019 x[61] = (ui>>1) + vi;
1021 cplxMultDiv2(&vi, &vr, x[47], x[46], fft32_w32[5]);
1025 x[15] = (ui>>1) + vi;
1027 x[47] = (ui>>1) - vi;
1029 cplxMultDiv2(&vr, &vi, x[63], x[62], fft32_w32[5]);
1033 x[31] = (ui>>1) - vi;
1035 x[63] = (ui>>1) + vi;