905     float x0r, x0i, x1r, x1i, x2r, x2i, x3r, x3i;  in cftfsub()  local
922             x0i = a[j + 1] + a[j1 + 1];  in cftfsub()
930             a[j + 1] = x0i + x2i;  in cftfsub()
932             a[j2 + 1] = x0i - x2i;  in cftfsub()
942             x0i = a[j + 1] - a[j1 + 1];  in cftfsub()
946             a[j1 + 1] = x0i;  in cftfsub()
955     float x0r, x0i, x1r, x1i, x2r, x2i, x3r, x3i;  in cftbsub()  local
972             x0i = -a[j + 1] - a[j1 + 1];  in cftbsub()
980             a[j + 1] = x0i - x2i;  in cftbsub()
982             a[j2 + 1] = x0i + x2i;  in cftbsub()
992             x0i = -a[j + 1] + a[j1 + 1];  in cftbsub()
996             a[j1 + 1] = x0i;  in cftbsub()
1006     float x0r, x0i, x1r, x1i, x2r, x2i, x3r, x3i;  in cft1st()  local
1009     x0i = a[1] + a[3];  in cft1st()
1017     a[1] = x0i + x2i;  in cft1st()
1019     a[5] = x0i - x2i;  in cft1st()
1026     x0i = a[9] + a[11];  in cft1st()
1034     a[9] = x0i + x2i;  in cft1st()
1035     a[12] = x2i - x0i;  in cft1st()
1038     x0i = x1i + x3r;  in cft1st()
1039     a[10] = wk1r * (x0r - x0i);  in cft1st()
1040     a[11] = wk1r * (x0r + x0i);  in cft1st()
1042     x0i = x3r - x1i;  in cft1st()
1043     a[14] = wk1r * (x0i - x0r);  in cft1st()
1044     a[15] = wk1r * (x0i + x0r);  in cft1st()
1056         x0i = a[j + 1] + a[j + 3];  in cft1st()
1064         a[j + 1] = x0i + x2i;  in cft1st()
1066         x0i -= x2i;  in cft1st()
1067         a[j + 4] = wk2r * x0r - wk2i * x0i;  in cft1st()
1068         a[j + 5] = wk2r * x0i + wk2i * x0r;  in cft1st()
1070         x0i = x1i + x3r;  in cft1st()
1071         a[j + 2] = wk1r * x0r - wk1i * x0i;  in cft1st()
1072         a[j + 3] = wk1r * x0i + wk1i * x0r;  in cft1st()
1074         x0i = x1i - x3r;  in cft1st()
1075         a[j + 6] = wk3r * x0r - wk3i * x0i;  in cft1st()
1076         a[j + 7] = wk3r * x0i + wk3i * x0r;  in cft1st()
1082         x0i = a[j + 9] + a[j + 11];  in cft1st()
1090         a[j + 9] = x0i + x2i;  in cft1st()
1092         x0i -= x2i;  in cft1st()
1093         a[j + 12] = -wk2i * x0r - wk2r * x0i;  in cft1st()
1094         a[j + 13] = -wk2i * x0i + wk2r * x0r;  in cft1st()
1096         x0i = x1i + x3r;  in cft1st()
1097         a[j + 10] = wk1r * x0r - wk1i * x0i;  in cft1st()
1098         a[j + 11] = wk1r * x0i + wk1i * x0r;  in cft1st()
1100         x0i = x1i - x3r;  in cft1st()
1101         a[j + 14] = wk3r * x0r - wk3i * x0i;  in cft1st()
1102         a[j + 15] = wk3r * x0i + wk3i * x0r;  in cft1st()
1111     float x0r, x0i, x1r, x1i, x2r, x2i, x3r, x3i;  in cftmdl()  local
1119         x0i = a[j + 1] + a[j1 + 1];  in cftmdl()
1127         a[j + 1] = x0i + x2i;  in cftmdl()
1129         a[j2 + 1] = x0i - x2i;  in cftmdl()
1141         x0i = a[j + 1] + a[j1 + 1];  in cftmdl()
1149         a[j + 1] = x0i + x2i;  in cftmdl()
1150         a[j2] = x2i - x0i;  in cftmdl()
1153         x0i = x1i + x3r;  in cftmdl()
1154         a[j1] = wk1r * (x0r - x0i);  in cftmdl()
1155         a[j1 + 1] = wk1r * (x0r + x0i);  in cftmdl()
1157         x0i = x3r - x1i;  in cftmdl()
1158         a[j3] = wk1r * (x0i - x0r);  in cftmdl()
1159         a[j3 + 1] = wk1r * (x0i + x0r);  in cftmdl()
1177             x0i = a[j + 1] + a[j1 + 1];  in cftmdl()
1185             a[j + 1] = x0i + x2i;  in cftmdl()
1187             x0i -= x2i;  in cftmdl()
1188             a[j2] = wk2r * x0r - wk2i * x0i;  in cftmdl()
1189             a[j2 + 1] = wk2r * x0i + wk2i * x0r;  in cftmdl()
1191             x0i = x1i + x3r;  in cftmdl()
1192             a[j1] = wk1r * x0r - wk1i * x0i;  in cftmdl()
1193             a[j1 + 1] = wk1r * x0i + wk1i * x0r;  in cftmdl()
1195             x0i = x1i - x3r;  in cftmdl()
1196             a[j3] = wk3r * x0r - wk3i * x0i;  in cftmdl()
1197             a[j3 + 1] = wk3r * x0i + wk3i * x0r;  in cftmdl()
1208             x0i = a[j + 1] + a[j1 + 1];  in cftmdl()
1216             a[j + 1] = x0i + x2i;  in cftmdl()
1218             x0i -= x2i;  in cftmdl()
1219             a[j2] = -wk2i * x0r - wk2r * x0i;  in cftmdl()
1220             a[j2 + 1] = -wk2i * x0i + wk2r * x0r;  in cftmdl()
1222             x0i = x1i + x3r;  in cftmdl()
1223             a[j1] = wk1r * x0r - wk1i * x0i;  in cftmdl()
1224             a[j1 + 1] = wk1r * x0i + wk1i * x0r;  in cftmdl()
1226             x0i = x1i - x3r;  in cftmdl()
1227             a[j3] = wk3r * x0r - wk3i * x0i;  in cftmdl()
1228             a[j3 + 1] = wk3r * x0i + wk3i * x0r;  in cftmdl()