libs/mapm/mapm_fft.c

388     int j0, k0, j1, k1, l, m, i, j, k;
410                 j1 = j + m;
412                 xr = a[j1];
413                 xi = a[j1 + 1];
416                 a[j1] = yr;
417                 a[j1 + 1] = yi;
420                 j1 += m;
422                 xr = a[j1];
423                 xi = a[j1 + 1];
426                 a[j1] = yr;
427                 a[j1 + 1] = yi;
430                 j1 += m;
432                 xr = a[j1];
433                 xi = a[j1 + 1];
436                 a[j1] = yr;
437                 a[j1 + 1] = yi;
442             j1 = j0 + k0 + m;
443             k1 = j1 + m;
444             xr = a[j1];
445             xi = a[j1 + 1];
448             a[j1] = yr;
449             a[j1 + 1] = yi;
468                 j1 = j + m;
470                 xr = a[j1];
471                 xi = a[j1 + 1];
474                 a[j1] = yr;
475                 a[j1 + 1] = yi;
488     int j, j1, j2, j3, l;
502             j1 = j + l;
503             j2 = j1 + l;
505             x0r = a[j] + a[j1];
506             x0i = a[j + 1] + a[j1 + 1];
507             x1r = a[j] - a[j1];
508             x1i = a[j + 1] - a[j1 + 1];
517             a[j1] = x1r - x3i;
518             a[j1 + 1] = x1i + x3r;
524             j1 = j + l;
525             x0r = a[j] - a[j1];
526             x0i = a[j + 1] - a[j1 + 1];
527             a[j] += a[j1];
528             a[j + 1] += a[j1 + 1];
529             a[j1] = x0r;
530             a[j1 + 1] = x0i;
539     int j, j1, j2, j3, l;
553             j1 = j + l;
554             j2 = j1 + l;
556             x0r = a[j] + a[j1];
557             x0i = -a[j + 1] - a[j1 + 1];
558             x1r = a[j] - a[j1];
559             x1i = -a[j + 1] + a[j1 + 1];
568             a[j1] = x1r - x3i;
569             a[j1 + 1] = x1i - x3r;
575             j1 = j + l;
576             x0r = a[j] - a[j1];
577             x0i = -a[j + 1] + a[j1 + 1];
578             a[j] += a[j1];
579             a[j + 1] = -a[j + 1] - a[j1 + 1];
580             a[j1] = x0r;
581             a[j1 + 1] = x0i;
697     int j, j1, j2, j3, k, kj, kr, m, m2;
703         j1 = j + l;
704         j2 = j1 + l;
706         x0r = a[j] + a[j1];
707         x0i = a[j + 1] + a[j1 + 1];
708         x1r = a[j] - a[j1];
709         x1i = a[j + 1] - a[j1 + 1];
718         a[j1] = x1r - x3i;
719         a[j1 + 1] = x1i + x3r;
725         j1 = j + l;
726         j2 = j1 + l;
728         x0r = a[j] + a[j1];
729         x0i = a[j + 1] + a[j1 + 1];
730         x1r = a[j] - a[j1];
731         x1i = a[j + 1] - a[j1 + 1];
742         a[j1] = wn4r * (x0r - x0i);
743         a[j1 + 1] = wn4r * (x0r + x0i);
761             j1 = j + l;
762             j2 = j1 + l;
764             x0r = a[j] + a[j1];
765             x0i = a[j + 1] + a[j1 + 1];
766             x1r = a[j] - a[j1];
767             x1i = a[j + 1] - a[j1 + 1];
780             a[j1] = wk1r * x0r - wk1i * x0i;
781             a[j1 + 1] = wk1r * x0i + wk1i * x0r;
793             j1 = j + l;
794             j2 = j1 + l;
796             x0r = a[j] + a[j1];
797             x0i = a[j + 1] + a[j1 + 1];
798             x1r = a[j] - a[j1];
799             x1i = a[j + 1] - a[j1 + 1];
812             a[j1] = wk1r * x0r - wk1i * x0i;
813             a[j1 + 1] = wk1r * x0i + wk1i * x0r;