940     int j, j1, j2, j3, l;  in cftfsub()  local
955             j2 = j1 + l;  in cftfsub()
956             j3 = j2 + l;  in cftfsub()
961             x2r = a[j2] + a[j3];  in cftfsub()
962             x2i = a[j2 + 1] + a[j3 + 1];  in cftfsub()
963             x3r = a[j2] - a[j3];  in cftfsub()
964             x3i = a[j2 + 1] - a[j3 + 1];  in cftfsub()
967             a[j2] = x0r - x2r;  in cftfsub()
968             a[j2 + 1] = x0i - x2i;  in cftfsub()
990     int j, j1, j2, j3, l;  in cftbsub()  local
1005             j2 = j1 + l;  in cftbsub()
1006             j3 = j2 + l;  in cftbsub()
1011             x2r = a[j2] + a[j3];  in cftbsub()
1012             x2i = a[j2 + 1] + a[j3 + 1];  in cftbsub()
1013             x3r = a[j2] - a[j3];  in cftbsub()
1014             x3i = a[j2 + 1] - a[j3 + 1];  in cftbsub()
1017             a[j2] = x0r - x2r;  in cftbsub()
1018             a[j2 + 1] = x0i + x2i;  in cftbsub()
1145     int j, j1, j2, j3, k, k1, k2, m, m2;  in cftmdl()  local
1152         j2 = j1 + l;  in cftmdl()
1153         j3 = j2 + l;  in cftmdl()
1158         x2r = a[j2] + a[j3];  in cftmdl()
1159         x2i = a[j2 + 1] + a[j3 + 1];  in cftmdl()
1160         x3r = a[j2] - a[j3];  in cftmdl()
1161         x3i = a[j2 + 1] - a[j3 + 1];  in cftmdl()
1164         a[j2] = x0r - x2r;  in cftmdl()
1165         a[j2 + 1] = x0i - x2i;  in cftmdl()
1174         j2 = j1 + l;  in cftmdl()
1175         j3 = j2 + l;  in cftmdl()
1180         x2r = a[j2] + a[j3];  in cftmdl()
1181         x2i = a[j2 + 1] + a[j3 + 1];  in cftmdl()
1182         x3r = a[j2] - a[j3];  in cftmdl()
1183         x3i = a[j2 + 1] - a[j3 + 1];  in cftmdl()
1186         a[j2] = x2i - x0i;  in cftmdl()
1187         a[j2 + 1] = x0r - x2r;  in cftmdl()
1210             j2 = j1 + l;  in cftmdl()
1211             j3 = j2 + l;  in cftmdl()
1216             x2r = a[j2] + a[j3];  in cftmdl()
1217             x2i = a[j2 + 1] + a[j3 + 1];  in cftmdl()
1218             x3r = a[j2] - a[j3];  in cftmdl()
1219             x3i = a[j2 + 1] - a[j3 + 1];  in cftmdl()
1224             a[j2] = wk2r * x0r - wk2i * x0i;  in cftmdl()
1225             a[j2 + 1] = wk2r * x0i + wk2i * x0r;  in cftmdl()
1241             j2 = j1 + l;  in cftmdl()
1242             j3 = j2 + l;  in cftmdl()
1247             x2r = a[j2] + a[j3];  in cftmdl()
1248             x2i = a[j2 + 1] + a[j3 + 1];  in cftmdl()
1249             x3r = a[j2] - a[j3];  in cftmdl()
1250             x3i = a[j2 + 1] - a[j3 + 1];  in cftmdl()
1255             a[j2] = -wk2i * x0r - wk2r * x0i;  in cftmdl()
1256             a[j2 + 1] = -wk2i * x0i + wk2r * x0r;  in cftmdl()