904     int j, j1, j2, j3, l;  in cftfsub()  local
919             j2 = j1 + l;  in cftfsub()
920             j3 = j2 + l;  in cftfsub()
925             x2r = a[j2] + a[j3];  in cftfsub()
926             x2i = a[j2 + 1] + a[j3 + 1];  in cftfsub()
927             x3r = a[j2] - a[j3];  in cftfsub()
928             x3i = a[j2 + 1] - a[j3 + 1];  in cftfsub()
931             a[j2] = x0r - x2r;  in cftfsub()
932             a[j2 + 1] = x0i - x2i;  in cftfsub()
954     int j, j1, j2, j3, l;  in cftbsub()  local
969             j2 = j1 + l;  in cftbsub()
970             j3 = j2 + l;  in cftbsub()
975             x2r = a[j2] + a[j3];  in cftbsub()
976             x2i = a[j2 + 1] + a[j3 + 1];  in cftbsub()
977             x3r = a[j2] - a[j3];  in cftbsub()
978             x3i = a[j2 + 1] - a[j3 + 1];  in cftbsub()
981             a[j2] = x0r - x2r;  in cftbsub()
982             a[j2 + 1] = x0i + x2i;  in cftbsub()
1109     int j, j1, j2, j3, k, k1, k2, m, m2;  in cftmdl()  local
1116         j2 = j1 + l;  in cftmdl()
1117         j3 = j2 + l;  in cftmdl()
1122         x2r = a[j2] + a[j3];  in cftmdl()
1123         x2i = a[j2 + 1] + a[j3 + 1];  in cftmdl()
1124         x3r = a[j2] - a[j3];  in cftmdl()
1125         x3i = a[j2 + 1] - a[j3 + 1];  in cftmdl()
1128         a[j2] = x0r - x2r;  in cftmdl()
1129         a[j2 + 1] = x0i - x2i;  in cftmdl()
1138         j2 = j1 + l;  in cftmdl()
1139         j3 = j2 + l;  in cftmdl()
1144         x2r = a[j2] + a[j3];  in cftmdl()
1145         x2i = a[j2 + 1] + a[j3 + 1];  in cftmdl()
1146         x3r = a[j2] - a[j3];  in cftmdl()
1147         x3i = a[j2 + 1] - a[j3 + 1];  in cftmdl()
1150         a[j2] = x2i - x0i;  in cftmdl()
1151         a[j2 + 1] = x0r - x2r;  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()
1188             a[j2] = wk2r * x0r - wk2i * x0i;  in cftmdl()
1189             a[j2 + 1] = wk2r * x0i + wk2i * x0r;  in cftmdl()
1205             j2 = j1 + l;  in cftmdl()
1206             j3 = j2 + l;  in cftmdl()
1211             x2r = a[j2] + a[j3];  in cftmdl()
1212             x2i = a[j2 + 1] + a[j3 + 1];  in cftmdl()
1213             x3r = a[j2] - a[j3];  in cftmdl()
1214             x3i = a[j2 + 1] - a[j3 + 1];  in cftmdl()
1219             a[j2] = -wk2i * x0r - wk2r * x0i;  in cftmdl()
1220             a[j2 + 1] = -wk2i * x0i + wk2r * x0r;  in cftmdl()