941     double x0r, x0i, x1r, x1i, x2r, x2i, x3r, x3i;  in cftfsub()  local
958             x0i = a[j + 1] + a[j1 + 1];  in cftfsub()
966             a[j + 1] = x0i + x2i;  in cftfsub()
968             a[j2 + 1] = x0i - x2i;  in cftfsub()
978             x0i = a[j + 1] - a[j1 + 1];  in cftfsub()
982             a[j1 + 1] = x0i;  in cftfsub()
991     double x0r, x0i, x1r, x1i, x2r, x2i, x3r, x3i;  in cftbsub()  local
1008             x0i = -a[j + 1] - a[j1 + 1];  in cftbsub()
1016             a[j + 1] = x0i - x2i;  in cftbsub()
1018             a[j2 + 1] = x0i + x2i;  in cftbsub()
1028             x0i = -a[j + 1] + a[j1 + 1];  in cftbsub()
1032             a[j1 + 1] = x0i;  in cftbsub()
1042     double x0r, x0i, x1r, x1i, x2r, x2i, x3r, x3i;  in cft1st()  local
1045     x0i = a[1] + a[3];  in cft1st()
1053     a[1] = x0i + x2i;  in cft1st()
1055     a[5] = x0i - x2i;  in cft1st()
1062     x0i = a[9] + a[11];  in cft1st()
1070     a[9] = x0i + x2i;  in cft1st()
1071     a[12] = x2i - x0i;  in cft1st()
1074     x0i = x1i + x3r;  in cft1st()
1075     a[10] = wk1r * (x0r - x0i);  in cft1st()
1076     a[11] = wk1r * (x0r + x0i);  in cft1st()
1078     x0i = x3r - x1i;  in cft1st()
1079     a[14] = wk1r * (x0i - x0r);  in cft1st()
1080     a[15] = wk1r * (x0i + x0r);  in cft1st()
1092         x0i = a[j + 1] + a[j + 3];  in cft1st()
1100         a[j + 1] = x0i + x2i;  in cft1st()
1102         x0i -= x2i;  in cft1st()
1103         a[j + 4] = wk2r * x0r - wk2i * x0i;  in cft1st()
1104         a[j + 5] = wk2r * x0i + wk2i * x0r;  in cft1st()
1106         x0i = x1i + x3r;  in cft1st()
1107         a[j + 2] = wk1r * x0r - wk1i * x0i;  in cft1st()
1108         a[j + 3] = wk1r * x0i + wk1i * x0r;  in cft1st()
1110         x0i = x1i - x3r;  in cft1st()
1111         a[j + 6] = wk3r * x0r - wk3i * x0i;  in cft1st()
1112         a[j + 7] = wk3r * x0i + wk3i * x0r;  in cft1st()
1118         x0i = a[j + 9] + a[j + 11];  in cft1st()
1126         a[j + 9] = x0i + x2i;  in cft1st()
1128         x0i -= x2i;  in cft1st()
1129         a[j + 12] = -wk2i * x0r - wk2r * x0i;  in cft1st()
1130         a[j + 13] = -wk2i * x0i + wk2r * x0r;  in cft1st()
1132         x0i = x1i + x3r;  in cft1st()
1133         a[j + 10] = wk1r * x0r - wk1i * x0i;  in cft1st()
1134         a[j + 11] = wk1r * x0i + wk1i * x0r;  in cft1st()
1136         x0i = x1i - x3r;  in cft1st()
1137         a[j + 14] = wk3r * x0r - wk3i * x0i;  in cft1st()
1138         a[j + 15] = wk3r * x0i + wk3i * x0r;  in cft1st()
1147     double x0r, x0i, x1r, x1i, x2r, x2i, x3r, x3i;  in cftmdl()  local
1155         x0i = a[j + 1] + a[j1 + 1];  in cftmdl()
1163         a[j + 1] = x0i + x2i;  in cftmdl()
1165         a[j2 + 1] = x0i - x2i;  in cftmdl()
1177         x0i = a[j + 1] + a[j1 + 1];  in cftmdl()
1185         a[j + 1] = x0i + x2i;  in cftmdl()
1186         a[j2] = x2i - x0i;  in cftmdl()
1189         x0i = x1i + x3r;  in cftmdl()
1190         a[j1] = wk1r * (x0r - x0i);  in cftmdl()
1191         a[j1 + 1] = wk1r * (x0r + x0i);  in cftmdl()
1193         x0i = x3r - x1i;  in cftmdl()
1194         a[j3] = wk1r * (x0i - x0r);  in cftmdl()
1195         a[j3 + 1] = wk1r * (x0i + x0r);  in cftmdl()
1213             x0i = a[j + 1] + a[j1 + 1];  in cftmdl()
1221             a[j + 1] = x0i + x2i;  in cftmdl()
1223             x0i -= x2i;  in cftmdl()
1224             a[j2] = wk2r * x0r - wk2i * x0i;  in cftmdl()
1225             a[j2 + 1] = wk2r * x0i + wk2i * x0r;  in cftmdl()
1227             x0i = x1i + x3r;  in cftmdl()
1228             a[j1] = wk1r * x0r - wk1i * x0i;  in cftmdl()
1229             a[j1 + 1] = wk1r * x0i + wk1i * x0r;  in cftmdl()
1231             x0i = x1i - x3r;  in cftmdl()
1232             a[j3] = wk3r * x0r - wk3i * x0i;  in cftmdl()
1233             a[j3 + 1] = wk3r * x0i + wk3i * x0r;  in cftmdl()
1244             x0i = a[j + 1] + a[j1 + 1];  in cftmdl()
1252             a[j + 1] = x0i + x2i;  in cftmdl()
1254             x0i -= x2i;  in cftmdl()
1255             a[j2] = -wk2i * x0r - wk2r * x0i;  in cftmdl()
1256             a[j2 + 1] = -wk2i * x0i + wk2r * x0r;  in cftmdl()
1258             x0i = x1i + x3r;  in cftmdl()
1259             a[j1] = wk1r * x0r - wk1i * x0i;  in cftmdl()
1260             a[j1 + 1] = wk1r * x0i + wk1i * x0r;  in cftmdl()
1262             x0i = x1i - x3r;  in cftmdl()
1263             a[j3] = wk3r * x0r - wk3i * x0i;  in cftmdl()
1264             a[j3 + 1] = wk3r * x0i + wk3i * x0r;  in cftmdl()