[Libre-soc-dev] FFT, DCT, REMAP
Luke Kenneth Casson Leighton
lkcl at lkcl.net
Tue Jun 22 12:13:20 BST 2021
https://www.nayuki.io/res/free-small-fft-in-multiple-languages/fft.py
rright.
so let us imagine that QTY2 halves of say length 32 FFT have already
been carried out.
if there was space for 64
:
* the bitreverse would have shifted up each address by one more bit
and ORed in 0 0 0 0... into the lower half of the indices and 1 1 1 1
.. into the upper.
* an extra while loop would have be run, multiplying and
add/subtracting between the two
however, because the offset between them is a fixed quantity when only
running the outer loop we may:
* load the first part of the first sub-FFT
* load the first part of the 2nd sub-FFT
* perform mul-add-sub
* save
* load the next chunks
* repeat
the question is, then, how to work out which elements from the batches
need to be part-loaded?
l.
More information about the Libre-soc-dev
mailing list