| 2353 | } |
| 2354 | |
| 2355 | void cftrec4(int n, float *a) |
| 2356 | { |
| 2357 | int isplt, j, k, m; |
| 2358 | |
| 2359 | m = n; |
| 2360 | while (m > 512) { |
| 2361 | m >>= 2; |
| 2362 | cftmdl1(m, &a[n - m]); |
| 2363 | } |
| 2364 | cftleaf(m, 1, &a[n - m]); |
| 2365 | k = 0; |
| 2366 | for (j = n - m; j > 0; j -= m) { |
| 2367 | k++; |
| 2368 | isplt = cfttree(m, j, k, a); |
| 2369 | cftleaf(m, isplt, &a[j - m]); |
| 2370 | } |
| 2371 | } |
| 2372 | |
| 2373 | void cftfsub(int n, float *a) |
| 2374 | { |