Hi, I am trying to ponder out what dct32 is actually doing. The comment in matrixing() reads:
x[i] = x'[i + 16] i = 0..15 x[i + 17] = -x'[31 - i] i = 0..15 x[i + 32] = -x'[16 - i] i = 0..15 x[i + 48] = -x'[i] i = 0..15 x[16] = 0
the hi and lo are only 16 entries, filter is 32. 'i + 32' and definately 'i + 48' go beyond these bounds.
Now, I understand matrix math, but the numbers just don't seem right here.