498     int j, k, l, m, mh, nw, nc;  in dfct()  local
523             k = m - j;  in dfct()
526             yr = a[k] - a[n - k];  in dfct()
527             yi = a[k] + a[n - k];  in dfct()
529             a[k] = yr;  in dfct()
531             t[k] = xi + yi;  in dfct()
562             k = 0;  in dfct()
564                 k += l << 2;  in dfct()
565                 a[k - l] = t[j] - t[j + 1];  in dfct()
566                 a[k + l] = t[j] + t[j + 1];  in dfct()
571                 k = m - j;  in dfct()
572                 t[j] = t[m + k] - t[m + j];  in dfct()
573                 t[k] = t[m + k] + t[m + j];  in dfct()
591     int j, k, l, m, mh, nw, nc;  in dfst()  local
611             k = m - j;  in dfst()
614             yr = a[k] + a[n - k];  in dfst()
615             yi = a[k] - a[n - k];  in dfst()
617             a[k] = yr;  in dfst()
619             t[k] = xi - yi;  in dfst()
651             k = 0;  in dfst()
653                 k += l << 2;  in dfst()
654                 a[k - l] = -t[j] - t[j + 1];  in dfst()
655                 a[k + l] = t[j] - t[j + 1];  in dfst()
660                 k = m - j;  in dfst()
661                 t[j] = t[m + k] + t[m + j];  in dfst()
662                 t[k] = t[m + k] - t[m + j];  in dfst()
729     int j, j1, k, k1, l, m, m2, ip[256];  in bitrv2()  local
745         for (k = 0; k < m; k++) {  in bitrv2()
746             for (j = 0; j < k; j++) {  in bitrv2()
747                 j1 = 2 * j + ip[k];  in bitrv2()
748                 k1 = 2 * k + ip[j];  in bitrv2()
788             j1 = 2 * k + m2 + ip[k];  in bitrv2()
800         for (k = 1; k < m; k++) {  in bitrv2()
801             for (j = 0; j < k; j++) {  in bitrv2()
802                 j1 = 2 * j + ip[k];  in bitrv2()
803                 k1 = 2 * k + ip[j];  in bitrv2()
830     int j, j1, k, k1, l, m, m2, ip[256];  in bitrv2conj()  local
846         for (k = 0; k < m; k++) {  in bitrv2conj()
847             for (j = 0; j < k; j++) {  in bitrv2conj()
848                 j1 = 2 * j + ip[k];  in bitrv2conj()
849                 k1 = 2 * k + ip[j];  in bitrv2conj()
889             k1 = 2 * k + ip[k];  in bitrv2conj()
907         for (k = 1; k < m; k++) {  in bitrv2conj()
908             for (j = 0; j < k; j++) {  in bitrv2conj()
909                 j1 = 2 * j + ip[k];  in bitrv2conj()
910                 k1 = 2 * k + ip[j];  in bitrv2conj()
930             k1 = 2 * k + ip[k];  in bitrv2conj()
1145     int j, j1, j2, j3, k, k1, k2, m, m2;  in cftmdl()  local
1199     for (k = m2; k < n; k += m2) {  in cftmdl()
1208         for (j = k; j < l + k; j += 2) {  in cftmdl()
1239         for (j = k + m; j < l + (k + m); j += 2) {  in cftmdl()
1272     int j, k, kk, ks, m;  in rftfsub()  local
1279         k = n - j;  in rftfsub()
1283         xr = a[j] - a[k];  in rftfsub()
1284         xi = a[j + 1] + a[k + 1];  in rftfsub()
1289         a[k] += yr;  in rftfsub()
1290         a[k + 1] -= yi;  in rftfsub()
1297     int j, k, kk, ks, m;  in rftbsub()  local
1305         k = n - j;  in rftbsub()
1309         xr = a[j] - a[k];  in rftbsub()
1310         xi = a[j + 1] + a[k + 1];  in rftbsub()
1315         a[k] += yr;  in rftbsub()
1316         a[k + 1] = yi - a[k + 1];  in rftbsub()
1324     int j, k, kk, ks, m;  in dctsub()  local
1331         k = n - j;  in dctsub()
1335         xr = wki * a[j] - wkr * a[k];  in dctsub()
1336         a[j] = wkr * a[j] + wki * a[k];  in dctsub()
1337         a[k] = xr;  in dctsub()
1345     int j, k, kk, ks, m;  in dstsub()  local
1352         k = n - j;  in dstsub()
1356         xr = wki * a[k] - wkr * a[j];  in dstsub()
1357         a[k] = wkr * a[k] + wki * a[j];  in dstsub()