Chris@82: TODO before FFTW-$2\pi$: Chris@82: Chris@82: * figure out how to autodetect NEON at runtime Chris@82: Chris@82: * figure out the arm cycle counter business Chris@82: Chris@82: * Wisdom: make it clear that it is specific to the exact fftw version Chris@82: and configuration. Report error codes when reading wisdom. Maybe Chris@82: have multiple system wisdom files, one per version? Chris@82: Chris@82: * DCT/DST codelets? which kinds? Chris@82: Chris@82: * investigate the addition-chain trig computation Chris@82: Chris@82: * I can't believe that there isn't a closed form for the omega Chris@82: array in Rader. Chris@82: Chris@82: * convolution problem type(s) Chris@82: Chris@82: * Explore the idea of having n < 0 in tensors, possibly to mean Chris@82: inverse DFT. Chris@82: Chris@82: * better estimator: possibly, let "other" cost be coef * n, where Chris@82: coef is a per-solver constant determined via some big numerical Chris@82: optimization/fit. Chris@82: Chris@82: * vector radix, multidimensional codelets Chris@82: Chris@82: * it may be a good idea to unify all those little loops that do Chris@82: copying, (X[i], X[n-i]) <- (X[i] + X[n-i], X[i] - X[n-i]), Chris@82: and multiplication of vectors by twiddle factors. Chris@82: Chris@82: * Pruned FFTs (basically, a vecloop that skips zeros). Chris@82: Chris@82: * Try FFTPACK-style back-and-forth (Stockham) FFT. (We tried this a Chris@82: few years ago and it was slower, but perhaps matters have changed.) Chris@82: Chris@82: * Generate assembly directly for more processors, or maybe fork gcc. =) Chris@82: Chris@82: * ensure that threaded solvers generate (block_size % 4 == 0) Chris@82: to allow SIMD to be used. Chris@82: Chris@82: * memoize triggen.